博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
bzoj3504: [Cqoi2014]危桥
阅读量:6158 次
发布时间:2019-06-21

本文共 2901 字,大约阅读时间需要 9 分钟。

题意:给出一个图,有的边可以无限走,有的只能走两次(从一头到另一头为一次),给定两个起点以及对应的终点以及对应要走几个来回,求判断是否能完成。

先来一个NAIVE的建图:直接限制边建为容量1,无限制为INF,按照原图连,然后跑最大流就可以了。

可惜这样还不够,因为有可能有一部分流量不是对应的起点流过来的,即两条路有流量交换,这样就不一定可以满足题意了。

解决方法是:再跑一遍网络流,但是建图要改变一下,即将a路线的起点终点调换一下(当然b也可以),再接着跑,如果仍然满足则是真的有解。

证明看了网上的,都说的不太清楚,然后与ihopenot大爷讨论了下,有了一个证明(如果错了请指出)。

我们现在要证如果一开始的第一遍网络流跑出来满足条件并且第二遍也满足,则说明一定有真正满足题意的方案存在。

如果有交叉流,我们画一张图。

(因为使用的画图软件不能用减号,所以用下划线代替)其中As为起点,At为终点(B同理)

我们第二次建图后,如果能没有交叉流,那很好,直接说明有这种方案(因为倒过来是一样的嘛),但我们纠结的是,如果第二次仍然有交叉流怎么办?

我们画出来,发现,Bs到As有x的流量,第一张图又有As到Bt的x流量,那么这就是另一条x流量,与第一张图合起来不就是Bn吗,同理,A也是可以这样证明有满足条件且不交叉的流,那么就说明,我们如果两遍跑出来都满足,那就真的是有解的了。

送上代码

1 #include
2 using namespace std; 3 #define N 5005 4 #define INF 1e9 5 inline int read(){ 6 int x=0,f=1; char a=getchar(); 7 while(a>'9' || a<'0') {
if(a=='-') f=-1; a=getchar();} 8 while(a>='0' && a<='9') x=x*10+a-'0',a=getchar(); 9 return x*f;10 }11 int n,d[N],head[N],cur[N],cnt,S,T,a1,a2,an,b1,b2,bn,ans;12 bool vis[N],g[55][55],b[55][55];13 queue
q;14 char st[55][55];15 struct edges{16 int to,cap,flow,next;17 }e[2*N];18 inline void insert(int u,int v,int c){19 e[cnt]=(edges){v,c,0,head[u]};head[u]=cnt++;20 e[cnt]=(edges){u,0,0,head[v]};head[v]=cnt++;21 }22 inline bool bfs(){23 memset(vis,0,sizeof(vis));24 d[S]=0; vis[S]=1; q.push(S);25 while(!q.empty()){26 int x=q.front(); q.pop();27 for(int i=head[x];i>=0;i=e[i].next)28 if(!vis[e[i].to] && e[i].cap>e[i].flow)29 d[e[i].to]=d[x]+1,vis[e[i].to]=1,q.push(e[i].to);30 }31 return vis[T];32 }33 int dfs(int x,int a){34 if(x==T || !a) return a;35 int f,flow=0;36 for(int& i=cur[x];i>=0;i=e[i].next){37 if(d[e[i].to]==d[x]+1 && (f=dfs(e[i].to,min(a,e[i].cap-e[i].flow)))>0)38 flow+=f,e[i].flow+=f,e[i^1].flow-=f,a-=f;39 if(!a) break;40 }41 return flow;42 }43 inline int maxflow(){44 int flow=0;45 while(bfs()){46 for(int i=S;i<=T;i++) cur[i]=head[i];47 flow+=dfs(S,INF);48 }49 return flow;50 }51 inline void build(){52 memset(head,-1,sizeof(head)); cnt=0;53 for(int i=1;i<=n;i++)54 for(int j=1;j<=n;j++)55 if(g[i][j]) insert(i,j,INF);56 else if(b[i][j]) insert(i,j,1);57 }58 int main(){59 S=0; T=5000;60 while(scanf("%d",&n)!=EOF){61 memset(g,0,sizeof(g));62 memset(b,0,sizeof(b));63 a1=read()+1; a2=read()+1; an=read();64 b1=read()+1; b2=read()+1; bn=read();65 for(int i=1;i<=n;i++) scanf("%s",st[i]+1);66 for(int i=1;i<=n;i++)67 for(int j=1;j<=n;j++)68 if(st[i][j]=='O') b[i][j]=1;69 else if(st[i][j]=='N') g[i][j]=1;70 bool flag=0;71 build(); insert(S,a1,an); insert(a2,T,an); insert(S,b1,bn); insert(b2,T,bn);72 if(maxflow()

 

转载于:https://www.cnblogs.com/enigma-aw/p/6230579.html

你可能感兴趣的文章
C++类的继承一
查看>>
数据库分库分表(sharding)系列(五) 一种支持自由规划无须数据迁移和修改路由代码的Sharding扩容方案...
查看>>
巧用VMware Workstation的clone来制作虚拟机模板
查看>>
Spring-Mybatis MapperScannerConfigurer 取不到PropertyPlaceholderConfigurer里的值
查看>>
HP DL380G4服务器前面板指示灯的含义
查看>>
数据结构_树结构
查看>>
常用URL地址
查看>>
每天一个linux命令(19):find 命令概览
查看>>
MySQL kill操作
查看>>
windows下看端口占用
查看>>
Decommissioning a Domain Controller 降域控
查看>>
Character中的奇葩
查看>>
c++书籍推荐
查看>>
互联网通用架构技术----缓存雪崩
查看>>
Dell R710服务器磁盘恢复数据库一例(记录)
查看>>
我的友情链接
查看>>
Ionic3 通讯录索引的实现
查看>>
轻松监听Azure service health 状态
查看>>
Matlab 进行FFT
查看>>
Eclipse 工作台用户指导>视图和编辑器
查看>>