HDU 3395 Special Fish[最大流最大費用]
阿新 • • 發佈:2018-12-16
題意:有n個點,每個點可以選擇攻擊其他的點,且攻擊次數不限,每次攻擊可獲得的價值為v[i]^v[j]。每個點只能被攻擊一次,求可獲得的最大價值
思路:
設影子節點(n+1)~2n,如下建邊
(1,i,1,0) i∈[1,n]
(i,j+n,1,-(v[i]^v[j])) if(mp[i][j]==1)
(i,sink,1,0) 每條魚可以不攻擊其他魚,因為可能把攻擊的機會讓給別人反而獲益更多
(j,sink,1,0)
#include <cstdio> #include <cstring> #include <queue> #include <stack> #include <vector> #include <cmath> #include <cstdlib> #include <algorithm> #include <iostream> #define mAXn 500+50 #define mAXm 20000+5 #define InF 0x3f3f3f3f using namespace std; struct Edge { int from, to, cap, flow, cost, next; }; Edge edge[mAXm]; int head[mAXn], edgenum; int pre[mAXn];//記錄增廣路徑上 到達點i的邊的編號 int dist[mAXn]; bool vis[mAXn]; int n, m;//點數 邊數 int source, sink;//超級源點 超級匯點 void init() { edgenum = 0; memset(head, -1, sizeof(head)); } void addEdge(int u, int v, int w, int c) { Edge E1 = {u, v, w, 0, c, head[u]}; edge[edgenum] = E1; head[u] = edgenum++; Edge E2 = {v, u, 0, 0, -c, head[v]}; edge[edgenum] = E2; head[v] = edgenum++; } bool SPFA(int s, int t)//尋找花銷最少的路徑 { //跑一遍SPFA 找s——t的最少花銷路徑 且該路徑上每一條邊不能滿流 //若存在 說明可以繼續增廣,反之不能 queue<int> Q; memset(dist, InF, sizeof(dist)); memset(vis, false, sizeof(vis)); memset(pre, -1, sizeof(pre)); dist[s] = 0; vis[s] = true; Q.push(s); while(!Q.empty()) { int u = Q.front(); Q.pop(); vis[u] = false; for(int i = head[u]; i != -1; i = edge[i].next) { Edge E = edge[i]; if(dist[E.to] > dist[u] + E.cost && E.cap > E.flow)//可以鬆弛 且 沒有滿流 { dist[E.to] = dist[u] + E.cost; pre[E.to] = i;//記錄前驅邊 的編號 if(!vis[E.to]) { vis[E.to] = true; Q.push(E.to); } } } } return pre[t] != -1;//可達返回true } void mCmF(int s, int t, int &cost, int &flow) { flow = 0;//總流量 cost = 0;//總費用 while(SPFA(s, t))//每次尋找花銷最小的路徑 { int Min = InF; //通過反向弧 在源點到匯點的最少花費路徑 找最小增廣流 for(int i = pre[t]; i != -1; i = pre[edge[i^1].to]) { Edge E = edge[i]; Min = min(Min, E.cap - E.flow); } //增廣 for(int i = pre[t]; i != -1; i = pre[edge[i^1].to]) { edge[i].flow += Min; edge[i^1].flow -= Min; cost += edge[i].cost * Min;//增廣流的花銷 } flow += Min;//總流量累加 } } char mp[105][105]; int a[105][105]; int v[105]; void getmap(){ for(int i=1;i<=n;i++) scanf("%d",&v[i]); for(int i=1;i<=n;i++) scanf("%s",mp[i]+1); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) a[i][j]=mp[i][j]-'0'; source=0,sink=2*n+1; for(int i=1;i<=n;i++) addEdge(source,i,1,0); for(int i=1;i<=n;i++) addEdge(i,sink,1,0); for(int i=1;i<=n;i++){ for(int j=1;j<=n;j++){ if(!a[i][j]||i==j) continue; addEdge(i,j+n,1,-(v[i]^v[j])); } } for(int i=n+1;i<=2*n;i++) addEdge(i,sink,1,0); return ; } int main() { while(scanf("%d", &n)==1&&n){ init(); getmap();//建圖 int cost, flow;//最小費用 最大流 mCmF(source, sink, cost, flow); printf("%d\n", -cost);//最小費用 最大流 // printf("%d\n", flow);//最小費用 最大流 } return 0; } /* 5 22 33 44 55 99 10001 11111 11110 01011 11110 */