#include #include #include #include #include #define repeat(i,n) for (int i = 0; (i) < int(n); ++(i)) #define whole(f,x,...) ([&](decltype((x)) whole) { return (f)(begin(whole), end(whole), ## __VA_ARGS__); })(x) using namespace std; int main() { // input int n; cin >> n; vector l(n); vector > g(n); vector g_inv(n); repeat (i,n) { int s; cin >> l[i] >> s; -- s; g[s].push_back(i); g_inv[i] = s; } // solve double ans = whole(accumulate, l, 0) / 2.0; vector used(n); function use = [&](int i) { if (used[i]) return; used[i] = true; for (int j : g[i]) if (not used[j]) use(j); }; function &)> roots = [&](int i, vector & acc) { auto it = whole(find, acc, i); if (it != acc.end()) { acc.erase(acc.begin(), it); } else { acc.push_back(i); roots(g_inv[i], acc); } }; repeat (i,n) if (not used[i]) { vector acc; roots(i, acc); int j = *whole(min_element, acc, [&](int i, int j) { return l[i] < l[j]; }); use(j); ans += l[j] / 2.0; } // output printf("%.1f\n", ans); return 0; }