結果

問題 No.1069 電柱 / Pole (Hard)
ユーザー mugen_1337mugen_1337
提出日時 2021-03-05 18:02:10
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 196 ms / 2,000 ms
コード長 6,269 bytes
コンパイル時間 3,049 ms
コンパイル使用メモリ 239,056 KB
実行使用メモリ 10,992 KB
最終ジャッジ日時 2023-09-30 21:21:29
合計ジャッジ時間 6,374 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 196 ms
6,396 KB
testcase_05 AC 176 ms
10,992 KB
testcase_06 AC 7 ms
4,380 KB
testcase_07 AC 8 ms
4,376 KB
testcase_08 AC 7 ms
4,380 KB
testcase_09 AC 7 ms
4,380 KB
testcase_10 AC 5 ms
4,380 KB
testcase_11 AC 2 ms
4,380 KB
testcase_12 AC 3 ms
4,376 KB
testcase_13 AC 6 ms
4,380 KB
testcase_14 AC 7 ms
4,380 KB
testcase_15 AC 5 ms
4,376 KB
testcase_16 AC 5 ms
4,376 KB
testcase_17 AC 4 ms
4,376 KB
testcase_18 AC 6 ms
4,380 KB
testcase_19 AC 8 ms
4,380 KB
testcase_20 AC 4 ms
4,380 KB
testcase_21 AC 8 ms
4,384 KB
testcase_22 AC 4 ms
4,376 KB
testcase_23 AC 18 ms
4,376 KB
testcase_24 AC 7 ms
4,376 KB
testcase_25 AC 8 ms
4,376 KB
testcase_26 AC 6 ms
4,376 KB
testcase_27 AC 6 ms
4,380 KB
testcase_28 AC 6 ms
4,376 KB
testcase_29 AC 3 ms
4,380 KB
testcase_30 AC 2 ms
4,376 KB
testcase_31 AC 1 ms
4,376 KB
testcase_32 AC 2 ms
4,376 KB
testcase_33 AC 1 ms
4,380 KB
testcase_34 AC 2 ms
4,376 KB
testcase_35 AC 2 ms
4,380 KB
testcase_36 AC 2 ms
4,376 KB
testcase_37 AC 3 ms
4,376 KB
testcase_38 AC 3 ms
4,380 KB
testcase_39 AC 2 ms
4,376 KB
testcase_40 AC 3 ms
4,376 KB
testcase_41 AC 4 ms
4,376 KB
testcase_42 AC 3 ms
4,380 KB
testcase_43 AC 4 ms
4,376 KB
testcase_44 AC 5 ms
4,380 KB
testcase_45 AC 6 ms
4,380 KB
testcase_46 AC 7 ms
4,380 KB
testcase_47 AC 6 ms
4,376 KB
testcase_48 AC 6 ms
4,380 KB
testcase_49 AC 6 ms
4,376 KB
testcase_50 AC 3 ms
4,380 KB
testcase_51 AC 6 ms
4,376 KB
testcase_52 AC 3 ms
4,376 KB
testcase_53 AC 7 ms
4,376 KB
testcase_54 AC 2 ms
4,376 KB
testcase_55 AC 2 ms
4,376 KB
testcase_56 AC 2 ms
4,380 KB
testcase_57 AC 2 ms
4,380 KB
testcase_58 AC 1 ms
4,500 KB
testcase_59 AC 2 ms
4,376 KB
testcase_60 AC 2 ms
4,376 KB
testcase_61 AC 2 ms
4,376 KB
testcase_62 AC 1 ms
4,376 KB
testcase_63 AC 1 ms
4,380 KB
testcase_64 AC 2 ms
4,376 KB
testcase_65 AC 2 ms
4,384 KB
testcase_66 AC 2 ms
4,376 KB
testcase_67 AC 2 ms
4,376 KB
testcase_68 AC 2 ms
4,376 KB
testcase_69 AC 4 ms
4,380 KB
testcase_70 AC 3 ms
4,380 KB
testcase_71 AC 4 ms
4,376 KB
testcase_72 AC 3 ms
4,376 KB
testcase_73 AC 3 ms
4,380 KB
testcase_74 AC 4 ms
4,380 KB
testcase_75 AC 4 ms
4,376 KB
testcase_76 AC 7 ms
4,376 KB
testcase_77 AC 7 ms
4,376 KB
testcase_78 AC 2 ms
4,380 KB
testcase_79 AC 8 ms
4,376 KB
testcase_80 AC 6 ms
4,380 KB
testcase_81 AC 9 ms
4,380 KB
testcase_82 AC 5 ms
4,380 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;
}

// graph template
// ref : https://ei1333.github.io/library/graph/graph-template.cpp
template<typename T=int>
struct Edge{
    int from,to;
    T w;
    int idx;
    Edge()=default;
    Edge(int from,int to,T w=1,int idx=-1):from(from),to(to),w(w),idx(idx){}
    operator int() const{return to;}
};

template<typename T=int>
struct Graph{
    vector<vector<Edge<T>>> g;
    int V,E;
    Graph()=default;
    Graph(int n):g(n),V(n),E(0){}

    size_t size(){
        return g.size();
    }
    void resize(int k){
        g.resize(k);
    }
    inline const vector<Edge<T>> &operator[](int k)const{
        return (g.at(k));
    }
    inline vector<Edge<T>> &operator[](int k){
        return (g.at(k));
    }
    void add_directed_edge(int from,int to,T cost=1){
        g[from].emplace_back(from,to,cost,E++);
    }
    void add_edge(int from,int to,T cost=1){
        g[from].emplace_back(from,to,cost,E);
        g[to].emplace_back(to,from,cost,E++);
    }
    void read(int m,int pad=-1,bool weighted=false,bool directed=false){
        for(int i=0;i<m;i++){
            int u,v;cin>>u>>v;
            u+=pad,v+=pad;
            T w=T(1);
            if(weighted) cin>>w;
            if(directed) add_directed_edge(u,v,w);
            else         add_edge(u,v,w);
        }
    }
};

// 有向グラフ向け
// 無向グラフでやりたいなら辺のidxが独立になるように振りなおす
// O(KN * (N+M)logN)のはずだけどどこか悪化してるかも
/*
ref: https://qiita.com/nariaki3551/items/821dc6ffdc552d3d5f22
     https://ei1333.github.io/library/graph/shortest-path/k-shortest-path.cpp
*/
template<typename T>
vector<pair<T,vector<int>>> Yen(Graph<T> g,int s,int t,int k){
    T inf=numeric_limits<T>::max();

    vector<Edge<T>> edge(g.E);
    for(int i=0;i<g.V;i++)for(auto &e:g[i]) edge[e.idx]=e;

    unordered_set<int> removed_edge;
    using P=pair<T,int>;
    auto dijkstra=[&](vector<T> &d,vector<int> &from_v,vector<int> &from_e,int st){
        priority_queue<P,vector<P>,greater<P>> que;
        que.emplace(d[st],st);
        while(!que.empty()){
            auto p=que.top();que.pop();
            int cur=p.second;
            if(d[cur]<p.first) continue;
            for(auto &e:g[cur])if(!removed_edge.count(e.idx)){
                if(chmin(d[e.to],d[cur]+e.w)){
                    from_v[e.to]=cur;
                    from_e[e.to]=e.idx;
                    que.emplace(d[e.to],e.to);
                }
            }
        }
    };

    // (s, t)へのパスの辺番号のvector
    auto get_path=[&](vector<int> &from_v,vector<int> &from_e,int st,int gl){
        vector<int> e_idx;
        for(;gl!=st;gl=from_v[gl]) e_idx.push_back(from_e[gl]);
        reverse(begin(e_idx),end(e_idx));
        return e_idx;
    };

    vector<pair<T,vector<int>>> A;
    set<pair<T,vector<int>>> B;

    vector<T> d;
    vector<int> from_v,from_e;
    d.assign(g.V,inf);
    d[s]=0;
    from_v.assign(g.V,-1);
    from_e.assign(g.V,-1);
    dijkstra(d,from_v,from_e,s);

    if(d[t]==inf) return A;

    A.emplace_back(d[t],get_path(from_v,from_e,s,t));

    for(int i=1;i<k;i++){
        vector<int> &last_path=A.back().second;
        d.assign(g.V,inf);
        d[s]=0;
        from_v.assign(g.V,-1);
        from_e.assign(g.V,-1);

        vector<int> same_path(i); // last pathに沿った道のリスト
        iota(begin(same_path),end(same_path),0);

        // O(N)
        for(int j=0;j<(int)last_path.size();j++){
            int spur=edge[last_path[j]].from;// spur node

            unordered_set<int> removed_edge_new;
            for(auto &l:same_path)if(j<(int)A[l].second.size()) removed_edge_new.insert(A[l].second[j]);
            swap(removed_edge,removed_edge_new);
            
            vector<T> d_t{d};
            vector<int> from_v_t{from_v},from_e_t{from_e};

            dijkstra(d_t,from_v_t,from_e_t,spur);
            if(d_t[t]<inf){
                auto path=get_path(from_v_t,from_e_t,s,t);
                bool ok=true;
                for(auto l:same_path){
                    if(path==A[l].second){
                        ok=false;
                        break;
                    }
                }
                if(ok) B.emplace(d_t[t],path);
            }

            // last pathに沿って1つ進める
            auto &spur_root=edge[last_path[j]];
            d[spur_root.to]=d[spur_root.from]+spur_root.w;
            from_e[spur_root.to]=spur_root.idx;
            from_v[spur_root.to]=spur_root.from;

            vector<int> tmp;
            for(auto l:same_path)if(j<(int)A[l].second.size() and A[l].second[j]==last_path[j])
                tmp.push_back(l);
            swap(tmp,same_path);
        }
        if(B.empty()) return A;
        
        A.push_back(*begin(B));
        B.erase(begin(B));
    }
    return A;
}

signed main(){
    int n,m,k;cin>>n>>m>>k;
    int s,t;cin>>s>>t;s--,t--;
    vector<double> p(n),q(n);
    rep(i,n) cin>>p[i]>>q[i];

    Graph<double> g(n);
    rep(i,m){
        int u,v;cin>>u>>v;u--,v--;
        double c=hypot(p[u]-p[v],q[u]-q[v]);
        g.add_directed_edge(u,v,c);
        g.add_directed_edge(v,u,c);
    }
    auto res=Yen(g,s,t,k);
    rep(i,k){
        if(i<(int)res.size()) cout<<res[i].first<<endl;
        else cout<<-1<<endl;
    }

    return 0;
}
0