#include<bits/stdc++.h>
using namespace std;
using f128=long double;
int n,m,a[200010];
f128 pres[200010];
struct N_{
   f128 v;int c;
   N_ operator+(const N_&o)const{
      return{v+o.v,c+o.c};
   }
   bool operator<(const N_&o)const{
      return v<o.v;
   }
}f[200010];
namespace ds{
int t[200010<<2];
N_ eval(int i,int k){
   return k<=i?N_{(f128)1e18,0}:f[i]+N_{sqrt((k-i)*(pres[k]-pres[i])),0};
}
void insert(int u,int l,int r,int x){
   if(t[u]==-1){t[u]=x;return;}
   if(l==r)return;
   int mid=l+r>>1;
   if(eval(x,mid)<eval(t[u],mid))swap(t[u],x);
   if(x<t[u])insert(u<<1,l,mid,x);
   else insert(u<<1|1,mid+1,r,x);
}
N_ qry(int u,int l,int r,int p){
   if(l==r)return eval(t[u],p);
   int mid=l+r>>1;
   return min(eval(t[u],p),p<=mid?qry(u<<1,l,mid,p):qry(u<<1|1,mid+1,r,p));
}
N_ work(f128 w){
   memset(t,-1,sizeof(t));
   t[1]=0;
   for(int i=1;i<=n;i++){
      f[i]=qry(1,1,n,i)+N_{-w,1},insert(1,1,n,i);
   }
   return f[n];
}
}
int main(){
   scanf("%d%d",&n,&m);
   for(int i=1;i<=n;i++)scanf("%d",&a[i]);
   sort(a+1,a+n+1);
   for(int i=1;i<=n;i++)pres[i]=pres[i-1]+a[i];
   f128 l=-1e10,r=0;
   for(int t=200;t--;){
      f128 mid=(l+r)/2;
      if(ds::work(mid).c<m)l=mid;
      else r=mid;
   }
   printf("%.20Lf",ds::work(l).v+l*m);
   return 0;
}