JZOJ4813. 【NOIP2016提高A组五校联考2】

Description

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

Input

第一行两个整数n,m。
接下来一行有m 个正整数，代表a1; a2……am

Output

输出一个整数，代表永远不会被同学跑到的格子的数目。

Sample Input

6 1
2

Sample Output

3

Data Constraint

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

分析

对于每一个人i，他可以走到的格子一定是$gcd(a_i,n)$的倍数。
所以我们枚举n的约数d，如果有一个i，$gcd(a_i,n)|d$,说明所有gcd(j,n)=d 的格子都能被到达，答案加上$φ(\dfrac{n}{d})$即可。

code（c++）

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h> 
using namespace std;
int n,m,a[53],ans;
int euler(int n)
{ 
     int res=n,a=n;  
     for(int i=2;i*i<=a;i++){  
         if(a%i==0){  
             res=res*(i-1)/i;
             while(a%i==0) a/=i;  
         }  
     }  
     if(a>1) res=res/a*(a-1);  
     return res;  
}  
int gcd(int x,int y){if(x%y==0)return y;else return gcd(y,x%y);}
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]);
    for(int i=1;i<=trunc(sqrt(n));i++)
    {
        if(n%i==0)
        {
            for(int j=1;j<=m;j++)
                if(i%gcd(a[j],n)==0)
                {
                    ans+=euler(n/i);
                    break;
                }
            if(n/i!=i)
            {
                for(int j=1;j<=m;j++)
                    if((n/i)%gcd(a[j],n)==0)
                    {
                        ans+=euler(i);
                        break;
                    }
            }
        }
    }
    printf("%d\n",n-ans);
}