(HDU) 6185 Covering (矩阵快速幂)

传送门

题意：给你一个4*n的矩阵，然后用1*2和2*1的砖去把他铺满，求方案数。

解题思路：这个无疑是推出一个递推式，但是由于n很大，所以必须要用矩阵快速幂，但是菜啊，不是多算就是少算了，这个大哥图文解说，还是自己的思维太差了啊。

#include<bits/stdc++.h>
#define il inline
#define pb push_back
#define fi first
#define se second
#define ms(_data,v) memset(_data,v,sizeof(_data))
#define sc(n) scanf("%d",&n)
#define SC(n,m) scanf("%d %d",&n,&m)
#define SZ(a) int((a).size())
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define drep(i,a,b)	for(int i=a;i>=b;--i)
using namespace std;
typedef long long ll;
const ll inf=0x3f3f3f3f;
const double PI=acos(-1.0);
const double eps=1e-9;
const int mod=1e9+7;
ll n;
struct Ma{
	ll a[5][5];
	Ma(){ms(a,0);}
	Ma operator *(const Ma &tp)const{
		Ma res;
		rep(i,1,4)	rep(j,1,4)	rep(k,1,4){
			res.a[i][j]=(res.a[i][j]+a[i][k]*tp.a[k][j]%mod)%mod;
		}
		return res;		
	}
};
il Ma mpow(Ma x,ll n){
	Ma res;
	rep(i,1,4)	res.a[i][i]=1;
	while(n){
		if(n&1)	res=res*x;
		x=x*x;
		n>>=1;
	}
	return res;
}
int main(){
	std::ios::sync_with_stdio(0);
	while(cin>>n){
		if(n==1)	cout<<1<<endl;
		else if(n==2)	cout<<5<<endl;
		else if(n==3)	cout<<11<<endl;
		else if(n==4)	cout<<36<<endl; 
		else{
			Ma x;
			x.a[1][1]=1,x.a[1][2]=5;
			x.a[1][3]=1,x.a[1][4]=-1;
			x.a[2][1]=1,x.a[3][2]=1;
			x.a[4][3]=1;
			x=mpow(x,n-4);
			ll ans=(x.a[1][1]*36%mod+x.a[1][2]*11%mod+mod)%mod;
			ans=(ans+x.a[1][3]*5%mod+mod)%mod;
			ans=(ans+x.a[1][4]%mod+mod)%mod;
			cout<<ans<<endl;;
		}
	}
	return 0;
}