POJ-2992 Divisors(数学知识)

Divisors

Time Limit: 1000MS		Memory Limit: 65536K
Total Submissions: 12085		Accepted: 3600

Description

Your task in this problem is to determine the number of divisors of C_n^k. Just for fun -- or do you need any special reason for such a useful computation?

Input

The input consists of several instances. Each instance consists of a single line containing two integers n and k (0 ≤ k ≤ n ≤ 431), separated by a single space.

Output

For each instance, output a line containing exactly one integer -- the number of distinct divisors of C_n^k. For the input instances, this number does not exceed 2⁶³ - 1.

Sample Input

5 1
6 3
10 4

Sample Output

2
6
16

Source

CTU Open 2005

求一个数的因子个数 可以将它分解质因数分成若干个质数的幂的积，设他们的指数为a1,a2,a3……an 则这个数的因子个数为(a1+1)*(a2+1)*(a3+1)*……*(an+1)  当求a!中出现的多少个p的时候 a!可以表示成1*2*3*……*(p-1)*p*(p+1)*……*(2*p-1)*(2*p)*(2*p+1)*……其中p, 2*p , 3*p……会有因子p所以先加上a/p 但是p², 2*p²……这些数中有两个p    类似的p³ , 2*p³ , 3*p³ ……会有三个p  所以在处理的时候要考虑全 具体看代码

#include "bits/stdc++.h"
using namespace std;
typedef long long LL;
const int MAX=435;
int n,m,len,pri[MAX];LL ans[MAX][MAX];
bool t[MAX];
void prime(){
    register int i,j;
    memset(t,true,sizeof(t));
    for (i=2;i<MAX;i++){
        if (t[i]) pri[++len]=i;
        for (j=1;j<=len && pri[j]*i<MAX;j++) {t[pri[j]*i]=false; if (i%pri[j]==0) break;}
    }
}
inline int calc(register int x,register int y){
    register int an=0;
    while (x){an+=x/y;x/=y;}
    return an;
}
int main(){
    freopen ("divisors.in","r",stdin);freopen ("divisors.out","w",stdout);
    register int i,j,k;
    prime();ans[0][0]=1;
    for (i=1;i<MAX;i++){
        ans[i][i]=ans[i][0]=1;
        for (j=1;j<i;j++){
            ans[i][j]=1;
            for (k=1;k<=len;k++)
                ans[i][j]*=calc(i,pri[k])-calc(j,pri[k])-calc(i-j,pri[k])+1;
        }
    }
    while (~scanf("%d%d",&n,&m)) printf("%lld
",ans[n][m]);
    return 0;
}