51nod 1013【快速幂+逆元】

等比式子:

Sn=(a1-an*q)/(1-q)

n很大,搞一发快速幂,除法不适用于取膜,逆元一下(利用费马小定理) 假如p是质数,且gcd(a,p)=1,那么 a^(p-1)≡1(mod p)。刚好在本道题目一样适用,mod=1e9+7就是质数,那么gcd也就是=1,OK,那么b*k=1 这个逆元就等于 a^(mod-2);

#include <cstdio>
#include <stack>
#include <iostream>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long LL;

const LL mod=1e9+7;

LL cal(LL x,LL g)
{
    LL ans=1;
    while(g)
    {
        if(g&1)
            ans=(ans*x)%mod;
        x=(x*x)%mod;
        g>>=1;
    }
    return ans%mod;
}

LL solve(LL n)
{
    LL ans;
    ans=(cal(3,n)*3%mod-1)%mod;
    ans=(ans*cal(2,mod-2))%mod;
    return ans;
}

int main()
{
    LL n;
    scanf("%lld",&n);
    if(!n)
        puts("1");
    else
        printf("%lld\n",solve(n));
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值