結果

問題 No.1308 ジャンプビーコン
コンテスト
ユーザー 沙耶花
提出日時 2020-12-05 13:47:20
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 3,849 bytes
コンパイル時間 3,592 ms
コンパイル使用メモリ 234,680 KB
最終ジャッジ日時 2025-01-16 17:48:27
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18 WA * 1 TLE * 18
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:13:33: warning: use of ‘auto’ in parameter declaration only available with ‘-std=c++20’ or ‘-fconcepts’
   13 | vector<long long> get(int start,auto &E){
      |                                 ^~~~
main.cpp: In function ‘int main()’:
main.cpp:128:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
  128 |                 scanf("%d %d %d",&u,&v,&l);
      |                 ~~~~~^~~~~~~~~~~~~~~~~~~~~
main.cpp:140:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
  140 |                 scanf("%d",&x[i]);
      |                 ~~~~~^~~~~~~~~~~~

ソースコード

diff # 

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

#include <stdio.h>
#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000000000000



vector<long long> get(int start,auto &E){
	vector<long long> dis(E.size(),Inf);
	dis[start] = 0;
	queue<int> Q;
	Q.push(start);
	
	while(Q.size()>0){
		int u = Q.front();
		Q.pop();
		rep(i,E[u].size()){
			int v = E[u][i].first;
			if(dis[v]!=Inf)continue;
			dis[v] = dis[u] + E[u][i].second;
			Q.push(v);
		}
	}
	
	return dis;
}


struct HLD{
	vector<int> sz,parent,depth,root,pos;
	vector<int> arr;
	HLD(vector<vector<int>> &E){
		sz.resize(E.size(),1);
		parent.resize(E.size(),0);
		depth.resize(E.size(),0);
		root.resize(E.size(),0);
		pos.resize(E.size(),0);
		
		dfs(0,-1,E);
		dfs2(0,-1,E,0);
	}
	
	void dfs(int now,int p,vector<vector<int>> &E){
		parent[now] = p;
		if(p==-1){
			depth[now] = 0;
		}
		else{
			depth[now] = depth[p]+1;
		}
		for(int i=0;i<E[now].size();i++){
			int to = E[now][i];
			if(to==p)continue;
			dfs(to,now,E);
			sz[now] += sz[to];
		}
	}
	
	void dfs2(int now,int p,vector<vector<int>> &E,int r){
		pos[now] = arr.size();
		arr.push_back(now);
		root[now] = r;
		int maxi = 0;
		int ind = -1;
		for(int i=0;i<E[now].size();i++){
			int to = E[now][i];
			if(to==p)continue;
			if(maxi<sz[to]){
				maxi = sz[to];
				ind = to;
			}
		}
		if(ind==-1)return;
		dfs2(ind,now,E,r);
		for(int i=0;i<E[now].size();i++){
			int to = E[now][i];
			if(to==p||to==ind)continue;
			dfs2(to,now,E,to);
		}
	}
	
	vector<pair<int,int>> query(int u,int v){
		vector<pair<int,int>> ret;
		int t = 0;
		while(root[u]!=root[v]){
			if(depth[root[u]] <= depth[root[v]]){
				ret.insert(ret.begin()+t,{pos[root[v]], pos[v]});
				v = parent[root[v]];
			}
			else{
				ret.insert(ret.begin()+t,{pos[u],pos[root[u]]});
				u = parent[root[u]];
				t++;
			}
		}
		ret.insert(ret.begin()+t,{pos[u],pos[v]});
		return ret;
	}
	
	int lca(int u,int v){
		for(;;v=parent[root[v]]){
			if(pos[u]>pos[v])swap(u,v);
			if(root[u]==root[v])return u;
		}
	}
	
	int get_distance(int u,int v){
		return depth[u] + depth[v] - 2 * depth[lca(u,v)];
	}
	
};

int main(){
	
	int N,Q;
	long long C;
	cin>>N>>Q>>C;
	
	vector E(N,vector<pair<int,long long>>());
	vector EE(N,vector<int>());
	rep(i,N-1){
		int u,v,l;
		scanf("%d %d %d",&u,&v,&l);
		
		u--;v--;
		E[u].emplace_back(v,l);
		E[v].emplace_back(u,l);
		
		EE[u].push_back(v);
		EE[v].push_back(u);
	}
	
	vector<int> x(Q);
	rep(i,Q){
		scanf("%d",&x[i]);
		x[i]--;
	}
	
	vector dis(N,vector<long long>());
	
	rep(i,N){
		dis[i] = get(i,E);
	}
	
	HLD H(EE);
	
	vector<long long> dp(N,Inf);
	
	dp[H.pos[x[0]]] = 0LL;
	
	rep(i,Q-1){
		vector<long long> ndp(N,Inf);

		vector<pair<long long,int>> V;
		set<int> S;
		rep(j,N){
			if(j==x[i]){
				long long t = dp[H.pos[j]]+dis[j][x[i+1]];
				if(t<Inf)V.emplace_back(t,j);
			}
			else{
				long long t = dp[H.pos[j]]+dis[j][x[i+1]]+C;
				if(t<Inf)V.emplace_back(t,j);
			}
			
			long long t = dp[H.pos[j]] + min(dis[x[i]][x[i+1]],C+dis[j][x[i+1]]);
			
			if(t<Inf)V.emplace_back(t,-j-1);
			S.emplace_hint(S.end(),j);
			
		}
		
		sort(V.begin(),V.end());
		
		rep(j,V.size()){		
			if(S.size()==0)break;
			long long t = V[j].first;
			if(t>=Inf)break;
			
			
			int ind = V[j].second;
			if(ind<0){
				ind++;
				ind *= -1;
				
				if(S.count(H.pos[ind])){
					ndp[H.pos[ind]] = t;
					S.erase(H.pos[ind]);
				}
				continue;
			}
			auto ret = H.query(ind,x[i+1]);
			for(auto &x:ret){
				if(x.first>x.second)swap(x.first,x.second);
				auto it = S.lower_bound(x.first);
				while(it!=S.end() && (*it)<=x.second){
					ndp[*it] = t;
					it = S.erase(it);
				}
			}
			
			
		}
		
		swap(dp,ndp);
		
		
	}
	
	long long ans = Inf;
	rep(i,N)ans = min(ans,dp[i]);
	
	cout<<ans<<endl;
	
	
    return 0;
}
0