hdu5072 容斥+枚举

这题说的是给了 n 个数字 每个数值大于1 小于100000，n小于100000 ，找出满足下面要求的三人组有多少种 比如abc ( (ab)==(bc)==(ac) ==1 )||( (ab)!=1&&(bc)!=1&&(ac)!=1 )

(()----表示gcd )计算出这样的三元组的个数。

这样考虑 三元组必须满足上面的要求 ， 那么我们减去不满足的可能情况 abc 计算与 a互质的个数 和与a不互质的个数 得到了b在内的不可行的方案数，然后这样发现会多算一次 如果abc

(ab)==1 (ac)!=1 在计算c的时候又会计算一次。 然后我们得到了所有的 不肯能方案除以2 得到了不可能的方案数  总方案数减去不可能方案数就得到了（因为数比较小用容斥便可以计算出答案）

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include<string.h>
#include <vector>
using namespace std;
typedef pair<int,int> PF;
bool vis[];
int prime[],num;
void pri_p(){
   num=;
   memset(vis,false,sizeof(vis));
   for(long long i=; i<=; ++i)
    if(vis[i]==false){
       prime[num++]=i;
       for(long long j=i*i; j<=; j+=i) vis[j]=true;
   }
}
int A[];
int P[];
int sizea[];
int D[];
vector<PF>F[];
void fenjie(int loc,int len){
    for(int s=; s<(<<len); ++s){
       int ss=,num=;
       for(int i=;i<len ;++i){
           if(s&(<<i)){
             ss*=P[i]; num++;
           }
       }
       F[loc].push_back(make_pair(ss,num));
   }
   sizea[loc]=F[loc].size();
}
void inti(){
    pri_p();
    for(int i=; i<=; ++i){
            int len=;
            int s=i,num=;
            while(s>&&prime[num]!=){
                 if(s%prime[num]==){
                    P[len++]=prime[num];
                    while(s>&&s%prime[num]==){
                         s/=prime[num];
                    }
                 }
                 num++;
            }
        fenjie(i,len);
    }
}
void sum(int loc){
      int sz=sizea[loc];
      for(int i=; i<sz; ++i){
           PF t=F[loc][i];
           D[t.first]++;
      }
}
int jud(int loc){
    int ans=;
    for(int i=; i<sizea[loc]; ++i){
        PF f=F[loc][i];
        if(f.second%==)ans+=D[f.first];
        else ans-=D[f.first];
    }
    return ans;
}
int main()
{
    inti();
    int cas;
    scanf("%d",&cas);
    for(int cc=; cc<=cas; ++cc){
        long long n;
        scanf("%I64d",&n);
        memset(D,,sizeof(D));
        for(int i=; i<n; ++i){
            scanf("%d",&A[i]);
            sum(A[i]);
        }
        long long ans=;
        for(int i=; i<n; ++i){
         long long    S= jud(A[i]);
             long long noeven = n-S;
             if(S<=||noeven==) continue;
             ans = ans+noeven*(S-);
        }
        printf("%I64d\n",n*(n-)*(n-)/-ans/);
    }
    return ;
}