51Nod - 1013 3的幂的和

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

Input

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

Output

输出:计算结果

Sample Input

3

Sample Output

40
题解:
等比数列前(n+1)项和,公比为3首项为1所以t=( 3^(n+1)-1 ) / 2 % 1000000007。
3^(n+1)用快速幂解决,而除法取模只需将除法变乘法即可,也就是用被除数去乘以除数的乘法逆元。
求乘法逆元用费马小定理:当模为素数,a的逆元为myPow(a,mod-2)。

代码:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;
#define mod 1000000007

long long myPow(long long t,long long N){

	long long mid = 1;

	while(N){
		if(N&1){
			mid *= t;
			mid %= mod;
		}
		t *= t;
		t %= mod;
		N >>= 1;
	}

	return mid;
}

int main(){

	long long N;
	cin>>N;
	
	long long t = myPow(3,N+1)-1;
	t = (t*myPow(2,mod-2))%mod;
	
	cout<<t<<endl;

	return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值