#zkw费用流，最小费用最大流#洛谷 4012 codevs 1917 ssl 2620 深海机器人问题

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/sugar_free_mint/article/details/85011326

题目大意

在一个平面直角坐标系中，机器人只能往右和上采集标本，每个格点都有不同的价值，现在若干个机器人从某点出发目的地为某点，问采集到的最大价值

---

分析

其实这道题类比于K取方格数，容易建出这样一张图

然后跑一遍最大费用最大流就可以了，但是我把费用取反，跑的是最小费用最大流

---

代码

#include <cstdio>
#include <deque>
#include <cstring>
#define rr register
#define id(x,y) ((x-1)*m+y)
using namespace std;
const int inf=707406378;
struct node{
    int y,w,f,next;
}e[3001];
int ss,tt,n,m,s,t,ans,k=1,ls[301],dis[301]; bool v[301];
inline void add(int x,int y,int w,int f){
    e[++k]=(node){y,w,f,ls[x]}; ls[x]=k;
    e[++k]=(node){x,0,-f,ls[y]}; ls[y]=k;
}
inline signed spfa(){
    memset(v,0,sizeof(v)); memset(dis,127/3,sizeof(dis));
    dis[t]=0; v[t]=1; rr deque<int>q; q.push_back(t);
    while (q.size()){
        rr int x=q.front(); q.pop_front();
        for (rr int i=ls[x];i;i=e[i].next)
        if (e[i^1].w&&dis[e[i].y]>dis[x]-e[i].f){
            dis[e[i].y]=dis[x]-e[i].f;
            if (!v[e[i].y]){
                v[e[i].y]=1;
                if (q.size()&&dis[e[i].y]<dis[q.front()])
                q.push_front(e[i].y); else q.push_back(e[i].y);
            }
        }
        v[x]=0;
    }
    return dis[s]<707406378;
}
inline signed dfs(int x,int now){
    if (x==t) {v[t]=1; return now;}
    rr int rest=0,f; v[x]=1;
    for (rr int i=ls[x];i;i=e[i].next)
    if (!v[e[i].y]&&e[i].w&&dis[e[i].y]+e[i].f==dis[x]){
        rest+=(f=dfs(e[i].y,min(e[i].w,now-rest)));
        if (f) ans-=f*e[i].f,e[i].w-=f,e[i^1].w+=f;
        if (rest==now) break;
    }
    return rest;
}
inline void answ(){
    while (spfa()){
        v[t]=1;
        while (v[t]){
            memset(v,0,sizeof(v));
            dfs(s,1e9);
        }
    }
}
signed main(){
    scanf("%d%d%d%d",&ss,&tt,&n,&m); s=++n*++m+1,t=s+1;
    for (rr int i=1;i<=n;++i)
    for (rr int j=1,x;j<m;++j){
        scanf("%d",&x);
        add(id(i,j),id(i,j+1),1,-x);
        add(id(i,j),id(i,j+1),inf,0);
    }
    for (rr int j=1;j<=m;++j)
    for (rr int i=1,x;i<n;++i){
        scanf("%d",&x);
        add(id(i,j),id(i+1,j),1,-x);
        add(id(i,j),id(i+1,j),inf,0);		
    }
    for (rr int i=1,x,y,w;i<=ss;++i){
    	scanf("%d%d%d",&w,&x,&y);
    	add(s,id(++x,++y),w,0);
    }
    for (rr int i=1,x,y,w;i<=tt;++i){
    	scanf("%d%d%d",&w,&x,&y);
    	add(id(++x,++y),t,w,0);
    }
    answ();
    printf("%d",ans);
    return 0;
}