#include<bits/stdc++.h>
using namespace std;
#define int long long
#define LD long double
const LD eps = 1e-8;
const LD INF = 1e18;
const int MOD = ((1LL<<31)-1);

int sgn(LD x){return fabs(x)<=eps ? 0 : (x>eps ? 1 : -1);}
int getNext(int x){ return (x*233811181LL+1)%MOD;}

const int N = 2005;
struct Pt{
	LD x, y;
	Pt operator-(Pt b)const{return Pt{x-b.x, y-b.y};}
	Pt operator+(Pt b)const{return Pt{x+b.x, y+b.y};}
	Pt operator*(LD b)const{return Pt{x*b, y*b};}
	LD crs(Pt b)const{return x*b.y-y*b.x;}
	
	LD len2()const{return x*x+y*y;}
	LD len()const{return sqrt(x*x+y*y);}
	Pt rot90()const{return Pt{-y, x};}
}pt[N];

struct Cir{
	Pt c; LD r;	
	bool cover(Pt p){return sgn(r - (p-c).len())>=0;}
};

Pt pt_l_l(Pt p1, Pt v1, Pt p2, Pt v2){
	return p1 + v1*(1.0L*v2.crs(p1-p2)/v1.crs(v2));																												
}

Cir getCir(Pt a, Pt b, Pt c){
	Pt o = pt_l_l((b+a)*0.5L, (b-a).rot90(), (c+a)*0.5L, (c-a).rot90());
	LD r = (o-a).len();
	return Cir{o, r};
}

int n, k, xx[N], yy[N];
LD cost[N][N], dp[N];
vector<int> vec[N];

signed main(){
//	ios::sync_with_stdio(0); cin.tie(0);
//	cin >> n >> k >> xx[1];
	scanf("%lld%lld%lld", &n, &k, &xx[1]);
	yy[1] = getNext(xx[1]);
	for (int i=2; i<=n; ++i) xx[i]=getNext(yy[i-1]), yy[i]=getNext(xx[i]);
//	for (int i=1; i<=n; ++i) printf("(%lld %lld)\n", xx[i], yy[i]);
	for (int i=1; i<=n; ++i) pt[i]={1.0L*xx[i], 1.0L*yy[i]};
	
	for (int r=1; r<=n; ++r){
		Cir cir={pt[r], 0.0L};
		cost[r][r]=0.0L;
		for (int i=r-1; i>0; --i){
			if (!cir.cover(pt[i])){
				vec[r].push_back(i);
				cir={pt[i], 0.0L};
				for (int j=r; j>i; --j){
					if (cir.cover(pt[j])) continue;
					
					cir={(pt[i]+pt[j])*0.5L, (pt[i]-pt[j]).len()*0.5L};
					for (int k=r; k>j; --k){
						if (cir.cover(pt[k])) continue;
						cir=getCir(pt[i], pt[j], pt[k]);
					}
				}
//				printf("cost[%lld][%lld]=%Lf\n", i, r, cir.r);
			}
			cost[i][r]=cir.r;
		}
		vec[r].push_back(0);
	}
	
//	for (int r=1; r<=n; ++r){
//		printf("r=%lld:", r); for (int x : vec[r]) printf("%lld ", x); puts("");
//	}
//	
	for (int j=1; j<=n; ++j) dp[j]=INF;
	for (int i=1; i<=k; ++i){
		for (int j=n; j>0; --j){
			for (int x : vec[j]) dp[j] = min(dp[j], dp[x]+cost[x+1][j]);
		}
//		for (int j=1; j<=n; ++j) printf("dp[%lld][%lld]=%Lf\n", i, j, dp[j]);
	}
	printf("%.10Lf\n", dp[n]);
	return 0;	
}