51Nod 1013 3的幂的和 快速幂+逆元

24 篇文章 0 订阅

求:3^0 + 3^1 +...+ 3^(N) mod 1000000007

收起

输入

输入一个数N(0 <= N <= 10^9)

输出

输出:计算结果

输入样例

3

输出样例

40

 等比数列求和公式:

Sum=a1*(1-q^n)/(1-q),所以可以转化为:

Sum=(3^(n+1)-1)/2;

因为数太大,需要利用逆元来代替/2运算。

因为2与1000000007互质,所以可以利用公式2^(Mod-2)求出逆元。

所以最终公式为:

Sum=(3^(n+1)-1)*(2^(Mod-2))%Mod;

套用快速幂求出答案。

代码如下:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
typedef long long ll;
const ll Mod=1000000007;
ll n;
ll Fast (ll a,ll b)
{
    ll sum=1;
    while (b)
    {
        if(b&1)
        {
            sum=sum*a%Mod;
        }
        a=a*a%Mod;
        b/=2;
    }
    return sum;
}
int main()
{
    scanf("%lld",&n);
    ll re=(Fast(3,n+1)-1)*Fast(2,Mod-2)%Mod;
    printf("%lld\n",re);
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值