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

题目

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

样例输入：
第一行两个整数n,m。
接下来一行有m 个正整数，代表a1; a2……am
6 1
2

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

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

---

剖解题目

。。。。。。

---

思路

比赛时，很明显的想到了对于i这个人，所有是gcd(ai,n)的倍数的格子都会到达。
一下子就想到容斥原理，只需要将这m个数容斥处理一下就行了。
然而，发现容斥的理论复杂度是O(2^50)，这肯定会炸，然后就炸了。

---

解法

部分分就不说了。
解法一：
容斥原理，如上述思路，在容斥时进行优化。对于重复的，成倍数关系的数都删去。
然后dfs(v,1,lcm)−>dfs(v+1,1,lcm),dfs(v+1,−1,lcm∗a[v+1]/gcd(lcm,a[v+1]))这是互为相反的贡献，就不用继续往下弄了。（我不打这解法，是copy过来的）
解法二：
枚举n的约数d，如果存在一个ai使得gcd(ai,n)|d，那么所有gcd(i,n)=d的格子皆是能到达的。
那么就求出gcd(i,n)=d的数的个数。
也就是求gcd(i/d,n/d)=1的个数。
也就是求$\phi(n/d)$。
给answer减去$\phi(n/d)$。

---

代码

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)

using namespace std;

int a[50],n,k,num,m,ans;
int hg[9]={2,3,5,7,11,13,17,19,23};

int gcd(int x,int y)
{
    if (!y) return x; else return gcd(y,x%y);
}
int phi(int x)
{
    if (x==1) return 1;
    int xx=x;
    fo(i,0,8){
        bool bz=true;
        while (x%hg[i]==0){
            x/=hg[i];
            if (bz) {
                xx=xx*(hg[i]-1)/hg[i];
                bz=false;
            }
        }
    }
    if (x!=1) xx=xx*(x-1)/x;
    return xx;
}
void doing(int d)
{
    fo(j,1,m)
    if (d%gcd(a[j],n)==0){
        ans+=phi(n/d);
        break;
    }
}
int main()
{
    freopen("running.in","r",stdin);
    freopen("running.out","w",stdout);
    scanf("%d%d",&n,&m);
    int nn=n;
    fo(i,1,m) scanf("%d",&a[i]);
    n=nn;
    fo(i,1,sqrt(n)){
        if(n%i==0) {
            doing(i);
            if (n%i==i) continue;
            doing(n/i);
        }
    }
    printf("%d",n-ans);
}