poj 3308 （最大流 最小点权覆盖）

dinic递归法调了一下午最后TLE。。。

二分图最小点权覆盖：

从x或者y集合中选取一些点，使这些点覆盖所有的边，并且选出来的点的权值尽可能小。

行为一个集合，列为一个集合。因为要乘积，所以用对数来算，就变成加法了。也就转换成求最小割。用最大流来做。

#include<stdio.h>
#include<string.h>
#include<cmath>
#include<iostream>
#include<queue>
#define MX 100000
using namespace std;
double map[200][200],f[200][200];
int d[1000],m,n,l;
bool bfs(){
	queue <int> Q;
	Q.push(0);
	memset(d,-1,sizeof(d));
	d[0]=0;
	while(!Q.empty()){
		int v=Q.front();
		Q.pop();
		for(int i=0;i<=n+m+1;i++){
			if(d[i]==-1&&map[v][i]-f[v][i]>0){
				d[i]=v;
				Q.push(i);
				if(i==n+m+1)return true;
			}
		}
	}
	return false;
}
double dinic(){
	double min_flow=0,min_f;
	while(bfs()){
		min_f=MX;
		for(int i=n+m+1,j=d[i];j!=i;){
			if(map[j][i]-f[j][i]<min_f){
				min_f=map[j][i]-f[j][i];
			}
			i=j;
			j=d[j];
		}
//		printf("!! %lf\n",min_f);
		for(int i=n+m+1,j=d[i];j!=i;){
			f[j][i]+=min_f;
			f[i][j]-=min_f;
			i=j;
			j=d[j];
		}
		min_flow+=min_f;
//		printf("! %lf\n",min_flow);
	}
	return min_flow;
}
int main(){
	int cas;
	scanf("%d",&cas);
	while(cas--){
		memset(map,0,sizeof(map));
		memset(f,0,sizeof(f));
		scanf("%d%d%d",&m,&n,&l);
		for(int i=1;i<=m;i++){
			double a;
			scanf("%lf",&a);
			map[0][i]=log(a);
		}
		for(int i=1;i<=n;i++){
			double a;
			scanf("%lf",&a);
			map[i+m][m+n+1]=log(a);
		}
		for(int i=1;i<=l;i++){
			int x,y;
			scanf("%d%d",&x,&y);
			map[x][y+m]=MX;
		}
		double ans=0;
		ans=dinic();
		printf("%.4lf\n",exp(ans));
	}
	return 0;
}