P4159 [SCOI2009] 迷路

使用博客园食用更佳：https://www.cnblogs.com/sonnety-v0cali0d-kksk/p/17471440.html

题目链接

---

题目内容

[SCOI2009] 迷路

题目背景

windy 在有向图中迷路了。

题目描述

该有向图有 $n$ 个节点，节点从 $1$ 至 $n$ 编号，windy 从节点 $1$ 出发，他必须恰好在 $t$ 时刻到达节点 $n$。

现在给出该有向图，你能告诉 windy 总共有多少种不同的路径吗？

答案对 $2009$ 取模。

注意：windy 不能在某个节点逗留，且通过某有向边的时间严格为给定的时间。

输入格式

第一行包含两个整数，分别代表 $n$ 和 $t$。

第 $2$ 到第 $(n+1)$ 行，每行一个长度为 $n$ 的字符串，第 $(i+1)$ 行的第 $j$ 个字符 $c_{i,j}$ 是一个数字字符，若为 $0$，则代表节点 $i$ 到节点 $j$ 无边，否则代表节点 $i$ 到节点 $j$ 的边的长度为 $c_{i,j}$。

输出格式

输出一行一个整数代表答案对 $2009$ 取模的结果。

样例 #1

样例输入 #1

2 2
11
00

样例输出 #1

1

样例 #2

样例输入 #2

5 30
12045
07105
47805
12024
12345

样例输出 #2

852

提示

样例输入输出 1 解释

路径为 $1\to 1\to 2$。

数据规模与约定

• 对于 $30\%$ 的数据，保证 $n\le 5$，$t\le 30$。
• 对于 $100\%$ 的数据，保证 $2\le n\le 10$，$1\le t\le 10^9$。

---

前置知识：矩阵快速幂

题目链接

一个幂数k=$\sum 2^i$ $i$ $\in$ $Z$
所以快速幂板子就出来啦（）

#include<bits/stdc++.h>
using namespace std;

const int maxn=150;
const int m=1e9+7; 
typedef long long intx;

int n; 
intx k; 

struct matrix{
	intx num[maxn][maxn];
	matrix(){memset(num,0,sizeof(num));}
	void build(){for(int i=1;i<=n;++i)	num[i][i]=1;}	//单位矩阵 
};matrix a,ans;
matrix operator*(const matrix &x,const matrix &y){
	matrix z;
	for(int k=1;k<=n;++k){
		for(int i=1;i<=n;++i){
			for(int j=1;j<=n;++j){
				z.num[i][j]=(z.num[i][j]+x.num[i][k]*y.num[k][j]%m)%m;
			}
		}
	}
	return z;
}

void input(){
	scanf("%d%lld",&n,&k);
	for(int i=1;i<=n;++i){
		for(int j=1;j<=n;++j){
			scanf("%d",&a.num[i][j]);
		}
	}
}

int main(){
	input();
	ans.build();
	while(k){
		if(k&1)	ans=ans*a;
		a=a*a;
		k=k>>1;
	}
	for(int i=1;i<=n;++i){
		for(int j=1;j<=n;++j){
			printf("%d ",ans.num[i][j]);
		}
		printf("\n");
	}
	return 0;
}

---

私货：话题#虚拟歌姬也会迷路吗#

思路历程

1.我寻思这图里咋还有自环呢

对角线还有值这自己和自己连边啊。
干不会了。
然后就看这个专题名字叫“矩阵快速幂”
但是我还没有打过板子所以速速去打了一下（）

2.ok快乐的板子时光总是短暂的（）

我想了想觉得边权肯定是附加的属性，于是我先不管边权只建边或许可以的得到一些启发：
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-jk9WIO2T-1686387618437)(https://img2023.cnblogs.com/blog/3106747/202306/3106747-20230610153146268-1758510698.png)]

好丑的图（）
嗯所以和矩阵有什么关系呢（其实当我画出这么多边的时候我就知道没什么用了）
怎么还越扯越远了（）

3.额我们还是不看边权，但是不扯到图上去了。

如果，我是说如果，这是个邻接矩阵：
样例1：
1 1
0 0
问有几种到达方法。
那么就是说1要一个一个点到达2，第一行都要有数。
也就是说啊，i可以一个一个点的到达j，意味着$a_i{,}_2$ ———— $a_i{,}_j$要有数
那经过的每一个点又都有一个出度，出度就是第i行1的个数
是不是说从1起点能到达的点的出度相乘再除以2就是方案数呢？
那么样例2呢：
1 1 0 1 1
0 1 1 0 1
1 1 1 0 1
1 1 0 1 1
1 1 1 1 1
怎么感觉不太对啊……
话说我要是从点1到点5再回点1再到点5这算不算一种方案啊……

4.那我们现在加上边权吧

那该怎么算边权呢？
我们假设边权只有1
$f_t[i,j]$表示从i点到j点花费时间为t的方案数
就有了方程：
$f_t=\sum _{k=1}^n{f}_{t-1}[i,k]+f_1[k,j]$
那$f_t$也就是$f_1^t$。
那就出来了a。
$f_1$就是最初的矩阵啊。

5.回归本题

本题的区别在于权值是1-9啊
1-9，1-9……
我想了想，我他喵的$f_t$是$f_1^t$
那$f_9$不就是$f_1^9$吗！！！
那我把这个边多拆成几条边。
反正是求不小于t的方案数。
但是我如何建边才能保证没有多得方案数啊。
总不能搁节点1和节点2之间塞点吧
……
去看看题解（）
……
我们1个点建9个小点，只有第0小点可以跨越小点的集合。
那么我们建立第i小点到第i-1小点的w=1的单向边
而要连的两个点1到点2则将小点集合1中的小点0和集合2中的小点w-1相连。
如下图（只画出了点1到点2经过的路径）

那么我们就转换成了边权都为1的9n*9n大小的矩阵。

---

代码实现：

注：与以上解释不同，代码中a[i][j]表示j指向i的边，并且不再设置第0小点改为设置第1小点为唯一可以外接点。

Miku's Code

#include<bits/stdc++.h>
using namespace std;

const int maxn=12;
const int mod=2009;

char s[maxn];
int n,sn,t;

struct matrix{
	int num[maxn*9][maxn*9];
	void clear(){
		memset(num,0,sizeof(num));
	}
	 
};matrix a;
matrix operator *(matrix x,matrix y){
	matrix res;
	res.clear();
	for(int k=1;k<=n;++k){
		for(int i=1;i<=n;++i){
			for(int j=1;j<=n;++j){
				res.num[i][j]=(res.num[i][j]+x.num[i][k]*y.num[k][j]%mod)%mod;
			}
		}
	}
	return res;
}

matrix operator ^(matrix x,int y){
	matrix res;
	res.clear();
	for(int i=1;i<=n;++i){
		res.num[i][i]=1;		//单位矩阵 
	}
	while(y){
		if(y&1)	res=res*x;
		x=x*x;
		y=y>>1;
	} 
	return res;
}

void inner_link(){				//小点间建边 
	for(int i=1;i<=sn;++i){
		for(int j=1;j<=8;++j){
			a.num[9*(i-1)+j][9*(i-1)+j+1]=1;
		}
	} 
}

void outer_link(){				//集合间建边 
	for(int i=1;i<=sn;++i){
		scanf("%s",s+1);
		for(int j=1;j<=sn;++j){
			if(s[j]>'0'){
				a.num[9*(i-1)+s[j]-'0'][9*(j-1)+1]=1;
			}
		}
	}
} 

int main(){
	scanf("%d%d",&n,&t);
	sn=n;
	n*=9;
	inner_link();
	outer_link();
	a=a^t;
	printf("%d",a.num[1][sn*9-8]);
	return 0;
}