P4159 [SCOI2009]迷路

最新推荐文章于 2022-04-25 15:12:32 发布

_Sunstrider

最新推荐文章于 2022-04-25 15:12:32 发布

阅读量222

点赞数

CC 4.0 BY-SA版权

分类专栏：图论文章标签：邻接矩阵自乘

本文链接：https://blog.youkuaiyun.com/qq_37973923/article/details/87466934

图论专栏收录该内容

4 篇文章

订阅专栏

本文详细解析了SCOI2009迷路问题的算法解决方案，通过有向无权图邻接矩阵自乘获取方案数量，介绍了如何将有向有权图转换为有向无权图，并提供了具体的C++实现代码，关键在于正确的转化方式以确保算法复杂度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

P4159 [SCOI2009]迷路

模型总结

有向无权图邻接矩阵自乘获得方案数
有向有权图转化为有向无权图

关键点

注意转化方式，保证复杂度正确

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
const int inf=1e9;
int read(){
    int x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-') f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return f*x;
}
const int mod=2009;
int n,t,a[110][110],N,tmp[110][110];
char ch[15];
int cnt[110][110];
void mul1(){
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			for(int k=1;k<=N;k++){
				tmp[i][j]=(tmp[i][j]+cnt[i][k]*a[k][j])%mod;
			}
		}
	}
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			cnt[i][j]=tmp[i][j];
			tmp[i][j]=0;
		}
	}
}
void mul2(){
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			for(int k=1;k<=N;k++){
				tmp[i][j]=(tmp[i][j]+a[i][k]*a[k][j])%mod;
			}
		}
	}
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++){
			a[i][j]=tmp[i][j];
			tmp[i][j]=0;
		}
	}
}
void qpow(int b){
	for(int i=1;i<=N;i++){
		for(int j=1;j<=N;j++)
			cnt[i][j]=a[i][j];
	}
	b--;
	while(b){
		if(b&1) mul1();
		mul2();
		b>>=1;
	}
}
int main(){
	n=read(); t=read(); N=n*9;
	for(int i=1;i<=n;i++){
		scanf("%s",ch+1);
		for(int j=1;j<=n;j++){
			int w=ch[j]-'0';
			if(w) a[i*9][j*9-w+1]=1;
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=8;j++) a[(i-1)*9+j][(i-1)*9+j+1]=1; 
	}
	qpow(t);
	printf("%d\n",cnt[9][N]);
    return 0;
}