题意:
可转化为已知N+1个数,x1, x2, .. xn, m;使得gcd(x1, x2, ..., xn) = 1的解的个数(xi<=m)。
思路:m^n为所有不同序列个数,然后减去公因数>1且为m的素因子的解数,通过容斥原理思想,答案为
ans=m^n-(奇数个m的素因子相乘作为公因子的n元组)+(偶数个m的素因子相乘作为公因子的n元组)。
公因数为d的n元组,每个位置上有 (m/d)个选择(1 ~ m里面有m/d个d的倍数),根据乘法原理,可以得出有公因数d的n元组有 (m/d)^n 个.
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<stdio.h>
#include<math.h>
#include <string>
#include<string.h>
#include<map>
#include<queue>
#include<set>
#include<utility>
#include<vector>
#include<algorithm>
#include<stdlib.h>
using namespace std;
#define eps 1e-8
#define pii pair<int,int>
#define inf 0x3f3f3f3f
#define rd(x) scanf("%d",&x)
#define rd2(x,y) scanf("%d%d",&x,&y)
#define ll long long int
#define mod 1000000007
#define maxn 200100
#define maxm 10001005
int n,m;
int pr[maxn],pn;
ll pow(ll a,int b){
ll ans=1;
while(b){
if(b&1) ans*=a;
a*=a;
b>>=1;
}
return ans;
}
int main()
{
while(~rd2(n,m)){
ll sum=pow(m,n);
pn=0;
int mm=m;
for(int i=2;i*i<=mm;i++){//m的素因子
if(!(mm%i)) pr[pn++]=i;
while(mm%i==0) mm/=i;
}
if(mm!=1) pr[pn++]=mm;
for(int i=1;i<(1<<pn);i++){//数位枚举所有奇偶素因子情况
int k=1;
ll kk=1;
int nn=0;
for(int j=0;j<pn;j++){
if(k&i) {
kk*=pr[j];
nn++;
}
k<<=1;
}
if(nn%2) sum-=pow(m/kk,n);
else sum+=pow(m/kk,n);
}
printf("%I64d\n",sum);
}
return 0;
}