GXOI/GZOI2019 逼死强迫症 题解

题目传送门

题目大意： 有一个 $2\times n$ 的网格图，有 $n$ 个 $1\times 2$ 的方块，其中一个被劈开变成了两个 $1\times 1$ 大小的方块，要求两个小方块不相邻，问有多少种不同的用这些方块放满网格图的方案。

题解

有一个结论，设小方块所在的列分别为 $x,y$，那么 $x$ ~ $y$ 这段区间放方块的方案唯一，考虑在同行和异行的情况。

在同一行时，中间至少隔着一个横着的 $1\times 2$ 的方块，即 $y-x\ge 3$ 且 $(y-x)\mathrm{m}\mathrm{o}\mathrm{d}2=1$，然后另外一行中 $x$ 到 $y$ 列只能横着放满 $1\times 2$ 的方块，即方案唯一。

不在同一行时，两行中的 $x$ 到 $y$ 列这一段，依然只能选择横着放满 $1\times 2$，所以类似的满足 $y-x\ge 2$ 且 $(y-x)\mathrm{m}\mathrm{o}\mathrm{d}2=0$，且方案也唯一。

也就是说，$y-x>1$，并且能够随意放 $1\times 2$ 的方块的区间只有 $[1,x-1]$ 和 $[y+1,n]$，设 $f_i$ 表示一段长度为 $i$ 的区间随意放的方案数，可以发现，放方块只有两种方法：竖着放一个，横着放两个（上下各一个），所以有 $f_i=f_{i-1}+f_{i-2}$，即斐波那契数列，不过这里的初值为 $f_0=f_1=1$。

设 $g_i$ 表示 $n=i$ 时的答案，那么第 $i$ 位的放法也是类似的，竖着放一个，横着放两个，又或者放一个 $1\times 1$ 的方块。

如果放一个 $1\times 1$ 的方块，那么考虑枚举另一个方块的位置，方案数就是 $2\times {\sum }_{j=0}^{i-3}{f}_j$，乘 $2$ 是因为可以交换上下两行得到新的方案。

根据斐波那契数列的性质，可以知道它的前缀和满足 $sum(n)=f_{n+2}-1$，所以得到 $g$ 的递推式：$g_n=g_{n-1}+g_{n-2}+2sum(n-3)=g_{n-1}+g_{n-2}+2f_{n-1}-2$。

然后可以表示成矩阵的形式，用矩阵快速幂加速即可：
$$\left[\begin{array}{ccccc}1& 1& 2& 0& -2\\ 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}{g}_{n-1}\\ g_{n-2}\\ f_{n-1}\\ f_{n-2}\\ 1\end{array}\right]=\left[\begin{array}{c}{g}_n\\ g_{n-1}\\ f_n\\ f_{n-1}\\ 1\end{array}\right]$$

代码如下：

#include <cstdio>
#include <cstring>
#define mod 1000000007

int T,n;
void add(int &x,int y){x=(x+y>=mod?x+y-mod:x+y);}
struct matrix{
	int h,l,a[6][6];matrix(){memset(a,0,sizeof(a));}
	matrix operator *(const matrix &b)
	{
		matrix c;c.h=h;c.l=b.l;
		for(int i=1;i<=c.h;i++)
		for(int j=1;j<=c.l;j++)
		for(int k=1;k<=l;k++)
		add(c.a[i][j],1ll*a[i][k]*b.a[k][j]%mod);
		return c;
	}
}A,B,C,D;

int main()
{
	scanf("%d",&T);A.h=A.l=5;B.h=5;B.l=1;
	A.a[1][1]=1;A.a[1][2]=1;A.a[1][3]=2;A.a[1][5]=mod-2;
	A.a[2][1]=1;A.a[3][3]=1;A.a[3][4]=1;A.a[4][3]=1;A.a[5][5]=1;
	B.a[3][1]=2;B.a[4][1]=1;B.a[5][1]=1;//从f[2]开始，所以a[3][1]=f[2]=2
	while(T--)
	{
		scanf("%d",&n);n-=2;C=A;D=B;
		while(n>0){if(n&1)D=C*D;n>>=1;C=C*C;}
		printf("%d\n",D.a[1][1]);
	}
}