Bzoj1297:[SCOI2009]迷路:矩阵快速幂

题目链接1297:[SCOI2009]迷路

看完题就想到了沼泽鳄鱼那道题

但是这道题通过每条边的花费并不是1，所以单纯地邻接矩阵T次方并不是答案

但是我们可以拆点，每个点拆成9个，标号为i1--i9

从ij->ij-1连边

然后对于一条边(x,y,z)表示x到y花费z的时间，我们连边x->yz

这样x->y就会被拆成z步，然后就是邻接矩阵快速幂了

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define id(x,y) ((y-1)*n+x)
#define LL long long
using namespace std;
const int mod=2009;
const int maxid=101;
const int N=101;
int a[N][N],b[N][N],n;
LL T;
char str[N];

void Matrix_mul(int a[N][N],int b[N][N],int c[N][N]){
	int tmp[N][N]={0};
	for (int i=1;i<=n*9;++i)
	    for (int j=1;j<=n*9;++j)
	        for (int k=1;k<=n*9;++k)
	            tmp[i][j]=(tmp[i][j]+a[i][k]*b[k][j])%mod;
	for (int i=1;i<=n*9;++i)
	    for (int j=1;j<=n*9;++j)
	        c[i][j]=tmp[i][j];
}

void Matrix_pow(LL y){
	while (y){
		if (y&1) Matrix_mul(a,b,a);
		Matrix_mul(b,b,b); y>>=1;
	}
}

int main(){
	scanf("%d%lld",&n,&T);
	for (int i=1;i<=n;++i){
		scanf("%s",str+1);
		for (int j=1;j<=n;++j){
			int x=str[j]-'0';
			if (!x) continue;
			b[i][id(j,x)]=1;
		}
	}
	for (int i=1;i<=n;++i)
	    for (int j=2;j<=9;++j)
	        b[id(i,j)][id(i,j-1)]=1;
	for (int i=1;i<=n*9;++i) a[i][i]=1;
	Matrix_pow(T);
	printf("%d\n",a[1][n]);
}