Language: C++

#include<iostream>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<functional>
#include<iomanip>
#include<queue>
#include<ciso646>
#include<random>
#include<map>
#include<set>
#include<complex>
#include<bitset>
#include<stack>
#include<unordered_map>
#include<utility>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const ll mod = 1000000007;
const ll INF = (ll)1000000007 * 1000000007;
typedef pair<int, int> P;
#define stop char nyaa;cin>>nyaa;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
typedef long double ld;
typedef complex<ld> Point;
const ld eps = 1e-6;
const ld pi = acos(-1.0);
typedef pair<ll, ll> LP;
int n, m;
double x[14], y[14];
double d[14][14];
double rd[1 << 14];
bool can(int s,ld r) {
	vector<int> v;
	rep(i, n) {
		if (s&(1 << i))v.push_back(i);
	}
	if (v.size() <= 1)return true;
	int len = v.size();
	rep(i, len) {
		Rep(j, i + 1, len) {
			int id = v[i], jd = v[j];
			if (2*r < d[id][jd])return false;
			ld X = (x[jd] - x[id]) / 2.0;
			ld Y = (y[jd] - y[id]) / 2.0;
			ld dist = sqrt(pow(r, 2) - pow(d[id][jd] / 2.0, 2));
			ld sr = d[id][jd] / 2.0;
			ld x3 = Y * dist / sr; ld y3 = -X * dist / sr;
			ld midx = (x[id] + x[jd]) / 2.0, midy = (y[id] + y[jd]) / 2.0;
			ld cx = x3 + midx, cy = y3 + midy;
			bool valid = true;
			rep(k, len) {
				int kd = v[k];
				ld dis = sqrt(pow(cx - x[kd], 2) + pow(cy - y[kd], 2));
				if (dis > r+eps) {
					valid = false; break;
				}
			}
			if (valid)return true;
			valid = true;
			cx = -x3 + midx, cy = -y3 + midy;
			rep(k, len) {
				int kd = v[k];
				ld dis = sqrt(pow(cx - x[kd], 2) + pow(cy - y[kd], 2));
				if (dis > r+eps) {
					valid = false; break;
				}
			}
			if (valid)return true;
		}
	}
	return false;
}
ld dp[1 << 14], cop[1 << 14];
int main() {
	cin >> n;
    m = 1;
	rep(i, n) {
		cin >> x[i] >> y[i];
	}
	rep(i, n) {
		rep(j, n) {
			d[i][j] = sqrt(pow(x[i] - x[j], 2) + pow(y[i] - y[j], 2));
		}
	}
	rep(i, (1 << n)) {
		ld le = 0, ri = 200000;
		while (ri - le > eps) {
			ld r = (ri + le) / 2;
			if (can(i, r))ri = r;
			else le = r;
		}
		rd[i] = ri;
	}
	fill(dp, dp + (1 << n), mod);
	dp[0] = 0;
	rep(aa, m) {
		rep(j, (1 << n)) {
			cop[j] = dp[j];
		}
		rep(j, (1 << n)) {
			vector<int> v;
			rep(i, n) {
				if (j&(1 << i))continue;
				v.push_back(i);
			}
			int len = v.size();
			rep(i, (1 << len)) {
				int s = 0;
				rep(k, len) {
					if (i&(1 << k))s += (1 << v[k]);
				}
				cop[j + s] = min(cop[j + s], max(dp[j], rd[s]));
			}
		}
		rep(j, (1 << n)) {
			dp[j] = cop[j];
		}
	}
	cout << fixed << setprecision(5);
	// cout << dp[(1 << n) - 1] << endl;
    printf("%.7lf\n", dp[(1 << n) - 1]);
	//stop
	return 0;
}

Compiler output:
a.cc: In function 'int main()':
a.cc:121:65: error: no matching function for call to 'max(ld&, double&)'
  121 |                                 cop[j + s] = min(cop[j + s], max(dp[j], rd[s]));
      |                                                              ~~~^~~~~~~~~~~~~~
In file included from /usr/include/c++/14/string:51,
                 from /usr/include/c++/14/bits/locale_classes.h:40,
                 from /usr/include/c++/14/bits/ios_base.h:41,
                 from /usr/include/c++/14/ios:44,
                 from /usr/include/c++/14/ostream:40,
                 from /usr/include/c++/14/iostream:41,
                 from a.cc:1:
/usr/include/c++/14/bits/stl_algobase.h:257:5: note: candidate: 'template<class _Tp> constexpr const _Tp& std::max(const _Tp&, const _Tp&)'
  257 |     max(const _Tp& __a, const _Tp& __b)
      |     ^~~
/usr/include/c++/14/bits/stl_algobase.h:257:5: note:   template argument deduction/substitution failed:
a.cc:121:65: note:   deduced conflicting types for parameter 'const _Tp' ('long double' and 'double')
  121 |                                 cop[j + s] = min(cop[j + s], max(dp[j], rd[s]));
      |                                                              ~~~^~~~~~~~~~~~~~
/usr/include/c++/14/bits/stl_algobase.h:303:5: note: candidate: 'template<class _Tp, class _Compare> constexpr const _Tp& std::max(const _Tp&, const _Tp&, _Compare)'
  303 |     max(const _Tp& __a, const _Tp& __b, _Compare __comp)
      |     ^~~
/usr/include/c++/14/bits/stl_algobase.h:303:5: note:   candidate expects 3 arguments, 2 provided
In file included from /usr/include/c++/14/algorithm:61,
                 from a.cc:6:
/usr/include/c++/14/bits/stl_algo.h:5706:5: note: candidate: 'template<class _Tp> constexpr _Tp std::max(initializer_list<_Tp>)'
 5706 |     max(initializer_list<_Tp> __l)
      |     ^~~
/usr/include/c++/14/bits/stl_algo.h:5706:5: note:   candidate expects 1 argument, 2 provided
/usr/include/c++/14/bits/stl_algo.h:5716:5: note: candidate: 'template<class _Tp, class _Compare> constexpr _Tp std::max(initializer_list<_Tp>, _Compare)'
 5716 |     max(initializer_list<_Tp> __l, _Compare __comp)
      |     ^~~
/usr/include/c++/14/bits/stl_algo.h:5716:5: note:   template argument deduction/substitution failed:
a.cc:121:65: note:   mismatched types 'std::initializer_list<_Tp>' and 'long double'
  121 |                                 cop[j + s] = min(cop[j + s], max(dp[j], rd[s]));
      |                                                              ~~~^~~~~~~~~~~~~~