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;
}

猜你喜欢

转载自blog.csdn.net/vocaloid01/article/details/80369883