【数论】ural1091

最新推荐文章于 2015-08-16 22:08:09 发布

穆林幕

最新推荐文章于 2015-08-16 22:08:09 发布

阅读量491

点赞数

分类专栏：数论文章标签：数学容斥原理

本文链接：https://blog.csdn.net/u010582475/article/details/40478003

版权

数论专栏收录该内容

19 篇文章 0 订阅

订阅专栏

题目连接 http://acm.timus.ru/problem.aspx?space=1&num=1091

题意就是求不大与s的数中有k个数具有相同的公因子的数目。当大于10000的时候直接输出10000

刚开始想的是水题，然后直接按照C（n,r）来计算的，后来发现有重复的。。不晓得肿么办了，众多大神纷纷说是容斥定理，莫比乌斯反演。。然后我就晕了。渣渣都没学过啊。。。

最后这个题可以水过的。。

而且，这个题中应该看到，所有可能使因数大于2的素数因子不大与23，那么就简化了很多运算，还有就是pirme[i]中，当cnt=s/prime[i];那么数中拥有prime[i]因子的数的个数就是cnt,直接C(cnt,k)就好，（我都没发现。。~~~~(>_<)~~~~ ）

不用再遍历s,然后求出cnt了。。。。。囧。。

还有最关键的去重，直接找到因子的公倍数，然后除去公倍数的计算就好了，这也算是容斥定理的思想了吧，但是还未用到容斥。

//求不超过S的，K个具有相同因子的集合的个数。------数论
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int prime[]={2,3,5,7,11,13,17,19,23};

int c[25][25]; //dp可以很快求出系数，杨辉三角 

void init()
{
	memset(c,0,sizeof(c));
	int i,j;
	for(i=1;i<=25;i++)
	{
		for(j=0;j<=i;j++)
		{
			if(i==j||j==0)c[i][j]=1;
			else 
			c[i][j]=c[i-1][j-1]+c[i-1][j];
		}
	}
}
int main()
{
	int s,k;
    int i,j;
	scanf("%d%d",&k,&s);
	int cnt,res=0,tmp;
	init();
	for(i=0;i<9;i++)
	{
		cnt=s/prime[i];
		if(cnt<k)
		{
			tmp=i;
			break;
		}
		else res+=c[cnt][k];
	} 
	for(i=0;i<9;i++)
	{
		for(j=i+1;j<9;j++)
		{
			cnt=s/(prime[i]*prime[j]);
			if(cnt>=k)res-=c[cnt][k];
		}
	}

	if(res<10000)
	cout<<res<<endl;
	else cout<<"10000"<<endl;
}

大神们都有着不必坚强的内心，渣渣我也要成为一名坚强的渣渣~~~哈哈哈哈\(^o^)/~