POJ 3904 JZYZOJ 1202 Sky Code 莫比乌斯反演 组合数

http://poj.org/problem?id=3904

 

题意：给一些数，求在这些数中找出四个数互质的方案数。

 

莫比乌斯反演的式子有两种形式http://blog.csdn.net/outer_form/article/details/50588307

这里用的是第二种形式。

求出四个数的公约数为x的倍数的方案数，即可得到，四个数的公约数为x的方案数。

这里x为1。

代码

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<cmath>
 #include<queue>
 using namespace std;
 const int maxn=;
 const long long modn=;
 int n;
 int a[maxn*]={};
 long long cnt[maxn]={},f[maxn]={};
 long long mic[maxn]={},su[maxn]={},tot=;
 bool vis[maxn]={};
 int main(){
     mic[]=;
     for(int i=;i<=maxn;i++){
         if(!vis[i]){
             su[++tot]=i;
             mic[i]=-;
         }
         for(int j=;j<=tot;j++){
             int k=i*su[j];
             if(k>maxn)break;
             vis[k]=;
             if(i%su[j])mic[k]=-mic[i];
             else break;
         }
     }
     while(~scanf("%d",&n)){
         int ma=;memset(cnt,,sizeof(cnt));
         for(int i=;i<=n;i++){
             scanf("%d",&a[i]);
             ma=max(a[i],ma);
             int w=sqrt(double(a[i]));
             for(int j=;j<=w;j++){
                 if(a[i]%j==){cnt[j]+=;cnt[a[i]/j]+=;}
             }
             if(w*w==a[i])cnt[w]--;
         }
         long long ans=;
         for(int i=;i<=ma;i++){
             if(cnt[i]<||mic[i]==)continue;
             ans+=mic[i]*cnt[i]*(cnt[i]-)*(cnt[i]-)*(cnt[i]-)/;
         }
         printf("%I64d\n",ans);
     }
 }