poj 2195 Going Home

  1 /*
  2    做网络流的题建图真的是太重要了!
  3    本题是将人所在的位置和房子所在的位置建立边的联系,其中man到house这一条边的流量为 1, 费用为两者的距离
  4    而方向边的流量为 0, 费用为正向边的相反数(也就是沿着反向边进行增广时,费用要减少,更改先前错误的选择) 
  5    最后增加一个源点和一个汇点完毕(源点映射到man, house映射到汇点, 费用为0, 流量为1) 
    
6 */ 7 #include<iostream> 8 #include<cmath> 9 #include<cstdio> 10 #include<cstring> 11 #include<queue> 12 #define Max 0x3f3f3f3f 13 #define N 205 14 using namespace std; 15 16 class node 17 { 18 public: 19 int x, y; 20 }; 21 22 node xyM[N]; 23 node xyH[N]; 24 int cost[N][N], cap[N][N]; 25 int cntM, cntH; 26 int pre[N*2], dist[N*2], vis[N*2], n, m; 27 28 void addE(int i, int j, int ct, int cp) 29 { 30 cost[i][j]=ct; 31 cap[i][j]=cp; 32 cost[j][i]=-ct; 33 //cap[j][i]=0; 34 } 35 36 int s, t, minCost; 37 38 void buildMap() 39 { 40 int i, j; 41 memset(cap, 0, sizeof(cap)); 42 s=0; t=cntM+cntH+1; 43 for(i=0; i<cntM; ++i) 44 addE(0, i+1, 0, 1); 45 for(i=0; i<cntH; ++i) 46 addE(cntM+i+1, t, 0, 1); 47 for(i=0; i<cntM; ++i) 48 for(j=0; j<cntH; ++j) 49 addE(i+1, cntM+j+1, abs(xyM[i].x-xyH[j].x)+abs(xyM[i].y-xyH[j].y), 1); 50 } 51 52 queue<int>q; 53 54 int spfa() 55 { 56 int u, v; 57 memset(dist, 0x3f, sizeof(dist)); 58 dist[0]=0; 59 q.push(0); 60 vis[0]=1; 61 while(!q.empty()) 62 { 63 u=q.front(); 64 q.pop(); 65 vis[u]=0; 66 for(v=0; v<=t; ++v) 67 if(cap[u][v]>0 && dist[v]>dist[u]+cost[u][v]) 68 { 69 dist[v]=dist[u]+cost[u][v]; 70 pre[v]=u; 71 if(!vis[v]) 72 { 73 vis[v]=1; 74 q.push(v); 75 } 76 } 77 } 78 if(dist[t]==Max) 79 return 0; 80 return 1; 81 } 82 83 void updateEdge() 84 { 85 int u, minFlow=Max; 86 for(u=t; u!=s; u=pre[u])//通过最短路径寻找这条路径上的最小流量 87 if(cap[pre[u]][u]<minFlow) 88 minFlow=cap[pre[u]][u]; 89 for(u=t; u!=s; u=pre[u]) 90 { 91 cap[pre[u]][u]-=minFlow; 92 cap[u][pre[u]]+=minFlow; 93 minCost+=cost[pre[u]][u]; 94 } 95 } 96 97 int main() 98 { 99 int i, j; 100 char c; 101 while(scanf("%d%d", &n, &m) && (n||m)) 102 { 103 cntM=cntH=0; 104 minCost=0; 105 for(i=1; i<=n; ++i) 106 { 107 getchar(); 108 for(j=1; j<=m; ++j) 109 { 110 scanf("%c", &c); 111 if(c=='m') 112 { 113 xyM[cntM].x=i; 114 xyM[cntM++].y=j; 115 } 116 else if(c=='H') 117 { 118 xyH[cntH].x=i; 119 xyH[cntH++].y=j; 120 } 121 } 122 } 123 buildMap(); 124 while(spfa()) 125 updateEdge(); 126 printf("%d\n", minCost); 127 } 128 return 0; 129 }

 //邻接表

  1 #include<iostream>
  2 #include<queue>
  3 #include<cstring>
  4 #include<cstdio>
  5 #include<algorithm>
  6 #define INF 0x3f3f3f3f
  7 #define N 1000005
  8 using namespace std;
  9 
 10 int cntH, cntM;
 11 
 12 struct node{
 13     int x, y;
 14 };
 15 
 16 struct EDGE{
 17     int u, v, cap, cost, nt;
 18 };
 19 EDGE edge[N];
 20 
 21 queue<int>q;
 22 node man[105], house[105];
 23 int first[205];
 24 int dist[205];
 25 int pre[205], flow[205], vis[205];
 26 int cnt, t;
 27 int minCost;
 28 int n, m;
 29 
 30 void addEdge(int u, int v, int cap, int cost){
 31     edge[cnt].u=u;
 32     edge[cnt].v=v;
 33     edge[cnt].cap=cap;
 34     edge[cnt].nt=first[u];
 35     edge[cnt].cost=cost;
 36     first[u]=cnt++;
 37     
 38     edge[cnt].u=v;
 39     edge[cnt].v=u;
 40     edge[cnt].cap=0;
 41     edge[cnt].nt=first[v];
 42     edge[cnt].cost=-cost;
 43     first[v]=cnt++;
 44 }
 45 
 46 void buildMap(){
 47     memset(first, -1, sizeof(first));
 48     t=cntH+cntM+1;
 49     for(int i=1; i<=cntM; ++i)
 50         for(int j=1; j<=cntH; ++j)
 51             addEdge(i, cntM+j, 1, abs(man[i].x-house[j].x) + abs(man[i].y-house[j].y));
 52     for(int i=1; i<=cntM; ++i)
 53         addEdge(0, i, 1, 0);
 54     for(int i=1; i<=cntH; ++i)
 55         addEdge(cntM+i, t, 1, 0); 
 56 }
 57 
 58 bool MCMF(){
 59     memset(dist, 0x3f, sizeof(dist));
 60     memset(vis, 0, sizeof(vis));
 61     q.push(0);
 62     flow[0]=INF;
 63     dist[0]=0;
 64     vis[0]=1;
 65     while(!q.empty()){
 66         int u=q.front(); q.pop();
 67         vis[u]=0;
 68         for(int e=first[u]; ~e; e=edge[e].nt){
 69             int v=edge[e].v, cap=edge[e].cap, cost=edge[e].cost;
 70             if(cap>0 && dist[v]>dist[u]+cost){
 71                 dist[v]=dist[u]+cost;
 72                 flow[v]=min(flow[u], cap);
 73                 pre[v]=e;
 74                 if(!vis[v]){
 75                     vis[v]=1;
 76                     q.push(v);
 77                 }
 78             }
 79         }
 80     }
 81     if(dist[t]==INF) return false;
 82     minCost+=dist[t];
 83     int x=t;
 84     while(x!=0){
 85         edge[pre[x]].cap-=flow[t];
 86         edge[pre[x]^1].cap+=flow[t];
 87         x=edge[pre[x]].u;
 88     }
 89     return true;
 90 }
 91  
 92 int main(){
 93     while(scanf("%d%d", &n, &m) && (n||m)){
 94         cnt=cntH=cntM=0;
 95         for(int i=1; i<=n; ++i){
 96             getchar();
 97             for(int j=1; j<=m; ++j){
 98                 char ch;
 99                 scanf("%c", &ch);
100                 if(ch=='m'){
101                     man[++cntM].x=i;
102                     man[cntM].y=j;
103                 }
104                 else if(ch=='H'){
105                     house[++cntH].x=i;
106                     house[cntH].y=j;
107                 }
108             }
109         }
110         buildMap();
111         minCost=0;
112         while(MCMF());
113         printf("%d\n", minCost);
114     }
115     return 0;
116 } 

 

转载于:https://www.cnblogs.com/hujunzheng/p/3809620.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/531961.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

CardLayout布局练习(小的图片浏览器)

1 /*2 涉及Panel中的图片的加载&#xff0c;还有Frame的关闭的方法&#xff0c; CardLayout&#xff08;int hgap, int vgap&#xff09;就会决定卡片面板的大小3 匿名类的使用。。。4 */5 import java.awt.*;6 import java.awt.event.*;7 import javax.swing.*;8 public class…

python求逆矩阵的方法,Python 如何求矩阵的逆

我就废话不多说了&#xff0c;大家还是直接看代码吧~import numpy as npkernel np.array([1, 1, 1, 2]).reshape((2, 2))print(kernel)print(np.linalg.inv(kernel))注意&#xff0c;Singular matrix奇异矩阵不可求逆补充&#xff1a;pythonnumpy中矩阵的逆和伪逆的区别定义&a…

plsql存过声明游标_plsql编程学习之游标一

oralce plsql编程的游标游标分类1显示游标2隐式游标隐式游标&#xff0c;oracle自动管理&#xff0c;不用声明&#xff0c;打开和关闭&#xff0c;ORACLE自动处理&#xff0c;使用隐式游标%FOUND时&#xff0c;需要加上 SQL%FOUND显示游标&#xff0c;需要自己声明&#xff0c;…

用命令行编译java并生成可执行的jar包

用命令行编译java并生成可执行的jar包 1.编写源代码。 编写源文件&#xff1a;CardLayoutDemo.java并保存&#xff0c;例如&#xff1a;I:\myApp\CardLayoutDemo.java。程序结构如下&#xff1a;package test;import java.awt.*; import javax.swing.*; //更多包的导入...clas…

python计时器单位,python(计时器)

计时器要求&#xff1a;定制一个计时器的类start 和 stop方法代表启动计时和停止计时假设计时器对象 t1&#xff0c;print(t1)和直接调用t1 均显示结果当计时器未启动或已停止计时&#xff0c;调用stop方法能给予温馨提示两个计时器对象可以相加&#xff1a; t1 t2只能使用提供…

查询分析300万笔记录_给你100万条数据的一张表,你将如何查询优化?

1.两种查询引擎查询速度(myIsam 引擎)InnoDB 中不保存表的具体行数&#xff0c;也就是说&#xff0c;执行select count(*) from table时&#xff0c;InnoDB要扫描一遍整个表来计算有多少行。MyISAM只要简单的读出保存好的行数即可。注意的是&#xff0c;当count(*)语句包含 whe…

poj 3321 Apple Trie

/*poj 3321 Apple Trie 这道题的关键是如何将一个树建成一个一维数组利用树状数组来解题&#xff01;可以利用dfs&#xff08;&#xff09;来搞定&#xff0c;我们在对一个节点深搜后&#xff0c;所经过的节点的数目就是该节点的子树的数目所以我们利用start[i]数组来记录 i 节…

php美团项目分享,美团项目(纯代码)(示例代码)

一.框架搭建1.icon规格要求可从文档中查找,搜索app icon.2.因为很多界面重复利用,所以不用storyboarda.删除stroyboard,在设置中Info -> Main storyboard file base name 项直接去除b.创建ZXHomeViewController(UICollectionViewController)和ZXNavigationController(UINavi…

ioc spring 上机案例_Spring的IoC入门案例

1、创建工程&#xff0c;导入坐标1.1 创建工程1.2 导入坐标xmlns:xsi"http://www.w3.org/2001/XMLSchema-instance"xsi:schemaLocation"http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">4.0.0org.examplespring_01_io…

java中父类与子类, 不同的两个类中的因为构造函数由于递归调用导致栈溢出问题...

1 /*2 对于类中对成员变量的初始化和代码块中的代码全部都挪到了构造函数中&#xff0c;3 并且是按照java源文件的初始化顺序依次对成员变量进行初始化的&#xff0c;而原构造函数中的代码则移到了构造函数的最后执行4 */5 import static java.lang.System.out;6 7 public clas…

liunx php的项目地址,在 Linux 配置 PHP 项目

在 Linux 配置 PHP 项目一, 搭建测试环境软件环境:(PHP 项目)PHP5.4Apache(httpd2.4)mysql5.7二, 安装1挂载:1. 把 iso 的镜像文件放到虚拟机 Linux 的 CD/ROM(在右下角 (网络适配器 / 桥接模式) 旁有个光盘, 点击连接, 之后页面出现一个光盘)2. 使用挂载命令, 把 CD/ROM 设备里…

springwebflux 页面_【SpringBoot WEB系列】WebFlux静态资源配置与访问

上一篇博文介绍SpringMVC的静态资源访问&#xff0c;那么在WebFlux中&#xff0c;静态资源的访问姿势是否一致呢I. 默认配置与SpringBoot的默认配置一样&#xff0c;WebFlux同样是classpath:/META-INF/resources/,classpath:/resources/,classpath:/static/,classpath:/public/…

java中TreeSet集合如何实现元素的判重

1 /*2 看一下部分的TreeSet源码....3 public class TreeSet<E> extends AbstractSet<E>4 implements NavigableSet<E>, Cloneable, java.io.Serializable5 {6 private transient NavigableMap<E,Object> m;7 //NavigableMap继承SortedMap&…

php中改变函数路由,通过PHP重启路由器以更换IP(原创)

在采集大批量数据时常常会触发对方服务器的“自我保护”&#xff0c;请求过于频繁就限制访问。这时需要停留很长一段时间(十几分钟到几十分钟不等)才能恢复访问&#xff0c;这样采集数据的速度就受到非常大的限制。解决方法有两个&#xff1a;1 通过图片识别绕过验证码机制&…

axure 画小程序效果图_APP详情页如何用Axure画出来

详情页是App原型中比较复杂的页面类型&#xff0c;熟悉它的常用套路有助于快速画出。之前的文章已经讲解了APP常见功能中的页面模板、下导航、上导航、列表页怎么画出来&#xff0c;请继续关注浪子教你画APP原型后续的其他功能模块。APP详情页往往包含上导航&#xff0c;内容区…

HashSet中实现不插入重复的元素

/* 看一下部分的HashSet源码.... public class HashSet<E>extends AbstractSet<E>implements Set<E>, Cloneable, java.io.Serializable {static final long serialVersionUID -5024744406713321676L;private transient HashMap<E,Object> map;privat…

tuxedo错误码6_TUXEDE返回的所有错误代码

TUXEDE返回的所有错误代码tuxedo/include/atmi.h定于了TUXEDE返回的所有错误代码。/** tperrno values - error codes* The man pages explain the context in which the following error codes* can return.*/#define TPMINVAL 0 /* minimum error message */#define TPEABORT…

java中finally和return的执行顺序

注意&#xff1a;return的位置。。。从这几个例子中可以看到&#xff0c;如果try之前没有有条件的return&#xff0c;则try..catch..finally语句块中的语句都是顺序执行&#xff08;如果try中或者catch中 有return语句&#xff0c;那么先执行该return&#xff0c;然后执行final…

oracle如何设置权限,ORACLE的权限设置

创建用户create user abc identified by 123;----------------------------------------------------授权grant create session,create table to abcgrant create sysdba to database----------------------------------------------------然后conn abc密码&#xff1a;123----…

有关try..catch..finally处理异常的总结

//看一下下面的程序&#xff0c;你能正确的写出不同的testEx2()方法时&#xff0c;程序的最终打印出来的数据吗....先不要看下面的答案 public class ExceptionTest { public ExceptionTest() { } boolean testEx() throws Exception { boolean ret true; try { ret te…