【51Nod】1013 3的幂的和

题意

求：3^0 + 3^1 +…+ 3^(N) mod 1000000007

解题思路

由等比数列的前n相和，我们可以得到
$\frac{3^{n+1}-1}{2} \% 1000000007$
我们知道$\frac{A}{B} \% C \neq \frac{A\%C}{B\%C}\%C$和$(A*B)\%C =(A\%C * B\%C)\%C$
所以我们要将除法转为乘法，引入乘法逆元。
乘法逆元：如果x与y的积除以z所得的余数为1，即xy = 1 (mod z)，则称x和y对于模数z来说互为逆元。
假定$(B*X)\%C=1$则
$\frac{A}{B} \% C=\frac{A}{B} * 1 \% C = \frac{A}{B}\%C * 1 \% C= \frac{A}{B}\%C *(B*X)\%C=(\frac{A}{B}*B*X)\%C=(A*X)\%C$
在本题中，2*500000004 = 1 （mod 1000000007），所以500000004是2对于1000000007的逆元。则
$\frac{3^{n+1}-1}{2} \% 1000000007 = (3^{n+1}-1)*500000004 \% 1000000007$
对于$a^b\%c$，采用快速幂运算即可。

参考代码

#include <iostream>
using namespace std;
typedef long long ll;
ll MOD = 1000000007;
ll mod_pow(ll n){
    ll k=3,ans=1;
    while (n){
        if (n&1) ans=(ans*k)%MOD;
        n/=2;
        k=(k*k)%MOD;
    }
    return ans;
}
int main(){
    ll n;
    while (cin>>n){
        ll ans=((mod_pow(n+1)-1)*500000004)%MOD;
        cout<<ans<<endl;
    }
    return 0;
}