結果

問題 No.1301 Strange Graph Shortest Path
ユーザー momoharamomohara
提出日時 2020-11-27 22:52:13
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 206 ms / 3,000 ms
コード長 4,468 bytes
コンパイル時間 2,770 ms
コンパイル使用メモリ 216,668 KB
実行使用メモリ 34,992 KB
最終ジャッジ日時 2024-07-26 20:05:38
合計ジャッジ時間 10,067 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 161 ms
32,508 KB
testcase_03 AC 132 ms
28,928 KB
testcase_04 AC 201 ms
31,360 KB
testcase_05 AC 130 ms
32,000 KB
testcase_06 AC 177 ms
29,312 KB
testcase_07 AC 166 ms
30,336 KB
testcase_08 AC 132 ms
29,312 KB
testcase_09 AC 164 ms
27,904 KB
testcase_10 AC 132 ms
29,056 KB
testcase_11 AC 182 ms
30,208 KB
testcase_12 AC 186 ms
30,080 KB
testcase_13 AC 162 ms
32,808 KB
testcase_14 AC 163 ms
27,904 KB
testcase_15 AC 155 ms
28,800 KB
testcase_16 AC 203 ms
31,488 KB
testcase_17 AC 172 ms
32,908 KB
testcase_18 AC 151 ms
29,824 KB
testcase_19 AC 183 ms
29,312 KB
testcase_20 AC 177 ms
28,544 KB
testcase_21 AC 167 ms
31,484 KB
testcase_22 AC 185 ms
29,184 KB
testcase_23 AC 163 ms
32,384 KB
testcase_24 AC 182 ms
29,056 KB
testcase_25 AC 192 ms
31,488 KB
testcase_26 AC 167 ms
30,208 KB
testcase_27 AC 176 ms
30,208 KB
testcase_28 AC 138 ms
31,872 KB
testcase_29 AC 206 ms
30,848 KB
testcase_30 AC 186 ms
31,104 KB
testcase_31 AC 191 ms
31,104 KB
testcase_32 AC 2 ms
6,944 KB
testcase_33 AC 84 ms
25,472 KB
testcase_34 AC 190 ms
34,992 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
/*
#include <atcoder/all>
using namespace atcoder;
*/
#define all(hoge) (hoge).begin(), (hoge).end()
#define en '\n'
using ll = long long;
using ull = unsigned long long;
#define rep(i, m, n) for(ll i = (ll)(m); i < (ll)(n); ++i)
#define rep2(i, m, n) for(ll i = (ll)(n)-1; i >= (ll)(m); --i)
#define REP(i, n) rep(i, 0, n)
#define REP2(i, n) rep2(i, 0, n)
template <class T>
using vec = vector<T>;
template <class T>
using vvec = vector<vec<T>>;
typedef pair<ll, ll> P;
using tp = tuple<ll, ll, ll>;

constexpr long long INF = 1LL << 60;
constexpr int INF_INT = 1 << 25;
//constexpr long long MOD = (ll)1e9 + 7;
constexpr long long MOD = 998244353LL;
using ld = long double;
static const ld pi = 3.141592653589793L;

using Array = vector<ll>;
using Matrix = vector<Array>;

/*
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
*/

template <class T>
inline bool chmin(T &a, T b) {
    if(a > b) {
        a = b;
        return true;
    }
    return false;
}
template <class T>
inline bool chmax(T &a, T b) {
    if(a < b) {
        a = b;
        return true;
    }
    return false;
}
template <typename flow_t, typename cost_t>
struct PrimalDual {
    const cost_t INF;

    struct edge {
        int to;
        flow_t cap;
        cost_t cost;
        int rev;
        bool isrev;
    };
    vector<vector<edge>> graph;
    vector<cost_t> potential, min_cost;
    vector<int> prevv, preve;

    PrimalDual(int V) : graph(V), INF(numeric_limits<cost_t>::max()) {}

    void add_edge(int from, int to, flow_t cap, cost_t cost) {
        graph[from].emplace_back((edge){to, cap, cost, (int)graph[to].size(), false});
        graph[to].emplace_back((edge){from, 0, -cost, (int)graph[from].size() - 1, true});
    }

    cost_t min_cost_flow(int s, int t, flow_t f) {
        int V = (int)graph.size();
        cost_t ret = 0;
        using Pi = pair<cost_t, int>;
        priority_queue<Pi, vector<Pi>, greater<Pi>> que;
        potential.assign(V, 0);
        preve.assign(V, -1);
        prevv.assign(V, -1);

        while(f > 0) {
            min_cost.assign(V, INF);
            que.emplace(0, s);
            min_cost[s] = 0;
            while(!que.empty()) {
                Pi p = que.top();
                que.pop();
                if(min_cost[p.second] < p.first)
                    continue;
                for(int i = 0; i < graph[p.second].size(); i++) {
                    edge &e = graph[p.second][i];
                    cost_t nextCost = min_cost[p.second] + e.cost + potential[p.second] - potential[e.to];
                    if(e.cap > 0 && min_cost[e.to] > nextCost) {
                        min_cost[e.to] = nextCost;
                        prevv[e.to] = p.second, preve[e.to] = i;
                        que.emplace(min_cost[e.to], e.to);
                    }
                }
            }
            if(min_cost[t] == INF)
                return -1;
            for(int v = 0; v < V; v++)
                potential[v] += min_cost[v];
            flow_t addflow = f;
            for(int v = t; v != s; v = prevv[v]) {
                addflow = min(addflow, graph[prevv[v]][preve[v]].cap);
            }
            f -= addflow;
            ret += addflow * potential[t];
            for(int v = t; v != s; v = prevv[v]) {
                edge &e = graph[prevv[v]][preve[v]];
                e.cap -= addflow;
                graph[v][e.rev].cap += addflow;
            }
        }
        return ret;
    }

    void output() {
        for(int i = 0; i < graph.size(); i++) {
            for(auto &e : graph[i]) {
                if(e.isrev)
                    continue;
                auto &rev_e = graph[e.to][e.rev];
                cout << i << "->" << e.to << " (flow: " << rev_e.cap << "/" << rev_e.cap + e.cap << ")" << endl;
            }
        }
    }
};

void solve() {
    int n, m;
    cin >> n >> m;
    PrimalDual<int, ll> pd(n);
    REP(i, m) {
        int u, v, c, d;
        cin >> u >> v >> c >> d;
        u--;
        v--;
        pd.add_edge(u, v, 1, c);
        pd.add_edge(v, u, 1, c);
        pd.add_edge(u, v, 1, d);
        pd.add_edge(v, u, 1, d);
    }

    auto ans = pd.min_cost_flow(0, n - 1, 2);
    cout << ans << en;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    /*
    ll t;
    cin >> t;
    REP(i, t - 1) {
        solve();
    }*/

    solve();

    return 0;
}
0