意甲冠軍:
一m*n該網絡的規模格。詳細地點稱為傘兵著陸(行和列)。
現在,在一排(或列)
安裝激光槍,激光槍可以殺死線(或塔)所有傘兵。在第一i安裝一排
費用是Ri。在第i列安裝的費用是Ci。
要安裝整個激光槍系統,總費用為這些
激光槍費用的乘積。
求殺死全部傘兵的最小費用。
構圖:
把傘兵視為邊,行與列視為頂點。添加源點和匯點,對于第i行。從源點向頂點i連接一條
容量為Ri的邊。對于第j列。從頂點j向匯點連接一條容量為Rj的邊。
假設某一點(i,j)有傘兵降落,則從頂點Ri向頂點Cj連接一條容量為無窮大的邊。
算法:
依據割的性質,源點和匯點必不連通。則割邊必然在S->R,R->C,C->T其一。
為了求得最小容量,
將R->C設為無窮大,則其不可能被選中。這樣割邊集為S-->R,C-->T的集合,也就是選中行或列。
此時求得的最小割為花費最小的方案。
因為花費為行和列的乘積。則通過對數運算把乘法轉化為加法。
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<queue>
#define maxm 15000
#define maxn 105
#define eps 1e-6
using namespace std;struct node
{int v,next;double val;
}e[maxm<<1];
int st,en,n,m,l,cnt;
int d[maxn];
int head[maxn],cur[maxn];
const double INF = 1000007;
queue<int> q;void init()
{st = 0,en = n+m+1;memset(head,-1,sizeof(head));cnt = 0;
}
void add(int x,int y,double z)
{e[cnt].v = y;e[cnt].val = z;e[cnt].next = head[x];head[x]=cnt++;e[cnt].v = x;e[cnt].val = 0;e[cnt].next = head[y];head[y]=cnt++;
}
bool bfs()
{while(!q.empty())q.pop();memset(d,-1,sizeof(d));int u;d[st] = 0;q.push(st);while(!q.empty()){u = q.front();q.pop();for(int i=head[u];i!=-1;i=e[i].next){int t = e[i].v;if(e[i].val>0 && d[t]==-1){d[t] = d[u]+1;q.push(t);if(t==en) return true;}}}return false;
}double dfs(int x,double flow)
{if(x==en || fabs(flow)<=eps) return flow;double ret = 0,dd;for(int& i=cur[x];i!=-1;i=e[i].next){int t = e[i].v;if(d[t] == d[x]+1 && (dd = dfs(t,min(flow,e[i].val)))>0){e[i].val-=dd;e[i^1].val+=dd;flow-=dd;ret+=dd;if (fabs(flow) <= eps) break;}}return ret;
}
double Dinic()
{double tmp = 0,maxflow = 0;while(bfs()){for(int i=0;i<=en;i++)cur[i] = head[i];maxflow+=dfs(st,INF);}return maxflow;
}int main()
{int T,a,b;double x;scanf("%d",&T);while(T--){scanf("%d%d%d",&m,&n,&l);init();for(int i=1;i<=m;i++){scanf("%lf",&x);add(st,i,log(x));}for(int i=1;i<=n;i++){scanf("%lf",&x);add(i+m,en,log(x));}for(int i=1;i<=l;i++){scanf("%d%d",&a,&b);add(a,b+m,INF);}printf("%.4f\n",exp(Dinic()));}return 0;
}
版權聲明:本文博主原創文章。博客,未經同意不得轉載。
總結
以上是生活随笔為你收集整理的zoj 2874 amp; poj 3308 Paratroopers (最小割)的全部內容,希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。