HDU 5793 A Boring Question(快速幂&求逆元)

这道题就是给你一个公式,然后再给数据范围给你,然后求答案。

这道题的公式推到最后就是一个等比数列的求和公式,也就是说,输入n和m,然后从m的0次方一直加到n次方的和就是答案。

但是由于这道题的数据太大了,有1e9,所以我们不可能用for循环一直加到最后,所以根据公式,我们需要用到快速幂来求m^n。用快速幂求出来之后,我们根据等比数列的求和公式,就成了(1-m^n)/(1-m)%mod。因为还要对1e9+7取模,所以mod=1e9+7。但是因为要取模,所以直接当做除法来做会出错,所以,要把(1-m)和mod的逆元求出来,当做乘法来做。所以求逆元,其实就是求倒数。所以(1-m)和mod就是其倒数。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <vector>
#include <set>
#include <list>
#include <queue>
#include <map>
#include <stack>
#include <ctime>
using namespace std;
const int mod=1e9+7;
const int maxn = 1000010;

long long n,m;
long long X,Y;
long long mi(long long x,long long a,int p)
{
    long long ans = 1;
    x = x % p;
    while(a)
    {
        if(a&1)
            ans = (ans * x) % p;
        x = (x * x) % p;
        a >>= 1;
    }
    return ans;
}
long long gcd(long long a,long long b)
{
    int t,d;
    if(b == 0)
    {
        X = 1;
        Y = 0;
        return a;
    }
    d = gcd(b,a%b);
    t = X;
    X = Y;
    Y = t - (a / b) * Y;
    return d;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&m);
        gcd(m-1,mod);
        if(X < 0)
            X += mod;
        long long ans = (mi(m,n+1,mod)-1)*X%mod;
        printf("%lld\n",ans);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值