POJ 2992 : Divisors - 求组合数的因子个数

题意：

输入为n，k，求出C(n,k)(n>=k）的约数个数。 n，k<= 2⁶³ - 1

分析：

k，n比较大，不可能直接求组合数。

现在要求约数的个数，首先进行素因子分解，组合数写成阶乘形式，先对阶乘进行素因子分解。

步骤：

（1）筛选431内素数；

（2）对每个阶乘进行素因子分解；

（3）求组合数的素因子分解；

  得到约数个数。

n=p1^e1 * p2^e2 * ... * pn^en

其中p1,,p2,,pn为素数，

则n的约数的个数为(e1+1) * (e2+1) * ... * (en+1)

利用筛法,,我们可以求出1~n中每个数的所有素约数的指数

C(n,k)=n!/k!/(n-k)!

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#define LL long long
using namespace std;
int main(){
	bool pri[440];
	memset(pri,0,sizeof(pri));
	int p[85],i,j,k=0;
	for(i=2;i<440;i++){
		if(!pri[i]) p[k++]=i;
		for(j=0;j<k && i*p[j]<440;j++){
			pri[i*p[j]]=1;
			if(!(i%p[j])) break;
		}
	}

	int jie[440][86];//阶乘分解为素数积
	memset(jie,0,sizeof(jie));
	for(i=0;i<k;i++)
		for(j=2;j<440;j++)
			jie[j][i]=j/p[i]+jie[j/p[i]][i];

	LL C[432][432];//组合数的二维数组
	for(i=2;i<432;i++)
		for(j=1;j<i;j++){
			C[i][j]=1;
			for(int ii=0;ii<k && jie[i][ii];ii++){
				int side=jie[i][ii]-jie[j][ii]-jie[i-j][ii];
				if(side)
					C[i][j]*=(side+1);//求因子个数
			}
		}

	int n,x;
	while(~scanf("%d%d",&n,&x)){
		if(x==0||x==n){
			printf("1\n");
			continue;
		}
		printf("%lld\n",C[n][x]);
	}

	return 0;
}