poj 3308

题意：一个矩阵，有的方格里有一些敌人，每行或每列都可以安放一架枪，每架枪都有一个花费，而且能消灭他所在的行或列的所有敌人，最后的花费为所有的枪花费的乘积。

乘积问题可以先取对数。转化成求和。

然后问题很像以前的一个二分匹配问题。把每一行加入到X集合。每一列加入到Y集合。则每一个外星人所在的（row,col）就连条边。求最小点覆盖。但是这个问题有权值。可以转化成KM或者最小割来求解。

最小割建图：

src 与每一行相连，权值为每行的花费。

en 与每一列相连,权值为每列花费。

假如外星人出现在(r,c)则连一条r -->c权值为inf的边。

// File Name: 3308.cpp
// Author: Missa
// Created Time: 2013/4/17 星期三 11:24:48

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<cstdlib>
#include<map>
#include<set>
using namespace std;
#define CL(x,v) memset(x,v,sizeof(x));
#define R(i,st,en) for(int i=st;i<en;++i)
#define LL long long

const double eps = 1e-6;
const int inf = 0x3f3f3f3f;
const int maxn = 100+5;
const int maxm = 15000;
struct Edge
{
    int v, next;
    double c;
}p[maxm << 1];
int head[maxn], e;
int d[maxn], cur[maxn];
int n, m, st, en;
void init()
{
    e = 0;
    //n = en;//记住n赋值为点的个数
    memset(head, -1, sizeof(head));
}
void addEdge(int u, int v, double c)
{
    p[e].v = v; p[e].c = c;
    p[e].next = head[u]; head[u] = e++;
    swap(u,v);
    p[e].v = v; p[e].c = 0;
    p[e].next = head[u]; head[u] = e++;
}
int bfs(int st, int en)
{
    queue <int > q;
    memset(d, 0, sizeof(d));
    d[st] = 1;
    q.push(st);
    while (!q.empty())
    {
        int u = q.front();q.pop();
        for (int i = head[u]; i != -1; i = p[i].next)
        {
            if (p[i].c > 0 && !d[p[i].v])
            {
                d[p[i].v] = d[u] + 1;
                q.push(p[i].v);
            }
        }
    }
    return d[en];
}
double dfs(int u, double a)
{
    if (u == en || fabs(a) <= eps) return a;
    double f, flow = 0;
    for (int& i = cur[u]; i != -1; i = p[i].next)
    {
        if (d[u] + 1 == d[p[i].v] && (f = dfs(p[i].v, min(a, p[i].c))) > 0)
        {
            p[i].c -= f;
            p[i^1].c += f;
            flow += f;
            a -= f;
            if (fabs(a) <= eps) break;
        }
    }
    return flow;
}
double dinic(int st, int en)
{
    double ret = 0, tmp;
    while (bfs(st, en))
    {
        for (int i = 0; i <= n; ++i)
            cur[i] = head[i];
        ret += dfs(st, inf);
    }
    return ret;
}
int row, col, le;
int main()
{
    int cas;
    double c;
    scanf("%d", &cas);
    while (cas--)
    {
        init();
        scanf("%d%d%d", &row, &col, &le);
        st = 0, en = row + col + 1;
        n = en;
        for (int i = 1; i <= row; ++i)
        {
            scanf("%lf", &c);
            addEdge(st, i, log(c));
        }
        for (int i = 1; i <= col; ++i)
        {
            scanf("%lf", &c);
            addEdge(i + row, en, log(c));
        }
        while (le--)
        {
            int u, v;
            scanf("%d%d",&u, &v);
            addEdge(u, v + row, inf);
        }
        printf("%.4f\n",exp(dinic(st, en)));
    }
    return 0;
}