結果
| 問題 |
No.654 Air E869120
|
| コンテスト | |
| ユーザー |
dsytk7
|
| 提出日時 | 2018-03-27 11:50:12 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 3,120 bytes |
| コンパイル時間 | 2,155 ms |
| コンパイル使用メモリ | 188,792 KB |
| 実行使用メモリ | 6,948 KB |
| 最終ジャッジ日時 | 2024-06-25 12:41:47 |
| 合計ジャッジ時間 | 6,698 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 WA * 1 RE * 3 |
| other | WA * 12 RE * 23 |
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:93:25: warning: overflow in conversion from 'long long int' to 'std::vector<int>::value_type' {aka 'int'} changes value from '1152921504606846976' to '0' [-Woverflow]
93 | t[N-1].push_back(1LL<<60);
| ~~~^~~~
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for (int (i)=(0);(i)<(int)(n);++(i))
using ll = long long;
using P = pair<int, int>;
using namespace std;
template<class T> void vin(vector<T>& v, int n) {
v.resize(n);
for (int i = 0; i < n; ++i) {
cin >> v[i];
}
}
template<class T>
class Dinic {
private:
const T INF;
struct edge {
T to, cap, rev;
edge(T to, T cap, T rev) :to(to), cap(cap), rev(rev) {}
};
T V;
vector<vector<edge>> g;
vector<T> level; // sourceからの距離
vector<T> itr; // どこまで調べ終わったか
public:
Dinic(T V) :INF(1e9+7), V(V), g(V, vector<edge>()) {}
void add_edge(T from, T to, T cap) {
g[from].emplace_back(to, cap, g[to].size());
g[to].emplace_back(from, 0, g[from].size()-1);
}
// sourceからの最短路を幅優先で探索
void bfs(T s) {
level = vector<T>(V, -1);
queue<T> que;
que.push(s);
level[s] = 0;
while (que.size()) {
int v = que.front(); que.pop();
for (edge &e : g[v]) {
if (e.cap > 0 and level[e.to] < 0) {
level[e.to] = level[v] + 1;
que.push(e.to);
}
}
}
}
// v->tの増加路を探す
T dfs(T v, T t, T f) {
if (v == t) return f;
for (T& i = itr[v]; i < g[v].size(); ++i) {
edge& e = g[v][i];
if (e.cap > 0 and level[v] < level[e.to]) {
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;
}
T max_flow(T s, T t) {
T flow = 0, f;
while (true) {
bfs(s);
if (level[t] < 0) return flow;
itr = vector<T>(V, 0);
while ((f = dfs(s, t, numeric_limits<T>::max())) > 0) {
flow += f;
}
}
}
};
int N, M, d;
int u[1010], v[1010], p[1010], q[1010], w[1010];
vector<int> t[1010];
int sz[1010];
int main() {
cin >> N >> M >> d;
rep(i, N) t[i].push_back(0);
t[N-1].push_back(1LL<<60);
rep(i, M) {
cin >> u[i] >> v[i] >> p[i] >> q[i] >> w[i];
u[i]--;
v[i]--;
q[i] += d;
t[u[i]].push_back(p[i]);
t[v[i]].push_back(q[i]);
}
rep(i, N) {
if (i) sz[i] += sz[i-1]+t[i-1].size();
}
Dinic<ll> flow(sz[N-1]+t[N-1].size()-1);
rep(i, N) {
sort(t[i].begin(), t[i].end());
t[i].erase(unique(t[i].begin(), t[i].end()));
rep(j, t[i].size()-1) flow.add_edge(sz[i]+j, sz[j]+j+1, 1LL<<60);
}
rep(i, N) {
ll x = sz[u[i]] + lower_bound(t[u[i]].begin(), t[u[i]].end(), p[i]) - t[u[i]].begin();
ll y = sz[v[i]] + lower_bound(t[v[i]].begin(), t[v[i]].end(), q[i]) - t[v[i]].begin();
flow.add_edge(x, y, w[i]);
}
cout << flow.max_flow(0, sz[N-1]+t[N-1].size()-1) << endl;
return 0;
}
dsytk7