#include #include #include #include #include #include #include #include #include #include #include #include namespace lib { templatestruct modnum; templateconstexpr T pow(T base, std::size_t p) { if (p == 0) { return T(1); } else if (p == 1) { return base; } else if (p == 2) { return base*base; } else if (p % 2 == 0) { return pow(pow(base, p / 2), 2); } else { return pow(pow(base, p / 2), 2)*base; } } templateconstexpr auto inverse(modnum const&); templatestruct modnum { static constexpr auto mod = Mod; std::uint64_t val; modnum() = default; constexpr modnum(std::uint64_t v) :val(v%Mod) { } constexpr modnum& operator+=(modnum const& v) { val += v.val; val %= mod; return *this; } constexpr modnum& operator-=(modnum const& v) { val += mod - v.val; val %= mod; return *this; } constexpr modnum& operator*=(modnum const& v) { val *= v.val; val %= mod; return *this; } constexpr modnum& operator/=(modnum const& v) { return operator*=(inverse(v)); } }; templateconstexpr auto operator+(modnum lhs, modnumconst& rhs) { return lhs += rhs; } templateconstexpr auto operator-(modnum lhs, modnumconst& rhs) { return lhs -= rhs; } templateconstexpr auto operator*(modnum lhs, modnumconst& rhs) { return lhs *= rhs; } templateconstexpr auto operator/(modnum lhs, modnumconst& rhs) { return lhs /= rhs; } templateconstexpr auto inverse(modnumconst& base) { return pow(base, Mod - 2); } templateconstexpr auto clamp(T v) { return std::max(v, T()); } templatevoid sort(std::vector& vec) { std::sort(vec.begin(), vec.end()); } templateauto lower_bound(std::vectorconst& vec, T v) { return std::distance(vec.begin(), std::lower_bound(vec.begin(), vec.end(), v)); } templateauto upper_bound(std::vectorconst& vec, T v) { return std::distance(vec.begin(), std::upper_bound(vec.begin(), vec.end(), v)); } templatestruct dmod { std::uint64_t val; static std::uint64_t mod; dmod(std::uint64_t v = 0) :val(v%mod) { } dmod& operator+=(dmod const& rhs) { val += rhs.val; val %= mod; return *this; } dmod& operator-=(dmod const& rhs) { val += mod - rhs.val; val %= mod; return *this; } dmod& operator*=(dmod const& rhs) { val *= rhs.val; val %= mod; return *this; } }; templatedmod operator+(dmod lhs, dmod const& rhs) { return lhs += rhs; } templatedmod operator-(dmod lhs, dmod const& rhs) { return lhs -= rhs; } templatedmod operator*(dmod lhs, dmod const& rhs) { return lhs *= rhs; } templatestd::uint64_t dmod::mod; struct scope_exit { std::function func; scope_exit(std::function f) :func(f) { } ~scope_exit() { func(); } }; templateusing int_tag = std::integral_constant; } void Main(); int main() { std::cin.tie(nullptr); std::cin.sync_with_stdio(false); Main(); } bool check_cycle(int index, int start, std::vector>const&dest) { for (auto v : dest[index]) { if (v == start) { return true; } if (check_cycle(v, start, dest)) { return true; } } return false; } void Main() { int N; std::cin >> N; std::vector data(N); std::vector> dest(N); std::vector cycle(N); std::set rest; for (int i = 0; i < N; ++i) { std::cin >> data[i]; int x; std::cin >> x; dest[x - 1].emplace(i); rest.emplace(i); } for (int i = 0; i < N; ++i) { cycle[i] = check_cycle(i, i, dest); } int ret = 0; std::vector sum(N); while (rest.size() > 0) { auto top = *rest.begin(); for (auto index : rest) { if (std::make_pair(!cycle[index], data[index]) < std::make_pair(!cycle[top], data[top])) { top = index; } } std::stack stack; stack.emplace(top); ret += data[top]; while (stack.size() > 0) { auto index = stack.top(); stack.pop(); if (!rest.count(index)) { continue; } rest.erase(index); ret += data[index]; for (auto v : dest[index]) { stack.emplace(v); } } } std::cout << ret / 2 << "." << (ret % 2 == 0 ? 0 : 5) << std::endl; }