题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=1796
题目大意:
给一个n,给m个数,求出1~n中能够被m中某个数整除的数的个数。
解题思路:
简单容斥,不解释。
先对m个除掉0,然后数排序,然后容斥。每次求出前面的最小公倍数。
注意:
10 2
2 11 keys:4
10 2
11 2 keys:4
代码:
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<stack>
#include<list>
#include<queue>
#define eps 1e-6
#define INF 0x1f1f1f1f
#define PI acos(-1.0)
#define ll __int64
#define lson l,m,(rt<<1)
#define rson m+1,r,(rt<<1)|1
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
//freopen("data.in","r",stdin);
//freopen("data.out","w",stdout);
ll pp[15],cnt,n;
ll ans;
ll gcd(ll a,ll b)
{
if(a%b==0)
return b;
return gcd(b,a%b);
}
ll lcm(ll a,ll b)
{
return a/gcd(a,b)*b;
}
void dfs(ll la,ll pos,ll flag)
{
if(pos>cnt)
return ;
for(int i=pos;i<=cnt;i++)
{
ll t=(n-1)/lcm(la,pp[i]);
if(flag)
ans+=t;
else
ans-=t;
dfs(lcm(la,pp[i]),i+1,flag^1);
}
}
int main()
{
while(~scanf("%I64d",&n))
{
cnt=0;
int m;
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
ll a;
scanf("%I64d",&a);
if(a)
pp[++cnt]=a;
}
sort(pp+1,pp+cnt+1); //排序 不排序是有问题的
ans=0;
dfs(1,1,1);
printf("%I64d\n",ans);
}
return 0;
}