# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll INF=1e18;
struct Info
{
	ll d;
	int a,b;
	Info(ll _d=0,int _a=0,int _b=0):d(_d),a(_a),b(_b){}
	bool operator<(const Info &x)const{return d<x.d;}
}val[2010][2010][4],sum[300010][4];
struct node2
{
	int v;ll w;
	node2(int _v=0,ll _w=0):v(_v),w(_w){}
	bool operator<(const node2 &a)const{return w>a.w;}
};
int n,m,cnt,X[100010],head[2010],to[4010],edge[4010],nxt[4010],tot=1;
ll dis[4010][4010];
vector<Info> vec[2010][2010];
vector<pair<int,int>> G[10010];
void addedge(int u,int v,int w)
{
	to[++tot]=v;
	edge[tot]=w;
	nxt[tot]=head[u];
	head[u]=tot;
}
void add(int u,int v,int w){G[u].emplace_back(v,w);}
void Dij(int s)
{
	static bool vis[10010];
	static priority_queue<node2> que;
	static ll _dis[10010];
	fill(_dis+1,_dis+cnt+1,INF);
	memset(vis+1,0,sizeof(bool)*cnt);
	_dis[s]=edge[s];que.emplace(s,_dis[s]);
	while(!que.empty())
	{
		node2 t=que.top();que.pop();
		if(vis[t.v]) continue;
		int u=t.v;vis[u]=1;
		for(auto p:G[u])
		{
			int v=p.first,w=p.second;
			if(_dis[v]>_dis[u]+w) _dis[v]=_dis[u]+w,que.emplace(v,_dis[v]);
		}
	}
	memcpy(dis[s]+2,_dis+2,sizeof(ll)*2*m);
}
void calc(vector<Info> vec,Info *val)
{
	sort(vec.begin(),vec.end());
	fill(val,val+4,Info(INF,0,0));
	int sz=vec.size();
	if(!sz) return;
	val[0]=vec[0];
	for(int i=1;i<sz;i++)
		if(vec[i].a!=val[0].a && vec[i].b!=val[0].b)
		{
			val[1]=vec[i];
			break;
		}
	for(int i=1;i<sz;i++)
		if(vec[i].a!=val[0].a && vec[i].b!=val[1].b)
		{
			val[2]=vec[i];
			break;
		}
	for(int i=1;i<sz;i++)
		if(vec[i].a!=val[1].a && vec[i].b!=val[0].b)
		{
			val[3]=vec[i];
			break;
		}
}
void merge(Info *v,Info *v1,Info *v2)
{
	vector<Info> vec;
	for(int i=0;i<4;i++)if(v1[i].a)
		for(int j=0;j<4;j++)if(v2[j].a)
			if(v1[i].b!=v2[j].a)
				vec.emplace_back(v1[i].d+v2[j].d,v1[i].a,v2[j].b);
	calc(vec,v);
}
void build(int rt,int l,int r)
{
	if(l==r)
	{
		for(int i=0;i<4;i++) sum[rt][i]=val[X[l]][X[l+1]][i];
		return;
	}
	int mid=(l+r)/2;
	build(rt*2,l,mid);build(rt*2+1,mid+1,r);
	merge(sum[rt],sum[rt*2],sum[rt*2+1]);
}
void update(int rt,int l,int r,int x)
{
	if(l==r)
	{
		for(int i=0;i<4;i++) sum[rt][i]=val[X[l]][X[l+1]][i];
		return;
	}
	int mid=(l+r)/2;
	if(x<=mid) update(rt*2,l,mid,x);
	else update(rt*2+1,mid+1,r,x);
	merge(sum[rt],sum[rt*2],sum[rt*2+1]);
}
int main()
{
	int q,L,u,v,w;
	cin>>n>>m>>q>>L;cnt=2*m+1;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&u,&v,&w);
		addedge(u,v,w);addedge(v,u,w);
	}
	for(int u=1;u<=n;u++)
	{
		int in=0,out=0;
		for(int i=head[u];i;i=nxt[i])
		{
			int w=edge[i];
			if(out) add(out,i,w),add(i^1,in,0);
			add(i^1,++cnt,0);
			if(out) add(out,cnt,0);
			out=cnt;
			add(++cnt,i,w);
			if(in) add(cnt,in,0);
			in=cnt;
		}
	}
	for(int i=2;i<=2*m+1;i++) Dij(i);
	for(int i=2;i<=2*m+1;i++)
		for(int j=2;j<=2*m+1;j++)
			if(dis[i][j]<INF)
				vec[to[i^1]][to[j]].emplace_back(dis[i][j],to[i],to[j^1]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			calc(vec[i][j],val[i][j]);
	for(int i=1;i<=L;i++) scanf("%d",&X[i]);
	build(1,1,L-1);
	while(q--)
	{
		scanf("%d%d",&u,&v);X[u]=v;
		if(u<L) update(1,1,L-1,u);
		if(u>1) update(1,1,L-1,u-1);
		ll ans=sum[1][0].d;
		if(ans<INF) printf("%lld\n",ans);
		else puts("-1");
	}
	return 0;
}