51nod 计算 1556

题意：

1556 计算
1.0 秒 524,288.0 KB 80 分 5级题
有一个1*n的矩阵 固定第一个数为1 其他填正整数 且相邻数的差不能超过1 求方案数%1e9+7的结果

输入
一个数n 表示1*n的矩阵（n<=10^6）
输出
一个数 表示方案数%1e9+7的结果
输入样例
3
输出样例
5

思路：

（1）题意很短，很简洁，很清晰，很明了。不说了hhh

（2）这道题用到了默慈金数，那么什么是默慈金数呢？有一种说法就是从（0，0）开始走，每次只能向右走，横向，右上，右下都可以，但是不可以走到y = 0下面，问走到（n，0）的方案数。 方法的数量就是默慈金数。

（3）默慈金数的求解

（4）答案的计算：走的第一步只有两种方案，剩下的都有三种方案。。

ans[i] = ans[i - 1] * 3 - M[i - 1]

答案 = 总的方案数 - 不可行的方案数

代码实现：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1e6 + 5;
const ll mod = 1e9 + 7;

int n;
ll M[maxn];

//快速幂
ll qpow(ll a,ll b){
	ll res = 1;
	while(b){
		if(b & 1) res = (res * a) % mod;
		b >>= 1;
		a = (a * a) % mod;
	}
	return res % mod;
}

//求解默慈金数
void getM(){
	M[0] = 1; M[1] = 1; M[2] = 2;
	for(int i = 3;i <= n;i++){
		M[i]=(((((2*(i-1)+3)*M[i-1])%mod+((3*(i-1))*M[i-2])%mod)%mod)*qpow(i-1+3,mod-2))%mod;
	}
}

int main(){
	scanf("%d",&n);
	getM();
	ll ans = 2;
	for(int i = 2;i < n;i++){
		ans = (ans * 3 - M[i - 1] + mod) % mod;
	}
	printf("%lld\n",ans);
	return 0;
}