海底机器人问题（最小费用流）

Description

    深海资源考察探险队的潜艇将到达深海的海底进行科学考察。潜艇内有多个深海机器人。潜艇到达深海海底后，深海机器人将离开潜艇向预定目标移动。深海机器人在移动中还必须沿途采集海底生物标本。沿途生物标本由最先遇到它的深海机器人完成采集。每条预定路径上的生物标本的价值是已知的，而且生物标本只能被采集一次。本题限定深海机器人只能从其出发位置沿着向北或向东的方向移动，而且多个深海机器人可以在同一时间占据同一
位置。
    用一个P′Q 网格表示深海机器人的可移动位置。西南角的坐标为（0,0），东北角的坐标为 (Q,P)。

    给定每个深海机器人的出发位置和目标位置，以及每条网格边上生物标本的价值。计算深海机器人的最优移动方案，使深海机器人到达目的地后，采集到的生物标本的总价值最高。

Input

多组数据输入.
每组输入第1 行为深海机器人的出发位置数a，和目的地数b，第2 行为P和Q 的值。接下来的P+1 行，每行有Q 个正整数，表示向东移动路径上生物标本的价值，行数据依从南到北方向排列。再接下来的Q+1 行，每行有P 个正整数，表示向北移动路径上生物标本的价值，行数据依从西到东方向排列。接下来的a行，每行有
3 个正整数k,x,y，表示有k个深海机器人从(x,y)位置坐标出发。再接下来的b行，每行有3个正整数r,x,y，表示有r个深海机器人可选择(x,y)位置坐标作为目的地。

Output

每组输出采集到的生物标本的最高总价值

Sample Input

1 1
2 2
1 2
3 4
5 6
7 2
8 10
9 3
2 0 0
2 2 2

Sample Output

42

题目出自nefu494

思路：这题乍一看以为是最大流，直接套模板一顿建图，直接调试数据都过不了，然后想想发现的确不能用最大流直接建图，因为源点和汇点的边根本没法设置权值，然后改用最小费用最大流来建图了，每条矩阵上的边流量为1，表示只能有一次流过去（我太天真），费用是输入的数据的相反数，然后源点连起点，流量是输入数据，费用0，汇点和源点连法相同，结果wa了……

看了看题解发现原先矩阵的每条边少连了一条边，流量为oo，费用是0的边……

【建模方法】

把网格中每个位置抽象成网络中一个节点，建立附加源S汇T。

1、对于每个顶点i，j为i东边或南边相邻的一个节点，连接节点i与节点j一条容量为1，费用为该边价值的有向边。
2、对于每个顶点i，j为i东边或南边相邻的一个节点，连接节点i与节点j一条容量为无穷大，费用为0的有向边。
3、从S到每个出发点i连接一条容量为该点出发的机器人数量，费用为0的有向边。
4、从每个目标点i到T连接一条容量为可以到达该点的机器人数量，费用为0的有向边。

求最大费用最大流，最大费用流值就采集到的生物标本的最高总价值。

【建模分析】

这个问题可以看做是多出发点和目的地的网络运输问题。每条边的价值只能计算一次，容量限制要设为1。同时还将要连接上容量不限，费用为0的重边。由于“多个深海机器人可以在同一时间占据同一位置”，所以不需限制点的流量，直接求费用流即可。

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
using namespace std;
typedef long long ll;
const   int oo=1e9;
const   int mm=11111111;
const   int mn=888888;
int node,src,dest,edge;
int ver[mm],flow[mm],cost[mm],nex[mm];
int head[mn],dis[mn],p[mn],q[mn],vis[mn];
/**这些变量基本与最大流相同，增加了
 cost 表示边的费用，
 p 记录可行流上节点对应的反向边
 */
void prepare(int _node,int _src,int _dest)
{
    node=_node,src=_src,dest=_dest;
    for(int i=0; i<node; i++)head[i]=-1,vis[i]=0;
    edge=0;
}
void addedge(int u,int v,int f,int c)
{
    ver[edge]=v,flow[edge]=f,cost[edge]=c,nex[edge]=head[u],head[u]=edge++;
    ver[edge]=u,flow[edge]=0,cost[edge]=-c,nex[edge]=head[v],head[v]=edge++;
}
/**以上同最大流*/
/**spfa 求最短路，并用 p 记录最短路上的边*/
bool spfa()
{
    int i,u,v,l,r=0,tmp;
    for(i=0; i<node; ++i)dis[i]=oo;
    dis[q[r++]=src]=0;
    p[src]=p[dest]=-1;
    for(l=0; l!=r; (++l>=mn)?l=0:l)
        for(i=head[u=q[l]],vis[u]=0; i>=0; i=nex[i])
            if(flow[i]&&dis[v=ver[i]]>(tmp=dis[u]+cost[i]))
            {
                dis[v]=tmp;
                p[v]=i^1;
                if(vis[v]) continue;
                vis[q[r++]=v]=1;
                if(r>=mn)r=0;
            }
    return p[dest]>-1;
}
/**源点到汇点的一条最短路即可行流，不断的找这样的可行流*/
int SpfaFlow()
{
    int i,ret=0,delta;
    while(spfa())
    {
        for(i=p[dest],delta=oo; i>=0; i=p[ver[i]])
            if(flow[i^1]<delta)delta=flow[i^1];
        for(i=p[dest]; i>=0; i=p[ver[i]])
            flow[i]+=delta,flow[i^1]-=delta;
        ret+=delta*dis[dest];
    }
    return ret;
}

int main()
{
    int a,b,p,q;
    while(~scanf("%d%d",&a,&b))
    {
        scanf("%d%d",&p,&q);
        int sum=(p+1)*(q+1);
        prepare(sum+2,sum,sum+1);
        for(int i=0; i<=p; i++)
            for(int j=0; j<q; j++)
            {
                int tt;
                scanf("%d",&tt);
                addedge(i*(q+1)+j,i*(q+1)+j+1,1,-tt);
                addedge(i*(q+1)+j,i*(q+1)+j+1,oo,0);
            }
        for(int i=0; i<=q; i++)
            for(int j=0; j<p; j++)
            {
                int tt;
                scanf("%d",&tt);
                addedge(j*(q+1)+i,(j+1)*(q+1)+i,1,-tt);
                addedge(j*(q+1)+i,(j+1)*(q+1)+i,oo,0);
            }
        for(int i=0; i<a; i++)
        {
            int r,x,y;
            scanf("%d%d%d",&r,&x,&y);
            addedge(sum,x*(q+1)+y,r,0);
        }
        for(int i=0; i<b; i++)
        {
            int r,x,y;
            scanf("%d%d%d",&r,&x,&y);
            addedge(x*(q+1)+y,sum+1,r,0);
        }
        printf("%d\n",-SpfaFlow());
    }
    return 0;
}