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

题目就是求\(G^{\sum_{d|n}C^d_n}\%p\),p为质数。那么考虑使用欧拉定理对\(C^d_n\)取模,但是φ(p)即p-1就不是质数了,得用exlucas。
考虑使用lucas的方法,p-1是个有特殊性质的质数(这不看题解很难想出来..)它可以被分解为2,3,4679,35617几个质数,并且都只有一次幂。
考虑分别对每个质数计算组合数,那么就构成了一个线性同余方程组,并且满足模数互质,于是crt合并出来一个答案就行了(注意crt时的模数是φ(p)即p-1)。
还有一个很难察觉的坑点就是,如果G为p的整数倍的话答案是0,但是因为常用的快速幂模板,所以mod出来会很奇怪会是1.所以必须特判一下(但是不确定到底是不是这个原因,望在评论中指出)

#include <bits/stdc++.h>
using namespace std;

#define Mod 999911659
#define ll long long

ll n;
ll G;
ll fac[35617+10], inv[35617+10];
ll mod[4] = {2, 3, 4679, 35617}, ans[4];

ll power(ll a, ll b, ll p) {
    ll ans = 1ll;
    while(b) {
        if(b & 1) ans = ans * a % p;
        a = a * a % p; b >>= 1; 
    } return ans % p;
}
//C(n, m) % p 
ll C(ll n, ll m, ll p) {
    if(m < n) return 0;
    return fac[m] * inv[n] % p * inv[m - n] % p; 
}

ll lucas(ll n, ll m, ll p) {
    if(m < n) return 0; if(!m) return 1;
    return lucas(n / p, m / p, p) * C(n % p, m % p, p) % p; 
}

ll md[4];
ll crt() {
    for(int i = 0; i < 4; ++i) md[i] = (Mod - 1) / mod[i];
    ll Ans = 0;
    for(int i = 0; i < 4; ++i) {
        Ans += power(md[i], mod[i] - 2, mod[i]) * ans[i] % (Mod - 1) * md[i] % (Mod - 1);
        Ans %= (Mod - 1);
    }
    Ans += (Mod - 1); Ans %= (Mod - 1);
    return power(G, Ans, Mod);
}

void init(ll p) { fac[0] = 1;
    for(ll i = 1; i <= p; ++i) fac[i] = fac[i - 1] * i % p;
    for(ll i = 0; i <= p; ++i) inv[i] = power(fac[i], p - 2, p);
}

int main() {
    scanf("%lld%lld", &n, &G);
    if(G % Mod == 0) return puts("0"), 0;
    for(int j = 0; j < 4; ++j) {
        init(mod[j]);
        for(ll i = 1; i * i <= n; ++i) {
            if(n % i == 0) {
                ans[j] = (ans[j] + lucas(i, n, mod[j]) % mod[j]) % mod[j];
                if((n / i) != i) ans[j] = (ans[j] + lucas(n / i, n, mod[j]) % mod[j]) % mod[j];
            }
        }
        ans[j] %= mod[j];
    }
    printf("%lld\n", crt());
}