POJ 2992 Divisors

最新推荐文章于 2020-04-30 08:36:37 发布

Napoleon2004

最新推荐文章于 2020-04-30 08:36:37 发布

阅读量161

点赞数

CC 4.0 BY-SA版权

分类专栏：数论

本文链接：https://blog.youkuaiyun.com/Napoleon2004/article/details/80297134

数论专栏收录该内容

21 篇文章

订阅专栏

本文介绍了一种利用素数筛法进行快速求解的算法，并通过组合计数的方法来解决特定数学问题。该算法首先使用筛法找出一定范围内的所有素数，然后通过递归方式计算组合数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#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;
}