#include<bits/stdc++.h>

#define ll long long
#define pi pair<int,int>
#define vi vector<int>
#define cpy(x,y,s) memcpy(x,y,sizeof(x[0])*(s))
#define mset(x,v,s) memset(x,v,sizeof(x[0])*(s))
#define all(x) begin(x),end(x)
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define ary array

using namespace std;

#define N 500005

int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9')f=(ch=='-'?-1:f),ch=getchar();
    while(ch>='0' && ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    return x*f;
}
void write(int x){
    if(x<0)x=-x,putchar('-');
    if(x/10)write(x/10);
    putchar(x%10+'0');
}

int n,m;
int a[N],s[N];
struct tup{
    int x,l,r;
    tup(int x=0,int l=0,int r=0):x(x),l(l),r(r){}
};

deque<tup> q;
long double f[N];
int g[N];

long double w(int pos,int r){
    return sqrt((r-pos)*(s[r]-s[pos]));
}
long double calc(long double k){
    while(q.size())q.pop_front();
    q.push_front(tup(0,1,n));
    for(int i=1;i<=n;i++){
        int pos=q.front().x;
        f[i]=f[pos]+w(pos,i)+k;g[i]=g[pos]+1;
        if(q.front().l==q.front().r)q.pop_front();
        else{
            auto [x,l,r]=q.front();q.pop_front();
            q.push_front(tup(x,l+1,r));
        }
        while(q.size() && f[i]+w(i,q.back().l)<=f[q.back().x]+w(q.back().x,q.back().l))q.pop_back();
        if(q.size()){
            if(f[i]+w(i,q.back().r)<=f[q.back().x]+w(q.back().x,q.back().r)){
                pos=q.back().x;int l=q.back().l,r=q.back().r;
                while(l<r){
                    int mid=(l+r)>>1;
                    if(f[i]+w(i,mid)<=f[pos]+w(pos,mid))r=mid;
                    else l=mid+1;
                }
                auto [X,L,R]=q.back();q.pop_back();
                q.push_back(tup(X,L,l-1));
            }
        }
        
        if(q.size() && q.back().r!=n)q.push_back(tup(i,q.back().r+1,n));
        if(!q.size())q.push_back(tup(i,i+1,n));
    }
    return f[n];
}

int main(){
    #ifdef EAST_CLOUD
    freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    #endif

    n=read(),m=read();
    for(int i=1;i<=n;i++)a[i]=read();
    sort(a+1,a+n+1);
    for(int i=1;i<=n;i++)s[i]=s[i-1]+a[i];
    long double L=-0x3f3f3f3f,R=0x3f3f3f3f,eps=1e-12;
    for(int i=1;i<=n;i++)
    while(R-L>eps){
        long double mid=(L+R)/2;calc(mid);
        //cerr<<mid<<' '<<g[n]<<endl;
        if(g[n]<=m)R=mid;
        else L=mid;
    }
    printf("%.12Lf",calc(L)-L*m);
    return 0;
}