結果
| 問題 |
No.957 植林
|
| コンテスト | |
| ユーザー |
shin3110
|
| 提出日時 | 2023-11-12 18:17:04 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,955 bytes |
| コンパイル時間 | 1,905 ms |
| コンパイル使用メモリ | 182,636 KB |
| 実行使用メモリ | 13,184 KB |
| 最終ジャッジ日時 | 2024-09-26 03:02:39 |
| 合計ジャッジ時間 | 5,544 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | TLE * 1 -- * 44 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define ll long long
struct edge {
int to, seq;
ll cup;
};
class max_flow {
public:
int siz;
vector<vector<edge>> side;
vector<bool> used;
max_flow (int n) {
siz = n;
side.assign(siz, {});
used.assign(siz, false);
return ;
}
void set(int pos, int to, ll x) {
int a = side[to].size();
int b = side[pos].size();
side[pos].push_back(edge{to, a, x});
side[to].push_back(edge{pos, b, 0});
return ;
}
ll dfs(int pos, int goal, ll F) {
if(pos == goal) return F;
used[pos] = true;
for(int i=0; i<side[pos].size(); i++) {
if(side[pos][i].cup == 0) continue;
if(used[side[pos][i].to]) continue;
ll G = dfs(side[pos][i].to, goal, min(F, side[pos][i].cup));
if(G > 0) {
side[pos][i].cup -= G;
side[side[pos][i].to][side[pos][i].seq].cup += G;
return G;
}
}
return 0;
}
ll len(int s, int t) {
ll res = 0;
while(1) {
for(int i=0; i<siz; i++) used[i] = false;
ll F = dfs(s, t, 1e18);
if(F == 0) break;
res += F;
}
return res;
}
};
signed main() {
int h, w; cin >> h >> w;
int ax = h+w+1;
max_flow F(ax+1);
vector<ll> c(h+1, 0);
for(int i=0; i<h*w; i++) {
ll g; cin >> g;
int x = i/w;
c[x+1] += g;
int y = i%w;
F.set(x+1, h+y+1, g);
}
ll ans = 0;
for(int j=1; j<=h; j++) {
ll r; cin >> r;
r -= c[j];
if(r < 0) F.set(0, j, -r);
else {
ans += r;
F.set(j, ax, r);
}
}
for(int i=1; i<=w; i++) {
ll r; cin >> r;
F.set(h+i, ax, r);
ans += r;
}
ans -= F.len(0, ax);
cout << ans << endl;
}
shin3110