浅谈Floyd的三种用法 By cellur925
Floyd大家可能第一时间想到的是他求多源最短路的n³算法。其实它还有另外两种算法的嘛qwq。写一发总结好了qwq。
一、多源最短路
放段代码跑,注意枚举顺序,用邻接矩阵存图。本质是一种动规。
复杂度O(n³)。
1 for(int k=1;k<=n;k++) 2 for(int i=1;i<=n;i++) 3 for(int j=1;j<=n;j++) 4 f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
View Code
放个例题跑。
灾后重建
二、传递闭包
在交际网络中,给定若干个元素,若干个二元关系,关系有传递性。传递闭包就是一种“通过传递性推导出尽量多的元素之间关系的问题”,求出可确定排名的元素个数。
实现用一个布尔型的邻接矩阵,f[i][j]=1表示i与j有关系,否则则没有关系。
我们每次可以枚举k点,来解决那些间接相关的关系处理。
1 for(int k=1;k<=n;k++) 2 for(int i=1;i<=n;i++) 3 for(int j=1;j<=n;j++) 4 f[i][j]|=f[i][k]&f[k][j];
View Code
例题 [USACO08JAN]牛大赛Cow Contest
对于奶牛的编程能力,用f[i][j]=1表示i比j强,之后就是一个裸的传递闭包。跑一遍后n²统计每只牛它与其他牛的关系是否已经确定,意思就是说只要有f[i]j]=1或f[j][i]=1其中一个就行,来统计答案。
Code
1 #include<cstdio> 2 #include<algorithm> 3 4 using namespace std; 5 6 int n,m,ans; 7 int f[200][200]; 8 9 int main() 10 { 11 scanf("%d%d",&n,&m); 12 for(int i=1;i<=m;i++) 13 { 14 int x=0,y=0; 15 scanf("%d%d",&x,&y); 16 f[x][y]=1; 17 } 18 for(int k=1;k<=n;k++) 19 for(int i=1;i<=n;i++) 20 for(int j=1;j<=n;j++) 21 f[i][j]|=f[i][k]&f[k][j]; 22 for(int i=1;i<=n;i++) 23 { 24 int j; 25 for(j=1;j<=n;j++) 26 { 27 if(i==j) continue; 28 if(f[i][j]==0&&f[j][i]==0) break; 29 } 30 if(j>n) ans++; 31 } 32 printf("%d",ans); 33 return 0; 34 }
View Code
三、求无向图最小环
例题1 USACO4.1篱笆回路
这道题难在建图,图建好以后就是裸的跑floyd找最小环了。
(瞎说一句,这题竟然有个数组开了1000的空间,但是越界了呀qwq)
Code
1 /* 2 ID:cellur_2 3 TASK:fence6 4 LANG:C++ 5 */ 6 #include<cstdio> 7 #include<algorithm> 8 #include<cstring> 9 10 using namespace std; 11 const int inf=0x3f3f3f3f; 12 13 int n,num,ans=inf; 14 int dis[300][300],mapp[300][300]; 15 struct node{ 16 int len; 17 int lcnt,rcnt,lid,rid,id; 18 int l[300],r[300]; 19 }edge[300]; 20 21 int main() 22 { 23 scanf("%d",&n); 24 for(int i=1;i<=n;i++) 25 { 26 scanf("%d",&edge[i].id); 27 int x=edge[i].id; 28 scanf("%d",&edge[x].len); 29 scanf("%d%d",&edge[x].lcnt,&edge[x].rcnt); 30 for(int j=1;j<=edge[x].lcnt;j++) 31 scanf("%d",&edge[x].l[j]); 32 for(int j=1;j<=edge[x].rcnt;j++) 33 scanf("%d",&edge[x].r[j]); 34 } 35 for(int i=1;i<=n;i++) 36 {// lid 这条边左端点的点编号 37 // rid 这条边右端点的点编号 38 if(!edge[i].lid) edge[i].lid=++num; 39 for(int j=1;j<=edge[i].lcnt;j++) 40 { 41 int x=edge[i].l[j]; 42 bool flag=0; 43 for(int k=1;k<=edge[x].lcnt;k++) 44 if(edge[x].l[k]==i) 45 { 46 flag=1; 47 break; 48 } 49 if(flag) edge[x].lid=edge[i].lid; 50 else edge[x].rid=edge[i].lid; 51 } 52 if(!edge[i].rid) edge[i].rid=++num; 53 for(int j=1;j<=edge[i].rcnt;j++) 54 { 55 int x=edge[i].r[j]; 56 bool flag=0; 57 for(int k=1;k<=edge[x].lcnt;k++) 58 if(edge[x].l[k]==i) 59 { 60 flag=1; 61 break; 62 } 63 if(flag) edge[x].lid=edge[i].rid; 64 else edge[x].rid=edge[i].rid; 65 } 66 } 67 memset(mapp,0x3f,sizeof(mapp)); 68 memset(dis,0x3f,sizeof(dis)); 69 ans=dis[2][33]; 70 for(int i=1;i<=n;i++) mapp[i][i]=0,dis[i][i]=0; 71 for(int i=1;i<=n;i++) 72 { 73 int lid=edge[i].lid; 74 int rid=edge[i].rid; 75 int len=edge[i].len; 76 mapp[rid][lid]=mapp[lid][rid]=len; 77 dis[rid][lid]=dis[lid][rid]=len; 78 } 79 //floyd找最小环 80 for(int k=1;k<=num;k++) 81 { 82 for(int i=1;i<k;i++) 83 for(int j=i+1;j<k;j++) 84 ans=min(ans,dis[i][j]+mapp[i][k]+mapp[k][j]); 85 for(int i=1;i<=num;i++) 86 for(int j=1;j<=num;j++) 87 dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]); 88 } 89 printf("%d\n",ans); 90 return 0; 91 }
View Code
例题2 POJ 1734 Sightseeing Trip
其实是floyd找最小环的板子题,但是由于题目要求输出一种合法的方案,所以我们只要再开一个vector就行了。
Code
1 #include<cstdio> 2 #include<algorithm> 3 #include<vector> 4 #include<cstring> 5 6 using namespace std; 7 typedef long long ll; 8 9 int n,m; 10 int ans=0x3f3f3f3f; 11 int dis[200][200],mapp[200][200],pos[200][200]; 12 vector<int>path; 13 14 void get_path(int x,int y) 15 { 16 if(pos[x][y]==0) return ; 17 get_path(x,pos[x][y]); 18 path.push_back(pos[x][y]); 19 get_path(pos[x][y],y); 20 } 21 22 int main() 23 { 24 scanf("%d%d",&n,&m); 25 memset(dis,0x3f,sizeof(dis)); 26 for(int i=1;i<=n;i++) dis[i][i]=0; 27 for(int i=1;i<=m;i++) 28 { 29 int x=0,y=0,z=0; 30 scanf("%d%d%d",&x,&y,&z); 31 dis[x][y]=dis[y][x]=min(dis[x][y],z); 32 } 33 memcpy(mapp,dis,sizeof(dis)); 34 for(int k=1;k<=n;k++) 35 { 36 for(int i=1;i<k;i++) 37 for(int j=i+1;j<k;j++) 38 if((ll)mapp[i][j]+dis[j][k]+dis[i][k]<ans) 39 { 40 ans=mapp[i][j]+dis[i][k]+dis[k][j]; 41 path.clear(); 42 path.push_back(i); 43 get_path(i,j); 44 path.push_back(j); 45 path.push_back(k); 46 } 47 for(int i=1;i<=n;i++) 48 for(int j=1;j<=n;j++) 49 if(mapp[i][j]>mapp[i][k]+mapp[k][j]) 50 { 51 mapp[i][j]=mapp[i][k]+mapp[k][j]; 52 pos[i][j]=k; 53 } 54 } 55 if(ans==0x3f3f3f3f) 56 { 57 printf("No solution."); 58 return 0; 59 } 60 for(int i=0;i<path.size();i++) 61 printf("%d ",path[i]); 62 return 0; 63 }
View Code
转载于:https://www.cnblogs.com/nopartyfoucaodong/p/9590907.html
浅谈Floyd的三种用法 By cellur925相关推荐
- 浅谈js函数三种定义方式 四种调用方式 调用顺序
在Javascript定义一个函数一般有如下三种方式: 函数关键字(function)语句: function fnMethodName(x){alert(x);} 函数字面量(Function Li ...
- 从petshop中一实例谈using 的三种用法
在petshop4.0中有这样一段代码: Code 1 public IList<CategoryInfo> GetCategories() { 2 3 IList ...
- 浅谈STM32的三种Boot模式
文章目录 一.关于BOOT模式的介绍 二.实际例子进行分析 三.分析上电启动流程 四.小结 五.参考资料 一.关于BOOT模式的介绍 所谓启动,一般来说就是指我们下好程序后,重启芯片时,SYSCLK的 ...
- 浅谈STM32的三种Boot模式的差异以及用汇编语言设计一个LED灯程序
STM32的三种Boot模式的差异 1.三种方式的比较 用汇编程序完成LED的程序 1.三种方式的比较 所谓启动,一般来说就是指我们下好程序后,重启芯片时,SYSCLK的第4个上升沿,BOOT引脚的值 ...
- java布尔类型比较器_浅谈Java中几种常见的比较器的实现方法
在java中经常会涉及到对象数组的排序问题,那么就涉及到对象之间的比较问题. 通常对象之间的比较可以从两个方面去看: 第一个方面:对象的地址是否一样,也就是是否引用自同一个对象.这种方式可以直接使用& ...
- 浅谈常用的几种web攻击方式以及解决办法
身在互联网的时候,web在给我们带来便利的同时,有些人也在盯着这些便利,因此出现了攻击网站的现象.所以我们在开发的时候,要注意这些容易被攻击的地方,以及做好防御的措施,下面将介绍一些这些 常见的攻击手 ...
- using 的三种用法
using 有哪三种用法? 1)引入命名空间. 2)给命名空间或者类型起别名. 3)划定作用域.自动释放资源,使用该方法的类型必须实现了 System.IDisposable接口,当对象脱离作用域之后 ...
- static的三种用法,定义静态变量,静态函数,静态代码块!
static的三种用法,定义静态变量,静态函数,静态代码块! 1.定义静态变量 class PersonStatic { //静态变量的特点 //1.静态变量无需生成对象就可被调用,可以使用类名和对象 ...
- Vue——基本的代码结构和插值表达式、v-cloak||Vue指令之v-text和v-html||v-bind的三种用法||Vue指令之v-on
Vue的基本代码 <!DOCTYPE html> <html lang="en"><head><meta charset="UT ...
最新文章
- hihocoder offer收割编程练习赛11 B 物品价值
- webapi 异步返回
- 通过使用Byte Buddy,便捷地创建Java Agent
- config之安全(用户认证)
- Nginx的rewrite之if指令(一)
- MySQL使用二进制日志恢复数据库
- 二十、PHP框架Laravel学习笔记——模型的作用域
- 配置 Syslog 守护程序
- C语言二叉树字符统计,C语言实现二叉树-利用二叉树统计单词数目
- 万达辟谣王健林去世:已报警!造谣账号已封禁
- XMind2TestCase自定义测试用例设计模板
- 鸿蒙系统可以上外网吗,【图片】华为鸿蒙系统的厉害之处在于 你可能非用不可
!【手机吧】_百度贴吧...
- velodyne16驱动调试记录
- 在网页中创建自己的调试控制台
- 关于React报Too many re-renders. React limits the number of renders to prevent an infinite错误的解决方案
- [附源码]java毕业设计医疗预约系统
- 交易日九点到九点半的挂单撤单以及价格的一点心得
- 豆瓣电台WP7客户端 开发记录5
- 【学习笔记】《如何构建敏捷项目管理团队》第四章 改变自己的风格
- VS2019运行OpenGL时出现的常见错误及解决办法
热门文章
- __asm__ __volatile__ GCC的内嵌汇编语法 ATT汇编语言语法(Z)
- python实现最长公共子序列的求解
- api数据库管理_API管理平台如何增强您的数据科学项目
- 客户细分_客户细分:K-Means聚类和A / B测试
- 64位Ubuntu kylin 16.04使用fastboot下载内核到tiny4412开发板
- 大数据职业理解_到底什么是大数据,大数据职业的具体工作内容是什么?
- java命令_JAVA与模式之命令模式
- 计算机机房用户不规则行为,网络及网管机房管理理论练习
- matlab数据游标属性设置,将节点属性添加到图论图数据提示
- python -pass的用法