HDU 5514 Frogs 容斥

http://acm.hdu.edu.cn/showproblem.php?pid=5514

没有想到怎么利用gcd的性质来改进指数级别的容斥，还是没有理解容斥的思想，只知道基础的。

已经知道结果就在因数之中，只要计算因数可以算几次就可以了。

代码：

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 10005;
int yinshu[maxn];
int a[maxn];
int gcd(int a,int b){
	return b==0?a:gcd(b,a%b);
}
int used[maxn],vis[maxn];
int main()
{
	int t,i,j;
	int n,m;
	scanf("%d",&t);
	int icase=0;
	while(t--){
		memset(used,0,sizeof(used));
		memset(vis,0,sizeof(vis));
		scanf("%d%d",&n,&m);
		for(i=0;i<n;i++)scanf("%d",&a[i]);
		for(i=0;i<n;i++)a[i]=gcd(a[i],m);
		int cnt=0;
		for(i=1;i*i<=m;i++){
			if(m%i==0){
				yinshu[cnt++]=i;
				if(i*i!=m)yinshu[cnt++]=m/i;
			}
		}
		sort(yinshu,yinshu+cnt);
		sort(a,a+n);
		for(i=0;i<cnt;i++){
			for(j=0;j<n;j++)
			if(yinshu[i]%a[j]==0){
				vis[i]=1;
				break;
			}
		} 
		ll ans=0;
		for(i=0;i<cnt;i++){
			if((vis[i]-used[i])==0)continue;
			ll tmp=(m-1)/yinshu[i];
			ans+=(tmp+1)*tmp/2*yinshu[i]*(vis[i]-used[i]);
			for(j=i+1;j<cnt;j++){
				if(yinshu[j]%yinshu[i])continue;
				used[j]+=vis[i]-used[i];
			}
		}
		printf("Case #%d: %lld\n",++icase,ans);
	}

   // cout << "Hello world!" << endl;
    return 0;
}