結果

問題 No.1301 Strange Graph Shortest Path
ユーザー carrot46carrot46
提出日時 2020-11-28 12:40:23
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 285 ms / 3,000 ms
コード長 3,699 bytes
コンパイル時間 2,039 ms
コンパイル使用メモリ 179,400 KB
実行使用メモリ 34,516 KB
最終ジャッジ日時 2023-10-09 23:52:40
合計ジャッジ時間 11,432 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,372 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 231 ms
32,792 KB
testcase_03 AC 191 ms
29,072 KB
testcase_04 AC 285 ms
31,204 KB
testcase_05 AC 201 ms
32,104 KB
testcase_06 AC 242 ms
29,308 KB
testcase_07 AC 233 ms
31,192 KB
testcase_08 AC 194 ms
29,760 KB
testcase_09 AC 244 ms
27,848 KB
testcase_10 AC 196 ms
29,332 KB
testcase_11 AC 257 ms
30,016 KB
testcase_12 AC 250 ms
30,156 KB
testcase_13 AC 231 ms
33,016 KB
testcase_14 AC 224 ms
28,144 KB
testcase_15 AC 223 ms
28,952 KB
testcase_16 AC 273 ms
31,236 KB
testcase_17 AC 245 ms
33,316 KB
testcase_18 AC 222 ms
30,264 KB
testcase_19 AC 254 ms
29,416 KB
testcase_20 AC 265 ms
28,272 KB
testcase_21 AC 238 ms
31,560 KB
testcase_22 AC 256 ms
29,052 KB
testcase_23 AC 239 ms
32,532 KB
testcase_24 AC 256 ms
28,980 KB
testcase_25 AC 267 ms
31,296 KB
testcase_26 AC 241 ms
30,304 KB
testcase_27 AC 249 ms
30,256 KB
testcase_28 AC 208 ms
32,048 KB
testcase_29 AC 280 ms
30,796 KB
testcase_30 AC 263 ms
31,280 KB
testcase_31 AC 269 ms
30,772 KB
testcase_32 AC 1 ms
4,372 KB
testcase_33 AC 159 ms
26,068 KB
testcase_34 AC 247 ms
34,516 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
//#include <chrono>
//#pragma GCC optimize("O3")
using namespace std;
#define reps(i,s,n) for(int i = s; i < n; i++)
#define rep(i,n) reps(i,0,n)
#define Rreps(i,n,e) for(int i = n - 1; i >= e; --i)
#define Rrep(i,n) Rreps(i,n,0)
#define ALL(a) a.begin(), a.end()

using ll = long long;
using vec = vector<ll>;
using mat = vector<vec>;

ll N,M,H,W,Q,K,A,B;
string S;
using P = pair<ll, ll>;
const ll INF = (1LL<<60);

template<class T> bool chmin(T &a, const T &b){
    if(a > b) {a = b; return true;}
    else return false;
}
template<class T> bool chmax(T &a, const T &b){
    if(a < b) {a = b; return true;}
    else return false;
}

struct edge{
    int to, cap, rev;
    ll cost;
    edge(){}
    edge(int a, int b, ll c, int d){
        to = a;
        cap = b;
        cost = c;
        rev = d;
    }
};
using ve = vector<edge>;

void add_edge(int from, int to, int cap, ll cost, vector<ve> &G, bool directed = true){
    rep(_, (directed ? 1 : 2)) {
        G[from].emplace_back(to, cap, cost, (int) G[to].size());
        G[to].emplace_back(from, 0, -cost, (int) G[from].size() - 1);
        swap(to, from);
    }
}

bool bellman_ford_for_flow(int s, vector<ve> &G, vec &dist){
    //O(|V|^2 + |V||E|)
    dist[s] = 0;
    int num(0), sz = G.size();
    while(num++ < sz){
        bool update = false;
        rep(v, sz){
            if(dist[v] != INF) {
                for (edge &e : G[v]) {
                    if (e.cap > 0 && chmin(dist[e.to], dist[v] + e.cost)) update = true;
                }
            }
        }
        if(!update) return true;
    }

    return false; //failed
}

void dijkstra_for_flow(int s, vector<ve> &G, vec &h, vec &dist, vector<int> &prev_v, vector<int> &prev_e){
    typedef struct _comp{
        ll cost;
        int v;
        _comp(){}
        _comp(ll a, int b){
            cost = a;
            v = b;
        }
        bool operator > (const _comp a){
            return cost > a.cost;
        }
    } mp;
    priority_queue<mp, vector<mp>, greater<> > pque;
    dist[s] = 0;
    pque.emplace(0, s);
    while(!pque.empty()){
        mp p = pque.top(); pque.pop();
        if(dist[p.v] < p.cost) continue;
        rep(i, (int)G[p.v].size()){
            edge &e = G[p.v][i];
            if(e.cap > 0 && chmin(dist[e.to], p.cost + e.cost + h[p.v] - h[e.to])){
                prev_v[e.to] = p.v;
                prev_e[e.to] = i;
                pque.emplace(dist[e.to], e.to);
            }
        }
    }
}

//全然Validateできていないので注意
ll min_cost_flow(int s, int t, int f, vector<ve> &G, bool minus_exist = true){
    //minus_exist : if (edge.cost < 0)
    //for graph without negative cycles
    int sz = G.size();
    ll res(0);
    vector<int> prev_v(sz), prev_e(sz);
    vec h(sz, minus_exist ? INF : 0);
    if(minus_exist && (!bellman_ford_for_flow(s, G, h) || h[t] == INF)) return -1;
    while(f > 0){
        vec dist(sz, INF);
        dijkstra_for_flow(s, G, h, dist, prev_v, prev_e);
        if(dist[t] == INF) return -1;
        rep(v, sz) h[v] += dist[v];
        int d = f;
        for(int v = t; v != s; v = prev_v[v]){
            chmin(d, G[prev_v[v]][prev_e[v]].cap);
        }
        f -= d;
        res += h[t] * d;
        for(int v = t; v != s; v = prev_v[v]){
            edge &e = G[prev_v[v]][prev_e[v]];
            e.cap -= d;
            G[v][e.rev].cap += d;
        }
    }
    return res;
}

int main(){
    cin>>N>>M;
    vector<ve> G(N);
    rep(i, M){
        int u, v, c, d;
        cin>>u>>v>>c>>d;
        --u; --v;
        add_edge(u, v, 1, c, G, false);
        add_edge(u, v, 1, d, G, false);
    }
    cout<<min_cost_flow(0, N - 1, 2, G, false)<<endl;
}
0