网络流二十四题 ————(二十）、P4012 深海机器人问题 费用流经典建模

20分钟AC，这速度还可以。

做了20道网络流题目后，基础建模已经没啥问题了。

首先还是分析题目，每条路径可以多次访问，但只有一次获得奖励。

多起点，多终点，显然网络流，S连向起点，终点连向T。

而且每次只能往东北走，也就是图是个DAG。则可以跑网络流。

由于每条边只有一次获得奖励的机会，则连一条流量1，花费奖励的边，再连一条流量inf，花费0的边，表示后面再访问这条边无奖励。

建完后跑最大费用流即可。

#include<bits/stdc++.h> 
using namespace std;
typedef long long ll;
const int N=1000005;
const int M=1000005;
const ll inf=0x7fffffff;
struct MCMF{
	struct EDGE {
		int nxt, to, cost, flow;
	}ee[M << 1];
	int n,s,t;
	int cur[N],head[N],vis[N];
	ll dis[N];
	ll mincost=0,maxflow=0;
	inline void AD(int from, int to, int flow, int cost) {
		ee[++cnt].nxt = head[from];
		ee[cnt].to = to;
		ee[cnt].cost = cost;
		ee[cnt].flow = flow;
		head[from] = cnt;
	}
	int cnt=1;
	inline void add(int u,int v,int flow,int cost){
	    AD(u,v,flow,cost);
	    AD(v,u,0,-cost);
	 //   cout<<u<<"  ->  "<<v<<endl;
	}
	inline int spfa()
	{
	    for(int i=1;i<=n;++i)dis[i]=inf,vis[i]=0;
	    queue<int>q;q.push(s);
	    dis[s]=0;vis[s]=1;
	    while(!q.empty()){
	        int u=q.front();
	        vis[u]=0;q.pop();
	        for(int i=head[u];i;i=ee[i].nxt)
	        {
	            int v=ee[i].to;
	            if(!ee[i].flow)continue;
	            if(ee[i].flow&&dis[v]>dis[u]+ee[i].cost){
	                dis[v]=dis[u]+ee[i].cost;
	                if(!vis[v]){
	                    vis[v]=1;
	                    q.push(v);
	                }
	            }
	        }
	    }
	    return dis[t]!=inf;
	}
	inline ll dfs(int u,ll flow)
	{
	    if(u==t)return flow;
	    ll rest=flow;
	    vis[u]=1;
	    for(int i=cur[u];i&&rest;i=ee[i].nxt){
	        cur[u]=i;
	        int v=ee[i].to;
	        ll w=ee[i].cost,flow=ee[i].flow;
	        if(ee[i].flow&&(dis[v]==dis[u]+w)&&!vis[v]){
	            ll k=dfs(v,min(flow,rest));
	            if(k){
	                mincost+=k*w;
	                ee[i].flow-=k;
	                ee[i^1].flow+=k;
	                rest-=k;
	            }
	        }
	    }
	    vis[u]=0;
	    return flow-rest;
	}
	void gao(){
		while(spfa()){
	        for(int i=1;i<=n;++i)cur[i]=head[i];
	        maxflow+=dfs(s,inf);
	    }
	}
	void init(int nn,int ss,int tt){
		n=nn,s=ss,t=tt;
		maxflow=mincost=0;
		for(int i=0;i<=n;i++)head[i]=0;
	}
}mc;
void path(int p,int t){
	while(p!=t){
		cout<<p<<" -> ";
		for(int i=mc.head[p];i;i=mc.ee[i].nxt){
			cout<<p<<"  -   "<<mc.ee[i].to<<"     "<<mc.ee[i].flow<<endl;
			if(i&1)continue;
			if(mc.ee[i^1].flow){
				mc.ee[i^1].flow--;
				p=mc.ee[i].to;
		//		cout<<p<<" == ";
				break;
			}
		}
	}
	cout<<endl;
}
int a,b,s,t,nm,n,m,k,x,y;
int id(int x,int y){
	return x*(m+1)+y+1;
}
int main()
{
	ios::sync_with_stdio(false);
  	cin.tie(0);
  	
  	cin>>a>>b;
  	cin>>n>>m;
  	nm=(n+1)*(m+1)+2;
  	s=nm-1,t=nm;
  	mc.init(nm,s,t);
  	for(int i=0;i<=n;i++){//n行
  		for(int j=1;j<=m;j++){
  			cin>>x;
  			int u=id(i,j-1),v=id(i,j);
  			mc.add(u,v,1,-x);
  			mc.add(u,v,inf,0);
		  }
	}
	for(int j=0;j<=m;j++){//m列 
  		for(int i=1;i<=n;i++){
  			cin>>x;
  			int u=id(i-1,j),v=id(i,j);
  			mc.add(u,v,1,-x);
  			mc.add(u,v,inf,0);
		  }
	}
	for(int i=1;i<=a;i++){
		cin>>k>>x>>y;
		mc.add(s,id(x,y),k,0);
	}
	for(int i=1;i<=b;i++){
		cin>>k>>x>>y;
		mc.add(id(x,y),t,k,0);
	}
//	cout<<"ok"<<endl;
	mc.gao();
//	path(s,t);
	cout<<-mc.mincost<<endl;
	return 0;
}