题意：n个数，挑4个数互质 的对数；
莫比乌斯反演
和这个差不多：
63MS
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<iostream>
#define ll long long
using namespace std;
#define N 10000
int mib[N+10];
int f[N+10];
bool vis[N+10];
void mobius()
{
int i,j;
for(i=1;i<=N;i++) mib[i]=1,vis[i]=false;
for(i=2;i<=N;i++)
{
if(vis[i]) continue;
for(j=i;j<=N;j+=i)
{
vis[j]=true;
if((j/i)%i==0)
{
mib[j]=0;
continue;
}
mib[j]=-mib[j];
}
}
f[0]=mib[0];
for(i=1;i<=N;i++)
f[i]=f[i-1]+mib[i];
}
int p[N+10];
int cnt[N+10],num[N+10];
ll h[N+10];
int main()
{
mobius();
int i,j,k;
int n;
h[0]=h[1]=h[2]=h[3]=0;
h[4]=1;
for(i=5;i<=N;i++)
h[i]=h[i-1]*i/(i-4);
while(scanf("%d",&n)!=EOF)
{
memset(num,0,sizeof(num));
memset(cnt,0,sizeof(cnt));
int mx=0;
for(i=1;i<=n;i++)
{
scanf("%d",&p[i]);
num[p[i]]++;
mx=max(mx,p[i]);
}
for(i=1;i<=N;i++)
for(j=i;j<=N;j+=i)
cnt[i]+=num[j];
ll ans=0;
for(i=1;i<=mx;i++)
{
if(!cnt[i]) continue;
ans+=(ll)mib[i]*h[cnt[i]];
}
cout<<ans<<endl;
}
}
容斥原理：
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define ll long long
#define N 10010
int p[N],g[N];
int t[20];
ll h[N];
void solve(int x)
{
int i,j,k,cnt=0;
for(i=2;i*i<=x;i++)
{
if(x%i==0)
{
t[cnt++]=i;
while(x%i==0)
x/=i;
}
}
if(x!=1)
t[cnt++]=x;
int tt=(1<<cnt);
for(i=1;i<tt;i++)
{
int tsum=1,tflag=0;
for(j=0;j<cnt;j++)
{
if((i>>j)&1)
{
tsum*=t[j];
tflag++;
}
}
p[tsum]++;
g[tsum]=tflag;
}
}
int main()
{
int n;
ll i,j,k;
h[0]=h[1]=h[2]=h[3]=0;
for(i=4;i<=10000;i++)
{
h[i]=i*(i-1)*(i-2)*(i-3)/24;
}
while(scanf("%d",&n)!=EOF)
{
memset(p,0,sizeof(p));
memset(g,0,sizeof(g));
for(i=0;i<n;i++)
{
scanf("%d",&k);
solve(k);
}
ll res=h[n];
for(i=1;i<=10000;i++)
{
if(p[i]!=0)
{
if(g[i]&1)
res-=h[p[i]];
else
res+=h[p[i]];
}
}
cout<<res<<endl;
}
}
