【51Nod】1284 - 2 3 5 7的倍数（容斥原理 & 二进制优化）

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

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

10

1

#include <cstdio>
#include <stack>
#include <queue>
#include <cmath>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
#define CLR(a,b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define LL long long
int main()
{
	LL n;
	int num[4] = {2,3,5,7};
	scanf ("%lld",&n);
	LL ans = n;
	int i = 1;
	while (i < (1 << 4))
	{
		int mul = 1;
		int ant = 0;
		int t = i;
		for (int j = 0 ; j < 4 ; j++)
		{
			if ((1 << j) & i)
			{
				ant++;
				mul *= num[j];
			}
		}
		if (ant & 1)		//奇数
			ans -= n / mul;
		else
			ans += n / mul;
		i++;
	}
	printf ("%lld\n",ans);
	return 0;
}