POJ 3308 Paratroopers

 

类比二分图匹配的最大流解法, 加入源点s, 汇点t, 将匹配以一条条s->u->v->t形式的流路径"串联"起来. 匹配的限制是在点上, 即要求选择一些边, 限制你边两边的点要满足某某条件, 利用流的容量限制, 以满足这些条件(比如一个点只能连在一条边上), 类似的, [PKU][2699][The Maximum Number of Strong Kings]中"一场比赛有两个人参加, 有且只有一个人获胜"的限制也是利用了流的限制. 而点覆盖集的限制在边上, 最小割是最大流的对偶问题, 对偶往往是将问题的性质从点转边, 从边转点.

基于以上的动机, 从s向每个横坐标连边, 从每个纵坐标向t连边, 权值为花费的对数, 以每个伞兵位置的横纵坐标值连一条从横坐标到纵坐标的有向边, 权值为inf. 则在此有向图上任意一条从s到t的路径一定具有s->u->v->t的形式. 割的性质是"不存在一条从s到t的路径", 故路径上的三条边<s,u>, <u,v>, <v,t>中至少有一条边在割中, 正好与点覆盖集限制条件形式相符(对于每条边(u,v), u,v至少有一个被选中). 而目标是最小化点权之和, 恰好也是最小割的优化目标.

详细解释参见Amber的论文[最小割模型在信息学竞赛中的应用].

注意本题求的是最小化乘积, 要先把边权取对数.

 

转自：http://www.answeror.com/archives/27505

 

代码:

#include<iostream> #include<cmath> using namespace std; const double inf=500.0; const int MAX=105; struct node { int v,next; double c; }g[100000]; int num[MAX],dis[MAX],cur[MAX],adj[MAX],pre[MAX]; int s,t,e,vn,r,c,l; void add(int u,int v,double c) { g[e].v=v; g[e].c=c; g[e].next=adj[u]; adj[u]=e++; g[e].v=u; g[e].c=0; g[e].next=adj[v]; adj[v]=e++; } double sap() { int i,u,v,flag; double aug=inf,flow=0.0; for(i=0;i<=vn;i++) { dis[i]=num[i]=0; cur[i]=adj[i]; } u=s; num[0]=vn; while(dis[s]<vn) { flag=0; for(i=cur[u];i!=-1;i=g[i].next) { v=g[i].v; if(g[i].c>0.0&&dis[u]==dis[v]+1) { flag=1; cur[u]=i; pre[v]=u; aug=min(aug,g[i].c); u=v; if(u==t) { flow+=aug; //cout<<flow<<endl; while(u!=s) { u=pre[u]; g[cur[u]].c-=aug; g[cur[u]^1].c+=aug; } aug=inf; } break; } } if(flag) continue; if(--num[dis[u]]==0) return flow; for(dis[u]=vn,i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(g[i].c>0.0&&dis[v]<dis[u]) { dis[u]=dis[v]; cur[u]=i; } } dis[u]++; num[dis[u]]++; if(u!=s) u=pre[u]; } return flow; } int main() { int i,j,test,a,b; double cost; scanf("%d",&test); while(test--) { memset(adj,-1,sizeof(adj)); e=0; scanf("%d %d %d",&r,&c,&l); s=0; t=r+c+1; vn=t+1; for(i=1;i<=r;i++) { scanf("%lf",&cost); add(s,i,log(cost)); } for(i=1;i<=c;i++) { scanf("%lf",&cost); add(i+r,t,log(cost)); } for(j=1;j<=l;j++) { scanf("%d %d",&a,&b); add(a,r+b,inf); } // cout<<sap()<<endl; printf("%.4lf/n",exp(sap())); } return 0; }