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)
所以
aix≡gcd(ai,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;
}