loj #3086. 「GXOI / GZOI2019」逼死强迫症

背景：

好像也没什么。

题目传送门：

https://loj.ac/problem/3086

题意：

给一个$2\ast n$的矩阵，现在要用$2\ast (n-1)$的矩形和$2$块$1\ast 1$的正方形铺满这个矩阵，要求这两个正方形不能相邻（边不能靠在一起）。

思路$1$：

显然要用$dp$。
设$f_{i,j}$表示在前$i$列已经用了$j$块$1\ast 1$的正方形的方案。
我们知道若没有这两块正方形，其方案数是$\text{Fibonacci}$数列。
用类似的思想。
$f_{i,0}=f_{i-1,0}+f_{i-2,0}$。
$f_{i,1}=f_{i-1,1}+2\ast f_{i-1,0}$。
$f_{i,2}=f_{i-2,1}+f_{i-1,2}+f_{i-2,2}$。
解释一下：
$[1]$：由于没有没有这两块正方形，其方案数是$\text{Fibonacci}$数列。
$[2]$：用了一块正方形，因此有两种情况：
$1.$在原来的列就用了，因此方案数为$f_{i-1,1}$。
$2.$在当前仅剩的一个这一列用，那么前面肯定$i-1$列没有用过，那么当前第$i$列可以放在上面或下面，即$2$种可能，因此方案数为$2\ast f_{i-1,0}$。
$[3]$：用了两块正方形，因此有两种情况：
$1.$在原来的列就用了，也就早就铺满了，剩下的肯定就是类似$\text{Fibonacci}$数列的东西了，因此方案数为$f_{i-1,2}+f_{i-2,2}$。
$2.$在当前第$i$列用，那么可以发现放置方法是固定的，不会产生新的贡献，直接继承$f_{i-2,1}$（不会出现两个正方形在相邻两列的情况）即可。（可以自己在草稿纸上画画。）
最后矩阵乘法优化即可。

思路$2$：

发现上面的式子用矩阵乘法真的不好做。
因此探求一种更优秀的做法。
$f_i=f_{i-1}+f_{i-2}+2\ast {\sum }_{j=1}^{i-3}fi{b}_j$
其中$fi{b}_j$表示$\text{Fibonacci}$数列第$j$项。
为什么可以这样呢？
对于思路$1$来说，由于我们只用求$f_{i,2}$，因此可以将$[1],[2]$带入$[3]$中，化简即可得到上面的式子。
那么$\sum$怎么消除呢？
$\text{Fibonacci}$数列有一个性质：${\sum }_{j=1}^{i-3}fi{b}_j=fi{b}_{i-1}-1$，因此有：
$$f_i=f_{i-1}+f_{i-2}+2\ast (fi{b}_{i-1}-1)$$

$$f_i=f_{i-1}+f_{i-2}+2fi{b}_{i-1}-2$$
好像好了许多。
矩阵乘法如下：
$$\left[\begin{array}{c}{f}_i\\ f_{i-1}\\ fi{b}_{i-1}\\ fi{b}_{i-2}\\ -2\end{array}\right]\ast \left[\begin{array}{ccccc}1& 1& 2& 2& 1\\ 1& 0& 0& 0& 0\\ 0& 0& 1& 1& 0\\ 0& 0& 1& 0& 0\\ 0& 0& 0& 0& 1\end{array}\right]=\left[\begin{array}{c}{f}_{i+1}\\ f_i\\ fi{b}_i\\ fi{b}_{i-1}\\ -2\end{array}\right]$$

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define mod 1000000007
#define LL long long
using namespace std;
struct node
{
	LL a[5][5];
	node clear()
	{
		memset(a,0,sizeof(a));
	}
	friend node operator * (const node a,const node b)
	{
		node c;
		c.clear();
		for(int i=0;i<5;i++)
			for(int j=0;j<5;j++)
				for(int k=0;k<5;k++)
					c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j]%mod)%mod;
		return c;
	}
} base,ans;
	int n;
LL work(int x)
{
	base.a[0][0]=1,base.a[0][1]=1,base.a[0][2]=2,base.a[0][3]=2,base.a[0][4]=1;
	base.a[1][0]=1,base.a[1][1]=0,base.a[1][2]=0,base.a[1][3]=0,base.a[1][4]=0;
	base.a[2][0]=0,base.a[2][1]=0,base.a[2][2]=1,base.a[2][3]=1,base.a[2][4]=0;
	base.a[3][0]=0,base.a[3][1]=0,base.a[3][2]=1,base.a[3][3]=0,base.a[3][4]=0;
	base.a[4][0]=0,base.a[4][1]=0,base.a[4][2]=0,base.a[4][3]=0,base.a[4][4]=1;
	ans.clear();
	for(int i=0;i<5;i++)
		ans.a[i][i]=1;
	while(x)
	{
		if(x&1) ans=ans*base;
		base=base*base;
		x>>=1;
	}
	return ((ans.a[1][2]+ans.a[1][3])%mod-ans.a[1][4]*2%mod+mod)%mod;
}
int main()
{
	int T,x;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&x);
		printf("%lld\n",work(x-1));
	}
}