POJ 2992 Divisors

#include<cstdio>
using namespace std;
const int M=431;
bool b[M+5];
int l,p[M+5],jie[M+5][M+5];
void jc(){
	int i,j;
	for(i=2;i<=M;i++)
		if(b[i]==0){
			for(j=i*i;j<=M;j+=i)
				b[j]=1;
			p[++l]=i;
		}
	for(i=1;i<=M;i++)
		for(j=1;j<=l;j++)
			jie[i][j]=i/p[j]+jie[i/p[j]][j];
}
long long cal(int n,int k){
	int i,x;
	long long ans=1;
	for(i=1;i<=l;i++){
		x=jie[n][i]-jie[k][i]-jie[n-k][i];
		if(x)ans*=(x+1);
	}
	return ans;
}
int main(){
	int n,k;
	jc();
	while(scanf("%d%d",&n,&k)==2)
		printf("%I64d\n",cal(n,k));
	return 0;
}