# include <bits/stdc++.h>
using namespace std;
using ll=long long;
constexpr ll INF=1e18;
int a[100010],W[8010],pre[4010],nxt[4010];
ll f[4010][4010],dis[8010];
vector<int> G[2010],rG[2010],H[8010];
struct Node
{
	int v;ll w;
	Node(int v,ll w):v(v),w(w){}
	bool operator<(const Node &t)const{return w>t.w;}
};
struct Path
{
	int s,t;ll d;
	Path(){}
	Path(int s,int t,ll d):s(s),t(t),d(d){}
	bool operator<(const Path &t)const{return d<t.d;}
	Path operator+(const Path &t)const{return Path(s,t.t,d+t.d);}
};
struct Info
{
	Path a[4];
	Info(){fill(a,a+4,Path(0,0,INF));}
	Info(vector<Path> &vec)
	{
		fill(a,a+4,Path(0,0,INF));
		if(vec.empty()) return;
		int sz=vec.size();
		static bool vis[4000010];
		sort(vec.begin(),vec.end());
		vis[0]=1;a[0]=vec[0];
		for(int i=1;i<sz;i++)
			if(!vis[i] && vec[i].s!=a[0].s && vec[i].t!=a[0].t)
			{
				vis[i]=1;a[1]=vec[i];
				break;
			}
		for(int i=1;i<sz;i++)
			if(!vis[i] && vec[i].s!=a[0].s && vec[i].t!=a[1].t)
			{
				vis[i]=1;a[2]=vec[i];
				break;
			}
		for(int i=1;i<sz;i++)
			if(!vis[i] && vec[i].s!=a[1].s && vec[i].t!=a[0].t)
			{
				vis[i]=1;a[3]=vec[i];
				break;
			}
		memset(vis,0,sizeof(bool)*sz);
	}
	Info operator+(const Info &t)const
	{
		vector<Path> vec;
		for(int i:{0,1,2,3})
			if(a[i].d<INF)
				for(int j:{0,1,2,3})
					if(t.a[j].d<INF && (a[i].t-1)/2!=(t.a[j].s-1)/2)
						vec.push_back(a[i]+t.a[j]);
		return Info(vec);
	}
}g[2010][2010],val[100010],sum[300010];
void Dijk(int n,int s)
{
	static bool vis[8010];
	static priority_queue<Node,vector<Node>> pq;
	fill(dis+1,dis+n+1,INF);
	memset(vis+1,0,sizeof(bool)*n);
	pq.emplace(s,dis[s]=W[s]);
	while(!pq.empty())
	{
		int u=pq.top().v;pq.pop();
		if(vis[u]) continue;
		vis[u]=1;
		for(int v:H[u])
			if(dis[u]+W[v]<dis[v])
				pq.emplace(v,dis[v]=dis[u]+W[v]);
	}
}
inline void pushup(int rt){sum[rt]=sum[rt*2]+sum[rt*2+1];}
void build(int rt,int l,int r)
{
	if(l==r) return sum[rt]=val[l],void();
	int mid=(l+r)/2;
	build(rt*2,l,mid);
	build(rt*2+1,mid+1,r);
	pushup(rt);
}
void update(int rt,int l,int r,int x,const Info &y)
{
	if(l==r) return sum[rt]=y,void();
	int mid=(l+r)/2;
	if(x<=mid) update(rt*2,l,mid,x,y);
	else update(rt*2+1,mid+1,r,x,y);
	pushup(rt);
}
int main()
{
	int n,m,q,L,u,v,w,tot=0;
	cin>>n>>m>>q>>L;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&u,&v,&w);
		W[++tot]=w;
		G[u].push_back(tot);
		rG[v].push_back(tot);
		W[++tot]=w;
		G[v].push_back(tot);
		rG[u].push_back(tot);
	}
	for(int i=1;i<=n;i++)
	{
		int sz=G[i].size();
		for(int j=0;j<sz;j++) H[tot+j+1].push_back(G[i][j]);
		for(int j=1;j<sz;j++) H[tot+j+1].push_back(tot+j),pre[G[i][j]]=tot+j;
		tot+=sz;
		for(int j=0;j<sz;j++) H[tot+j+1].push_back(G[i][j]);
		for(int j=1;j<sz;j++) H[tot+j].push_back(tot+j+1),nxt[G[i][j-1]]=tot+j+1;
		tot+=sz;
	}
	for(int i=1;i<=2*m;i+=2)
	{
		if(pre[i+1]) H[i].push_back(pre[i+1]);
		if(nxt[i+1]) H[i].push_back(nxt[i+1]);
		if(pre[i]) H[i+1].push_back(pre[i]);
		if(nxt[i]) H[i+1].push_back(nxt[i]);
	}
	for(int i=1;i<=2*m;i++)
	{
		Dijk(tot,i);
		memcpy(f[i]+1,dis+1,sizeof(ll)*2*m);
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			if(i==j) continue;
			vector<Path> vec;
			for(int p:G[i])
				for(int q:rG[j])
					vec.emplace_back(p,q,f[p][q]);
			g[i][j]=Info(vec);
		}
	for(int i=1;i<=L;i++) scanf("%d",&a[i]);
	for(int i=1;i<L;i++) val[i]=g[a[i]][a[i+1]];
	build(1,1,L-1);
	while(q--)
	{
		scanf("%d%d",&u,&v);a[u]=v;
		if(u>1) update(1,1,L-1,u-1,g[a[u-1]][v]);
		if(u<L) update(1,1,L-1,u,g[v][a[u+1]]);
		if(sum[1].a[0].d<INF) printf("%lld\n",sum[1].a[0].d);
		else puts("-1");
	}
	return 0;
}