古代猪文--(全是数论QwQ）

题目链接

题意：

Sol：

首先当$q$为$999911659$时，答案为$0$.
否则的话，因为$999911659$是质数，根据欧拉定理的推论有：
$$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=2\times 3\times 4679\times 35617$四个质数的指数都是$1$.
那么我们枚举 $n$的约数$d$，利用$lucas$定理求组合数。
$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)$$

分别求出${\sum }_{d|n}{C}_n^d$对$2,3,4679,35617$四个质数取模的结果，记作：$a_1,a_2,a_3,a_4$。
计算过程中，可以先预处理$p$以内的阶乘，然后利用$lucas$快速求出组合数。
最后可以用中国剩余定理求解线性同余方程组：
$$\{\begin{array}{l}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{p}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{p}_2)\\ x\equiv a_3(\mathrm{m}\mathrm{o}\mathrm{d}{p}_3)\\ x\equiv a_4(\mathrm{m}\mathrm{o}\mathrm{d}{p}_4)\end{array}$$
就可以得到${\sum }_{d|n}{C}_n^d\mathrm{m}\mathrm{o}\mathrm{d}999911658$的最小非负正整数解$x$，然后利用快速幂求解$q^x\mathrm{m}\mathrm{o}\mathrm{d}999911659$就得到了问题的答案

Code:

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int mod = 999911658;

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

int inv(int x, int k){
    return qpow(x, k - 2, k);
}

int a[5], b[5] = {0, 2, 3, 4679, 35617};
const int N = 40010;
int jc[N];

void init(int k)
{
    jc[0] = 1;
    for(int i = 1; i <= k; ++i)
        jc[i] = (jc[i - 1] % k * i) % k;
}

int C(int n, int m, int p)
{
    if(n < m) return 0;
    return jc[n] % p * inv(jc[m], p) % p * inv(jc[n - m], p) % p;
}

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

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

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

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

最后：
我直呼好家伙，这几天少学一个（逆元，扩展欧几里得，欧拉定理，中国剩余定理，组合数，lucas ）就不会写这个了。我超 ！