古代猪文题解

原题传送门

本题前置知识：

1.欧拉定理：对于正整数$a,n$, 若其互质，则对于任何正整数$b$，有$a^b\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (n)}(\mathrm{m}\mathrm{o}\mathrm{d}n)$

2.$Lucas$定理：若$p$为质数，则对于任意整数$1\le m\le n$ 有：$C_n^m\equiv C_{n\mathrm{m}\mathrm{o}\mathrm{d}p}^{m\mathrm{m}\mathrm{o}\mathrm{d}p}\ast C_{n/p}^{m/p}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

3.若$n$为质数，则$\phi (n)=n-1$

证明：这样，然后这样，再这样，最后这样。

证毕

回到题目，给定整数$q$,$n$,计算
$$q^{{\sum }_{d|n}{C}_n^d}\mathrm{m}\mathrm{o}\mathrm{d}999911659$$
若$q=999911659$上式显然为$0$。由于$999911659$为质数，所以$999911659$和$q$互质。

由欧拉定理得：
$$q^{{\sum }_{d|n}{C}_n^d}\equiv q^{{\sum }_{d|n}{C}_n^d\mathrm{m}\mathrm{o}\mathrm{d}999911658}(\mathrm{m}\mathrm{o}\mathrm{d}999911659)$$
然后就是要求出
$$\sum _{d|n}{C}_n^d\mathrm{m}\mathrm{o}\mathrm{d}999911658$$
但是$999911658$并不是一个质数，这时$Lucas$定理就没法直接用了。

所以直接搞肯定不行了。(不然怎么是SDOI的题)

原以为这题搞不了了。。。。。

然后突然知道了个方法：分解质因数

对于$999911658$，对其进行质因数分解，$999911658=2\times 3\times 4679\times 35617$

然后枚举$n$的约数$d$，运用$Lucas$定理求出$C_n^d$的值，分别求出${\sum }_{d|n}{C}_n^d$对$2,3,4679,35617$四个质数取模之后的值

将其记为$a_1,a_2,a_3,a_4$

对于质数$p$，我们可以求出$p$以内的所有阶乘以及阶乘的模$p$的乘法逆元。用于表示$p$进制下的组合数。($Lucas$定理中的)

最后就是用中国剩余定理求出线性同余的解：
$$\{\begin{array}{l}x\mathrm{m}\mathrm{o}\mathrm{d}2=a_1\\ x\mathrm{m}\mathrm{o}\mathrm{d}3=a_2\\ x\mathrm{m}\mathrm{o}\mathrm{d}4679=a_3\\ x\mathrm{m}\mathrm{o}\mathrm{d}35617=a_4\end{array}$$

最终就能求得${\sum }_{d|n}{C}_n^d\mathrm{m}\mathrm{o}\mathrm{d}999911658$的最小正整数解$x$。再用快速幂可以求出最终的解。

这道题呢是一道综合性很强的题目，所需要的知识也相当的多。是一道不可多得的好题。

最后附上三个前置知识的无法完全保证正确性的证明：

1.欧拉定理：

设$b=q\ast \phi (n)+r$,其中$0\le r<\phi (n)$,即$r=b\mathrm{m}\mathrm{o}\mathrm{d}\phi (n)$

那么$a^b\equiv a^{q\ast \phi (n)+r}\equiv (a^{\phi (n)}{)}^q\ast 1^q\ast a^r\equiv a^r\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (n)}(\mathrm{m}\mathrm{o}\mathrm{d}n)$

2.$Lucas$定理：

$C_m^n\equiv C_{\lfloor \frac{m}{p}\rfloor }^{\lfloor \frac{n}{p}\rfloor }\ast C_{m\mathrm{m}\mathrm{o}\mathrm{d}p}^{n\mathrm{m}\mathrm{o}\mathrm{d}p}\mathrm{m}\mathrm{o}\mathrm{d}p$

对于任意的质数$p$:

$C_p^n\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p(n\ne 0且n\ne p)$

对于任意的一个实数$x$有：
$$(x+1{)}^p=\sum _{i=0}^p{C}_p^i{x}_i$$
在模$p$意义下有：

$(x+1{)}^p\equiv (x^p+1)\mathrm{m}\mathrm{o}\mathrm{d}p$

对于任意一个整数$m$:

$(x+1{)}^m=(x+1{)}^{\lfloor \frac{m}{p}\rfloor p}\ast (x+1{)}^{m\mathrm{m}\mathrm{o}\mathrm{d}p}$

$(x+1{)}^m=(x^p+1{)}^{\lfloor \frac{m}{p}\rfloor }\ast (x+1{)}^{m\mathrm{m}\mathrm{o}\mathrm{d}p}$

再用二项式定理展开：
$$\sum _{i=0}^m{C}_m^i{x}^i=\left(\sum _{i=0}^{\lfloor \frac{m}{p}\rfloor }{C}_{\lfloor \frac{m}{p}\rfloor }^i{x}^{p_i}\right)\left(\sum _{i=0}^{m\mathrm{m}\mathrm{o}\mathrm{d}p}{C}_{m\mathrm{m}\mathrm{o}\mathrm{d}p}^i{x}^i\right)$$
等号左边$i=n$时，等号右边唯一能组合出来$x^n$的就是$x^{\lfloor \frac{n}{p}\rfloor \ast p}$和$x^{n\mathrm{m}\mathrm{o}\mathrm{d}p}$

3.显然成立

最后奉上又臭又长的code：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int twx=+100;
const int inf=0x3f3f3f3f;
int MOD=999911658;
ll read()
{
    ll sum=0;
    ll flag=1;
    char c=getchar();
    while(c<'0'||c>'9')
    {
        if(c=='-')
        {
            flag=-1;
        }
        c=getchar();
    }
    while(c>='0'&&c<='9')
    {
        sum=((sum*10)+c-'0');
        c=getchar();
    }
    return sum*flag;
}
const int a[4]={2,3,4679,35617};
int p[36000];
int b[4];
int n;
int g;
int ans;
ll x,y;
ll POW(ll a,ll b,ll MOD)
{
    ll res=1;
    a%=MOD;
    for(;b;b>>=1)
    {
        if(b&1)
        {
            res=res*a%MOD;
        }
		a=a*a%MOD;
    }
    return res%MOD;
}
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
    ll res=exgcd(b,a%b,x,y);
	ll t=x;
	x=y;
	y=t-(a/b)*y;
	return res;
}
ll gcd(ll x,ll y)
{
    return !x?y:gcd(y%x,x);
}
ll inv(ll a,ll b)
{
    ll x,y;
    exgcd(a,b,x,y);
    return (x%b+b)%b;
}
ll calc(ll x,ll mod)
{
	ll ans=1;
    ll y;
    ll a,b;
	for(y=n;x;x/=mod,y/=mod)
	{
		a=x%mod;
        b=y%mod;
		ans=(ll)ans*p[b]%mod*inv(p[a],mod)%mod*inv(b<a?0:p[b-a],mod)%mod;
	}
	return ans;
}
void init()
{
	n=read();
    g=read();
}
void work()
{
	g%=MOD+1;
	if(!g)
    {
        puts("0");
        exit(0);
    }
    int i;
    for(p[0]=i=1;i<=a[3];i++)
    {
        p[i]=(ll)p[i-1]*i%MOD;
    }
    for(int i=1;i<=sqrt(n);i++)
    {
        if(n%i==0)
        {
            for(int j=0;j<=3;j++)
            {
                b[j]=(b[j]+calc(i,a[j]))%a[j];
            }
            if(i*i!=n)
            {
                for(int j=0;j<=3;j++)
                {
                    b[j]=(b[j]+calc(n/i,a[j]))%a[j]; 
                }
            }
        }
    }
    for(int i=0;i<=3;i++)
    {
        exgcd(MOD/a[i],a[i],x,y);
		ans=(ans+(ll)x*(MOD/a[i])%MOD*b[i])%MOD;
    }
	ans=(ans+MOD)%MOD;
    MOD++;
	ans=POW(g,ans,MOD);
	printf("%d\n",ans);
}
int main()
{
	//freopen("test.in","r",stdin);
	//freopen("test.out","w",stdout);
    init();
    work();
	return 0;
}