pku 3308 Paratroopers(带权二分图最小点覆盖)

把伞兵看成边，行列看成节点，转化为了带权二分图最小点覆盖。

加入超级源点和超级汇点，源点和所有行节点相连，所有列节点和汇点相连，如果a行b列有敌人，则把节点a和节点b相连。则问题又可以转化求最小割。

现在求源点和汇点之间的最小割。因为对任一敌人<a,b>，必然有source-->a-->b-->sink， 割的性质是"不存在一条从source到sink的路径", 故路径上的三条边<source,a>, <a,b>, <b,sink>中至少有一条边在割中,我们把<a,b>的权值设置为无限大，则其不可能被选中。于是割边集中必然有<source,a>和<b,sink>中的至少一条，也即对应选择了相应的行或列，我们把这些边的权值设置为花费，则最小割即是总花费的最小方案。

还有一个需要注意的地方是，这里问题是要求cost的乘积，可以通过使用log()把乘法先转换为加法，最后输出的时候再转换回去。

#include <iostream> #include "float.h" #include "math.h" using namespace std; #define MinN(x,y) ((x)<(y)?(x):(y)) #define Zero(x) (fabs(x)<EPS) const int maxn=105; const double INF=DBL_MAX; const double EPS=1e-10; double map[maxn][maxn]; double flow[maxn][maxn]; int pre[maxn]; double min_flow[maxn]; int stack[maxn]; bool visited[maxn]; bool Next_Path(int point_num,int source,int sink)//DFS { memset(pre,-1,sizeof(pre)); memset(visited,0,sizeof(visited)); int top=-1; stack[++top]=source; visited[source]=true; min_flow[source]=INF; while(top>=0) { int temp=stack[top--]; for(int i=0;i<point_num;i++) { if(!visited[i]&&!Zero(map[temp][i]-flow[temp][i])) { min_flow[i]=MinN(min_flow[temp],map[temp][i]-flow[temp][i]); pre[i]=temp; if(i==sink) return true; stack[++top]=i; visited[i]=true; } } } return false; } double Max_Flow(int point_num,int source,int sink) { double max_flow=0; memset(flow,0,sizeof(flow)); while(Next_Path(point_num,source,sink)) { int t=sink; while(pre[t]>=0) { flow[pre[t]][t]+=min_flow[sink]; flow[t][pre[t]]-=min_flow[sink]; t=pre[t]; } max_flow+=min_flow[sink]; } return max_flow; } int main() { int T,row,column,n; double cost; scanf("%d",&T); while(T--) { memset(map,0,sizeof(map)); scanf("%d%d%d",&row,&column,&n); for(int i=0;i<row;i++) { scanf("%lf",&cost); map[row+column][i]=log(cost); } for(int i=0;i<column;i++) { scanf("%lf",&cost); map[i+row][row+column+1]=log(cost); } for(int i=0;i<n;i++) { int a,b; scanf("%d%d",&a,&b); map[a-1][row+b-1]=INF; } double ans=Max_Flow(row+column+2,row+column,row+column+1); printf("%.4f/n",exp(ans)); } return 0; }