結果
問題 | No.2604 Initial Motion |
ユーザー | srjywrdnprkt |
提出日時 | 2024-01-14 23:47:55 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 681 ms / 3,000 ms |
コード長 | 3,225 bytes |
コンパイル時間 | 2,235 ms |
コンパイル使用メモリ | 223,608 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-09-29 10:38:48 |
合計ジャッジ時間 | 14,576 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 1 ms
5,248 KB |
testcase_03 | AC | 19 ms
5,248 KB |
testcase_04 | AC | 19 ms
5,248 KB |
testcase_05 | AC | 18 ms
5,248 KB |
testcase_06 | AC | 18 ms
5,248 KB |
testcase_07 | AC | 19 ms
5,248 KB |
testcase_08 | AC | 18 ms
5,248 KB |
testcase_09 | AC | 18 ms
5,248 KB |
testcase_10 | AC | 18 ms
5,248 KB |
testcase_11 | AC | 19 ms
5,248 KB |
testcase_12 | AC | 18 ms
5,248 KB |
testcase_13 | AC | 522 ms
5,248 KB |
testcase_14 | AC | 296 ms
5,248 KB |
testcase_15 | AC | 216 ms
5,248 KB |
testcase_16 | AC | 452 ms
5,248 KB |
testcase_17 | AC | 681 ms
5,248 KB |
testcase_18 | AC | 627 ms
5,248 KB |
testcase_19 | AC | 616 ms
5,248 KB |
testcase_20 | AC | 437 ms
5,248 KB |
testcase_21 | AC | 319 ms
5,248 KB |
testcase_22 | AC | 593 ms
6,816 KB |
testcase_23 | AC | 348 ms
6,820 KB |
testcase_24 | AC | 496 ms
6,816 KB |
testcase_25 | AC | 634 ms
6,820 KB |
testcase_26 | AC | 409 ms
6,816 KB |
testcase_27 | AC | 241 ms
6,820 KB |
testcase_28 | AC | 368 ms
6,820 KB |
testcase_29 | AC | 515 ms
6,816 KB |
testcase_30 | AC | 275 ms
6,820 KB |
testcase_31 | AC | 402 ms
6,816 KB |
testcase_32 | AC | 383 ms
6,820 KB |
testcase_33 | AC | 3 ms
6,820 KB |
testcase_34 | AC | 374 ms
6,820 KB |
testcase_35 | AC | 370 ms
6,816 KB |
testcase_36 | AC | 357 ms
6,820 KB |
testcase_37 | AC | 3 ms
6,816 KB |
testcase_38 | AC | 2 ms
6,816 KB |
testcase_39 | AC | 2 ms
6,816 KB |
testcase_40 | AC | 258 ms
6,816 KB |
testcase_41 | AC | 258 ms
6,820 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; template<typename T> using pq = priority_queue<T, vector<T>, greater<T>>; struct edge{ int to; ll cap, cost; int rev; }; struct MinCostFlow{ private: static const ll INF=9e18; vector<ll> dist, h; vector<int> prevv, preve; vector<pair<int, int>> idx; void dijkstra(int s){ int from, c; pq<pair<ll, int>> que; for (int i=0; i<N; i++) dist[i] = INF; que.push({0, s}); dist[s] = 0; while(!que.empty()){ tie(c, from) = que.top(); que.pop(); if (dist[from] < c) continue; for (int i=0; i<E[from].size(); i++){ edge &e = E[from][i]; if (e.cap > 0 && dist[e.to] > dist[from]+e.cost+h[from]-h[e.to]){ dist[e.to] = dist[from]+e.cost+h[from]-h[e.to]; prevv[e.to] = from; preve[e.to] = i; que.push({dist[e.to], e.to}); } } } } public: int N, M=0; vector<vector<edge>> E; MinCostFlow(int n) : N(n) { E.resize(N); dist.resize(N); h.resize(N); prevv.resize(N); preve.resize(N); } void add_edge(int from, int to, ll cap, ll cost){ M++; idx.push_back({from, (int)E[from].size()}); E[from].push_back({to, cap, cost, (int)E[to].size()}); E[to].push_back({from, 0, -cost, (int)E[from].size()-1}); } ll solve(int s, int t, ll f){ ll res = 0, flow; fill(h.begin(), h.end(), 0); while(f){ dijkstra(s); if (dist[t] == INF) return -1; //no solution for (int i=0; i<N; i++) h[i] += dist[i]; flow = f; for (int i=t; i!=s; i=prevv[i]) flow = min(flow, E[prevv[i]][preve[i]].cap); f -= flow; res += flow*h[t]; for (int i=t; i!=s; i=prevv[i]){ edge &e = E[prevv[i]][preve[i]]; e.cap -= flow; E[i][e.rev].cap += flow; } } return res; } //i番目の辺の(from, to, flow)を求める。 tuple<int, int, ll> get_flow(int i){ assert (i < M); edge &e = E[idx[i].first][idx[i].second]; return {idx[i].first, e.to, E[e.to][e.rev].cap}; } }; int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); ll K, N, M, u, v, d, s, t; cin >> K >> N >> M; MinCostFlow mcf(N+2); s = 0; t = N+1; vector<ll> A(N+1), B(N+1); for (int i=1; i<=K; i++){ cin >> u; A[u]++; } for (int i=1; i<=N; i++) cin >> B[i]; for (int i=1; i<=N; i++) mcf.add_edge(s, i, A[i], 0); for (int i=1; i<=N; i++) mcf.add_edge(i, t, B[i], 0); for (int i=0; i<M; i++){ cin >> u >> v >> d; mcf.add_edge(u, v, K, d); mcf.add_edge(v, u, K, d); } cout << mcf.solve(s, t, K) << endl; return 0; }