【HDU】2157 - How many ways??(矩阵快速幂)(从i到j经过k个点有多少种走法)
来源:互联网 发布:淘宝无线端其他来源 编辑:程序博客网 时间:2024/05/21 07:22
点击打开题目
How many ways??
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2721 Accepted Submission(s): 1022
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的时候输入结束
接下来的T行, 每行有三个整数 A, B, k, 表示问你从A 点到 B点恰好经过k个点的方案数 (k < 20), 可以走重复边。如果不存在这样的走法, 则输出0
当n, m都为0的时候输入结束
Output
计算每次询问的方案数, 由于走法很多, 输出其对1000取模的结果
Sample Input
4 40 10 21 32 320 3 20 3 33 60 11 00 22 01 22 121 2 10 1 30 0
Sample Output
2013
Author
小黑
Source
2008信息工程学院集训队——选拔赛
这里说一个知识点,好像是离散数学的内容,然而老师并没有讲。
这种题先构造矩阵A,(i , j)表示从 i 到 j 有无路径(有为1,否则为0)(注意此时的路径是单向的)。
然后从 i 点到 j 点经过 k 个点(终点也算),的走法数量为:求出A^k,然后A(i , j)的数值就是结果。
这道题就是这种方法,直接跑一个快速幂就出来了。
代码如下:
#include <cstdio>#include <cstring>#include <algorithm>using namespace std;#define MAX 20#define CLR(a,b) memset(a,b,sizeof(a))#define MOD 1000struct Matrix{__int64 m[MAX+4][MAX+4];int w,h;}ori,pr,res;void init(int n){CLR(res.m,0);res.h = res.w = n;for (int i = 0 ; i < n ; i++)//单位矩阵 res.m[i][i] = 1;pr = ori;}Matrix Matrix_multiply(Matrix x,Matrix y){Matrix t;CLR(t.m,0);t.h = x.h;t.w = y.w;for (int i = 0 ; i < x.h ; i++){for (int j = 0 ; j < x.w ; j++){if (x.m[i][j] == 0)continue;for (int k = 0 ; k < y.w ; k++)t.m[i][k] = (t.m[i][k] + x.m[i][j] * y.m[j][k] % MOD) % MOD;}}return t;}void Matrix_mod(int n){while (n){if (n & 1)res = Matrix_multiply(res , pr);pr = Matrix_multiply(pr , pr);n >>= 1;}}int main(){int n,m;while (~scanf ("%d %d",&n,&m) && (n || m)){CLR(ori.m,0);ori.h = ori.w = n;for (int i = 1 ; i <= m ; i++){int t1,t2;scanf ("%d %d",&t1,&t2);ori.m[t1][t2] = 1;}int Q;scanf ("%d",&Q);while (Q--){int st,endd,k;scanf ("%d %d %d",&st,&endd,&k);init(n);Matrix_mod(k);printf ("%I64d\n",res.m[st][endd]);}}return 0;}
0 0
- 【HDU】2157 - How many ways??(矩阵快速幂)(从i到j经过k个点有多少种走法)
- HDU 2157 How many ways?? 矩阵快速幂求A经过K个点到B方案数
- HDU 2157 How many ways??(经典矩阵快速幂)
- HDU 2157 How many ways??(矩阵快速幂)
- HDU 2157 How many ways??(矩阵快速幂)
- HDU 2157 How many ways??(矩阵的快速幂)
- HDU 2157 - How many ways??(矩阵快速幂)
- HDU 2157 How many ways??(矩阵快速幂)
- HDU 2157How many ways??(矩阵快速幂)
- 【矩阵快速幂】HDU 2157 How many ways??(矩阵快速幂经典问题)
- HDU 2157 How many ways??(简单线性DP | | 矩阵快速幂)
- HDU-2157 How many ways?? 矩阵快速幂
- hdu 2157 How many ways?? (矩阵快速幂)
- hdu 2157 How many ways?? 矩阵快速幂
- HDU 2157 How many ways?? 矩阵快速幂
- HDU 2157 How many ways??(矩阵快速幂)
- hdu 2157 How many ways?? 矩阵快速幂
- HDU 2157 How many ways?? (矩阵快速幂)
- hd1242 Rescue (bfs+队列优先)
- HDU 1849-Rabbit and Grass(尼姆博弈)
- 访问虚拟机上的本地网站
- windows安装python的mysql驱动mysql-connector-python
- 最全颜色表
- 【HDU】2157 - How many ways??(矩阵快速幂)(从i到j经过k个点有多少种走法)
- 解决android帮助文档打开慢
- js实现页面埋点
- Mysql远程不允许连接(Windows)
- Android Telephony分析(三) ---- RILJ详解
- 粒子滤波直方图
- 二叉树以及 二叉树的C++实现
- ios获取当前时间年月日时分秒
- 智能指针scoped_refptr