#include <bits/stdc++.h>
#define rep(i,j,k) for(int i=j;i<=(k);++i)
#define MAXN 1000005
typedef long double ld;
typedef long long ll;
int n,m,st,ed;
ll a[MAXN],sum[MAXN];
std::pair<ld,int> f[MAXN];
struct data{int pos,l,r;}seq[MAXN];
inline std::pair<ld,int> getval(int x,int y,ld k){
    return std::make_pair(f[y].first+sqrtl((x-y)*(sum[x]-sum[y]))+k,f[y].second-1);
}
inline std::pair<ld,int> calc(ld k){
    st=1,ed=0;
    seq[++ed]=((data){0,1,n});
    rep(i,1,n){
        while(seq[st].r<i) ++st;
        f[i]=getval(i,seq[st].pos,k);
        while(st<=ed&&getval(seq[ed].l,i,k)<=getval(seq[ed].l,seq[ed].pos,k)) --ed;
        if(st<=ed){
            int tl=seq[ed].l,tr=seq[ed].r;
            while(tl<tr){
                int tmid=tl+tr+1>>1;
                if(getval(tmid,i,k)<=getval(tmid,seq[ed].pos,k)) tr=tmid-1;
                else tl=tmid;
            }
            seq[ed].r=tl;
            if(tl!=n) seq[++ed]=(data){i,tl+1,n};
        }
        else seq[++ed]=(data){1,i+1,n};
    }
    return std::make_pair(f[n].first,-f[n].second);
}
inline ld solve(){
    ld l=0,r=1e10;
    while(r-l>1e-15){
        ld mid=(l+r)/2;
        if(calc(mid).second>=m) l=mid;
        else r=mid;
    }
    return calc((l+r)/2).first+f[n].second*(l+r)/2;
}
int main(){
    std::ios::sync_with_stdio(0);
    std::cin.tie(0),std::cout.tie(0);
    std::cin>>n>>m;
    rep(i,1,n) std::cin>>a[i];
    std::sort(a+1,a+n+1);
    rep(i,1,n) sum[i]=sum[i-1]+a[i];
    std::cout<<std::fixed<<std::setprecision(10)<<solve()<<'\n';
}