#include <bits/stdc++.h>

using namespace std;
using uint = unsigned int;
using ll = long long;
using ull = unsigned long long;
constexpr ll TEN(int n) { return (n==0) ? 1 : 10*TEN(n-1); }
template<class T> using V = vector<T>;
template<class T> using VV = V<V<T>>;

using D = double;
using P = complex<D>;

const D PI = acos(D(-1)), EPS = 1e-10;

int sgn(D a) {
    if (abs(a) <= EPS) return 0;
    return (a < 0) ? -1 : 1;
}
bool near(P a, P b) { return !sgn(abs(a-b)); }
D mabs(P p) { return abs(p.real()) + abs(p.imag()); }

bool lessP(P l, P r) {
    if (sgn(r.real()-l.real())) return l.real() < r.real();
    if (sgn(r.imag()-l.imag())) return l.imag() < r.imag();
    return false;
}

D cross(P a, P b) { return a.real()*b.imag() - a.imag()*b.real(); }
D dot(P a, P b) { return a.real()*b.real() + a.imag()*b.imag(); }

// -2, -1, 0, 1, 2 : front, clock, on, cclock, back
int ccw(P b, P c) {
    int s = sgn(cross(b, c) / mabs(b));
    if (s) return s;
    if (!sgn(mabs(c)) || !sgn(mabs(c-b))) return 0;
    if (dot(b, c) < 0) return 2;
    if (dot(-b, c-b) < 0) return -2;
}
int ccw(P a, P b, P c) { return ccw(b-a, c-a); }

struct L {
    P x, y;
    L() {}
    L(P x, P y) :x(x), y(y) {}
    P vec() const { return y-x; }
};
D abs(const L &l) { return abs(l.vec()); }

D distLP(const L &l, const P &p) {
    return abs(cross(l.vec(), p-l.x)/abs(l.vec()));
}
D distSP(const L &s, const P &p) {
    P s2 = s.vec()*P(0, 1);
    if (ccw(s.x, s.x+s2, p) == 1) return abs(s.x-p);
    if (ccw(s.y, s.y+s2, p) == -1) return abs(s.y-p);
    return min(min(abs(s.x-p), abs(s.y-p)), distLP(s, p));
}



using Pol = V<P>;
using Tri = array<P, 3>;

struct C {
    P p; D r;
    C() {}
    C(P p, D r) : p(p), r(r) {}
};

C circumCircle(P a, P b, P c) {
    b -= a; c -= a;
    D s = 2*cross(b, c);
    D x = norm(b-c), y = norm(c), z = norm(b);
    D S = x+y+z;
    P r = (y*(S-2*y)*b+z*(S-2*z)*c)/(s*s);
    return C(r + a, abs(r));
}

C smallestEnclosingCircle(Pol &p, int ps = 0, Tri q = Tri(), int qs = 0) {
    if (ps == p.size()) {
        switch (qs) {
            case 0:
                return C(P(0, 0), -1);
            case 1:
                return C(q[0], 0);
            case 2:
                return C((q[0]+q[1])/D(2.0), abs(q[0]-q[1])/D(2.0));
            case 3:
                return circumCircle(q[0], q[1], q[2]);
        }
        assert(false);
    }
    C c = smallestEnclosingCircle(p, ps+1, q, qs);
    if (sgn(abs(p[ps]-c.p) - c.r) == 1) {
        q[qs] = p[ps];
        return smallestEnclosingCircle(p, ps+1, q, qs+1);
    }
    return c;
}





const int MN = 14;
const int MM = 14;
int n, m;
P p[MN];
D pre[1<<MN];
D dp[MM+1][1<<MN];
bool vis[MM+1][1<<MN];

D solve(int f, int m) {
    if (f == 0) return 0;
    if (m == 0) return 1e100;
    if (vis[m][f]) return dp[m][f];
    vis[m][f] = true;
    D ans = 1e100;
    int g = 0;
    do {
        ans = min(ans, max(pre[g], solve(f^g, m-1)));
        g = ((g|~f)+1) & f;
    } while (g);
    return dp[m][f] = ans;
}
int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << setprecision(20);

    cin >> n >> m;
    for (int i = 0; i < n; i++) {
        D x, y;
        cin >> x >> y;
        p[i] = P(x, y);
    }
    for (int f = 0; f < (1<<n); f++) {
        V<P> v;
        for (int i = 0; i < n; i++) {
            if (!(f & (1<<i))) continue;
            v.push_back(p[i]);
        }
        pre[f] = smallestEnclosingCircle(v).r;

//        cout << f << " " << pre[f] << endl;
    }

    cout << solve((1<<n) - 1, m) << endl;
    return 0;
}



### Prompt

Please formulate a CPP solution to the following problem:

Problem

Great Demon King Megumi wants to defeat the $ N $ hero who lives on the ground.
Megumi can cast explosion magic up to $ M $ times.
Explosion magic is magic that extinguishes heroes who exist within a radius of $ r $ around arbitrary coordinates.
The hero is very thin, so you don't have to consider the size.
All $ M $ explosion spells shall be cast with the same radius.
I decided to use the magic of the minimum radius necessary to annihilate the hero.
Minimize the size of the explosive magic radius.

Constraints

The input satisfies the following conditions.

* $ 1 \ leq M \ leq N \ leq 14 $
* $ 0 \ leq x_i $, $ y_i \ leq 10 ^ 5 $

Input


$ N $ $ M $
$ x_1 $ $ y_1 $
...
$ x_N $ $ y_N $


The input is given in the following format.
The number of heroes $ N $ and the number of explosions to cast $ M $ are given as integers on the first line.
From the second line onward, the coordinates $ x_i $ $, $ $ y_i $ of each hero $ i $ are given as integers.

Output

Output the minimum value of the explosive magic radius as a real number.
Must not contain absolute error greater than $ 10 ^ {-3} $.

Examples

Input

5 2
0 0
5 5
10 10
100 100
200 200


Output

70.710678118755


Input

10 5
321 675
4312 6534
312 532
412 6543
21 43
654 321
543 0
32 5
41 76
5 1


Output

169.824909833728


Input

14 3
312 342
4893 432
321 4
389 4
23 543
0 0
1 1
2 2
432 12
435 32
1 5
2 10
100 100
20 50


Output

218.087711712613


Input

5 2
0 0
0 0
0 0
0 0
0 0


Output

0.000000000001

### Response

```cpp
#include <bits/stdc++.h>

using namespace std;
using uint = unsigned int;
using ll = long long;
using ull = unsigned long long;
constexpr ll TEN(int n) { return (n==0) ? 1 : 10*TEN(n-1); }
template<class T> using V = vector<T>;
template<class T> using VV = V<V<T>>;

using D = double;
using P = complex<D>;

const D PI = acos(D(-1)), EPS = 1e-10;

int sgn(D a) {
    if (abs(a) <= EPS) return 0;
    return (a < 0) ? -1 : 1;
}
bool near(P a, P b) { return !sgn(abs(a-b)); }
D mabs(P p) { return abs(p.real()) + abs(p.imag()); }

bool lessP(P l, P r) {
    if (sgn(r.real()-l.real())) return l.real() < r.real();
    if (sgn(r.imag()-l.imag())) return l.imag() < r.imag();
    return false;
}

D cross(P a, P b) { return a.real()*b.imag() - a.imag()*b.real(); }
D dot(P a, P b) { return a.real()*b.real() + a.imag()*b.imag(); }

// -2, -1, 0, 1, 2 : front, clock, on, cclock, back
int ccw(P b, P c) {
    int s = sgn(cross(b, c) / mabs(b));
    if (s) return s;
    if (!sgn(mabs(c)) || !sgn(mabs(c-b))) return 0;
    if (dot(b, c) < 0) return 2;
    if (dot(-b, c-b) < 0) return -2;
}
int ccw(P a, P b, P c) { return ccw(b-a, c-a); }

struct L {
    P x, y;
    L() {}
    L(P x, P y) :x(x), y(y) {}
    P vec() const { return y-x; }
};
D abs(const L &l) { return abs(l.vec()); }

D distLP(const L &l, const P &p) {
    return abs(cross(l.vec(), p-l.x)/abs(l.vec()));
}
D distSP(const L &s, const P &p) {
    P s2 = s.vec()*P(0, 1);
    if (ccw(s.x, s.x+s2, p) == 1) return abs(s.x-p);
    if (ccw(s.y, s.y+s2, p) == -1) return abs(s.y-p);
    return min(min(abs(s.x-p), abs(s.y-p)), distLP(s, p));
}



using Pol = V<P>;
using Tri = array<P, 3>;

struct C {
    P p; D r;
    C() {}
    C(P p, D r) : p(p), r(r) {}
};

C circumCircle(P a, P b, P c) {
    b -= a; c -= a;
    D s = 2*cross(b, c);
    D x = norm(b-c), y = norm(c), z = norm(b);
    D S = x+y+z;
    P r = (y*(S-2*y)*b+z*(S-2*z)*c)/(s*s);
    return C(r + a, abs(r));
}

C smallestEnclosingCircle(Pol &p, int ps = 0, Tri q = Tri(), int qs = 0) {
    if (ps == p.size()) {
        switch (qs) {
            case 0:
                return C(P(0, 0), -1);
            case 1:
                return C(q[0], 0);
            case 2:
                return C((q[0]+q[1])/D(2.0), abs(q[0]-q[1])/D(2.0));
            case 3:
                return circumCircle(q[0], q[1], q[2]);
        }
        assert(false);
    }
    C c = smallestEnclosingCircle(p, ps+1, q, qs);
    if (sgn(abs(p[ps]-c.p) - c.r) == 1) {
        q[qs] = p[ps];
        return smallestEnclosingCircle(p, ps+1, q, qs+1);
    }
    return c;
}





const int MN = 14;
const int MM = 14;
int n, m;
P p[MN];
D pre[1<<MN];
D dp[MM+1][1<<MN];
bool vis[MM+1][1<<MN];

D solve(int f, int m) {
    if (f == 0) return 0;
    if (m == 0) return 1e100;
    if (vis[m][f]) return dp[m][f];
    vis[m][f] = true;
    D ans = 1e100;
    int g = 0;
    do {
        ans = min(ans, max(pre[g], solve(f^g, m-1)));
        g = ((g|~f)+1) & f;
    } while (g);
    return dp[m][f] = ans;
}
int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << setprecision(20);

    cin >> n >> m;
    for (int i = 0; i < n; i++) {
        D x, y;
        cin >> x >> y;
        p[i] = P(x, y);
    }
    for (int f = 0; f < (1<<n); f++) {
        V<P> v;
        for (int i = 0; i < n; i++) {
            if (!(f & (1<<i))) continue;
            v.push_back(p[i]);
        }
        pre[f] = smallestEnclosingCircle(v).r;

//        cout << f << " " << pre[f] << endl;
    }

    cout << solve((1<<n) - 1, m) << endl;
    return 0;
}



```