hdu5793 A Boring Question

最新推荐文章于 2020-01-26 23:23:47 发布

i是一只程序猿

最新推荐文章于 2020-01-26 23:23:47 发布

阅读量415

点赞数

分类专栏： ACM 文章标签： ACM HDU5793 HDU 快速幂取余求逆元

本文链接：https://blog.csdn.net/zh9406/article/details/52133193

版权

ACM 专栏收录该内容

92 篇文章 0 订阅

订阅专栏

这个题的传送地址：http://acm.hdu.edu.cn/showproblem.php?pid=5793

题意很简单，大家看过即明白。这个题是关于二项式定理和推理的应用，可能是因为数学是对规律的把握，很多人都是直接打表找规律了。

\sum 0 \leq k 1, k 2, \dots k m \leq n \prod 1 \leq j < m (k j + 1 k j) = \sum 0 \leq k 1 \leq k 2 \leq \dots \leq k m \leq n \prod 1 \leq j < m (k j + 1 k j) = \sum k m = 0 n \sum k m - 1 = 0 k m \dots \sum k 1 = 0 k 2 \prod 1 \leq j < m (k j + 1 k j) = \sum k m = 0 n ⎧ ⎩ ⎨ (k m k m - 1) \sum k m - 1 = 0 k m ⎧ ⎩ ⎨ (k m - 1 k m - 2) \dots \sum k 1 = 0 k 2 (k 2 k 1) ⎫ ⎭ ⎬ ⎫ ⎭ ⎬ = \sum k m = 0 n ⎧ ⎩ ⎨ \sum k m - 1 = 0 k m ⎧ ⎩ ⎨ (k m k m - 1) \dots \sum k 2 = 0 k 3 (k 3 k 2) 2 k 2 ⎫ ⎭ ⎬ ⎫ ⎭ ⎬ = \sum k m = 0 n m k m = m n + 1 - 1 m - 1

$\sum_{0\leq k_{1},k_{2},\cdots k_{m}\leq n}\prod_{1\leq j< m}\binom{k_{j+1}}{k_{j}} \\ =\sum_{0\leq k_{1}\leq k_{2}\leq\cdots \leq k_{m}\leq n}\prod_{1\leq j< m}\binom{k_{j+1}}{k_{j}} \\ =\sum_{k_{m}=0}^{n}\sum_{k_{m-1}=0}^{k_{m}}\cdots \sum_{k_{1}=0}^{k_{2}}\prod_{1\leq j< m}\binom{k_{j+1}}{k_{j}} \\ =\sum_{k_{m}=0}^{n}\left \{ \binom{k_{m}}{k_{m-1}} \sum_{k_{m-1}=0}^{k_{m}} \left \{ \binom{k_{m-1}}{k_{m-2}} \cdots \sum_{k_{1}=0}^{k_{2}}\binom{k_{2}}{k_{1}} \right \} \right \} \\ =\sum_{k_{m}=0}^{n}\left \{ \sum_{k_{m-1}=0}^{k_{m}} \left \{ \binom{k_{m}}{k_{m-1}} \cdots \sum_{k_{2}=0}^{k_{3}}\binom{k_{3}}{k_{2}}2^{k_{2}} \right \} \right \} \\ =\sum_{k_{m}=0}^{n}m^{k_{m}} \\ =\frac{m^{n+1} - 1}{m - 1} \\$
这是我得到最后答案的步骤，其中第三个等号到第四个等号的转换用到的是二项式定理

2 n = \sum 0 \leq k \leq n C k n

$2^n = \sum_{0\leq k\leq n}{C_n^k}$
第三个等号到第四个等号的转换用到二项式推论

(1 + x) n = \sum k = 0 n C k n x k

$(1+x)^n = \sum_{k=0}^{n} {C_n^k}x^k$
最后就得到最终的答案了。后面就是快速幂取余以及求逆元了。
以上是对此题的运算过程，如有问题欢迎指教。

#include<stdio.h>
#include<string.h>
#include<algorithm>

using namespace std;
typedef long long LL;
const LL mod = 1000000007;
LL pow(LL a,LL n)
{
    LL ans = 1;
    while(n)
    {
        if(n&1)
            ans = ans * a % mod;
        a = a * a % mod;
        n>>=1;
    }
    return ans;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        LL n,m;
        scanf("%lld%lld",&n,&m);
        LL ans =  pow(m,n+1);
        LL mod_ = pow(m-1,mod-2);
        printf("%lld\n",((ans-1)*mod_%mod+mod)%mod);
    }
    return 0;
}