P2480 [SDOI2010] 古代猪文

P2480 [SDOI2010] 古代猪文

题目大意：求G∑d∣n C(n, d) mod 999911659

一开始想出题目大意，求n中因子的C(n, d)的和，那么直接Lucas开求！！！

结果当然是秒挂qwq

那么我们想一下为什么会挂，Lucas的适用条件是模数比较小，那么本题就又有了新的解题思路了

如何将模数缩小呢？

欧拉函数推论！！！

G ∑d∣n C(n,d)mod 999911659=G (∑d∣n C(n,d)mod 999911658)mod 999911659

手打的比较烂，具体可以去看欧拉函数

既然我们都把模数变成了999911658，那后面该怎么继续缩小呢？

将999911658分解发现它可以变成2, 3, 4679, 35617四个质数

那么取模四个两两互质的数来求最终那个答案该怎么做呢？

中国剩余定理！！！

而中国剩余定理需要扩展欧几里得来求每项M的逆元，所以这道题真的是一道特别好的数论全家桶！！！

(因为洛谷的两道中国剩余定理的板子都会让你的结果超longlong，所以在这里我的乘积用了效率比较高的浮点运算快速乘mul(模数不超过浮点运算的范围就可以用)。当然这题直接乘也是可以过的)

那么接下来就是愉快的总代码了hhh

#include<bits/stdc++.h>

using namespace std;

#define endl "\n"
#define ll long long
#define int long long

const int N = 4e4 + 10, mod = 999911658;
int a[5], b[5] = {0, 2, 3, 4679, 35617};
int v[N], inv[N], n, q;

int ksm(int a, int b, int p) {
    int ans = 1;
    while(b) {
        if(b & 1) ans = ans * a % p;
        a = a * a % p;
        b >>= 1;
    }
    return ans;
}

void init(int p) {
    v[0] = inv[0] = 1;
    for(int i = 1; i <= p; i++) {
        v[i] = v[i - 1] * i % p;
        inv[i] = inv[i - 1] * ksm(i, p - 2, p) % p;
    }
}

int C(int a, int b, int p) {
    if(a < b) return 0;
    return v[a] % p * inv[b] % p * inv[a - b] % p;
}

int lucas(int a, int b, int p) {
    if(a < p && b < p) return C(a, b, p);
    return C(a % p, b % p, p) % p * lucas(a / p, b / p, p) % p;
}

int exgcd(int a, int b, int &x, int &y) {
    if(!b) {
        x = 1, y = 0;
        return a;
    }
    
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    
    return d;
}

int crt(int a[], int m[], int n) {
    int M = 1, ans = 0;
    for(int i = 1; i <= n; i++) M *= m[i];
    
    for(int i = 1; i <= n; i++) {
        int mi = M / m[i];
        int x, y;
        int d = exgcd(mi, m[i], x, y);
        
        x = (x % m[i] + m[i]) % m[i];
        ans = (ans + a[i] % M * x % M * mi % M);
    }
    
    return ans;
}

signed main() {
    cin >> n >> q;
    
    if(q % (mod + 1) == 0) {
        cout << 0 << endl;
        return 0;
    }
    
    for(int p = 1; p <= 4; p++) {
        init(b[p]);
        for(int i = 1; i * i <= n; i++) {
            if(n % i == 0) {
                a[p] = (a[p] + lucas(n, i, b[p])) % b[p];
                if(i * i != n) a[p] = (a[p] + lucas(n, n / i, b[p])) % b[p];
            }
        }
    }
    
    int x = crt(a, b, 4);
    cout << ksm(q, x, mod + 1);
    
    return 0;
}