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);
}
}