shuoj 72 容斥定理

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm> 
using namespace std;
int a[1000],num;
bool vis[25];
int Q[10000];
int wait[25];
int gcd(int a,int b)
{
	if(b==0) return a;
	else return gcd(b,a%b);
}
//int solve(int n)  //容斥定理，求
//{
//    int k,t,ans;
//    t = ans = 0;
//    Q[t++] = 1;
//    for(int i = 0; i < num; ++i)
//    {
//        k = t;
//        for(int j = 0; j < k; ++j)
//        {
//        	int temp=1;
//        	if(Q[j]<0)temp=-1;
//        	Q[t++] = -1*fabs(Q[j])/gcd(fabs(Q[j]),a[i])*a[i]*temp;
//		}
//    }
//    //A∪B∪C = A+B+C - A∩B - B∩C - C∩A + A∩B∩C
//    //Q数组存放的就是A∪B∪C右边边各项的因子数以及正负情况。
//    for(int i = 0; i < t; ++i)
//    {
//    	ans += n/Q[i];
//	}
//    //n/Q[i]累加起来就是A∪B∪C
//    return ans;
//}
int solve(int m)
{
    int que[10000],i,j,k,t=0,sum=0;
    que[t++]=-1;
    for(i=0;i<num;i++)
    {
        k=t;
        for(j=0;j<k;j++)
        {
        	que[t++]=que[j]*a[i]*(-1);
        	cout << que[t-1]<<'*'<<endl;
		}
    }
    for(i=1;i<t;i++)
        sum=sum+m/que[i];
    return sum;
}
int main()
{
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		memset(Q,0,sizeof(Q));
		memset(a,0,sizeof(a));
		for(int i=0;i<m;i++)
		{
			scanf("%d",&a[i]);
		}
		num=m;
		printf("%d\n",solve(n));
	} 
    return 0;
}