NOIP提高组【JZOJ4813】running

Description

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

Data Constraint

对于30% 的数据，1<=n<=100
对于60% 的数据，1<=n<=10^6
对于100% 的数据，1<=n<=10^9; 1<=m<=50; 1<=ai<=n

Solution

显然第i个同学跳的格子可以看成是a$_i$*x+b=n *y,那么a$_i$*x-n *y=b,设$gcd(a_i,n)$=p,那么有$p*(x*(a[i]/p)-n*(y/p)=b$。所以b一定是p的倍数。所以第i个同学跳的格子即gcd(a_i,n)的倍数。

那么怎样去掉所有被覆盖的格子呢？我们枚举n的约数d，假设存在一个i使$gcd(a_i,n)|d$,那么对于所有的$gcd(x,n)=d$都会被i覆盖，移一下就变成$gcd(x/d,n/d)=1$，数量即$\varphi{n/d}$。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=55,maxn1=1000005;
ll b[maxn],a[maxn],n,i,t,j,k,l,m,ans,x,d[maxn1],p;
bool bz,bz1[maxn1];
ll gcd(ll x,ll y){
    ll r=x%y,p=x,q=y;
    while (r) p=q,q=r,r=p%q;
    return q;
}
void pan(int i){
    bool bz=true;int j,t,p;
    for (j=1;j<=m;j++){
        if (i%a[j]) continue;
        bz=false;
        break;
    }
    if (bz) return;
    j=1;t=n/i;p=n/i;
    while (d[j]<=t){
        if (!(t%d[j])){
            p=p*(d[j]-1)/d[j];
            while (!(t%d[j])) t=t/d[j];
        }
        j++;
    }
    ans+=p;
}
int main(){
    freopen("running.in","r",stdin);freopen("running.out","w",stdout);
//  freopen("data.in","r",stdin);
    scanf("%d%d",&n,&m);bz1[1]=true;
    for (i=2;i<maxn1;i++){
        if (!bz1[i]) d[++d[0]]=i;
        for (j=1;j<=d[0];j++){
            if (i*d[j]>=maxn1) break;
            bz1[i*d[j]]=true;
            if (!(i%d[j]))break;
        }
    }
    for (i=1;i<=m;i++)
        scanf("%d",&a[i]),a[i]=gcd(a[i],n);
    t=sqrt(n);
    for (i=1;i<=t;i++){
        if (n%i) continue;
        bz=true;
        pan(i);
        pan(n/i);
    }
    ans=n-ans;
    printf("%d\n",ans);
}