結果
問題 | No.1812 Uribo Road |
ユーザー | PCTprobability |
提出日時 | 2022-01-14 21:51:47 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 107 ms / 5,000 ms |
コード長 | 7,952 bytes |
コンパイル時間 | 5,301 ms |
コンパイル使用メモリ | 291,844 KB |
実行使用メモリ | 7,680 KB |
最終ジャッジ日時 | 2024-11-20 09:36:37 |
合計ジャッジ時間 | 7,039 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 3 ms
5,248 KB |
testcase_02 | AC | 3 ms
5,248 KB |
testcase_03 | AC | 10 ms
5,248 KB |
testcase_04 | AC | 3 ms
5,248 KB |
testcase_05 | AC | 3 ms
5,248 KB |
testcase_06 | AC | 2 ms
5,248 KB |
testcase_07 | AC | 3 ms
5,248 KB |
testcase_08 | AC | 4 ms
5,248 KB |
testcase_09 | AC | 5 ms
5,248 KB |
testcase_10 | AC | 4 ms
5,248 KB |
testcase_11 | AC | 4 ms
5,248 KB |
testcase_12 | AC | 46 ms
7,040 KB |
testcase_13 | AC | 25 ms
7,168 KB |
testcase_14 | AC | 24 ms
7,040 KB |
testcase_15 | AC | 34 ms
5,760 KB |
testcase_16 | AC | 30 ms
6,272 KB |
testcase_17 | AC | 88 ms
7,296 KB |
testcase_18 | AC | 19 ms
5,248 KB |
testcase_19 | AC | 106 ms
7,680 KB |
testcase_20 | AC | 107 ms
7,680 KB |
testcase_21 | AC | 85 ms
7,552 KB |
testcase_22 | AC | 66 ms
7,680 KB |
testcase_23 | AC | 5 ms
5,248 KB |
testcase_24 | AC | 2 ms
5,248 KB |
testcase_25 | AC | 15 ms
5,472 KB |
testcase_26 | AC | 3 ms
5,248 KB |
testcase_27 | AC | 14 ms
5,248 KB |
testcase_28 | AC | 28 ms
6,016 KB |
testcase_29 | AC | 2 ms
5,248 KB |
testcase_30 | AC | 5 ms
5,248 KB |
testcase_31 | AC | 55 ms
5,760 KB |
testcase_32 | AC | 3 ms
5,248 KB |
testcase_33 | AC | 35 ms
5,888 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #if __has_include(<atcoder/all>) #include <atcoder/all> using namespace atcoder; #endif using ll = long long; using ld = long double; using ull = unsigned long long; #define endl "\n" typedef pair<int,int> Pii; #define REP(i, n) for (int i = 0; i < (n); ++i) #define REP3(i, m, n) for (int i = (m); (i) < int(n); ++ (i)) #define FOR(i,a,b) for(ll i=a;i<=(ll)(b);i++) #define rep(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define ALL(x) begin(x), end(x) #define PB push_back #define rrep(i,a,b) for(int i=a;i>=b;i--) #define fore(i,a) for(auto &i:a) #define all(s) (s).begin(),(s).end() #define drep2(i, m, n) for (int i = (m)-1; i >= (n); --i) #define drep(i, n) drep2(i, n, 0) #define rever(vec) reverse(vec.begin(), vec.end()) #define sor(vec) sort(vec.begin(), vec.end()) #define fi first #define se second #define pb push_back #define P pair<ll,ll> #define PQminll priority_queue<ll, vector<ll>, greater<ll>> #define PQmaxll priority_queue<ll,vector<ll>,less<ll>> #define PQminP priority_queue<P, vector<P>, greater<P>> #define PQmaxP priority_queue<P,vector<P>,less<P>> #define NP next_permutation //typedef string::const_iterator State; //class ParseError {}; //const ll mod = 1000000009; const ll mod = 998244353; //const ll mod = 1000000007; const ll inf = 4100000000000000000ll; const ld eps = ld(0.000000000001); const long double pi = 3.141592653589793; template<class T>void vcin(vector<T> &n){for(int i=0;i<int(n.size());i++) cin>>n[i];} template<class T,class K>void vcin(vector<T> &n,vector<K> &m){for(int i=0;i<int(n.size());i++) cin>>n[i]>>m[i];} template<class T>void vcout(vector<T> &n){for(int i=0;i<int(n.size());i++){cout<<n[i]<<" ";}cout<<endl;} template<class T>void vcin(vector<vector<T>> &n){for(int i=0;i<int(n.size());i++){for(int j=0;j<int(n[i].size());j++){cin>>n[i][j];}}} template<class T>void vcout(vector<vector<T>> &n){for(int i=0;i<int(n.size());i++){for(int j=0;j<int(n[i].size());j++){cout<<n[i][j]<<" ";}cout<<endl;}cout<<endl;} void yes(bool a){cout<<(a?"yes":"no")<<endl;} void YES(bool a){cout<<(a?"YES":"NO")<<endl;} void Yes(bool a){cout<<(a?"Yes":"No")<<endl;} void possible(bool a){ cout<<(a?"possible":"impossible")<<endl; } void Possible(bool a){ cout<<(a?"Possible":"Impossible")<<endl; } void POSSIBLE(bool a){ cout<<(a?"POSSIBLE":"IMPOSSIBLE")<<endl; } template<class T>void print(T a){cout<<a<<endl;} template<class T,class F>void print(pair<T,F> a){cout<<a.fi<<" "<<a.se<<endl;} template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0;} template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0;} template<class T> void ifmin(T t,T u){if(t>u){cout<<-1<<endl;}else{cout<<t<<endl;}} template<class T> void ifmax(T t,T u){if(t>u){cout<<-1<<endl;}else{cout<<t<<endl;}} ll fastgcd(ll u,ll v){ll shl=0;while(u&&v&&u!=v){bool eu=!(u&1);bool ev=!(v&1);if(eu&&ev){++shl;u>>=1;v>>=1;}else if(eu&&!ev){u>>=1;}else if(!eu&&ev){v>>=1;}else if(u>=v){u=(u-v)>>1;}else{ll tmp=u;u=(v-u)>>1;v=tmp;}}return !u?v<<shl:u<<shl;} ll modPow(ll a, ll n, ll mod) { if(mod==1) return 0;ll ret = 1; ll p = a % mod; while (n) { if (n & 1) ret = ret * p % mod; p = p * p % mod; n >>= 1; } return ret; } vector<ll> divisor(ll x){ vector<ll> ans; for(ll i = 1; i * i <= x; i++){ if(x % i == 0) {ans.push_back(i); if(i*i!=x){ ans.push_back(x / ans[i]);}}}sor(ans); return ans; } ll pop(ll x){return __builtin_popcountll(x);} ll poplong(ll x){ll y=-1;while(x){x/=2;y++;}return y;} void cincout() { ios::sync_with_stdio(false); std::cin.tie(nullptr); cout<< fixed << setprecision(10); } using mint = modint1000000007; struct graph{ struct edge{ ll to,cost; }; ll v; vector<vector<edge>> g; vector<ll> d; vector<bool> negative; vector<bool> diameter; vector<ll> topological_sort; ll diametercost; bool bipartitecheck; vector<ll> bipartite; graph(ll n){ init(n); } void init(ll n){ v=n; g.resize(n); d.resize(n); negative.resize(n); diameter.resize(n); bipartite.resize(n); for(int i=0;i<v;i++){ d[i]=inf; bipartite[i]=-1; negative[i]=false; diameter[i]=false; } } void addedge(ll s,ll t,ll cost){ edge e; e.to=t; e.cost=cost; g[s].push_back(e); } void dijkstra(ll s){ for(int i=0;i<v;i++){ d[i]=inf; } d[s]=0; priority_queue<P,vector<P>,greater<P>> que; que.push(P(0,s)); while(!que.empty()){ P p=que.top(); que.pop(); ll V=p.second; if(d[V]<p.first) continue; for(auto e:g[V]){ if(d[e.to]>d[V]+e.cost){ d[e.to]=d[V]+e.cost; que.push(P(d[e.to],e.to)); } } } } void BellmanFord(ll s){ for(int i=0;i<v;i++){ d[i]=inf; negative[i]=false; } d[s]=0; for(int i=0;i<v;i++){ for(int V=0;V<v;V++){ if(d[V]==inf){ continue; } for(auto e:g[V]){ if(d[e.to]>d[V]+e.cost){ d[e.to]=d[V]+e.cost; if(i==v-1){ negative[e.to]=true; negative[V]=true; } } } } } } void bfs(ll s){ for(int i=0;i<v;i++){ d[i]=inf; } d[s]=0; queue<ll> q; q.push(s); while(q.size()){ ll v=q.front(); q.pop(); for(auto e:g[v]){ assert(e.cost==1); if(chmin(d[e.to],d[v]+1)) q.push(e.to); } } } void treediameter(){ bfs(0); ll p=0; ll q=0; for(int i=0;i<v;i++){ if(q<d[i]){ q=d[i]; p=i; } } diameter[p]=true; bfs(p); ll p2=0; ll q2=0; for(int i=0;i<v;i++){ if(q2<d[i]){ q2=d[i]; p2=i; } } diameter[p2]=true; diametercost=d[p2]; } void Bipartite(){ for(int i=0;i<v;i++){ if(bipartite[i]==-1){ Bipartitedfs(i); } } } void Bipartitedfs(ll s,ll cur=0){ bipartite[s]=cur; for(auto e:g[s]){ if(bipartite[e.to]!=-1){ if((bipartite[e.to]==bipartite[s])^(!e.cost%2)){ bipartitecheck=false; } } else{ if(e.cost%2){ Bipartitedfs(e.to,1-cur); } else{ Bipartitedfs(e.to,cur); } } } } void topologicalsort(){ for(int i=0;i<v;i++){ d[i]=0; } for(int i=0;i<v;i++){ if(d[i]) continue; topologicaldfs(i); } rever(topological_sort); } void topologicaldfs(ll a){ d[a]=1; for(auto e:g[a]){ if(d[e.to]) continue; topologicaldfs(e.to); } topological_sort.push_back(a); } }; ll dp[1<<12][26]; int main() { cincout(); ll n,m,k; cin>>n>>m>>k; vector<ll> x(k); vcin(x); vector<ll> p(k),q(k); graph g(n); vector<ll> y(k); for(int i=0;i<m;i++){ ll a,b,c; cin>>a>>b>>c; a--; b--; ll ind=-1; for(int j=0;j<k;j++){ if(i+1==x[j]) ind=j; } if(ind>=0){ p[ind]=a; q[ind]=b; y[ind]=c; } g.addedge(a,b,c); g.addedge(b,a,c); } map<ll,ll> ind,rev; for(int i=0;i<k;i++){ ind[p[i]]=2*i; ind[q[i]]=2*i+1; rev[2*i]=p[i]; rev[2*i+1]=q[i]; } ind[0]=2*k; ind[n-1]=2*k+1; rev[2*k]=0; rev[2*k+1]=n-1; vector<vector<ll>> d(2*k+2); for(int i=0;i<2*k+2;i++){ g.dijkstra(rev[i]); d[i]=g.d; } for(int i=0;i<(1<<12);i++){ for(int j=0;j<26;j++){ dp[i][j]=inf; } } dp[0][2*k]=0; for(int i=0;i<(1<<k)-1;i++){ for(int j=0;j<2*k+2;j++){ if(dp[i][j]>inf/2) continue; for(int v=0;v<k;v++){ if(!((i>>v)&1)){ chmin(dp[i+(1<<v)][2*v],dp[i][j]+d[j][rev[2*v+1]]+y[v]); chmin(dp[i+(1<<v)][2*v+1],dp[i][j]+d[j][rev[2*v]]+y[v]); } } } } ll ans=inf; for(int i=0;i<2*k+2;i++){ //cout<<dp[(1<<k)-1][i]<<" "<<d[i][n-1]<<" "<<rev[i]<<endl; chmin(ans,dp[(1<<k)-1][i]+d[i][n-1]); } cout<<ans<<endl; }