[POJ3613] Cow Relays(Floyd+矩阵快速幂)
解题报告
感觉这道题gyz大佬以前好像讲过一道差不多的?然鹅我这个蒟蒻发现矩阵快速幂已经全被我还给老师了...又恶补了一遍,真是恶臭啊。
题意
给定一个T(2 <= T <= 100)
条边的无向图,求S
到E
恰好经过N(2 <= N <= 1000000)
条边的最短路。
Idea
用Floyd和矩阵快速幂优化的产物。具体等咱下了课再好好说...
用一个矩阵a(i, j)来表示i到j经过若干条边的最短路,初始化a为i到j边的长度,没有则是正无穷。
然后重载*运算符,比如a矩阵表示经过n条边,b矩阵表示经过m条边,那么a * b得到的矩阵表示经过m + n条边,采用Floyd的思想进行更新。
为了省时间,用了快速幂,这样时间复杂度达到了O(n^3logk),可以通过。
#include<cstdio>
#include<iostream>
#include<map>
#include<algorithm>
#include<cstring>
#define ll long long
#define Dio ios::sync_with_stdio(0);
using namespace std;
const int maxn=200+1;
ll n,m,x,y,z,K;
map<int ,int> mapp;//数据中最多有100条边,说明最多有200个点,然鹅点的编号可达1000,所以离散化一下
struct Matrix{//矩阵的类
int a[maxn][maxn];
Matrix operator *(const Matrix& r){//按floyd重载矩阵乘法
Matrix c;
memset(c.a,0x3f,sizeof(c.a));
for(int k=1;k<=n;k++)//floyd模板
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++){
c.a[i][j]=min(c.a[i][j],a[i][k]+r.a[k][j]);
}
return c;
}
}st,ans;
void Pow(){//快速幂模板
ans=st;
K--;
while(K){
if(K&1) ans=ans*st;
st=st*st;
K>>=1;
}
}
int main(){
Dio
ll s,t;
cin>>K>>m>>s>>t;
memset(st.a,0x3f,sizeof(st.a));
for(int i=1;i<=m;i++){
int x,y,z;
cin>>z>>x>>y;
if(mapp[x]) x=mapp[x];
else x=mapp[x]=++n;
if(mapp[y]) y=mapp[y];
else y=mapp[y]=++n;
st.a[x][y]=st.a[y][x]=z;
}
Pow();
cout<<ans.a[mapp[s]][mapp[t]]<<endl;
return 0;
}
[POJ3613] Cow Relays(Floyd+矩阵快速幂)的更多相关文章
- POJ 3631 Cow Relays Floyd+矩阵快速幂
题目描述 For their physical fitness program, N (2 ≤ N ≤ 1,000,000) cows have decided to run a relay race ...
- POJ3613 Cow Relays(矩阵快速幂)
题目大概要求从起点到终点恰好经过k条边的最短路. 离散数学告诉我们邻接矩阵的k次幂就能得出恰好经过k条路的信息,比如POJ2778. 这题也一样,矩阵的幂运算定义成min,而min满足结合律,所以可以 ...
- poj 3613 Cow Relays【矩阵快速幂+Floyd】
!:自环也算一条路径 矩阵快速幂,把矩阵乘法的部分替换成Floyd(只用一个点扩张),这样每"乘"一次,就是经过增加一条边的最短路,用矩阵快速幂优化,然后因为边数是100级别的,所 ...
- poj 3613 经过k条边最短路 floyd+矩阵快速幂
http://poj.org/problem?id=3613 s->t上经过k条边的最短路 先把1000范围的点离散化到200中,然后使用最短路可以使用floyd,由于求的是经过k条路的最短路, ...
- foj 2173 floyd+矩阵快速幂
Problem 2173 Nostop Accept: 52 Submit: 210 Time Limit: 3000 mSec Memory Limit : 32768 KB Pro ...
- POJ 3613 floyd+矩阵快速幂
题意: 求s到e恰好经过n边的最短路 思路: 这题已经被我放了好长时间了. 原来是不会矩阵乘法,快速幂什么的也一知半解 现在终于稍微明白了点了 其实就是把矩阵乘法稍微改改 改成能够满足结合律的矩阵&q ...
- fzu 2173 floyd+矩阵快速幂
#include<stdio.h> #define inf 1000000000000000 #define N 100 long long tmp[N][N],ma[N][N]; int ...
- poj3613Cow Relays——k边最短路(矩阵快速幂)
题目:http://poj.org/problem?id=3613 题意就是求从起点到终点的一条恰好经过k条边的最短路: floyd+矩阵快速幂,矩阵中的第i行第j列表示从i到j的最短路,矩阵本身代表 ...
- 疯子的算法总结(九) 图论中的矩阵应用 Part 1+POJ3613 Cow Relays
图的存储有邻接矩阵,那么他就具备一些矩阵的性质,设有一个图的demo[100][100];那么demo[M][N]就是M—>N的距离,若经过一次松弛操作demo[M][N]=demo[M][K] ...
随机推荐
- PAT 反转链表
给定一个常数 K 以及一个单链表 L,请编写程序将 L 中每 K 个结点反转.例如:给定 L 为 1→2→3→4→5→6,K 为 3,则输出应该为 3→2→1→6→5→4:如果 K 为 4,则输出应该 ...
- Python内存管理机制-《源码解析》
Python内存管理机制 Python 内存管理分层架构 /* An object allocator for Python. Here is an introduction to the layer ...
- list基本运用
#include<iostream> #include<list> using namespace std; list<int>list1,list2; void ...
- 构造函数继承与class继承
构造函数继承 1.子类通过apply方法或者call方法把this指向父类 js代码 function Parent(name, age) { this.name = name this.age = ...
- react后台管理系统路由方案及react-router原理解析
最近做了一个后台管理系统主体框架是基于React进行开发的,因此系统的路由管理,选用了react-router(4.3.1)插件进行路由页面的管理配置. 实现原理剖析 1.hash的方式 ...
- 单例模式之懒汉模式,懒汉模式之高效模式,DLC双判断模式
import java.util.concurrent.LinkedBlockingQueue; import java.util.concurrent.ThreadPoolExecutor; imp ...
- v-else-if(v-show)
<div id="app"> <div v-if="type === 'A'"> A </div> <div v-el ...
- 【JMeter_06】JMeter逻辑控制器__If控制器<If Controller>
If控制器<If Controller> 业务逻辑: 根据表达式的结果来决定是否执行控制器下的脚本内容,与编程语言中的if判断逻辑大致相同,表达式结果为布尔值 true或false; 当表 ...
- c语言"##"的使用
#include<stdio.h> #define Operations(x) operation_ ## x // ## 是黏贴字符串 int Operations(sum)(int x ...
- Vue错误汇总
1.Vue导入js后没反应{{msg}}仍是{{msg}} 错误原因:js里写错了,或者没new一个Vue html页面: 解决方法:代码加入 new