#include<bits/stdc++.h>
using namespace std;
const int N=2005;
int n,m,v,e;
int c[N],d[N];
int g[305][305];
double k[N],dp[N][N][2];
int main()
{
	//freopen("class.in","r",stdin);
	//freopen("class.out","w",stdout);
	memset(g,0x3f3f3f3f,sizeof(g));
	int a,b,w;
	cin>>n>>m>>v>>e;
	for(int i=1;i<=n;i++)	cin>>c[i];
	for(int i=1;i<=n;i++)	cin>>d[i];
	for(int i=1;i<=n;i++)	cin>>k[i];
	for(int i=1;i<=e;i++)
	{
		cin>>a>>b>>w;
		g[b][a]=g[a][b]=min(g[a][b],w);
	}
	for(int i=1;i<=v;i++)
		for(int j=1;j<=v;j++)
			for(int k=1;k<=v;k++)
				g[j][k]=min(g[j][k],g[j][i]+g[k][i]);
	for(int i=1;i<=v;i++)	g[i][i]=g[i][0]=g[0][i]=0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=min(m,i);j++)
			dp[i][j][0]=dp[i][j][1]=99999999.0;
	dp[1][0][0]=dp[1][1][1]=0.0;
	for(int i=2;i<=n;i++)
	{
		dp[i][0][0]=dp[i-1][0][0]+g[c[i-1]][c[i]];
		for(int j=1;j<=min(i,m);j++)
		{
			int c1=c[i-1],c2=d[i-1],c3=c[i],c4=d[i];
			dp[i][j][0]=min(dp[i][j][0],min(dp[i-1][j-1][0]+g[c1][c3],dp[i-1][j-1][1]+g[c1][c3]*(1-k[i-1])+g[c2][c3]*k[i-1]));
			dp[i][j][1]=min(dp[i][j][1],min(dp[i-1][j-1][0]+g[c1][c3]*(1-k[i])+g[c1][c4]*k[i],dp[i-1][j-1][1]+g[c2][c4]*k[i-1]*k[i]+g[c1][c4]*(1-k[i-1])*k[i]+g[c1][c3]*(1-k[i-1])*(1-k[i])+g[c2][c3]*k[i-1]*(1-k[i])));
		//	cout<<dp[i][j][0]<<' '<<dp[i][j][1]<<endl;
		}
	}
	double ans=99999999;
	for(int i=1;i<=m;i++)	ans=min(ans,min(dp[n][i][0],dp[n][i][1]));
	cout<<fixed<<setprecision(2)<<ans;
	return 0;
}