PKU 1091 跳蚤 数论 容斥原理

#include<iostream>
#include<cmath>
using namespace std;
#define MAXN 20
#define MAXP 200000
__int64 p[MAXP],tp,ans,n,m;
void get_pri(__int64 x) //求素因子
{
int i;
tp=0;
for(i=2;i*i<=x;i++)
{
if (x%i==0)
{
while(x%i==0)
x/=i;
p[++tp]=i;
}
}
if (x!=1) p[++tp]=x;
}
__int64 por(__int64 x,__int64 k)
{
__int64 ans=1;
int i;
for(i=1;i<=k;i++)
ans*=x;
return ans;
}
void gao(int k,__int64 v,int z)
{
if (k>tp) {if (v>1) ans-=z*por(m/v,n);return ;}
gao(k+1,v*p[k],-z);
gao(k+1,v,z);
}

int main()
{
while(scanf("%I64d %I64d",&n,&m)!=EOF)
{
get_pri(m);
ans=por(m,n);
gao(1,1,-1);
printf("%I64d/n",ans);
}
}