結果
問題 | No.957 植林 |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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; }