hdu 4309 最大流 + DFS
生活随笔
收集整理的這篇文章主要介紹了
hdu 4309 最大流 + DFS
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題意:
? ? ?給以三種有向邊
? ? (1) 隧道,可以過無數人,也可以藏c個人。
? ? (2) 路,只能過人(流量INF)。
? ? (3)古橋,如果不修理可以過1個人,修理可以過無數個人,但是要花費c那么多錢?同時還給了每個城市的人數,要求是城市的人去隧道里躲避,問你在躲避人數最多的前提下費用最小,費用就是修理古橋所用的錢。
思路:
? ? ? ?先看古橋數量,最多12個,2^12 也就四千多,直接二分枚舉哪些古橋是修理的,然后對于每一種狀態,我們都重新建圖,用最大流去找最多藏多少人,然后更新最優就行了。下面說下建圖吧,這個題目的建圖不難,直接建就行。
?(0)建立超級遠點和匯點,s,t。
?(1)隧道,既然可以藏人,就可以虛擬出一個點來,比如給的隧道a(點),b(點),c(人數)
我們可以虛擬出來一個點d,add(a ,d ,INF),add(d ,b ,INF) ,add(d ,t ,c);
?(2)路,這個沒話說 add(a ,b ,INF);
? ? ?給以三種有向邊
? ? (1) 隧道,可以過無數人,也可以藏c個人。
? ? (2) 路,只能過人(流量INF)。
? ? (3)古橋,如果不修理可以過1個人,修理可以過無數個人,但是要花費c那么多錢?同時還給了每個城市的人數,要求是城市的人去隧道里躲避,問你在躲避人數最多的前提下費用最小,費用就是修理古橋所用的錢。
思路:
? ? ? ?先看古橋數量,最多12個,2^12 也就四千多,直接二分枚舉哪些古橋是修理的,然后對于每一種狀態,我們都重新建圖,用最大流去找最多藏多少人,然后更新最優就行了。下面說下建圖吧,這個題目的建圖不難,直接建就行。
?(0)建立超級遠點和匯點,s,t。
?(1)隧道,既然可以藏人,就可以虛擬出一個點來,比如給的隧道a(點),b(點),c(人數)
我們可以虛擬出來一個點d,add(a ,d ,INF),add(d ,b ,INF) ,add(d ,t ,c);
?(2)路,這個沒話說 add(a ,b ,INF);
?(3)古橋,對于深搜的當前狀態的當前這個古橋,如果是修理add(a ,b ,INF),如果不是修理add(a ,b ,1)這樣一邊重建,一邊深搜更新最大值就行了。
#include<stdio.h> #include<queue> #include<string.h>#define N_node 500 #define N_edge 10000 #define INF 1000000000 using namespace std;typedef struct {int to ,next ,cost; }STAR;typedef struct {int a ,b ,c; }NODE;typedef struct {int x ,t; }DEP;STAR E[N_edge]; NODE L[1100] ,G[15] ,S[25]; int n1 ,n2 ,n3; DEP xin ,tou; int list[N_node] ,tot; int list1[N_node]; int deep[N_node]; int mk[15]; int pp[1100];void add(int a ,int b ,int c) {E[++tot].to = b;E[tot].cost = c;E[tot].next = list[a];list[a] = tot;E[++tot].to = a;E[tot].cost = 0;E[tot].next = list[b];list[b] = tot; }int minn(int x ,int y) {return x < y ? x : y; }bool BFS_Deep(int s ,int n ,int t) {memset(deep ,255 ,sizeof(deep));deep[s] = 0;xin.x = s ,xin.t = 0;queue<DEP>q;q.push(xin);while(!q.empty()){tou = q.front();q.pop();for(int k = list[tou.x] ;k ;k = E[k].next){xin.x = E[k].to;xin.t = tou.t + 1;if(deep[xin.x] != -1 || !E[k].cost)continue;deep[xin.x] = xin.t;q.push(xin);}}for(int i = 0 ;i <= n ;i ++)list1[i] = list[i];return deep[t] != -1; }int DFS_flow(int s ,int t ,int flow) {if(s == t) return flow;int nowflow = 0;for(int k = list1[s] ;k ;k = E[k].next){list1[s] = k;int to = E[k].to;int c = E[k].cost;if(deep[to] != deep[s] + 1 || !c)continue;int tmp = DFS_flow(to ,t ,minn(c ,flow - nowflow));nowflow += tmp;E[k].cost -= tmp;E[k^1].cost += tmp;if(flow == nowflow) break; } if(!nowflow) deep[s] = 0;return nowflow; }int DINIC(int s ,int t ,int n) {int ans = 0;while(BFS_Deep(s ,t ,n)){ans += DFS_flow(s ,t ,INF);}return ans; }void Rebuid(int n) {memset(list ,0 ,sizeof(list)) ,tot = 1;int T = n + n1 + 1;for(int i = 1 ;i <= n1 ;i ++){int a = S[i].a ,b = S[i].b ,c = S[i].c;add(i + n ,T ,c);add(a ,i + n ,INF) ,add(i + n ,b ,INF);}for(int i = 1 ;i <= n2 ;i ++)add(L[i].a ,L[i].b ,INF);for(int i = 1 ;i <= n3 ;i ++)if(mk[i]) add(G[i].a ,G[i].b ,INF);else add(G[i].a ,G[i].b ,1);for(int i = 1 ;i <= n ;i ++)add(0 ,i ,pp[i]); }int Min_C ,Max_F; void DFS(int s ,int t ,int ii ,int nowcost ,int n) { if(ii == n3 + 1){Rebuid(n);int now = DINIC(s ,t ,t);if(Max_F < now || Max_F == now && Min_C > nowcost){Max_F = now;Min_C = nowcost;}return ;}mk[ii] = 1;DFS(s ,t ,ii + 1 ,nowcost + G[ii].c ,n);mk[ii] = 0;DFS(s ,t ,ii + 1 ,nowcost ,n); }int main () {int n ,m ,i ,q;while(~scanf("%d %d" ,&n ,&m)){for(i = 1 ;i <= n ;i ++)scanf("%d" ,&pp[i]);n1 = n2 = n3 = 0;NODE A;for(i = 1 ;i <= m ;i ++){scanf("%d %d %d %d" ,&A.a ,&A.b ,&A.c ,&q);if(q < 0) S[++n1] = A;if(!q) L[++n2] = A;if(q > 0) G[++n3] = A;}Min_C = Max_F = 0;DFS(0 ,n + n1 + 1 ,1 ,0 ,n);if(!Max_F) puts("Poor Heaven Empire");else printf("%d %d\n" ,Max_F ,Min_C);}return 0; }
總結
以上是生活随笔為你收集整理的hdu 4309 最大流 + DFS的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: hdu1816 + POJ 2723开锁
- 下一篇: hdu4115 2sat 石头剪刀布