Divisors——poj 1992

题目链接：

    http://poj.org/problem?id=2992

题目大意：

    题目描述很清晰，就是求组合数C(n, k)的因子个数

思路:

    (1)  一个正整数N可以进行质因数分解：N = p1^a1*p2^a2*...*pn^an

    (2)N的因子个数为(a1+1)*(a2+1)*...*(an+1)

    (由于组合数很大，还有多组输入所以这样直接做会TLE，所以要先预处理一下）

    (3)用一个数组prime[]存素数， prime[i]表示第i个素数，定义一个数组c[][]， c[n][k]表示n!中能分解出多少个prime[i]

    (4)众所周知，C(n, k) = n!/((n-k)!*k!),则，C(n, k)能分解出prime[i]的个数是c[n][i] - c[n-k][i] - c[k][i];

题目代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;
const int M = 435;
long long prime[M];
bool visit[M];
long long c[M][M];
void init()
{
	int cnt = 0;
	memset(visit, true, sizeof(visit));
	for(int i = 2; i < M; i++)
	{
		if(visit[i])
		{
			prime[cnt++] = i;
			for(int j = i*i; j < M; j += i)
			{
				visit[j] = false;
			}
		}
	}
	
	prime[cnt++] = 1e9 + 7;
	memset(c, 0, sizeof(c));
	for(int i = 0; i < M; i++)
	{
		for(int j = 0;prime[j] <= i; j++)
		{ 
		
			long long temp = i;
			long long n = 0;
			while(temp)
			{
				temp /= prime[j];
				n += temp;
			}
			c[i][prime[j]] = n; 
		}
	//	cout<<i<<endl;
	}
}
int main()
{
	init();
	long long  n, k;
	while(scanf("%lld %lld", &n, &k) != EOF)
	{
		long long  ans = 1;
		for(int i = 0; prime[i] <= n; i++)
		{
			//cout<<c[n][prime[i]]<<endl;
			ans *=  (c[n][prime[i]] - c[k][prime[i]] - c[n-k][prime[i]] + 1);
		}
		cout<<ans<<endl;
	}
}