51nod 1013 3的幂数和（快速幂+费马小定理）

题目链接：传送门
这题就是求公比为3首项为0的等比数列前n项和
公式是 (3^n-1)/2
上面直接用快速幂
上面会很大，这样除以2会有误差，所以除以2就要乘以1/2的逆元
2关于mod1000000007的逆元是 2^mod-2
顺便附大神对费马小定理的理解 传送门
代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
using namespace std;
#define N 1e9
#define mod 1000000007
long long int Pow(int a,long long int b){
    long long int r=1,base=a;
    while(b){
        if(b&1) r=(r*base)%mod;
        base=(base*base)%mod;
        b/=2;
    }
    return r;
}
int main(){
    int i,j,n,m;
    while(cin>>n){
        long long int res=Pow(3,n+1);
        long long int t=Pow(2,mod-2)%mod;
        cout<<(res-1)*t%mod<<endl;
    }
    return 0;
}

加一点今天学的一个推法：
去沈阳打比赛的时候C题是一道矩阵快速幂，是f(n)=f(n-1)+2f(n-2)+n^4
据说直接构造一个7维矩阵能直接构造出来，但不会，用二维的，就多出了个 i^4的求和
推了一整场硬是推不出来，心情就是rlg
然后赛后看了下，这种推都是从大一次方的开始推
比如 1^2+2^2+3^2+4^2+…+n^2
就有 (n+1)^3-n^3=3n^2+3n+1
n^3-(n-1)^3=-3n^2+3n-1
…
2^3-1^3=7
全部相加化简得
(n+1)^3-1^3=3(1^2+2^2+…+n^2)-1-n^2-n(n+1)/2
所以 1^2+2^2+..+n^2=(n/2)(n+1)(2n+1)
同理也可以用来证明4次方的用5次方来证
就必须知道 1次方，2次方 ，3次方的求和公式
最后得出 sum(i^4)(i=1….) =n(n+1)(2n+1)(3n^2+3n-1)/30
感觉知道了推法推出来也要一整场ORZ，感觉有些队伍也是用这个推出来的ORZ
还是老老实实写几道矩阵快速幂的题吧..