poj3421 X-factor Chains——分解质因数

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
int x,p[25],pc[25],cnt;
void divide(int n)
{
    cnt=0;
    memset(pc,0,sizeof pc);
    for(int i=2;i*i<=n;i++)
        if(n%i==0)
        {
            p[++cnt]=i;
            while(n%i==0)pc[cnt]++,n/=i;
        }
    if(n>1)p[++cnt]=n,pc[cnt]=1;
}
ll C(int n,int m)
{
    if(n<m)return 0;
    m=min(m,n-m);
    ll a=1,b=1;
    for(int i=n-m+1;i<=n;i++)a*=i;
    for(int i=2;i<=m;i++)b*=i;
    return a/b;
}
int main()
{
    while(~scanf("%d",&x))
    {
        divide(x);
        int sum=0; ll ans=1;
        for(int i=1;i<=cnt;i++)sum+=pc[i];
        int n=sum;
        for(int i=1;i<=cnt;i++)
        {
            ans*=C(n,pc[i]);
            n-=pc[i];
        }
        printf("%d %lld\n",sum,ans);
    }
    return 0;
}