a点 1 4 8 12 13 21 24

#include<bits/stdc++.h>
using namespace std;
int n,m,v,e;
int c[2005],d[2005],dis[305][305];
double k[2005],f[2005][2005][2];
int main(){
	scanf("%d%d%d%d",&n,&m,&v,&e);
	for(int i=1;i<=n;i++){
		scanf("%d",&c[i]);
	}	
	for(int i=1;i<=n;i++){
		scanf("%d",&d[i]);
	}
	for(int i=1;i<=n;i++){
		scanf("%lf",&k[i]);
	}
	for(int i=1;i<=v;i++){
		for(int j=1;j<i;j++){
			dis[j][i]=dis[i][j]=10000005;
		}
	}		
	for(int i=1;i<=e;i++){
		int a,b,w;
		scanf("%d%d%d",&a,&b,&w);
		dis[a][b]=min(dis[a][b],w);
		dis[b][a]=dis[a][b];
    }
	for(int h=1;h<=v;h++){//弗洛伊德 
		for(int i=1;i<=v;i++){
			for(int j=1;j<i;j++){
				if(dis[i][j]>dis[i][h]+dis[h][j]){
					dis[i][j]=dis[j][i]=dis[i][h]+dis[h][j];
				}
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=0;j<=m;j++){
			f[i][j][0]=f[i][j][1]=10000005;
		}
	}
	f[1][0][0]=f[1][1][1]=0;
	for(int i=2;i<=n;i++){
		for(int j=0;j<=m;j++){
			f[i][j][0]=min(f[i-1][j][1]+k[i-1]*dis[d[i-1]][c[i]]+1-(k[i-1])*dis[c[i-1]][c[i]],f[i-1][j][0]+dis[c[i-1]][c[i]]);
			if(j!=0){
				f[i][j][1]=min(f[i-1][j-1][1]+k[i-1]*k[i]*dis[d[i-1]][d[i]]+k[i-1]*(1-k[i])*dis[d[i-1]][c[i]]+(1-k[i-1])*k[i]*dis[c[i-1]][d[i]]+(1-k[i-1])*(1-k[i])*dis[c[i-1]][c[i]],f[i-1][j-1][0]+k[i]*dis[c[i-1]][d[i]]+(1-k[i])*dis[c[i-1]][c[i]]);
			}
		}
	}
	double ans=10000005;
	for(int i=0;i<=m;i++){
		for(int j=0;j<=1;j++){
			ans=min(ans,f[n][i][j]);
		}
	}
	printf("%.2lf",ans);
	return 0;
}