bzoj3504: [Cqoi2014]危桥

Posted 2020-08-25 Izaya

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

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

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

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

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

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

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

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

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

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

送上代码

 1 #include<bits/stdc++.h>
 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<int>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()<an+bn) flag=1;
73         if(!flag){
74             build(); insert(S,a2,an); insert(a1,T,an); insert(S,b1,bn); insert(b2,T,bn);
75             if(maxflow()<an+bn) flag=1;
76         }
77         if(flag) puts("No");
78         else puts("Yes");
79     }
80     return 0;
81 }