51nod1556 默慈金数+乘法逆元

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

using namespace std;
const int MAXN = 1000005;
const int MOD = 1e9+7;
typedef long long LL;

int n,m[MAXN],f[MAXN];

int fast_pow(int x,int y){
    int ret=1;
    for(;y;y>>=1){
        if(y&1) ret=(LL)ret*x%MOD;
        x=(LL)x*x%MOD;      
    }
    return ret;
}

int main(){
    scanf("%d",&n);m[1]=1;m[2]=2;f[1]=2;
    for(int i=3;i<n;i++) m[i]=(LL)((LL)(2*i+1)*m[i-1]%MOD+(LL)3*(i-1)*m[i-2]%MOD)%MOD*fast_pow(i+2,MOD-2)%MOD;
    for(int i=2;i<n;i++) f[i]=((LL)f[i-1]*3-m[i-1]+MOD)%MOD;
    printf("%d\n",f[n-1]);
    return 0;   
}

看下这个