bzoj2126 排斥反应 矩阵乘法

       发现p很小而q很大。。矩阵乘法这么明显的提示。

       不过转化还是很巧妙的。将所有的数变成p行q列的矩阵，然后就变成求从中取出一些点使得没有点相邻（包括同一行的最前面和最后面也算相邻）。

       首先可以得到某一列没有相邻的状态（当p=10时有123种），然后就可以用矩阵乘法加速得到第一行转移到最后一行的方案数辣！！然后得到ans。

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define mod 19921107
#define N 155
#define ll long long
using namespace std;

struct matrix{ int p[N][N]; }a,b; int cnt,n,m,bin[15],f[N];
matrix tms(matrix x,matrix y){
	matrix z; int i,j,k;
	for (i=1; i<=cnt; i++)
		for (j=1; j<=cnt; j++){
			z.p[i][j]=0;
			for (k=1; k<=cnt; k++) z.p[i][j]=((ll)x.p[i][k]*y.p[k][j]+z.p[i][j])%mod;
		}
	return z;
}
void dfs(int k,int now){
	if (k>n){
		if (n==1 || !((now&1) && (now&bin[n-1]))) f[++cnt]=now;
		return;
	}
	dfs(k+1,now); dfs(k+2,now|bin[k-1]);
}
int main(){
	scanf("%d%d",&n,&m); int i,j;
	bin[0]=1; for (i=1; i<=n; i++) bin[i]=bin[i-1]<<1;
	dfs(1,0); if (m==1){ printf("%d\n",cnt); return 0; }
	for (i=1; i<=cnt; i++)
		for (j=1; j<=cnt; j++) if (!(f[i]&f[j])) a.p[i][j]=1;
	memcpy(b.p,a.p,sizeof(a.p));
	for (i=m-2; i; i>>=1,a=tms(a,a)) if (i&1) b=tms(b,a);
	int ans=0;
	for (i=1; i<=cnt; i++)
		for (j=1; j<=cnt; j++) if (!(f[i]&f[j])) ans=(ans+b.p[i][j])%mod;
	printf("%d\n",ans);
	return 0;
}<em>
</em>

by lych

2016.3.13