結果
| 問題 |
No.654 Air E869120
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-10-01 10:25:07 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 8 ms / 2,000 ms |
| コード長 | 1,830 bytes |
| コンパイル時間 | 3,638 ms |
| コンパイル使用メモリ | 84,848 KB |
| 最終ジャッジ日時 | 2025-01-24 18:48:01 |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 5 |
| other | AC * 35 |
ソースコード
#include <queue>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
struct edge { int to, rev; long long cap; };
const long long inf = 1LL << 55;
int N, M, V, d, a[1009], b[1009], p[1009], q[1009], w[1009], ps[1009]; bool vis[2009]; vector<int> t[1009]; vector<edge> g[2009];
void add_edge(int v1, int v2, long long cost) {
g[v1].push_back(edge{ v2, (int)g[v2].size(), cost });
g[v2].push_back(edge{ v1, (int)g[v1].size() - 1, 0LL });
}
long long find_augment(int pos, int tar, long long curflow) {
if (pos == tar) return curflow;
vis[pos] = true;
for (edge &e : g[pos]) {
if (!vis[e.to] && e.cap > 0) {
long long newflow = find_augment(e.to, tar, min(e.cap, curflow));
if (newflow > 0) {
e.cap -= newflow;
g[e.to][e.rev].cap += newflow;
return newflow;
}
}
}
return 0;
}
long long max_flow(int src, int dst) {
long long ret = 0;
while (true) {
fill(vis, vis + V, false);
long long res = find_augment(src, dst, inf);
if (res == 0) break;
ret += res;
}
return ret;
}
int main() {
cin >> N >> M >> d;
t[0] = { 0 }; t[N - 1] = { 1012345678 };
for (int i = 0; i < M; i++) {
cin >> a[i] >> b[i] >> p[i] >> q[i] >> w[i];
t[--a[i]].push_back(p[i]);
t[--b[i]].push_back(q[i] + d);
}
for (int i = 0; i < N; i++) {
sort(t[i].begin(), t[i].end());
t[i].erase(unique(t[i].begin(), t[i].end()), t[i].end());
ps[i + 1] = ps[i] + t[i].size();
for (int j = ps[i]; j < ps[i + 1] - 1; j++) {
add_edge(j, j + 1, inf);
}
}
for (int i = 0; i < M; i++) {
int pa = lower_bound(t[a[i]].begin(), t[a[i]].end(), p[i]) - t[a[i]].begin() + ps[a[i]];
int pb = lower_bound(t[b[i]].begin(), t[b[i]].end(), q[i] + d) - t[b[i]].begin() + ps[b[i]];
add_edge(pa, pb, w[i]);
}
V = ps[N];
long long ret = max_flow(0, V - 1);
cout << ret << '\n';
return 0;
}