#include <iostream>
#include <cmath>
#include <algorithm>
#include <iomanip>
using namespace std;
int v[200005];
//int opt[200005];
long long spar[200005];
pair<long double,int> dp[200005];
int n;
int check(long double plm)
{
    long double valus;
    dp[0]={0,0};
    int opt=0;
    for(int i=1;i<=n;i++)
    {
        dp[i]={1e9,0};
        while (opt+1<i && dp[opt].first+sqrtl((spar[i]-spar[opt])*(i-opt))>=dp[opt+1].first+sqrtl((spar[i]-spar[opt+1])*(i-opt-1)))
        {
            opt++;
        }
        dp[i].first=dp[opt].first+sqrtl((spar[i]-spar[opt])*(i-opt));
        dp[i].second=dp[opt].second+1;
        dp[i].first+=plm;
    }
    return dp[n].second;
}
int main()
{
    int i,j,k,l,m,val;
    long double st=-1e9,dr=1e9,mij;
    cin>>n>>m;
    for(i=1;i<=n;i++)
    {
        cin>>v[i];
    }
    sort(v+1,v+n+1);
    for(i=1;i<=n;i++)
    {
        spar[i]=spar[i-1]+v[i];
    }
    for(i=1;i<=65;i++)
    {
        mij=(st+dr)/2;
        val=check(mij);
        if(val==m)
        {
            break;
        }
        else if(val>m)
        {
            st=mij;
        }
        else
        {
            dr=mij;
        }
    }
    cout<<fixed<<setprecision(12)<<dp[n].first-m*mij;
    return 0;
}