9999

#include<bits/stdc++.h>
int N,n,m,p;
int h[8192],s[8192];//8192=2<<13
double f[16][8192];
int g[16][8192];
int a[400];
using namespace std;
#define P pair<int,double>
double operator%(const P&a,const P&b){
    return (b.second-a.second)/(b.first-a.first);
}
void div(int x){
    long long q=0;
    for(int i=0;i<=p;i++){
        q=q*1000000000+a[i];a[i]=q/x;q%=x;
    }
}
int main(){
    scanf("%d%d%d",&N,&m,&p);p=p/9+1;
    scanf("%d",&h[n=1]);
    for(int i=2;i<=N;i++){
        int t;scanf("%d",&t);if(h[1]<t)h[++n]=t;
    }
    sort(h+1,h+n+1);
    if(m>n)m=n;
    for(int i=1;i<=n;i++)f[0][i]=h[1],s[i]=s[i-1]+h[i];
    int l=14;if(l>m)l=m;
    for(int i=1;i<=l;i++){
        static P q[8024];
        for(int j=2,l=1,r=0;j<=n;j++){
            P c=P(j-2,s[j-1]-f[i-1][j-1]);
            for(;l<r&&(q[r-1]%q[r])>(q[r]%c);r--);
            q[++r]=c;
            P t=P(j,s[j]);
            for(;l<r&&(q[l]%t)<(q[l+1]%t);l++);
            g[i][j]=q[l].first+1;
            f[i][j]=(q[l]%t);
        }
    }
    int _[16];_[l]=n-(m-l);
    for(int i=l;i;i--)_[i-1]=g[i][_[i]];
    a[0]=h[1];
    for(int i=1;i<=l;i++)a[0]+=s[_[i]]-s[_[i-1]],div(_[i]-_[i-1]+1);
    for(int i=_[l]+1;i<=n;i++)a[0]+=h[i],div(2);
    printf("%d.",a[0]);
    for(int i=1;i<=p;i++)printf("%09d",a[i]);
}