古代猪文 [Bzoj 1951，Sdoi 2010]

题目地址请点击——

---

古代猪文

---

Description

“在那山的那边海的那边有一群小肥猪。他们活泼又聪明，他们调皮又灵敏。他们自由自在生活在那绿色的大草坪，他们善良勇敢相互都关心……” ——选自猪王国民歌
很久很久以前，在山的那边海的那边的某片风水宝地曾经存在过一个猪王国。
猪王国地理位置偏僻，实施的是适应当时社会的自给自足的庄园经济，很少与外界联系，商贸活动就更少了。
因此也很少有其他动物知道这样一个王国。
猪王国虽然不大，但是土地肥沃，屋舍俨然。
如果一定要拿什么与之相比的话，那就只能是东晋陶渊明笔下的大家想象中的桃花源了。
猪王勤政爱民，猪民安居乐业，邻里和睦相处，国家秩序井然，经济欣欣向荣，社会和谐稳定。
和谐的社会带给猪民们对工作火红的热情和对未来的粉色的憧憬。
小猪iPig是猪王国的一个很普通的公民。
小猪今年10岁了，在大肥猪学校上小学三年级。
和大多数猪一样，他不是很聪明，因此经常遇到很多或者稀奇古怪或者旁人看来轻而易举的事情令他大伤脑筋。
小猪后来参加了全猪信息学奥林匹克竞赛(Pig Olympiad in Informatics, POI)，取得了不错的名次，最终保送进入了猪王国大学(Pig Kingdom University, PKU)深造。
现在的小猪已经能用计算机解决简单的问题了，比如能用P++语言编写程序计算出A + B的值。
这个“成就”已经成为了他津津乐道的话题。
当然，不明真相的同学们也开始对他刮目相看啦~
小猪的故事就将从此展开，伴随大家两天时间，希望大家能够喜欢小猪。

猪王国的文明源远流长，博大精深。

iPig在大肥猪学校图书馆中查阅资料，得知远古时期猪文文字总个数为 $N$ 。
当然，一种语言如果字数很多，字典也相应会很大。
当时的猪王国国王考虑到如果修一本字典，规模有可能远远超过康熙字典，花费的猪力、物力将难以估量。
故考虑再三没有进行这一项劳猪伤财之举。
当然，猪王国的文字后来随着历史变迁逐渐进行了简化，去掉了一些不常用的字。
iPig打算研究古时某个朝代的猪文文字。
根据相关文献记载，那个朝代流传的猪文文字恰好为远古时期的 $k$ 分之一，其中 $k$ 是 $N$ 的一个正约数（可以是 $1$ 和 $N$）。
不过具体是哪 $k$ 分之一，以及 $k$ 是多少，由于历史过于久远，已经无从考证了。
iPig觉得只要符合文献，每一种能整除 $N$ 的 $k$ 都是有可能的。
他打算考虑到所有可能的 $k$。
显然当k等于某个定值时，该朝的猪文文字个数为 $N / k$。
然而从 $N$ 个文字中保留下 $N / k$ 个的情况也是相当多的。
iPig预计，如果所有可能的 $k$ 的所有情况数加起来为 $P$ 的话，那么他研究古代文字的代价将会是 $G$ 的 $P$ 次方。
现在他想知道猪王国研究古代文字的代价是多少。
由于iPig觉得这个数字可能是天文数字，所以你只需要告诉他答案除以 $999911659$ 的余数就可以了。

---

Input

有且仅有一行：两个数 $N$ 、$G$，用一个空格分开。

---

Output

有且仅有一行：一个数，表示答案除以999911659的余数。

---

Sample Input

4 2

---

Sample Output

2048

---

Hint

$10\%$ 的数据中，$1 <= N <= 50$；
$20\%$的数据中，$1 <= N <= 1000$；
$40\%$的数据中，$1 <= N <= 100000$；
$100\%$的数据中，$1 <= G <= 1000000000，1 <= N <= 1000000000$ 。

---

Solution

首先我们可以把畸形的题目简略一下：

ans=G∑i|nCin mod 999911659

$$ans=G^{\sum_{i|n}C^i_n}\ mod\ 999911659$$
根据某个不记得名字的定理， $$a^k\ mod\ p=(a^{k\ mod\ \varphi(p)+\varphi(p)})\ mod\ p$$
则有，
$$ans=(G^{\sum_{i|n}G^i_n\ mod\ 999911658+999911658})\ mod\ 999911659$$
由 $999911658=2×3×4679×35617$，所以我们可以运用孙子定理，将 $\sum_{i|n}G^i_n$ 对 $2， 3， 4679， 35617$ 分别取模的结果算出，然后合并成对 $999911658$ 取模的结果。

在进行组合数取模的时候，可以用 $Lucas$ 定理，将组合数的上、下标降到 $35617$ 以下。
将 $n!$ 与 $n!$ 的逆元初始化，然后就可以快速求出组合数了。

---

Code

#include <iostream>
#include <cstdio>
#include <cmath>

#define MAXN 100000
#define LL unsigned long long
#define p (999911659)
#define rp (999911658)

using namespace std;

LL n,g,ans,up;
LL Mods[]={0,2,3,4679,35617};
LL ta[100010],tb[100010],tc[100010],td[100010];
LL Ma[100010],Mb[100010],Mc[100010],Md[100010];

LL get_MOD(LL a,LL b,LL c,LL d){
    return (((a*499955829%rp+b*333303886%rp)%rp+c*213702%rp*1353%rp)%rp+d*28074%rp*31254%rp)%rp;
}

LL C(LL x,LL y,LL MOD,LL tot){
    if(x<y)return 0;
    switch(tot){
        case 1:return ta[x]*Ma[x-y]%MOD*Ma[y]%MOD;
        case 2:return tb[x]*Mb[x-y]%MOD*Mb[y]%MOD;
        case 3:return tc[x]*Mc[x-y]%MOD*Mc[y]%MOD;
        case 4:return td[x]*Md[x-y]%MOD*Md[y]%MOD;
        default:printf("I'm Wrong!!!!\n");
    }
}

LL Lucas(LL x,LL y,LL MOD,LL tot){
    if(y==0)return 1;
    return C(x%MOD,y%MOD,MOD,tot)*Lucas(x/MOD,y/MOD,MOD,tot)%MOD;
}

LL power(LL x,LL y,LL MOD){
    if(y==0)return 1;
    if(y==1)return x%MOD;
    LL tmp=power(x,y/2,MOD);
    tmp=tmp*tmp%MOD;
    if(y&1)tmp=tmp*(x%MOD)%MOD;
    return tmp;
}

int main(){

    ta[0]=tb[0]=tc[0]=td[0]=1;

    for(LL i=1;i<=MAXN;i++){
        ta[i]=ta[i-1]*i%Mods[1];
        tb[i]=tb[i-1]*i%Mods[2];
        tc[i]=tc[i-1]*i%Mods[3];
        td[i]=td[i-1]*i%Mods[4];
    }

    Ma[0]=Mb[0]=Mc[0]=Md[0]=1;
    Ma[Mods[1]-1]=power(ta[Mods[1]-1],Mods[1]-2,Mods[1]);
    Mb[Mods[2]-1]=power(tb[Mods[2]-1],Mods[2]-2,Mods[2]);
    Mc[Mods[3]-1]=power(tc[Mods[3]-1],Mods[3]-2,Mods[3]);
    Md[Mods[4]-1]=power(td[Mods[4]-1],Mods[4]-2,Mods[4]);

    for(LL i=MAXN-1;i>=1;i--){
        if(i<Mods[1]-1)Ma[i]=Ma[i+1]*(i+1)%Mods[1];
        if(i<Mods[2]-1)Mb[i]=Mb[i+1]*(i+1)%Mods[2];
        if(i<Mods[3]-1)Mc[i]=Mc[i+1]*(i+1)%Mods[3];
        if(i<Mods[4]-1)Md[i]=Md[i+1]*(i+1)%Mods[4];
    }

    scanf("%llu%llu",&n,&g);
    LL sqr=sqrt(n),x;

    for(LL i=1;i<=sqr;i++){
        if(n%i==0){
            up=(up+get_MOD(Lucas(n,i,Mods[1],1),Lucas(n,i,Mods[2],2),Lucas(n,i,Mods[3],3),Lucas(n,i,Mods[4],4)))%rp;
            if(n/i!=i)up=(up+get_MOD(Lucas(n,n/i,Mods[1],1),Lucas(n,n/i,Mods[2],2),Lucas(n,n/i,Mods[3],3),Lucas(n,n/i,Mods[4],4)))%rp;
        }
    }

    printf("%llu\n",power(g,up+rp,p));

    return 0;
}