2015 沈阳 区域赛 F FROGS 莫比乌斯反演

呃 等重现放出来再说吧。。。

#include <bits/stdc++.h>
inline int gcd(int a,int b)
{
    return a==0?b:gcd(b%a,a);
}
#define N 1e8
using namespace std;

int n,m;
long long xx;
long long a[1010];
int fai[10000];
int check[10000];
int p[10000];
vector<long long> vec;

void getfai()
{
    fai[1]=1;
    int tot=0;
    for(int i=2; i<=10000; i++)
    {
        if(!check[i])
        {
            p[tot++]=i;
            fai[i]=i-1;
        }
        for(int j=0; j<tot; j++)
        {
            if(i*p[j]>N) break;
            check[i*p[j]]=1;
            if(i%p[j]==0)
            {
                fai[i*p[j]]=fai[i]*p[j];
                break;

            }
            else
            {
                fai[i*p[j]]=fai[i]*(p[j]-1);
            }
        }
    }
}

int main()
{
    //freopen("in.txt","r",stdin);
    int t;
    scanf("%d",&t);
    getfai();
    while(t--)
    {
        vec.clear();
        scanf("%d %lld",&n,&m);
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&xx);
            a[i]=(long long)gcd(xx,m);
        }
        sort(a+1,a+n+1);
        n=unique(a+1,a+n+1)-a-1;
        for(int i=2; i*i<=m; i++)
        {
            if(m%i==0)
            {
                vec.push_back(i);
                if(i*i!=m) vec.push_back(m/i);
            }
        }
        long long re=0;
        vector<long long>::iterator it;
        for(it=vec.begin(); it!=vec.end(); it++)
          for(int i=1; i<=n; i++)
            {
                long long x=*it;
                if(x%a[i]==0)
                {
                    re+=m*fai[m/x]/2;
                    break;
                }
            }
        printf("%lld\n",re);
    }
    return 0;
}