题意：求 $A^B$ 的所有约数之和 $\bmod \ 9901$ 。数据范围： $0\le A,B\le 5\times 10^7$题解：由算术基本定理：$$A=\prod_{i=1}^k p_i^{a_i}$$所以：$$A^B=\prod_{i=1}^k p_i^{Ba_i}$$其约数和为：$$\prod_{i=1}^k \sum_{j=0}^{Ba_i} p_i^j$$其中 $\sum...

前言明日方舟的游戏内容实在太少了，闲的无聊就顺手把立绘提取了出来预先准备AssetStudioGalPhotoAuto在文末提供下载步骤用AssetStudio提取出立绘，文件位于Android/data/com.hypergryph.arknights/files/AB/Android下的charpack与avg中，若不全请解压安装包并从相同位置提取因为小人与立绘打包在一起，需要先筛选，按...

自适应辛普森法一、定积分定积分是积分的一种，是函数$f(x)$在区间$[a,b]$上积分和的极限定积分的几何意义是函数的曲线上$x$的一段区间与$x$轴围成的曲边梯形的带符号面积表示为$$\int_a^b f(x)\mathrm{d}x$$由定积分的定义，若函数$f(x)$在区间$[a,b]$上可积分，则有$n$等分的特殊分法$$\lim_{n \rightarrow +\infty} \s...

素数判断一、素数的定义质数（prime number）又称素数，有无限个。质数定义为在大于1的自然数中，除了1和它本身以外不再有其他因数。注意 $1$ 不是素数也不是合数。二、判断一个数是否为素数由定义可得，一个数 $n$ 为素数当且仅当它只含因数 $1$ 和它本身。于是可从 $2$ 枚举到 $n-1$ 判断其是否为 $n$ 的因数。inline bool isPrime(int x){ ...

卢卡斯定理一、卢卡斯（Lucas）定理求$$C_n^m \ (\bmod \ p)$$1.定理$$ \begin{aligned} & a=a_kp^k+a_{k-1}p^{k-1}+\dots+a_1p+a_0 \\ & b=b_kp^k+b_{k-1}p^{k-1}+\dots+b_1p+b_0 \end{aligned} $$其中 $0 \le a_i,b_i \le ...

乘法逆元一、逆元的定义若 $a \cdot x \equiv 1 \ (\bmod \ p)$ 且 $(a,p)=1$ ，则称 $a$ 关于模 $p$ 的乘法逆元为 $x$ ，记为 $a^{-1}$ 。二、逆元的求法1.费马（Fermat）小定理若 $p$ 为素数，且 $(a,p)=1$ ，则有 $a^{p-1} \equiv 1 \ (\bmod \ p)$ 。由费马小定理 $$a \cd...

欧拉定理 & 费马小定理一、欧拉（Euler）定理1.定理若 $(a,m)=1$ ，则$$a^{\varphi(m)} \equiv 1 \ (\bmod \ m)$$ 其中 $\varphi(m)$ 为欧拉函数，定义为小于或等于 $m$ 的正整数中与 $m$ 互质的数的数目。2.证明：记模 $m$ 意义下的缩系（小于 $m$ 且与 $m$ 互素的正整数集合）$$R=\{x_1,x_2,\d...