【bzoj1297】[SCOI2009]迷路

1297: [SCOI2009]迷路

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 1423   Solved: 1021
[ Submit][ Status][ Discuss]

Description

windy在有向图中迷路了。 该有向图有 N 个节点，windy从节点 0 出发，他必须恰好在 T 时刻到达节点 N-1。 现在给出该有向图，你能告诉windy总共有多少种不同的路径吗？ 注意：windy不能在某个节点逗留，且通过某有向边的时间严格为给定的时间。

Input

第一行包含两个整数，N T。 接下来有 N 行，每行一个长度为 N 的字符串。 第i行第j列为'0'表示从节点i到节点j没有边。 为'1'到'9'表示从节点i到节点j需要耗费的时间。

Output

包含一个整数，可能的路径数，这个数可能很大，只需输出这个数除以2009的余数。

Sample Input

【输入样例一】
2 2
11
00

【输入样例二】
5 30
12045
07105
47805
12024
12345

Sample Output

【输出样例一】
1

【样例解释一】
0->0->1

【输出样例二】
852

HINT

30%的数据，满足 2 <= N <= 5 ； 1 <= T <= 30 。 100%的数据，满足 2 <= N <= 10 ； 1 <= T <= 1000000000 。

Source

[ Submit][ Status][ Discuss]



矩乘 + DP

如果边权为1的话非常好稿，但是这里边权并不为1

继续发现边权的值很小，于是我们可以把一个点拆成10个点，使他们之间的权值都为1

例如，现在要从1向2连一条权值为3的边，那么就(1 , 0) --- > (1 , 1) --- > (1 , 2) --- > (2 , 0)

然后跑矩乘就好了

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;

typedef long long LL;

const int maxtot = 110;
const int maxn = 15;
const int crz = 2009;

struct matrix{
	int r,c,m[maxtot][maxtot];
	matrix(){r = c = 0; memset(m,0,sizeof(m));}
	matrix operator * (matrix b) const
	{
		matrix ret;
		ret.r = r; ret.c = b.c;
		for (int i = 1; i <= r; i++)
			for (int j = 1; j <= c; j++)
				for (int k = 1; k <= b.c; k++)
					ret.m[i][k] = (ret.m[i][k] + m[i][j] * b.m[j][k]) % crz;
		return ret;
	}
}o,frd;

int id[maxn][maxn],n,t,tot;

inline int getint()
{
	int ret = 0,f = 1;
	char c = getchar();
	while (c < '0' || c > '9')
	{
		if (c == '-') f = -1;
		c = getchar();
	}
	while (c >= '0' && c <= '9')
		ret = ret * 10 + c - '0',c = getchar();
	return ret * f;
}

inline int getcom()
{
	char c = getchar();
	while (c < '0' || c > '9') c = getchar();
	return c - '0';
}

inline matrix qpow(matrix a,int b)
{
	matrix t = a,ret = o;
	while (b)
	{
		if (b & 1) ret = ret * t;
		t = t * t;
		b >>= 1;
	}
	return ret;
}

int main()
{
	n = getint(); t = getint();
	for (int i = 1; i <= n; i++)
		for (int j = 0; j <= 8; j++)
			id[i][j] = ++tot;
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= 8; j++)
			frd.m[id[i][j]][id[i][j - 1]] = 1;
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= n; j++)
		{
			int x = getcom();
			if (x == 0) continue;
			frd.m[id[j][0]][id[i][x - 1]] = 1;
		}
	o.r = o.c = frd.r = frd.c = tot;
	for (int i = 1; i <= tot; i++) o.m[i][i] = 1;
	
	matrix ans;
	ans.r = tot; ans.c = 1; ans.m[id[1][0]][1] = 1;
	ans = qpow(frd,t) * ans;
	printf("%d",ans.m[id[n][0]][1]);
	return 0;
}