51nod 1284【容斥】

1284 2 3 5 7的倍数

基准时间限制：1 秒 空间限制：131072 KB 分值: 5  难度：1级算法题

 收藏

 关注

给出一个数N，求1至N中，有多少个数不是2 3 5 7的倍数。 例如N = 10，只有1不是2 3 5 7的倍数。

Input

输入1个数N(1 <= N <= 10^18)。

Output

输出不是2 3 5 7的倍数的数共有多少。

Input示例

10

Output示例

1

思路：理解了容斥就很好写了，主要是奇加偶减；

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define max_n 10010
typedef long long LL;
using namespace std;

LL p[]={2,3,5,7},k;

LL nop(LL m)
{
	LL top=0,sum=0;
	LL que[max_n];
	que[top++]=-1;
	for(int i=0;i<k;i++)
	{
		LL t=top;
		for(int j=0;j<t;j++)
			que[top++]=que[j]*p[i]*(-1);
		
	}
	for(int i=1;i<top;i++)
		sum+=m/que[i];
	return sum;
}

int main()
{
	LL n;
	k=4;
	scanf("%lld",&n);
	printf("%lld\n",n-nop(n));
	return 0;
}

#include<cstdio>
#include<cstring>
typedef long long LL;

LL p[1010],k;

LL gcd(LL a,LL b)
{
	return b?gcd(b,a%b):a;
}

LL lcm(LL a,LL b) 
{
	return a/gcd(a,b)*b;
}

LL nop(LL m)
{
	LL res,sum=0,use; 
	for(LL i=1;i<1<<k;i++)
	{
		res=1,use=0;
		for(LL j=0;j<k;j++)
		{
			if(i & (1<<j))
			{
				use++;
				res=lcm(res,p[j]);
			}
		}
		if(use & 1) sum+=m/res;
		else sum-=m/res;
	}
	return sum;
}

int main()
{
	LL n;
	scanf("%lld",&n);
	k=4;
	p[0]=2,p[1]=3,p[2]=5,p[3]=7;
	printf("%lld",n-nop(n));
	return 0;
}