【JZOJ 4813】【NOIP2016提高组 五校联考2】running

Description

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

Analysis

n好大，m好小
或许你数学好，或者你找规律发现ai可以替换成gcd(ai,n)
证明如下：
根据裴蜀定理，有x,y，满足
$ax+by=gcd(a,b)$
那么
$ax≡gcd(a,b) (mod\ b)$
所以
$a_ix≡gcd(a_i,n) (mod\ n)$

可能你想容斥，但是m=50会爆。
正解很神奇
考虑枚举n的因数d，若有a[i]满足gcd(a[i],n)|d则说明所有gcd(编号,n)=d的格子都会被走到。而且这样算不会重复。

Code

#include<cstdio>
#include<cmath>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int N=55;
ll m,n,ans,a[N];
ll gcd(ll x,ll y)
{
    if(x%y==0) return y;
    else return gcd(y,x%y);
}
ll phi(ll x)
{
    ll t=x;
    for(int i=2;i*i<=x;i++)
        if(x%i==0)
        {
            t=t*(i-1)/i;
            while(x%i==0) x/=i;
        }
    if(x>1) t=t*(x-1)/x;
    return t;
}
void calc(ll x)
{
    fo(i,1,n)
        if(x%gcd(a[i],m)==0)
        {
            ans+=phi(m/x);
            break;
        }
}
int main()
{
    freopen("running.in","r",stdin);
    freopen("running.out","w",stdout);
    scanf("%lld %lld",&m,&n);
    fo(i,1,n)
    {
        scanf("%lld",&a[i]);
        a[i]=gcd(a[i],m);
    }
    fo(i,1,sqrt(m))
        if(m%i==0)
        {
            calc(i);
            if(m/i!=i) calc(m/i);
        }
    printf("%lld",m-ans);
    return 0;
}