int n, m, x[14], y[14], p, q; double w[13], dp[1 << 13][13]; double dist(int i, int j) { return abs(x[i] - x[j]) + abs(y[i] - y[j]); } template void chmin(T& x, const T& v) { if (x > v) x = v; } { rd(p, q, n, (x, y, w)(n)); m = 1 << n; x[n] = p; y[n] = q; double w_sum = sum(w(n)), k = (w_sum + 100) / 120; rep(i, m) rep(j, n) dp[i][j] = ll_inf; rep(i, n) dp[1 << i][i] = dist(i, n) * k + w[i]; rep(s, m) { k = 0; rep(v, n) if (!(s & 1 << v)) k += w[v]; k = (k + 100) / 120; rep(u, n) if (!(s & 1 << u)) { rep(v, n) if (s & 1 << v) { chmin(dp[s | 1 << u][u], dp[s][v] + dist(v, u) * k + w[u]); } } } double ans = ll_inf; k = 5.0 / 6.0; rep(v, n) chmin(ans, dp[m - 1][v] + k * dist(v, n)); wt(ans); }