HDOJ How many ways?? 2157【矩阵快速幂】

How many ways??

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)

Total Submission(s): 2046    Accepted Submission(s): 758

Problem Description

春天到了, HDU校园里开满了花, 姹紫嫣红, 非常美丽. 葱头是个爱花的人, 看着校花校草竞相开放, 漫步校园, 心情也变得舒畅. 为了多看看这迷人的校园, 葱头决定, 每次上课都走不同的路线去教室, 但是由于时间问题, 每次只能经过k个地方, 比方说, 这次葱头决定经过2个地方, 那他可以先去问鼎广场看看喷泉, 再去教室, 也可以先到体育场跑几圈, 再到教室. 他非常想知道, 从A 点恰好经过k个点到达B点的方案数, 当然这个数有可能非常大, 所以你只要输出它模上1000的余数就可以了. 你能帮帮他么??
你可决定了葱头一天能看多少校花哦

Input

输入数据有多组, 每组的第一行是2个整数 n, m(0 < n <= 20, m <= 100) 表示校园内共有n个点, 为了方便起见, 点从0到n-1编号,接着有m行, 每行有两个整数 s, t (0<=s,t<n) 表示从s点能到t点, 注意图是有向的.接着的一行是两个整数T,表示有T组询问(1<=T<=100),

接下来的T行, 每行有三个整数 A, B, k, 表示问你从A 点到 B点恰好经过k个点的方案数 (k < 20), 可以走重复边。如果不存在这样的走法, 则输出0

当n, m都为0的时候输入结束

Output

计算每次询问的方案数, 由于走法很多, 输出其对1000取模的结果

Sample Input

4 4
0 1
0 2
1 3
2 3
2
0 3 2
0 3 3
3 6
0 1
1 0
0 2
2 0
1 2
2 1
2
1 2 1
0 1 3
0 0

Sample Output

2
0
1
3

Author

小黑

Source

2008信息工程学院集训队——选拔赛

Recommend

linle   |   We have carefully selected several similar problems for you:  2159 2154 2158 2156 2153

中文题~

解题思路:经典矩阵算法,把给定的图转为邻接矩阵,即A(i,j)=1当且仅当存在一条边i->j。令C=A*A,那么C(i,j)=ΣA(i,k)*A(k,j),实际上就等于从点i到点j恰好经过2条边的路径数(枚举k为中转点)。类似地,C*A的第i行第j列就表示从i到j经过3条边的路径数。同理,如果要求经过k步的路径数,我们只需要二分求出A^k即可。

第一道矩阵快速幂。写的比较乱,而且这种写法时间复杂度较高,没有优化。不过比较容易看懂。

矩阵快速幂预备知识:

①矩阵相乘规则:

矩阵与矩阵相乘 第一个矩阵的列数必须等于第二个矩阵的行数 假如第一个是m*n的矩阵 第二个是n*p的矩阵 则结果就是m*p的矩阵 且得出来的矩阵中元素具有以下特点:

第一行第一列元素为第一个矩阵的第一行的每个元素和第二个矩阵的第一列的每个元素乘积的和 以此类推 第i行第j列的元素就是第一个矩阵的第i行的每个元素与第二个矩阵第j列的每个元素的乘积的和。

②单位矩阵:

n*n的矩阵  mat ( i , i )=1;  任何一个矩阵乘以单位矩阵就是它本身 n*单位矩阵=n, 可以把单位矩阵等价为整数1.

③快速幂:

这里矩阵快速幂等价于整数的快速幂,这里不再详细讲述

上代码:

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <stdlib.h>
#include <algorithm>

int s[25][25];
int b[25][25];
int n,m;
int a[25][25];

void Mat(int x[25][25],int y[25][25],int modn)
{
	int c[25][25];
	memset(c,0,sizeof(c)); //记得初始化
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			for(int k=0;k<n;k++)
				c[i][j]=(c[i][j]+x[i][k]*y[k][j]%modn)%modn;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		x[i][j]=c[i][j];
}

int Matrix(int begin,int end,int k)
{
	for(int i=0;i<n;i++){ //初始化一个单位矩阵
		for(int j=0;j<n;j++){
			a[i][j]=(i==j);
		}
	}
	for(int i=0;i<n;i++){  //记得用s保存再赋给b,不然b值变了之后结果就不对了
		for(int j=0;j<n;j++){
			b[i][j]=s[i][j];
		}
	}
	while(k){
		if(k&1)Mat(a,b,1000);
		Mat(b,b,1000);
		k>>=1;
	}
	return a[begin][end];
}

int main()
{
	while(scanf("%d%d",&n,&m),n!=0||m!=0){
		int S,G;
		memset(b,0,sizeof(b));
		memset(s,0,sizeof(s));
		for(int i=0;i<m;i++){
			scanf("%d%d",&S,&G);
			s[S][G]=1;
		}
		int T;
		scanf("%d",&T);
		int B,E,k;
		while(T--){
			scanf("%d%d%d",&B,&E,&k);
			int res=Matrix(B,E,k);
			printf("%d\n",res);
		}
	}
	return 0;
}

版权声明:本文为博主原创文章,转载请注明出处。

时间: 2024-11-05 21:41:28

HDOJ How many ways?? 2157【矩阵快速幂】的相关文章

HDOJ 4686 Arc of Dream 矩阵快速幂

矩阵快速幂: 根据关系够建矩阵 , 快速幂解决. Arc of Dream Time Limit: 2000/2000 MS (Java/Others)    Memory Limit: 65535/65535 K (Java/Others) Total Submission(s): 2164    Accepted Submission(s): 680 Problem Description An Arc of Dream is a curve defined by following fun

HDOJ 5411 CRB and Puzzle 矩阵快速幂

直接构造矩阵,最上面一行加一排1.快速幂计算矩阵的m次方,统计第一行的和 CRB and Puzzle Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 133    Accepted Submission(s): 63 Problem Description CRB is now playing Jigsaw Puzzle. There

hdoj 1575 Tr A 【矩阵快速幂】

Tr A Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3020    Accepted Submission(s): 2251 Problem Description A为一个方阵,则Tr A表示A的迹(就是主对角线上各项的和),现要求Tr(A^k)%9973. Input 数据的第一行是一个T,表示有T组数据. 每组数据的第一行有

HDOJ 233 Matrix 5015【矩阵快速幂】

233 Matrix Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 1355    Accepted Submission(s): 806 Problem Description In our daily life we often use 233 to express our feelings. Actually, we may

hdoj 2157 How many ways?? 【矩阵快速幂】【求任意两点间的路径上 经过k个点的方案数】

How many ways?? Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2102    Accepted Submission(s): 771 Problem Description 春天到了, HDU校园里开满了花, 姹紫嫣红, 非常美丽. 葱头是个爱花的人, 看着校花校草竞相开放, 漫步校园, 心情也变得舒畅. 为了多看看这

HDU 2157 How many ways??:矩阵快速幂【i到j共经过k个节点的方法数】

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2157 题解: 给你一个有向图,n个节点m条边,问你从i到j共经过k个节点的方法数(不算i点). 题解: 先用邻接矩阵存图. 假设k = 2,那么从i到j的方法数 = ∑ way[i][x] * way[x][j] (0<=x<n && x!=i && x!=j) 诶?快看,那是矩阵乘法! 设邻接矩阵为A,若i到j有边则val[i][j] = 1. k = 2时答案矩

HDU 2157 How many ways?? 关系矩阵+快速幂

Problem Description 春天到了, HDU校园里开满了花, 姹紫嫣红, 非常美丽. 葱头是个爱花的人, 看着校花校草竞相开放, 漫步校园, 心情也变得舒畅. 为了多看看这迷人的校园, 葱头决定, 每次上课都走不同的路线去教室, 但是由于时间问题, 每次只能经过k个地方, 比方说, 这次葱头决定经过2个地方, 那他可以先去问鼎广场看看喷泉, 再去教室, 也可以先到体育场跑几圈, 再到教室. 他非常想知道, 从A 点恰好经过k个点到达B点的方案数, 当然这个数有可能非常大, 所以你只

Hdoj 1588 Gauss Fibonacci 【矩阵快速幂】

Gauss Fibonacci Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2584 Accepted Submission(s): 1078 Problem Description Without expecting, Angel replied quickly.She says: "I'v heard that you'r a ver

HDU 2157 How many ways??(经典矩阵快速幂)

题意:求A经过K个点到B方案数 方法一: 1个0 1 的矩阵 A a[i][j] = 1 表示i 到 j可达 或者说 i 到 j 有1条路 或者说i到j经过一个点的方案数 路可以重复走 而A2 = A* A a[i][j] 的含义是 从i到j经过2个点的方案数 A的k次方 A[i,j]代表 i到j走k步的方案有a[i][j] 矩阵乘法的定义居然和这个模型如此契合,佩服,所以要非常熟悉矩阵乘法的具体步骤才能在这个题目中抽象出矩阵乘法可以正好实现两个定点间的所有可能情况 方法二: 动态规划的思路,状