hdu2157 How many ways??
来源:互联网 发布:魔盒软件最新版 编辑:程序博客网 时间:2024/06/05 02:15
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2181 Accepted Submission(s): 816
Total Submission(s): 2181 Accepted Submission(s): 816
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
这题可以用矩阵快速幂做,先把矩阵用邻接矩阵表示出来,如果i,j之间有路,那么a[i][j]就为1,否则为0.那么把两个a相乘,得到c,c[i][j]=Σa(i,k)*a(k,j),就是从i经过一个地方k走到j,那么a^k就是从i经过k个地方走到j的路线条数。
#pragma comment(linker, "/STACK:102400000,102400000")#include<iostream>#include<stdio.h>#include<stdlib.h>#include<string.h>#include<math.h>#include<vector>#include<map>#include<set>#include<queue>#include<stack>#include<string>#include<algorithm>using namespace std;typedef long long ll;#define inf 99999999#define pi acos(-1.0)#define MOD 1000int gra[25][25];struct matrix{ ll n,m,i; ll data[99][99]; void init_danwei(){ for(i=0;i<n;i++){ data[i][i]=1; } }};matrix multi(matrix &a,matrix &b){ ll i,j,k; matrix temp; temp.n=a.n; temp.m=b.m; for(i=0;i<temp.n;i++){ for(j=0;j<temp.m;j++){ temp.data[i][j]=0; } } for(i=0;i<a.n;i++){ for(k=0;k<a.m;k++){ if(a.data[i][k]>0){ for(j=0;j<b.m;j++){ temp.data[i][j]=(temp.data[i][j]+(a.data[i][k]*b.data[k][j])%MOD )%MOD; } } } } return temp;}matrix fast_mod(matrix &a,ll n){ matrix ans; ans.n=a.n; ans.m=a.m; memset(ans.data,0,sizeof(ans.data)); ans.init_danwei(); while(n>0){ if(n&1)ans=multi(ans,a); a=multi(a,a); n>>=1; } return ans;}int main(){ ll n,m,i,j,c,d,k,T; while(scanf("%lld%lld",&n,&m)!=EOF) { if(n==0 && m==0)break; matrix a; memset(a.data,0,sizeof(a.data)); a.n=a.m=n; for(i=1;i<=m;i++){ scanf("%lld%lld",&c,&d); a.data[c][d]=1; } scanf("%lld",&T); while(T--) { scanf("%lld%lld%lld",&c,&d,&k); matrix b; memset(b.data,0,sizeof(b.data)); b.n=b.m=n; for(i=0;i<n;i++){ for(j=0;j<n;j++){ b.data[i][j]=a.data[i][j]; } } matrix ans; ans=fast_mod(b,k); printf("%lld\n",ans.data[c][d]%MOD); } } return 0;}
0 0
- hdu2157 How many ways??
- hdu2157 How many ways?
- hdu2157 how many ways 矩阵
- hdu2157 How many ways?? (矩阵快速幂)
- hdu2157 How many ways?? 矩阵快速幂
- hdu2157 How many ways?? (矩阵快速幂)
- [HDU2157]How many ways??(矩阵乘法)
- hdu2157 How many ways??(矩阵快速幂
- #HDU2157#How many ways?(矩阵基础)
- (矩阵快速幂)hdu2157 How many ways??
- 【矩阵快速幂】经典题 hdu2157 how many ways、woj642 Lost In WHU
- hdu1978 How many ways
- How many ways??
- How many ways
- HDU1978--How many ways
- hdu1978 How many ways
- How many ways-dp
- hdu how many ways
- svn的merge使用例子
- Win7系统下无法运行调试提示BORdbk70.dll未注册问题解决
- 苹果漏洞:Siri会泄露你的个人资料!
- Android_(自定义view_1之牛刀小试)
- Android支付平台集成调研
- hdu2157 How many ways??
- Xcode如何导出ipa包
- 关于EditText的一点深入的了解
- 融云-添加好友流程
- iOS 【UIKit-UINavigationItem 定制导航栏中元素的信息】
- 文章标题
- SpringMVC整合Swagger框架
- Entrence (.Lib .Dll)
- 每周小结