poj 1091（枚举中二进制的巧妙应用，最大公约数，容斥原理）

#include<iostream>
//https://www.cnblogs.com/shihuajie/archive/2013/03/25/2980709.html
using namespace std;
long long pow(long long x,int y){
    long long ans = 1;
    while(y--){
        ans *= x;
    }
    return ans;
}
int main(){
    int i,j,k,n,t,tot,cnt,prime[10];
    long long m,ans;
    scanf("%d%lld",&n,&m);
    ans = pow(m,n);
    tot = 0;
    t = m;
    for(i=2;i*i<=m;i++){
        if(t%i==0){
            prime[tot++] = i;
            while(t%i==0){
                t/=i;
            }
        }
    }
    if(t!=1){
        prime[tot++] = t;
    }
    for(i=1;i<(1<<tot);i++){
        k = 1;
        cnt = 0;
        for(j=0;j<tot;j++){
            if(i&(1<<j)){
                cnt++;
                k *= prime[j];
            }
        }
        if(cnt%2==1)ans -= pow(m/k,n);
        else ans += pow(m/k,n);
    }
    printf("%lld\n",ans);
    return 0;
}