深海机器人（cogs 742）

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

标为 (Q,P)。

给定每个深海机器人的出发位置和目标位置，以及每条网格边上生物标本的价值。计算
深海机器人的最优移动方案，使深海机器人到达目的地后，采集到的生物标本的总价值最高。
«数据输入：
由文件shinkai.in提供输入数据。文件的第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)位置坐标作为目的地。
«结果输出:
程序运行结束时，将采集到的生物标本的最高总价值输出到文件shinkai.out中。

shinkai.in
1 1
2 2
1 2
3 4
5 6
7 2
8 10
9 3
2 0 0

2 2 2

shinkai.out

42

1<=P,Q<=15 1<=a,b<=10

/*
    费用流。
    因为每条边可以随便走，并且只能获利一次，所以把每条边拆成两条边，
    一条容量为1，费用为w，另一条容量无限，费用为0。
    注意体面中的x,y坐标比较坑。
*/
#include<iostream>
#include<cstdio>
#include<queue>
#define N 1010
#define inf 1000000000
using namespace std;
int head[N],dis[N],inq[N],fa[N],a,b,P,Q,S,T,cnt=;
struct node{int v,f,w,pre;}e[N*];
queue<int> q;
void add(int u,int v,int f,int w){
    e[++cnt].v=v;e[cnt].f=f;e[cnt].w=w;e[cnt].pre=head[u];head[u]=cnt;
    e[++cnt].v=u;e[cnt].f=;e[cnt].w=-w;e[cnt].pre=head[v];head[v]=cnt;
}
bool spfa(){
    for(int i=;i<=T;i++) dis[i]=inf;
    dis[S]=;q.push(S);
    while(!q.empty()){
        int u=q.front();q.pop();inq[u]=;
        for(int i=head[u];i;i=e[i].pre)
            if(e[i].f&&dis[e[i].v]>dis[u]+e[i].w){
                dis[e[i].v]=dis[u]+e[i].w;
                fa[e[i].v]=i;
                if(!inq[e[i].v]){
                    inq[e[i].v]=;
                    q.push(e[i].v);
                }
            }
    }
    return dis[T]!=inf;
}
int updata(){
    int i=fa[T],x=inf;
    while(i){
        x=min(x,e[i].f);
        i=fa[e[i^].v];
    }
    i=fa[T];
    while(i){
        e[i].f-=x;
        e[i^].f+=x;
        i=fa[e[i^].v];
    }
    return x*dis[T];
}
int id(int x,int y){return (x-)*Q+y;}
int main(){
    freopen("shinkai.in","r",stdin);
    freopen("shinkai.out","w",stdout);
    scanf("%d%d%d%d",&a,&b,&P,&Q);P++;Q++;
    S=;T=P*Q+;
    for(int i=;i<=P;i++)
        for(int j=;j<Q;j++){
            int x;scanf("%d",&x);
            add(id(i,j),id(i,j+),,-x);
            add(id(i,j),id(i,j+),inf,);
        }
    for(int j=;j<=Q;j++)
        for(int i=;i<P;i++){
            int x;scanf("%d",&x);
            add(id(i,j),id(i+,j),,-x);
            add(id(i,j),id(i+,j),inf,);

        }
    for(int i=;i<=a;i++){
        int k,x,y;scanf("%d%d%d",&k,&x,&y);
        add(S,id(x+,y+),k,);
    }
    for(int i=;i<=b;i++){
        int k,x,y;scanf("%d%d%d",&k,&x,&y);
        add(id(x+,y+),T,k,);
    }
    int minv=;
    while(spfa())
        minv+=updata();
    printf("%d",-minv);
    return ;
}