[HDU1852] Beijing 2008 [唯一分解定理][除法取模]

[$\mathfrak{L}\mathfrak{i}\mathfrak{n}\mathfrak{k}$]

---

有多组测试数据。
每组给定$\mathrm{N},\mathrm{K}$，令$\mathrm{S}$为$2008^{\mathrm{N}}$的正因数和$\sigma (2008^{\mathrm{N}})$，$\mathrm{M}\equiv \mathrm{S}(\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{K})$。
求$2008^{\mathrm{M}}(\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{K})$。

$1\le \mathrm{N}\le 10^7$，$500\le \mathrm{K}\le 10^4$。

---

要求的是$\left[2008^{\sigma (2008^N)\%K}\right]\%K$。我们先设$k=\omega (2008^N)$。
数据范围很大，所以暴力应该是不行的，首先得对公式做一些简化。
提起因数，我就想到唯一分解定理
$$2008^N=\prod _{i=1}^k{p}_i^{a_i}\text{\hspace{1em}(1)}$$
那么，$\sigma (2008^N)$呢？
$$\sigma (2008^N)=\prod _{i=1}^k\sum _{j=1}^{a_i}{p}_i^j\text{\hspace{1em}(2)}$$
我们知道等比数列求和公式
$$\sum _{i=1}^n{a}_1{q}^{i-1}=a_1\cdot \frac{1-q^n}{1-q}\text{\hspace{1em}(3)}$$
由$(2)$、$(3)$，得到：
$$\sigma (2008^N)=\prod _{i=1}^k{p}_i\cdot \frac{1-p_i^{a_i}}{1-p_i}\text{\hspace{1em}(4)}$$
记$x$对模$K$的逆元为$x^{-1}$。那么：
$$\sigma (2008^N)\equiv \prod _{i=1}^k{p}_i(1-p_i^{a_i})(1-p_i{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}K)\text{\hspace{1em}(5)}$$
这是一种处理除法取模的方法。但是逆元有一个条件！要求逆元的那个应该和$K$互质！！


那我们就要思考一下除法取模还有什么办法了。
首先，很显然$(1-p_i)|(1-p_i^{a_i})$，那么，设$\eta =\sigma (2008^N)\mathrm{m}\mathrm{o}\mathrm{d}K$：
$$\eta =\left\{\prod _{i=1}^k(p_i\mathrm{m}\mathrm{o}\mathrm{d}K)\cdot \frac{(1-p_i^{a_i})\mathrm{m}\mathrm{o}\mathrm{d}[(1-p_i)\cdot K]}{(1-p_i)}\mathrm{m}\mathrm{o}\mathrm{d}K\right\}\mathrm{m}\mathrm{o}\mathrm{d}K\text{\hspace{1em}(6)}$$
这个应该还是蛮容易推的？
（推$\forall a,b\in {\mathbb{N}}^{\ast }$，$a|b$时$\frac{a}{b}\mathrm{m}\mathrm{o}\mathrm{d}c=\frac{a\mathrm{m}\mathrm{o}\mathrm{d}bc}{b}\mathrm{m}\mathrm{o}\mathrm{d}c$就可以了。也挺显然的。）
_{（当然你也可以上下分解质因数，不过我觉得挺麻烦的，万一爆常数了呢。）}


现在，我们只需要利用$(6)$求出$\eta$。

但是怎么对$2008^N$分解质因数？——设$2008=\prod _{i=1}^m{p}_i^{t_i}$，
$$2008^N=\prod _{i=1}^m{p}_i^{t_{i}N}\text{\hspace{1em}(7)}$$
利用$(7)$，又因为对$2008$做质因数分解不难，我们可以先唯一分解$2008$。
$$2008=251\ast 2^3\text{\hspace{1em}(8)}$$
考虑到并不用真的计算出$2008^N$：我们只需要分别记录$p_i$和$t_{i}N$就可以了。
$$2008^N=2^{3N}501^N\text{\hspace{1em}(9)}$$

然后求$2008^{\eta }(\mathrm{m}\mathrm{o}\mathrm{d}K)$。
这样就可以解决这道题目了。

---

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<cctype>
#include<algorithm>
using namespace std;
long long defpowm(long long a, long long b, long long m) {
    long long ans = 1;
    while(b) {
        if (b&1) {
            ans*=a;
            ans%=m;
        }
        a*=a;
        a%=m;
        b>>=1;
    }
    return ans;
}
int main() {
    long long n, k;
    while(~scanf("%lld%lld", &n, &k)) {
        if (n==0&&k==0) break;
        long long temp0 = k*250ll;
        long long temp1 = defpowm(251, n+1, temp0) - 1;
        long long temp2 = defpowm(2, 3*n+1, temp0) - 1;
        long long ans = temp1*temp2/250; ans %= k;
        printf("%lld\n", defpowm(2008, ans, k));
    }
    return 0;
}