結果

問題 No.1301 Strange Graph Shortest Path
ユーザー mugen_1337mugen_1337
提出日時 2020-11-28 19:30:34
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 222 ms / 3,000 ms
コード長 4,028 bytes
コンパイル時間 2,397 ms
コンパイル使用メモリ 213,096 KB
実行使用メモリ 43,976 KB
最終ジャッジ日時 2023-10-10 00:57:55
合計ジャッジ時間 10,471 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,352 KB
testcase_01 AC 1 ms
4,352 KB
testcase_02 AC 171 ms
40,780 KB
testcase_03 AC 141 ms
36,660 KB
testcase_04 AC 213 ms
38,776 KB
testcase_05 AC 142 ms
40,780 KB
testcase_06 AC 187 ms
36,096 KB
testcase_07 AC 179 ms
38,164 KB
testcase_08 AC 142 ms
37,160 KB
testcase_09 AC 181 ms
34,108 KB
testcase_10 AC 146 ms
36,568 KB
testcase_11 AC 191 ms
37,248 KB
testcase_12 AC 197 ms
37,056 KB
testcase_13 AC 173 ms
40,988 KB
testcase_14 AC 176 ms
34,472 KB
testcase_15 AC 172 ms
35,728 KB
testcase_16 AC 221 ms
38,596 KB
testcase_17 AC 188 ms
41,244 KB
testcase_18 AC 164 ms
37,692 KB
testcase_19 AC 193 ms
36,284 KB
testcase_20 AC 196 ms
34,792 KB
testcase_21 AC 182 ms
39,556 KB
testcase_22 AC 203 ms
35,736 KB
testcase_23 AC 175 ms
41,012 KB
testcase_24 AC 196 ms
35,636 KB
testcase_25 AC 221 ms
38,836 KB
testcase_26 AC 189 ms
37,564 KB
testcase_27 AC 195 ms
37,312 KB
testcase_28 AC 149 ms
40,612 KB
testcase_29 AC 222 ms
37,696 KB
testcase_30 AC 209 ms
38,692 KB
testcase_31 AC 212 ms
38,212 KB
testcase_32 AC 2 ms
4,348 KB
testcase_33 AC 90 ms
31,604 KB
testcase_34 AC 209 ms
43,976 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define ALL(x) begin(x),end(x)
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}

struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
template<typename T>
ostream &operator<<(ostream &os,const vector<T>&v){
    for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");
    return os;
}
template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}

template<typename flow_t, typename cost_t>
struct PrimalDual{
    struct edge{
        int to;
        flow_t cap;
        cost_t cost;
        int rev;//この辺の逆辺がg[from]の何番目にあるか
        bool isrev;
    };
 
    vector<vector<edge>> graph;
    vector<cost_t> potential,min_cost;
    vector<int> prevv,preve;//点，辺
    const cost_t TINF;
 
    PrimalDual(int V):graph(V),TINF(numeric_limits<cost_t>::max()){}
 
    void add_edge(int from,int to,flow_t cap,cost_t cost){
        graph[from].push_back((edge){to,cap,cost,(int)graph[to].size(),false});
        graph[to].push_back((edge){from,0,-cost,(int)graph[from].size()-1,true});
    }
 
    cost_t min_cost_flow(int s,int t,flow_t f,bool &ok){
        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,TINF);
            que.emplace(0,s);
            min_cost[s]=0;
            //dijkstraパート
            while(!que.empty()){
                Pi p=que.top();que.pop();
                if(min_cost[p.second]<p.first) continue;
                for(int i=0;i<(int)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 and 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]==TINF){
                ok=false;
                return ret;
            }
            // dijkstraの結果に応じてpotentialを調節
            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;
            }
        }
        ok=true;
        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;
            }
        }
    }
};

signed main(){
    int n,m;cin>>n>>m;
    PrimalDual<ll,ll> flow(n);
    rep(i,m){
        int u,v;ll c,d;cin>>u>>v>>c>>d;u--,v--;
        flow.add_edge(u,v,1,c);
        flow.add_edge(v,u,1,c);
        flow.add_edge(u,v,1,d);
        flow.add_edge(v,u,1,d);
    }
    bool ok;
    cout<<flow.min_cost_flow(0,n-1,2,ok)<<endl;
    return 0;
}
0