容斥原理1004 POJ 1091

题意:
Z城市居住着很多只跳蚤。在Z城市周六生活频道有一个娱乐节目。
一只跳蚤将被请上一个高空钢丝的正中央。钢丝很长，可以看作是无限长。
节目主持人会给该跳蚤发一张卡片。卡片上写有N+1个自然数。
其中最后一个是M，而前N个数都不超过M，卡片上允许有相同的数字。
跳蚤每次可以从卡片上任意选择一个自然数S，然后向左，或向右跳S个单位长度。
而他最终的任务是跳到距离他左边一个单位长度的地方，并捡起位于那里的礼物。
比如当N=2，M=18时，持有卡片(10, 15, 18)的跳蚤，
就可以完成任务：他可以先向左跳10个单位长度，然后再连向左跳3次，
每次15个单位长度，最后再向右连跳3次，每次18个单位长度。
而持有卡片(12, 15, 18)的跳蚤，则怎么也不可能跳到距他左边一个单位长度的地方。
当确定N和M后，显然一共有M^N张不同的卡片。现在的问题是，
在这所有的卡片中，有多少张可以完成任务。
思路:
因为POJ的数据相当的水…. 可以不用大数
根据题意,如果这张卡片可以
a1x1+a2x2+a3x3+….anxn=1必然有解
->gcd(a1,a2,a3,a4….an)=1
方案数->m^n-2|gcd(ai)-3|gcd(ai)-…pi|gcd(ai)+2*3|gcd(ai)+pi*pj|gcd(ai)….
奇加偶减就OK了

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 100005;
const int inf=(1<<28)-1;
LL quickpow(int a,int b)
{
    LL res=1,tmp=a;
    while(b)
    {
        if(b&1)
            res=(res*tmp);
        tmp*=tmp;
        b/=2;
    }
    return res;
}
int Fac[50];
void GetFac(int x)
{
    Fac[0]=0;
    for(int i=2;i*i<=x;++i)
    if(x%i==0)
    {
        Fac[++Fac[0]]=i;
        while(x%i==0) x/=i;
    }
    if(x>1) Fac[++Fac[0]]=x;
}
LL n,m;
LL dfs(int pos,LL Now,int step)
{
    LL Ans=0;
    if(pos==Fac[0]+1)
    {
        if(Now==1) return 0;
        Ans=quickpow(m/Now,n);
        if(step&1)
            return Ans;
        else
            return -Ans;
    }
    Ans+=dfs(pos+1,Now,step);
    Ans+=dfs(pos+1,Now*Fac[pos],step+1);
    return Ans;
}
int main()
{
    scanf("%lld%lld",&n,&m);
    LL Ans=quickpow(m,n);
    GetFac(m);
    Ans-=dfs(1,1,0);
    printf("%lld\n",Ans);
    return 0;
}