Problem:
JOI 国のお祭り事情(Festivals in JOI Kingdom)

JOI 国にはN 個の街があり，それらの間はM 本の双方向に通行可能な道路で結ばれている．国民は道路を通ってそれらの街を移動する．

JOI 国の国民には，お祭りが好きな人が多く，現在，K 個の街でお祭りが開催されており，非常に賑わっている．一方で，一部の国民はお祭りを騒がしいとして嫌い，お祭りにできるだけ近づきたくないと思っている．

そこで国王は，優秀なプログラマーであるあなたに，そのようなお祭りを嫌う国民のため，ある街からある街に移動するために，お祭りが開催されている街にどれだけ近づかずに移動することができるかを高速に調べることのできるプログラムの作成を依頼した．

課題

道路の情報とお祭りが開催されている街の情報，およびQ 個のクエリ（出発する街Si と行きたい街Tiの組）が与えられる．各クエリi に対し，街Si から街Ti へのすべての経路のうち，お祭りまでの距離が最大となる経路の，お祭りまでの距離を求めるプログラムを作成せよ．ただし，ある経路のお祭りまでの距離とは，経路上の街からお祭りが開催されている街までの移動距離の最小値のことである．

制限

2 ≤ N ≤ 100000 (= 105)     JOI 国の街の個数

1 ≤ M ≤ 200000 (= 2 × 105)     JOI 国の道路の本数

1 ≤ K ≤ N     お祭りが開催されている街の個数

1 ≤ Q ≤ 100000 (= 105)     クエリの個数

1 ≤ Li ≤ 1000     i 番目の道路の長さ

入力

標準入力から以下のデータを読み込め．

1 行目には整数N, M, K, Q が空白を区切りとして書かれている．N はJOI 国の街の個数を，M はJOI国の道路の本数を，K はお祭りが開催されている街の個数を，Q はクエリの個数をそれぞれ表す．街には1, 2, ... , N の番号がつけられている．

続くM 行は道路の情報を表す．i + 1 行目(1 ≤ i ≤ M) には整数Ai, Bi, Li (1 ≤ Ai ≤ N, 1 ≤ Bi ≤ N) が空白を区切りとして書かれている．これは，i 番目の道路が街Ai と街Bi を結んでおり，長さがLi であることを表す．道路の両端が同じ街であることはない．また，任意の2 つの街p, q に対し，p とqを結ぶ道路は2 本以上存在しない．どの街からどの街へもいくつかの道路をたどって行くことができる．

続くK 行はお祭りが開催されている街の情報を表す．i + M + 1 行目(1 ≤ i ≤ K) には1 つの整数Fi(1 ≤ Fi ≤ N) が書かれている．これは街Fi でお祭りが開催されていることを表す．F1, ..., FK の中に同じ値が2 回以上現れることはない．

続くQ 行はクエリを表す．i + M + K + 1 行目(1 ≤ i ≤ Q) には2 つの整数Si, Ti (1 ≤ Si ≤ N, 1 ≤ Ti ≤ N, Si ≠ Ti) が空白を区切りとして書かれている．これはi 番目のクエリの出発する街がSi であり行きたい街がTi であることを表す．

出力

標準出力に，全クエリへの答えをQ 行で出力せよ．すなわち，i 行目に，街Si から街Ti へのすべての経路のうち，お祭りまでの距離が最大となる経路の，お祭りまでの距離を表す整数を出力せよ．

採点基準

採点用データのうち，

配点の10% 分については，Q = 1 を満たす．

配点の20% 分については，N ≤ 5000, Q ≤ 5000 を満たす．

配点の30% 分については，これら2 つの条件の少なくとも一方を満たす．また，これら2 つの条件の両方を満たすような採点用データはない．

入出力例

入力例 1

6 6 2 3
1 2 5
2 3 4
2 4 6
3 5 9
4 5 3
5 6 7
1
6
3 4
5 2
1 4

出力例 1

7
5
0

6 つの街が6 本の道路で結ばれており，お祭りは街1, 6 の2 つの街で開催されている．クエリは以下の3 つである．

1 つ目のクエリは街3 から街4 へ行くというものである．街2 を経由する経路ではお祭りまでの距離は5，街5 を経由する経路ではお祭りまでの距離は7 となるため，答えは7 である．

2 つ目のクエリは街5 から街2 へ行くというものである．街3 と街4 のどちらを経由しても，街2 でお祭りまでの距離が最小となり，答えは5 である．

3 つ目のクエリは街1 から街4 へ行くというものである．街1 はお祭りが開催されている街であるため，答えは0 となる．

図1: 入力例1

入力例 2

12 17 2 5
1 3 6
1 6 7
2 3 8
2 4 4
2 8 11
2 12 2
3 6 3
3 7 8
3 11 2
4 12 2
5 10 3
6 10 5
8 9 6
8 12 7
9 10 6
11 9 10
12 9 5
8
7
2 6
5 2
1 10
8 9
9 4

出力例 2

8
8
11
0
6

図1: 入力例2

問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++ (judge verdict: Wrong Answer)

#include <cstdio>
#include <vector>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std;

#define MAX 100005
#define LMAX 20
#define DMAX 200000005
typedef pair<int,int> mp;
typedef pair<int,mp> mmp;
#define fir  first
#define sec  second
#define rep(i,s,n) for(int i=(s);i<(n);i++)

int n,m,k,q;

vector<mp> ver[MAX];
int dist[MAX];
//dist
vector<mp> tree[MAX];
int uni[MAX];
void uini(){
	rep(i,0,n+1){
		uni[i]=i;
	}
}
int upar(int p){
	if(uni[p]==p)return p;
	else return (uni[p]=upar(uni[p]));
}
bool usame(int p,int q){
	if(upar(p)==upar(q))return true;
	else return false;
}
void umer(int p,int q){
	uni[uni[p]]=uni[q];
}
//unifo

mp dob[MAX][LMAX];
int dep[MAX]={};
void dfs(int no,int pa,int di,int de){
	dep[no]=de;
	dob[no][0].sec=pa;
	dob[no][0].fir=di;
	rep(i,0,tree[no].size()){
		if(dep[tree[no][i].sec]==0){
			dfs(tree[no][i].sec,no,tree[no][i].fir,de+1);
		}
	}
}
//dob

int por(int fr,int by,int &ans){
	ans=min(ans,dob[fr][by].fir);
	return dob[fr][by].sec;
}
int lca(int p,int q){
	int ans=MAX;
	if(dep[p]<dep[q])swap(p,q);
	int a=LMAX;
	while(a>=0){
		if(dep[p]-dep[q]>=(1<<a))p=por(p,a,ans);
		a--;
	}
	if(p==q)return ans;
	int tr=0;
	for(int i=LMAX-1;i>=0;i--){
		if(por(p,i,tr)!=-1){
			p=por(p,i,ans);
			q=por(q,i,ans);
		}
	}
	return ans;
}
//lca&res

int main(void){
	memset(dist,-1,sizeof(dist));
	scanf("%d%d%d%d",&n,&m,&k,&q);
	rep(i,0,m){
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		ver[a].push_back(mp(c,b));
		ver[b].push_back(mp(c,a));
	}
	priority_queue<mp,vector<mp>,greater<mp> > que;
	rep(i,0,k){
		int a;
		scanf("%d",&a);
		que.push(mp(0,a));
	}
	
	while(!que.empty()){
		int to=que.top().sec,
			co=que.top().fir;
		que.pop();
		if(dist[to]!=-1)continue;
		dist[to]=co;
		rep(i,0,ver[to].size()){
			que.push( mp(co+ver[to][i].fir,ver[to][i].sec));
		}
	}
	//dist
	
	uini();
	priority_queue<mmp> krus;
	rep(i,1,n+1){
		rep(j,0,ver[i].size()){
			int to=ver[i][j].sec,
				co=ver[i][j].fir;
			krus.push( mmp(min(dist[i],dist[to]),mp(i,to)) );
		}
	}
	while(!krus.empty()){
		int a=krus.top().sec.fir,
		 	b=krus.top().sec.sec,
			c=krus.top().fir;
		krus.pop();
		if(!usame(a,b)){
			umer(a,b);
			tree[a].push_back(mp(c,b));
			tree[b].push_back(mp(c,a));
		}
	}
	//unifo
	
	dfs(1,-1,DMAX,1);
	rep(i,1,LMAX){
		rep(j,0,n){
			mp np=dob[j][i-1];
			if(np.sec==-1){
				dob[j][i]=np;
				continue;
			}
			mp pa=dob[np.sec][i-1];
			dob[j][i].sec=pa.sec;
			dob[j][i].fir=min(np.fir,pa.fir);
		}
	}
	//dob
	rep(i,0,q){
		int a,b;
		scanf("%d%d",&a,&b);
		printf("%d\n",lca(a,b));
	}
	//lca&res
	
	return 0;
}