poj 1091 跳蚤(分解质因数+容斥)

          首先满足要求的解要满足gcd(a1,a2,a3……an,M)=1，也就是求n个小于等于M的数与M的最大公约数的结果为1的组合的个数。找出M的所有质因子，通过容斥原理计算组合总数，所有的组合数为M^N，现在要减去所有组合的结果与M求最大公约数不为1的个数，这样的组合数必须N个数只要含有M的一个质因子，枚举质因子的组合情况，利用容斥原理，奇数个质因子减去，偶数个质因子加上。。。就能得到结果了。还有一点……看到网上有一个代码用了大数，因为M^N会超出long long的范围，但是网上的其他代码基本都是直接算的，没有考虑这一点，是因为这个题数据比较水？我认为不是的，其实在long long范围直接做的话，如果溢出相当于对2^64取模，只要最终的结果在2^64范围内，那么结果就是正确的。因此我尝试用了一个感觉最终结果不会超过的数实验了下，每一步都对这个数取模，然后发现也是过掉了。。。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
int d[100],a[100],cnt;
int n,m;
const ll mod=((1LL)<<61)+123456;
void getd(int m)
{
    cnt=0;
    for(ll i=2;i*i<=m;++i)
    {
        if(m%i==0)
        {
            d[cnt++]=i;
            while(m%i==0) m/=i;
        }
    }
    if(m>1) d[cnt++]=m;
}
ll pows(ll x,ll n)
{
    ll res=1;
    for(int i=0;i<n;++i)
        res=res*x%mod;
    return res;
}
ll dfs(int pos,int sp,int tot)
{
    if(pos==tot)
    {
        ll tmp=m;
        for(int i=0;i<tot;++i)
            tmp/=a[i];
        return pows(tmp,n);
    }
    ll sum=0;
    for(int i=sp;i<cnt;++i)
    {
        a[pos]=d[i];
        sum=(sum+dfs(pos+1,i+1,tot))%mod;
    }
    return sum;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    scanf("%d%d",&n,&m);
    getd(m);
    ll ans=pows(m,n),tmp;
    for(int i=1;i<=cnt;++i)
    {
        tmp=dfs(0,0,i);
        if(i&1) ans=(ans-tmp)%mod;
        else ans=(ans+tmp)%mod;
    }
    printf("%lld\n",ans%mod);
    return 0;
}