poj 1091 (容斥原理+数论

思路：
a1x1+a2x2…anxn+an+1 M=1有解，首先要知道最大公约数定理，**gcd(a1,…an+1)=d的话，则必然可以找到x1…xn+1使得a1x1+…an+1M=d，**所以假如要满足原式=1，则左边必定可被d整除，右边是1也必定可被d整除，能整除1的只有1，所以d=1，这样就好做了，但是互质很难判断，所以正难则反，我们分解M的因子，容斥一下就好了，这里并不需要高精（本来也想写的x)

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=1e5;
int n,m,prime[maxn],cnt,j;
ll num;
ll mypow(ll a,ll b)
{
   ll ans=1;
   while(b){ 
      if(b&1)ans=ans*a;
      a=a*a;
      b>>=1;
   } 
   return ans;
}
void dfs(int pos,int sum,int li)
{
   if(!li)
   {
      if(j&1)num-=mypow((m/sum),n);
      else num+=mypow((m/sum),n); 
      return;
   } 
   if(pos==cnt+1)return;
   if(prime[pos]*sum<=m)
       dfs(pos+1,prime[pos]*sum,li-1);
   dfs(pos+1,sum,li);
}
int main()
{
   cnt=0;
   scanf("%d%d",&n,&m);
   num=mypow(m,n);
   int k=m;
   for(int i=2;i*i<=k;++i)
   {
      if(k%i==0)
      { 
         prime[++cnt]=i; 
         while(k%i==0)k/=i;
      }
   }
   if(k>1)prime[++cnt]=k;
   for(j=1;j<=cnt;++j)
      dfs(1,1,j); 
   printf("%lld\n",num);
   return 0; 
}