結果

問題 No.957 植林
コンテスト
ユーザー chocorusk
提出日時 2019-12-19 22:24:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 760 ms / 2,000 ms
コード長 2,116 bytes
コンパイル時間 1,229 ms
コンパイル使用メモリ 118,848 KB
実行使用メモリ 15,724 KB
最終ジャッジ日時 2024-07-07 01:55:26
合計ジャッジ時間 19,770 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 45
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll INF=1e18;
struct edge {int to; long long int cap; int rev;} ;
 
vector<edge> G[1020];
int level[1020];
int iter[1020];
 
void add_edge(int from, int to, long long int cap){
	edge e;
	e.to=to, e.cap=cap, e.rev=G[to].size();
	G[from].push_back(e);
	e.to=from, e.cap=0, e.rev=G[from].size()-1;
	G[to].push_back(e);
}
 
void bfs(int s){
	memset(level, -1, sizeof(level));
	queue<int> que;
	level[s]=0;
	que.push(s);
	while(!que.empty()){
		int v=que.front();
		que.pop();
		for(int i=0; i<G[v].size(); i++){
			edge e=G[v][i];
			if(e.cap>0 && level[e.to]<0){
				level[e.to]=level[v]+1;
				que.push(e.to);
			}
		}
	}
}
 
long long int dfs(int v, int t, long long int f){
	if(v==t) return f;
	for(int &i=iter[v]; i<G[v].size(); i++){
		edge &e=G[v][i];
		if(e.cap>0 && level[v]<level[e.to]){
			long long int d=dfs(e.to, t, min(f, e.cap));
			if(d>0){
				e.cap-=d;
				G[e.to][e.rev].cap+=d;
				return d;
			}
		}
	}
	return 0;
}
 
long long int max_flow(int s, int t){
	long long int flow=0;
	while(1){
		bfs(s);
		if(level[t]<0) return flow;
		memset(iter, 0, sizeof(iter));
		long long int f;
		while((f=dfs(s, t, INF))>0){
			flow+=f;
		}
	}
}
int main()
{
	int h, w;
  cin>>h>>w;
  int s=h+w, t=s+1;
  ll sum=0;
  for(int i=0; i<h; i++){
    for(int j=0; j<w; j++){
      ll x; cin>>x;
      add_edge(i, t, x);
      add_edge(j+h, i, x);
    }
  }
  for(int i=0; i<h; i++){
    ll x; cin>>x;
    sum+=x;
    add_edge(s, i, x);
  }
  for(int i=0; i<w; i++){
    ll x; cin>>x;
    sum+=x;
    add_edge(s, i+h, x);
  }
  cout<<sum-max_flow(s, t)<<endl;
	return 0;
}
0