Nyoj 302 星际旅行[矩阵乘法求两点k步方案数]
来源:互联网 发布:21天学通java第7版pdf 编辑:程序博客网 时间:2024/04/30 06:35
题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=302
题目意思。。n*m的矩阵。。每个元素代表一个星球,每个元素的值为所属国家,也就是每个星球属于一个国家。。每个相邻(上下左右)的星球有一个航道。每个国家的任意的两个星球都有星际之门。。。。。问,从(1,1)星球到(n,m)星球的方案数。。。相同的方案为走过的星球顺序相同。。。。。
处理一下图就是很裸的问题。。。
处理图的方式有很多中。。。
我的方法为,把二维矩阵散列到一维上。。然后,两两判断是否有通路。。然后处理到邻接矩阵中。。。。直接矩阵乘法就好。。。
Code:
#include <iostream>#include <algorithm>#include <cstdio>#include <cstring>#include <cmath>using namespace std;#define LL long longconst int N = 105;const LL mod = 1000007;struct Matrix{ int n; LL a[N][N]; Matrix(){ memset(a, 0, sizeof(a)); }} ans, A;int node[N], n;int m, p;bool Judge(int x, int y){ if(x - n == y) return true; else if(x + n == y) return true; else if(x - 1 == y && y % n != 0) return true; else if(x + 1 == y && x % n != 0) return true; else return false;}Matrix operator * (Matrix a, Matrix b){ Matrix tmpans; tmpans.n = a.n; for(int i = 1; i <= a.n; i ++){ for(int j = 1; j <= a.n; j ++){ for(int k = 1; k <= a.n; k ++) tmpans.a[i][j] = (tmpans.a[i][j] + a.a[i][k] * b.a[k][j]) % mod;// printf("k = %d\n", tmpans.a[i][j]); } } return tmpans;}void power(int k){ while(k){ if(k & 1) ans = ans * A; A = A * A; k = k >> 1; }}int main(){// freopen("1.txt", "r", stdin); int T; cin >> T; while(T --){ cin >> n >> m >> p; for(int i = 0; i < m; i ++){ for(int j = 1; j <= n; j ++){ cin >> node[i * n + j]; } } for(int i = 1; i <= m * n; i ++){ for(int j = 1; j <= m * n; j ++){ A.a[i][j] = 0; ans.a[i][j] = 0; if(i == j) continue; if(node[i] == node[j]) A.a[i][j] = 1; else if(Judge(i, j)){ A.a[i][j] = 1; } } ans.a[i][i] = 1; } ans.n = n * m; A.n = n * m; power(p); printf("%lld\n", ans.a[1][n * m] % mod); } return 0;}
表示NYOJ上关于矩阵的题目数据的还是很强大的。。。。
0 0
- Nyoj 302 星际旅行[矩阵乘法求两点k步方案数]
- Nyoj 530 K steps[矩阵乘法求两点k步的方案数]
- 给定一个有向图,问从A点恰好走k步(允许重复经过边)到达B点的方案数---矩阵乘法
- hdoj 2157 How many ways?? 【矩阵快速幂】【求任意两点间的路径上 经过k个点的方案数】
- HDU 2157 How many ways?? 矩阵快速幂求A经过K个点到B方案数
- NYOJ题目460 项链(矩阵相乘最大步数)
- NYOJ 623 矩阵乘法
- 矩阵乘法 nyoj 623
- NYOJ 最少步数
- nyoj最小步数
- nyoj 最少步数
- NYOJ【最少步数】
- nyoj 最少步数
- NYOJ-最少步数
- NYOJ 最少步数
- nyoj-最少步数
- nyoj-最小步数
- NYOJ 最少步数
- iOS推送通知(Push Notification)的Erlang实现
- ACM基本输入输出
- 小白书39页sprintf与格式控制
- 最大子段和
- SQLSERVER存储过程分页
- Nyoj 302 星际旅行[矩阵乘法求两点k步方案数]
- nodejs调试问题(一)
- SSE _align_malloc 数据地址对齐
- 更新Xcode6.0.1之后,编译出现Use of undeclared identifier 'GL_BGRA_EXT'问题及解决方案
- POJ 2411 Mondriaan's Dream 状态压缩DP
- C++ STL Set 快速入门
- HDU 1025 Constructing Roads In JGShining's Kingdom
- Codeforces 468B Two Sets(二分图匹配)
- 实验一报告