poj3308 最大流最小割

本来说昨天把这道题给做出来的，结果昨天去听报告了，而且也没有听懂说的是什么，英语太差啦！！！我发现自己真的是越来越水了，做这道题的时候就一直是循环然后出不来，然后一直调试啊，就是出不来，后来参照了别人的。。。

大致题意：对于一个n*m的矩阵，可以放n+m个机关枪，有l个敌人将要分别降落在矩阵的一些格子中，在每行或者列放置机关枪的话费是不一样的，并且只要在某一行（列）放置一架机关枪，则这一行（列）的敌人都会死，要求最小的话费使得所有的敌人都被机关枪给杀死！

解题思路：最大流最小割，这点我就不解释了，因为我也是看别人的！

就是一个模版，下面是代码

#include<iostream>
#include<cstdio>
#include<string>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<stack>
#include<vector>
#include<climits>
using namespace std;
#define max(a,b) (a)>(b)?(a):(b)
#define min(a,b) (a)<(b)?(a):(b) 
#define rep(i,n) for(i=0; i<(n); i++)
#define reph(i,n,m) for(i=(n); i<=(m); i++)//正循环的
#define repd(i,n,m) for(i=(n); i>=(m); i--) //负循环的
#define max(a,b) (a)>(b)?(a):(b)
#define min(a,b) (a)<(b)?(a):(b)
#define fab(a) (a)>0?(a):0-(a)
#define arc(a) (a)*(a)
#define inf 10000000   //最大值的
#define eps 0.00000001     //浮点型的
#define N   150  //记录开的数组
typedef struct fun
{
	int x,y,pre;
	double e;
}rr;
int que[N];
fun a[10000];
int n,m,l,len,s,t;
int pre[N];
int stage[N];
int cur[N];
void add(int i,int j,double e)
{
	a[len].x=i;
	a[len].y=j;
	a[len].e=e;
	a[len].pre=pre[i];
	pre[i]=len++;


	a[len].x=j;
	a[len].y=i;
	a[len].e=0;
	a[len].pre=pre[j];
	pre[j]=len++;
}
int bfs(int n,int s,int t)
{
	queue<int>q;
	int i,x;
    memset(stage,-1,sizeof(stage));
	stage[s]=0;
	q.push(s);
	while(!q.empty())
	{
		x=q.front();
		q.pop();
		for(i=pre[x]; i!=-1; i=a[i].pre)
		{
			if(a[i].e>eps && stage[a[i].y]==-1)
			{
				stage[a[i].y]=stage[a[i].x]+1;
                q.push(a[i].y);
				if(a[i].y==t)//所有的路都找到了
					return 1;
			}
		}
	}
	return 0;
}
//double dfs(int N, int s, int t)
double dfs(int n,int s,int t)
{
      double maxflow=0,tmp;
	   int i;
	   int u,tail,e; 
      while(bfs(n,s,t)==1)
	   {
		   u=s;
		   tail=0; 
		   while(1)
		   {
		   if(t==u)
		   {
			   tmp=inf;
			   int mark;
			   rep(i,tail)//找到路中容量最小的管子，并且记录的
			    if(a[que[i]].e<tmp)
				{
					mark=i;
					tmp=a[que[i]].e;
				}
				rep(i,tail)
				{
					a[que[i]].e-=tmp;
					a[que[i]^1].e+=tmp;
				}
				maxflow+=tmp;
				tail=mark;
				u=a[que[tail]].x;
		   }
		   for(e=pre[u]; e!=-1; e=a[e].pre)
		   {
			   int y=a[e].y;
			   if(stage[y]==stage[u]+1 && a[e].e>eps)
			   {
				   break;
			   }
		   }
		   if(e!=-1)
		   {
			   que[tail++]=e;
			   u=a[e].y;
		   }
		   else
		   {
			   if(tail==0)//代表没有啦
				   break;
			   stage[a[que[--tail]].y]=-1;
			   u=a[que[tail]].x;
		   }
		   }
	   }
	   return maxflow;
}
int main()
{
	int T,i,j;
	double e;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d%d",&n,&m,&l);
        s=0;
    	t=n+m+1;
		memset(pre,-1,sizeof(pre));
		len=0;
		reph(i,1,n)
		{
			scanf("%lf",&e);
			add(s,i,log(e));
		}
		reph(i,n+1,n+m)
		{
			scanf("%lf",&e);
			add(i,t,log(e));
		}
		while(l--)
		{
			scanf("%d%d",&i,&j);
			add(i,j+n,inf);
		}
		double tmp=dfs(t+1,s,t);
		printf("%0.4lf\n",exp(tmp));
	}
    return 0;
}