[usOJ3384]幂和

题目

传送门 to usOJ

题目描述
对于给定的 $A,B,C,D,P$ ，求 ${\sum }_{i=A}^B{\sum }_{j=C}^D{j}^i$ 对 $P$ 取模的值。

保证 $P$ 是一个质数。

输入格式
输入一行五个整数 $A,B,C,D$ 和 $P$ 。

输出格式
输出一行，表示你的答案。

数据范围与约定
对于 $15\%$ 的数据，$A,B,C,D\le 1000$ 。

对于 $30\%$ 的数据，$C,D\le 500000$ 。

对于另外 $10\%$ 的数据，$A=B$ 。

对于 $50\%$ 的数据，$C,D\le 10^7$ 。

对于 $70\%$ 的数据，$P\le 10^7$ 。

对于 $100\%$ 的数据，$1\le C\le D\le 10^{18},1\le A\le B\le 3000,B+1＜P\le 10^9$ 。

思路

传说中的毒瘤数学题！

这道题显然是个数学结论题。考虑怎么处理 $S_k={\sum }_{i=0}^n{i}^k$ 。这样搞——

警告：下面是该死的真·数学变换

$$\begin{gathered} \because S_{k+1}=\sum _{i=0}^n{i}^{k+1} \\ \therefore S_{k+1}+(n+1{)}^{k+1}=\sum _{i=0}^n(i+1{)}^{k+1} \end{gathered}$$

二项式展开，得到

$$\begin{array}{rl}{S}_{k+1}+(n+1{)}^{k+1}& =\sum _{i=0}^n\sum _{j=0}^{k+1}{C}_{k+1}^j{i}^j\\ & =\sum _{j=0}^{k+1}{C}_{k+1}^j\sum _{i=0}^n{i}^j\end{array}$$

将 $j=k+1$ 的那一项提出来，得到

$$\begin{array}{rl}{S}_{k+1}+(n+1{)}^{k+1}& =C_{k+1}^{k+1}\sum _{i=0}^n{i}^{k+1}+\sum _{j=0}^k{C}_{k+1}^j\sum _{i=0}^n{i}^j\\ & =S_{k+1}+\sum _{j=0}^k{C}_{k+1}^j\sum _{i=0}^n{i}^j\end{array}$$

于是可以消除掉 $S_{k+1}$ ，得到

$$(n+1{)}^{k+1}=\sum _{j=0}^k{C}_{k+1}^j\sum _{i=0}^n{i}^j$$

于是我们把 $j=k$ 的那一项又拿出来，得到

$$\begin{array}{rl}(n+1{)}^{k+1}& =C_{k+1}^k\sum _{i=0}^n{i}^k+\sum _{j=0}^{k-1}{C}_{k+1}^j\sum _{i=0}^n{i}^j\\ & =(k+1)S_k+\sum _{j=0}^{k-1}{C}_{k+1}^j\sum _{i=0}^n{i}^j\end{array}$$

$$\therefore (k+1)S_k=(n+1{)}^{k+1}-\sum _{j=0}^{k-1}{C}_{k+1}^j\sum _{i=0}^n{i}^j$$

$$\therefore S_k=(k+1{)}^{-1}[(n+1{)}^{k+1}-\sum _{j=0}^{k-1}{C}_{k+1}^j{S}_j]$$

警报解除！

然后这就是一个 $\mathcal{O}(k^2)$ 的递推了。剩下的事情，就是前缀和做差。

代码

代码实现是比较简单的，这里就不给出了。