先给这道题的题目描述跪了(此处吐槽10分钟)
——————————————————————————
吐槽完毕,这道题首先看得出要算的结果是一个数的巨大无比次方,简单的快速幂显然不行了,但是最后的结果要mod一个质数,于是我们可以引入费马小定理:
a^(p-1)≡1(mod p)(当p为质数并且a不为p的倍数时成立),这能说明什么呢,就是说在mod p的意义下乘以一个数的p-1次方等于没乘!所以我们就可以将一个数巨大无比的次方mod(p-1)就行了。然而我们又发现了一个问题,这个次方是用组合数算出来的,由于组合数的递推公式无法使用(n,m过大),所以只能使用组合数的阶乘公式,然而就算我们用阶乘公式也要将阶乘从1一直乘到p-1,复杂度爆炸,所以我们又想到了另一个牛*的定理:中国剩余定理。
我们发现p-1恰好可以分解为35617*4679*3*2,那么我们只需算出这个次方分别在取余这4个质数之下的余数用中国剩余定理一合就可以了,中间还要穿插一下阶乘求逆(其实很简单)和鲁卡斯定理,鲁卡斯定理的内容就是从n个数中选出m个数的方案mod p==从n/p个数中选出m/p方案数*从n mod p个数中选出m mod p个数的方案,这样这道题就解决了。
鲁卡斯定理有毒…
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<iostream>
#include<iomanip>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
const long long maxn=999911659;
long long ni[4][40000];
long long jc[4][40000];
long long modd[4];
long long x[4];
long long Mi[4];
long long M;
void jisuan(long long num)
{
long long mod=modd[num];
long long mid=1;
jc[num][0]=1;
jc[num][1]=1;
for(int i=2;i<mod;i++)
{
mid=(mid*i)%mod;
jc[num][i]=(jc[num][i-1]*i)%mod;
}
ni[num][mod-1]=mod-1;
for(int i=mod-2;i>=0;i--) ni[num][i]=(ni[num][i+1]*(i+1))%mod;
}
long long zuhe(long long p,long long q,long long o)
{
long long sum=1;
while(p && q)
{
long long x=p%modd[o],y=q%modd[o];
if(y>x) return 0;
sum=(sum*(jc[o][x]*(ni[o][x-y]*ni[o][y])%modd[o])%modd[o])%modd[o];
p/=modd[o];
q/=modd[o];
}
return sum;
}
long long qiuni(long long yuan,long long p)
{
long long t=p-2;
long long mid=yuan;
long long ans=1;
while(t)
{
if(t&1) ans=(ans*mid)%p;
mid=(mid*mid)%p;
t>>=1;
}
return ans;
}
long long china()
{
long long t=0;
for(int i=0;i<=3;i++)
{
long long mid=(x[i]*Mi[i])%M;
mid=(mid*qiuni(Mi[i],modd[i]))%M;
t=(t+mid)%M;
}
return t;
}
int main()
{
long long n,g;
scanf("%lld%lld",&n,&g);
if(g%maxn==0) {cout<<0;return 0;}
modd[0]=35617;
modd[1]=4679;
modd[2]=3;
modd[3]=2;
jisuan(0);jisuan(1);jisuan(2);jisuan(3);
long long gen=sqrt(n);
for(int i=1;i<=gen;i++)
if(n%i==0)
{
for(int j=0;j<=3;j++)
{
x[j]=(x[j]+zuhe(n,i,j))%modd[j];
if(i*i<n) x[j]=(x[j]+zuhe(n,n/i,j))%modd[j];
}
}
M=modd[0]*modd[1]*modd[2]*modd[3];
for(int i=0;i<=3;i++) Mi[i]=M/modd[i];
long long goal=china();
long long t=goal;
long long ans=1;
long long mid=g;
while(t)
{
if(t&1) ans=(ans*mid)%maxn;
mid=(mid*mid)%maxn;
t>>=1;
}
cout<<ans<<endl;
return 0;
}