[JZOJ4813]running

题目大意

一个有$n$个位置的环，从$0$到$n-1$编号。有$m$个人，初始时都在位置$0$，第$i$个人不停向前走$a_i$个位置（到尽头循环）。
求有多少个位置是至始至终没人走过的。

$1\le n\le 10^9,1\le m\le 50,1\le a_i\le n$

---

题目分析

显然第$i$个人走过的位置一定是$\gcd(a_i,n)$的倍数。
正难则反，我们用总位置数减去一定会走过的位置数。现在问题在于如何不重不漏地计算出所有的位置。
枚举$n$的因数$d$，当$\exists a_i,\gcd(a_i,n)|d$的时候，$d$的倍数所在的位置一定会被走过。但是这里我们并不筛调$d$的所有倍数，因为这样会算重。我们只将$\gcd(x,n)=d$的位置$x$筛去。这样的位置显然有$\varphi(\frac{n}{d})$个，加进一定会走过的位置数量里面。
为什么这样就会不重不漏呢？首先位置$x$如果会被经过，只会被$\gcd(x,n)$筛掉，这保证了唯一性。其次，如果$x$会被经过，那么$\gcd(x,n)$一定会在枚举$n$的因子的时候枚举到。令$n=l\gcd(a_i,n)$，$\because\exists a_i,x=k\gcd(a_i,n)\therefore\gcd(x,n)=\gcd(a_i,n)\gcd(k,l)\therefore \gcd(a_i,n)|\gcd(x,n)$。这样保证了存在性。
然后就没有啦[]~(￣▽￣)~*

---

代码实现

#include <algorithm>
#include <iostream>
#include <cstdio>
#include <map>

using namespace std;

const int S=10000000;
const int M=55;

map<int,int> g;

int pri[S+5],phi[S+5],f[S+5];
int n,m,ans;
int a[M];

int Phi(int x)
{
    int ret=1;
    for (int i=1;i<=pri[0]&&x>S;i++)
        if (!(x%pri[i]))
        {
            int p=pri[i]-1,cnt=0;
            while (!(x%pri[i])) cnt++,x/=pri[i];
            for (int j=1;j<cnt;j++) p*=pri[i];
            ret*=p;
        }
    return ret*phi[x];
}

int gcd(int x,int y){return y?gcd(y,x%y):x;}

bool check(int x)
{
    for (int i=1;i<=m;i++) if (!(x%gcd(a[i],n))) return 1;
    return 0;
}

void solve()
{
    phi[1]=1;
    for (int i=1;i<=S;i++)
    {
        if (!phi[i]) phi[i]=i-1,pri[++pri[0]]=f[i]=i;
        for (int j=1;j<=pri[0];j++)
        {
            if (1ll*pri[j]*i>S) break;
            phi[pri[j]*i]=f[i]==pri[j]?phi[i]*pri[j]:phi[i]*(pri[j]-1);
            f[pri[j]*i]=pri[j];
            if (!(i%pri[j])) break;
        }
    }
    for (int d=1;1ll*d*d<=n;d++)
        if (!(n%d))
        {
            if (check(d))ans+=Phi(n/d);
            if (1ll*d*d!=n&&check(n/d)) ans+=Phi(d);
        }
    ans=n-ans;
}

int main()
{
    freopen("running.in","r",stdin),freopen("running.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++) scanf("%d",&a[i]);
    solve();
    printf("%d\n",ans);
    fclose(stdin),fclose(stdout);
    return 0;
}