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(ai,n)
的倍数。
所以我们枚举n的约数d,如果有一个i,
gcd(ai,n)|d
,说明所有gcd(j,n)=d 的格子都能被到达,答案加上
φ(nd)
即可。
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);
}