http://www.spoj.com/problems/NGM2/
容斥原理 第四题 水!
难题各种各样,水题都一个样,无语。
代码:
#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<set>
#include<map>
#include<cmath>
#include<algorithm>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
//#define ull unsigned long long
#define ll long long
using namespace std;
const int INF=0x3f3f3f3f;
const int MOD=1000000007;
const ll LMOD=1000000007;
const double eps=1e-6;
const int N=1050;
ll a[N];
ll gcd(ll x,ll y)
{
if(x%y==0)
return y;
return gcd(y,x%y);
}
ll lcm(ll x,ll y)
{
return x*y/gcd(x,y);
}
ll solve(ll n,ll *a,int m)
{
ll sum=0;
for(int k=1;k<(1<<m);++k)
{
int num=0;
ll LCM=1;
for(int i=0;i<m;++i)
{
if((k&(1<<i))>0)
{
++num;
LCM=lcm(LCM,(ll)a[i]);
if(LCM>n)
break;
}
}
if(LCM>n)
continue;
if((num&1)>0)
{
sum+=(n/LCM);
}else
{
sum-=(n/LCM);
}
}
return sum;
}
int main()
{
ll n;
int k;
while(cin>>n>>k)
{
for(int i=0;i<k;++i)
cin>>a[i];
cout<<(n-solve(n,a,k))<<endl;
}
return 0;
}