Poj3421 X-factor Chains 质因数分解 组合计数

一道组合计数题
题意:给定一个正整数X, 一个长度为m的X-因子链是由m+1个整数组成的。

其中 1 = X0, X1, X2, …, Xm = X 满足Xi < Xi+1 且 Xi 整除 Xi+1 。

现在要求X-因子链的最大长度和最大长度有多少条？

分析:一条最长链里面,要求是 $X_i+1$ % $X_i$ ==0 且 $X_i+1$ > $X_i$.
我们不难想到$X_i+1$是$X_i$乘上某个素数因子得到的.进一步地,如果想要让这个链的长度最大,那么唯一分解出来的每个素数因子和指数都要用上,用例子表达更加形象.
以$100$为例子,进行唯一分解,把它分解为若干个素数指数相乘的形式.
$$100=2^2\ast 5^2$$
最长链有 $1-2-4-20-100,1-2-10-20-100,1-5-10-50-100,1-5-25-50-100.共计四种,长度是4（不含1）$.
不难发现,答案的长度就是每个素数分解出来的指数相加，也就是
$$X={p1}^{a1}\ast {p2}^{a2}\ast p{3}^{a3}\ast \dots \ast {pn}^{an}$$.
该链的最长长度就是${\sum }_{i=1}^n{a}_i$.那么最长链有多少种呢,事实上就是高中的组合计数问题了,现在你有n种素数,每种素数有$a_i$个.假设最长链的长度是m.你的目标是往m个横线里入这全部的素数,有多少种组合方式.

如果现在有M条横线.现在我们想先放置第i个素数,把它放到$ai$个横线里面去，方案数是
C(M,ai).好了ai条横线用完,还剩余下M-ai条横线.
每个素数放置是乘法原理,应当相乘.
故答案为${\prod }_{i=1}^{n}C(M,ai)||M=m,m-=ai$
至于如何分解质因数,本题可以暴力打个素数表挨个找,因为数据弱.
我直接套用了Pollard-Rho.
AC代码:

/*
*/
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
const int maxn = 3e5+2;
const int INF = 1e9+7;
const ll MOD = 1e18;
	int prime[10]={2,3,5,7,11,13,19,61,29,37};
struct Prime{
	ll f_mul(ll x,ll y,ll mod){
    return ( x * y - (ll) ( (long double) x / mod*y )*mod + mod ) % mod;
	}
	ll f_pow(ll a,ll b,ll mod){
	ll ans = 1;
	while(b){
		if(b&1) ans= f_mul(ans,a,mod);
		b>>=1;a = f_mul(a,a,mod);
	}
	return ans%mod;
}
	ll gcd(ll a,ll b){
		if(b==0) return a;
		return gcd(b,a%b);
	}
	bool check(ll p,ll x){
		// to check  whether x is a prime.
		if(x % p ==0) return false;
		ll k = x - 1;
		ll t = f_pow(p%x,k,x);
		if(t!=1) return false;
		while(!(k&1)){ //一直到 k不是奇数为止 
			t = f_pow(p,k>>=1,x);
			if(t!=1&&t!=x-1) return false;
			if(t==x-1) return true;// t =x -1 而不是1的情况下 不能继续套用二次探索
			// 因为 x ^2  ≡1 是二次探索的右边 这时候右边是 x - 1 不能继续用了 认为通过测试. 
		}
		return true;
	}
	bool Miller_Rabin(ll x){
		for(int i=0;i<10;i++){
			if(prime[i]==x) return true;
			if(!check(prime[i],x)) return false;
		}
		return true;
	}
	vector<ll> div;//储存分解之后的质因数 是无序的  需要进行处理.
	ll ABS(ll x){
		return x>=0 ? x : (-x) ;
	}
    ll dfs(ll n){
    	ll x1 = 0 , x2 =0, c= rand()%(n-1) + 1;
    	ll val = 1;
    	for(ll i=1;;i<<=2,x1=x2,val=1){
    		for(ll j = 1;j<=i;++j){
    			x2 = (f_mul(x2,x2,n)+c)%n;
    			val =(f_mul(val,ABS(x1-x2),n))%n;
    			if(j%127==0){
    				ll d = gcd(val,n);
    				if(d>1) return d;
				}
			}
			ll d = gcd(val,n);
			if(d>1) return d;
		}
	}
    void fac(ll n){
    	if(n<2) return ;
    	if(Miller_Rabin(n)) {
    		div.push_back(n); return ;
		}
		ll p = n;
		while(p>=n) p = dfs(n);
		while(n%p==0) n/=p;
		fac(n);fac(p);
	}
		map<ll,int> cnt;
	void pollard(ll n,ll &ans){
		fac(n);
		sort(div.begin(),div.end());
		div.erase(unique(div.begin(),div.end()),div.end());
		for(vector<ll>::iterator it = div.begin();it!=div.end();it++){
			while(n%(*it)==0) n/=(*it),cnt[(*it)]++,ans++;
		}
	}
}p;
ll C(ll n,int k){
	//计算组合数C(n,k) 利用递推式子 C(n,k+1) = C(n,k) * (n-k) / (k+1)
	//令  k =0开始递推即可.用ans储存C(n,k) 
	ll ans = 1;
	for(ll i=1,j=n;i<=k;i++,j--){
		ans = ans * j / (i);
	}
	return ans;
}
int main(){
	ll x;
	while(cin>>x){
		p.div.clear();p.cnt.clear();
		ll ans1=0,ans2=1;
		//pollard - rho 分解数x 
		p.pollard(x,ans1);
		map<ll,int> &cnt = p.cnt;
		ll m = ans1;
		for(map<ll,int>::iterator it = cnt.begin();it!=cnt.end();it++){
			int k = (*it).second;
			ans2*=(C(m,k));
			m-=k;
		}
		cout<<ans1<<" "<<ans2<<"\n";
	}
return 0;}