BZOJ1951 [Sdoi2010]古代猪文【费马小定理 + Lucas定理 + 中国剩余定理 + 逆元递推 + 扩展欧几里得】

题目

“在那山的那边海的那边有一群小肥猪。他们活泼又聪明，他们调皮又灵敏。他们自由自在生活在那绿色的大草坪，他们善良勇敢相互都关心……” ——选自猪王国民歌很久很久以前，在山的那边海的那边的某片风水宝地曾经存在过一个猪王国。猪王国地理位置偏僻，实施的是适应当时社会的自给自足的庄园经济，很少与外界联系，商贸活动就更少了。因此也很少有其他动物知道这样一个王国。猪王国虽然不大，但是土地肥沃，屋舍俨然。如果一定要拿什么与之相比的话，那就只能是东晋陶渊明笔下的大家想象中的桃花源了。猪王勤政爱民，猪民安居乐业，邻里和睦相处，国家秩序井然，经济欣欣向荣，社会和谐稳定。和谐的社会带给猪民们对工作火红的热情和对未来的粉色的憧憬。小猪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的余数就可以了。

输入格式

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

输出格式

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

输入样例

4 2

输出样例

2048

提示

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

题解

一句话题意：

ans=G∑k|NC(N,k)mod999911659

真是神级题目，模数设计得十分巧妙，综合了很多知识，需要处理的细节也很多

由于999911659是一个质数，利用费马小定理
ans=G∑k|NC(N,k)mod99911658mod999911659

O(N−−√)枚举出N的所有因子，逐一算出其组合数相加
N很大，计算组合数便十分棘手，不能直接算，考虑Lucas定理

但是Lucas定理要求模为素数，999911658不满足要求，我们将其因式分解，得到2，3，4679，35617，有很好的性质就是每个质因子都不大而且指数都为1

我们就可以分开算模每个因子的情况下的结果，最后用中国剩余定理总和到一起算出最后的结果

Lucas算组合数时可以先预处理各个模意义下的阶乘和逆元阶乘

最后要注意当G可能等于P，这时就是0啦。
如果直接算可能会出一个数据使得指数为0，这样快速幂就会算出1来

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
using namespace std;
const int P = 999911659,PP = 999911658;
LL N,G,ans[5];
LL prime[4]={2,3,4679,35617},f[4][35618],inv[4][35618];
void init(){
    for (LL i = 0; i < 4; i++){
        LL p = prime[i];
        f[i][0] = 1;
        for (LL j = 1; j < p; j++) f[i][j] = f[i][j - 1] * j % p;
        inv[i][1] = 1;
        for (LL j = 2; j < p; j++) inv[i][j] = (p - p / j) * inv[i][p % j] % p;
        inv[i][0] = 1;
        for (LL j = 1; j <= p; j++) inv[i][j] = inv[i][j] * inv[i][j - 1] % p;
    }
}
LL qpow(LL a,LL b,LL Mod){
    LL ans = 1;
    for (; b; b >>= 1,a = (LL)a * a % Mod)
        if (b & 1) ans = (LL)ans * a % Mod;
    return ans;
}
void exgcd(LL a,LL b,LL& d,LL& x,LL& y){
    if (!b) {x = 1; y = 0; d = a;}
    else exgcd(b,a % b,d,y,x),y -= (a / b) * x;
}
LL Lucas(LL n,LL m,LL i){
    if (n < m) return 0;
    if (n < prime[i] && m < prime[i]) return f[i][n] * inv[i][m] % prime[i] * inv[i][n - m] % prime[i];
    return Lucas(n % prime[i],m % prime[i],i) * Lucas(n / prime[i],m / prime[i],i) % prime[i];
}
void cal(LL x){
    for (int i = 0; i < 4; i++)
        ans[i] = (ans[i] + Lucas(N,x,i)) % prime[i];
}
LL China(){
    LL a,b,d,x,y,Ans = 0;
    for (LL i = 0; i < 4; i++){
        exgcd(a = PP / prime[i],b = prime[i],d,x,y);
        x = (x % PP * (PP / prime[i]) % PP + PP) % PP;
        Ans = (Ans + ans[i] * x % PP) % PP;
    }
    return (Ans % PP + PP) % PP;
}
int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out1.txt","w",stdout);
    init();
    cin>>N>>G; LL i;
    if (G == P) {puts("0"); return 0;}
    for (i = 1; i * i < N; i++)
        if (N % i == 0) cal(i),cal(N / i);
    if (i * i == N) cal(i);
    ans[5] = qpow(G % P,China(),P);
    cout<<ans[5]<<endl;
    return 0;
}