bzoj4417 [Shoi2013]超级跳马 矩阵快速幂

Description

---

现有一个n行m列的棋盘，一只马欲从棋盘的左上角跳到右下角。每一步它向右跳奇数列，且跳到本行或相邻行。跳越期间，马不能离开棋盘。例如，当n = 3, m = 10时，下图是一种可行的跳法。

试求跳法种数mod 30011。

对于100%的数据，1 ≤ n ≤ 50，2 ≤ m ≤ 10^9

Solution

---

看起来就像是矩阵快速幂的题。。

计数嘛，f[i,j]表示走到（i，j）的方案数，然后写写转移可以发现本质是两个前缀和，用s1[]和s2[]表示
注意到我们转移是交替从s1转移到s2，再由s2转移到s1，一个技巧就是把s1和s2拼在一起，分别求出两个转移矩阵A和B，那么我们一次转移两步就行了，也就是令C=A*B。

我们最后求出的实际上是前缀和，因此要做两次

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

const int MOD=30011;

int n,m;

struct Mat {
	int r[102][102];

	Mat() {fill(r,0);}

	int* operator [](int x) {return r[x];}

	Mat operator *(Mat B) const {
		Mat A=*this,C;
		rep(i,1,n*2) rep(j,1,n*2) {
			rep(k,1,n*2) C[i][j]=(C[i][j]+1LL*A[i][k]*B[k][j]%MOD)%MOD;
		}
		return C;
	}

	Mat operator ^(int x) const {
		Mat A=*this,C=A; x--;
		for (;x;x>>=1) {
			if (x&1) C=C*A;
			A=A*A;
		}
		return C;
	}
} ;

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

int solve() {
	Mat A,B,C;
	rep(i,1,n+n) A[i][i]=1;
	rep(i,n+1,n+n) {
		if (i!=n+1) A[i-n-1][i]=1;
		A[i-n][i]=1;
		if (i!=n+n) A[i-n+1][i]=1;
	}
	rep(i,1,n+n) B[i][i]=1;
	rep(i,1,n) {
		if (i!=n) B[n+i+1][i]=1;
		B[n+i][i]=1;
		if (i!=1) B[n+i-1][i]=1;
	}
	C=A*B; m--;
	C=C^(m/2);
	if (m&1) {
		C=C*(A*B);
		return C[n+n][1];
	} else return C[n][1];
	// m++;
}

int main(void) {
	scanf("%d%d",&n,&m);
	int ans=solve(); m--;
	ans-=solve(); ans=(ans%MOD+MOD)%MOD;
	printf("%d\n", ans);
	return 0;
}