【JZOJ4813】running

Description

小胡同学是个热爱运动的好孩子。
每天晚上，小胡都会去操场上跑步，学校的操场可以看成一个由n 个格子排成的一个环形，格子按照顺时针顺序从0 到n-1 标号。
小胡观察到有m 个同学在跑步，最开始每个同学都在起点（即0 号格子），每个同学都有个步长ai，每跑一步，每个同学都会往顺时针方向前进ai 个格子。由于跑道是环形的，如果
一个同学站在n-1 这个格子上，如果他前进一个格子，他就会来到0。
他们就这样在跑道上上不知疲倦地跑呀跑呀。小胡同学惊奇地发现，似乎有些格子永远不会被同学跑到，他想知道这些永远不会被任何一个同学跑到的格子的数目，你能帮帮他
吗？（我们假定所有同学都跑到过0 号格子）。

Solution

经过找规律或证明得：所有$gcd(a_i,n)$都会被走到。

证明略。

正解是枚举一个$n$的约数$d$，若$gcd(a_i,n)|d$，那么所有编号$pos$满足$gcd(pos,n)|d$会被走到，答案加上$\varphi(\dfrac{n}{d})$即可。

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
#define ll long long
#define N 51
using namespace std;
int a[N];
int gcd(int x,int y)
{
    int z;
    while(x%y!=0)
    {
        z=x%y;
        x=y;
        y=z;
    }
    return y;
}
int phi(int x)
{
    if(x==1) return 1;
    int p=2;
    int s1=1,s2=1;
    int s;
    int xx=x;
    while(p*p<=x)
    {
        if(xx%p==0) s1*=p-1,s2*=p;
        while(xx%p==0) xx/=p;
        p++;
    }
    if(xx>1) s1*=xx-1,s2*=xx;
    s=x/s2*s1;
    return s;
}
int main()
{
    freopen("running.in","r",stdin);
    freopen("running.out","w",stdout);
    int n,m;
    cin>>n>>m;
    fo(i,1,m) scanf("%d",&a[i]);
    int ans=0;
    for(int i=1;i*i<=n;i++)
    if(n%i==0)
    {
        fo(j,1,m)
        {
            int t=gcd(a[j],n);
            if(i%t==0)
            {
                ans+=phi(n/i);
                break;
            }
        }
        fo(j,1,m)
        {
            int t=gcd(a[j],n);
            if(n/i!=i && (n/i)%t==0)
            {
                ans+=phi(i);
                break;
            }
        }
    }
    cout<<n-ans;
}