BZOJ 2095 [POI2010]Bridges (最大流、欧拉回路)
洛谷上有這題,但是輸出方案缺SPJ。。(而且我也懶得輸出方案了)
題目鏈接: https://www.lydsy.com/JudgeOnline/problem.php?id=2095
題解: 首先判掉度數有奇數的特殊情況,一眼能看出來二分答案(二分下界要設成每條邊較小權值的最大值),然后轉化成:
給定一張圖,有些邊有方向,有些邊無方向,問是否能給無向邊定向使得所有邊形成歐拉回路。
看著這個完全沒覺得這像網絡流啊……
好吧只要想到用網絡流來做就非常簡單了
考慮歐拉回路不就是每個點入度等于出度?現在入度不等于出度對吧?那我們考慮“入度等于出度”這個條件像什么?像流量平衡對吧。在網絡流中每個點流入的流量等于流出的流量。
那么我們的目標是對于每個點都要求其入度出度相等均為\(\frac{du[i]}{2}\) (\(du[i]\)為\(i\)的度數,\(dui[i],duo[i]\)分別為入度出度)
現在我們已有的有向邊入度比出度多\(dui[i]-duo[i]\), 那么它在無向邊定向之后就要出度比入度多\(dui[i]-duo[i]\).
所以如果不加任何改動拿無向邊跑最大流是出度等于入度,那么從\(S\)連一條\(dui[i]-duo[i]\)的邊跑最大流就相當于出度比入度多\(dui[i]-duo[i]\)了。同理,如果出度比入度多,這條邊就要連向T,權值為差的絕對值。
判斷是否滿流即可。
時間復雜度\(O(MaxFlow(n,m)\times \log W)\) (\(W\)為權值)
UPD: 為啥網上題解和我都不一樣??我的做法是錯的嗎??求教
網上的題解都是: 對于無向邊先隨便定個向,然后加這條邊的反向邊邊權為\(1\),對于每個點與\(S\)或者\(T\)連的邊權是出入度之差除以\(2\), 因為每反向一條邊相當于入度\(+1\)出度\(-1\), 對差的影響是\(2\).
代碼
#include<cstdio> #include<cstdlib> #include<cstring> #include<algorithm> using namespace std;const int N = 1002; const int M = 8000; const int INF = 1e7;namespace MaxFlow {struct Edge{int v,w,nxt,rev;} e[(M<<1)+3];int n,en,s,t;int dep[N+3];int que[N+3];int fe[N+3];int te[N+3];void clear(){for(int i=1; i<=n; i++) fe[i] = te[i] = 0;for(int i=1; i<=en; i++) e[i].v = e[i].w = e[i].nxt = e[i].rev = 0;n = en = s = t = 0;}void addedge(int u,int v,int w){en++; e[en].v = v; e[en].w = w;e[en].nxt = fe[u]; fe[u] = en; e[en].rev = en+1;en++; e[en].v = u; e[en].w = 0;e[en].nxt = fe[v]; fe[v] = en; e[en].rev = en-1;}bool bfs(){for(int i=1; i<=n; i++) dep[i] = 0;int head = 1,tail = 1; que[tail] = s; dep[s] = 1;while(head<=tail){int u = que[head]; head++;for(int i=fe[u]; i; i=e[i].nxt){if(dep[e[i].v]==0 && e[i].w>0){dep[e[i].v] = dep[u]+1;tail++; que[tail] = e[i].v;}}}return dep[t]!=0;}int dfs(int u,int cur){if(u==t) {return cur;}int rst = cur;for(int i=te[u]; i; i=e[i].nxt){if(dep[e[i].v]==dep[u]+1 && e[i].w>0 && rst>0){int flow = dfs(e[i].v,min(rst,e[i].w));if(flow>0){rst -= flow; e[i].w -= flow; e[e[i].rev].w += flow;if(e[i].w>0) te[u] = i;if(rst==0) {return cur;}}}}return cur-rst;}int dinic(int _n,int _s,int _t){n = _n,s = _s,t = _t;int ret = 0;while(bfs()){for(int i=1; i<=n; i++) te[i] = fe[i];ret += dfs(s,INF);}return ret;} } using MaxFlow::addedge; using MaxFlow::dinic;struct AEdge {int u,v,w1,w2; } ae[M+3]; int dui[N+3],duo[N+3],du[N+3]; int n,m;void clear() {for(int i=1; i<=n; i++) dui[i] = duo[i] = 0;MaxFlow::clear(); }int main() {scanf("%d%d",&n,&m); int left = 0,right = 0;for(int i=1; i<=m; i++){scanf("%d%d%d%d",&ae[i].u,&ae[i].v,&ae[i].w1,&ae[i].w2);left = max(left,min(ae[i].w1,ae[i].w2));right = max(right,max(ae[i].w1,ae[i].w2));du[ae[i].u]++; du[ae[i].v]++;}for(int i=1; i<=n; i++){if(du[i]&1) {printf("NIE"); return 0;}}while(left<right){int mid = left+((right-left)>>1);int std = 0;for(int i=1; i<=m; i++){if(ae[i].w1<=mid && ae[i].w2<=mid){addedge(ae[i].u+2,ae[i].v+2,1);addedge(ae[i].v+2,ae[i].u+2,1);}else if(ae[i].w1<=mid){duo[ae[i].u]++; dui[ae[i].v]++;}else if(ae[i].w2<=mid){duo[ae[i].v]++; dui[ae[i].u]++;}}for(int i=1; i<=n; i++){if(dui[i]>duo[i]){addedge(i+2,2,dui[i]-duo[i]);std += dui[i]-duo[i];}else if(duo[i]>dui[i]){addedge(1,i+2,duo[i]-dui[i]);}}int ans = dinic(n+2,1,2);if(ans==std) {right = mid;}else {left = mid+1;}clear();}printf("%d\n",right);return 0; }總結
以上是生活随笔為你收集整理的BZOJ 2095 [POI2010]Bridges (最大流、欧拉回路)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: BZOJ 2127 happiness
- 下一篇: BZOJ 4006 Luogu P326