#include<bits/stdc++.h>
using namespace std;
int n,m,k,s;
double w[2005],sw[2005],Sw[2005],dp[2005][105];
int main(){
    cin>>n>>m>>k>>s;n+=m;
    for(int i=1;i<=n;i++) cin>>w[i];sort(w+1,w+n+1);
    for(int i=0;i<=n;i++) for(int j=0;j<=k;j++) dp[i][j]=1e18;
    for(int i=1;i<=n;i++) sw[i]=w[i]+sw[i-1],Sw[i]=w[i]*w[i]+Sw[i-1];
    dp[0][0]=0;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=k;j++){
        auto &w=dp[i][j];
        for(int k=i-1;k>=0;k--) w=min(w,dp[k][j-1]+Sw[i]-Sw[k]-(sw[i]-sw[k])*(sw[i]-sw[k])/(i-k));
    }
    cout<<fixed<<setprecision(10)<<dp[n][k]+n*s*s/4.0<<'\n';
}