Submission #971798

# Submission time Handle Problem Language Result Execution time Memory
971798 2024-04-29T10:14:15 Z Cookie Winter Driving (CCO19_day1problem3) C++14
25 / 25
104 ms 59860 KB
#include<bits/stdc++.h>
#include<fstream>
using namespace std;
#define sz(a) (int)a.size()
#define ALL(v) v.begin(), v.end()
#define ALLR(v) v.rbegin(), v.rend()
#define ll long long
#define pb push_back
#define forr(i, a, b) for(int i = a; i < b; i++)
#define dorr(i, a, b) for(int i = a; i >= b; i--)
#define ld long double
#define vt vector
#include<fstream>
#define fi first
#define se second
#define pll pair<ll, ll>
#define pii pair<int, int>
#define mpp make_pair
const ld PI = 3.14159265359, prec = 1e-9;;
//using u128 = __uint128_t;
//const int x[4] = {1, 0, -1, 0};
//const int y[4] = {0, -1, 0, 1};
const ll mod =1e9 + 7, pr = 31;
const int mxn = 1e6 + 5, mxq = 1e5 + 5, sq = 500, mxv = 5e4 + 1;
//const int base = (1 <<18);
const ll inf = 1e13 + 5, neg = -69420, inf2 = 1e14;
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
// have fun!
int n;
ll a[mxn + 1], dp[mxn + 1], siz[mxn + 1];
int p[mxn + 1];
vt<int>adj[mxn + 1];
void dfs(int s, int pre, ll dep = 0){
    siz[s] = a[s];
    for(auto i: adj[s]){
        if(i != pre){
            dfs(i, s, dep + 1); siz[s] += siz[i];
        }
    }
    dp[1] += a[s] * (siz[s] - a[s]);
}
void dfs2(int s, int pre){
    for(auto i: adj[s]){
        if(i != pre){
            dp[i] = dp[s] + a[i] * (siz[1] - siz[i]) - siz[i] * a[s];
            dfs2(i, s);
        }
    }
}
void get(int l, int r, ll tot, vt<ll>&v, vt<ll>&cand){
    if(l == r){
        v.pb(tot); return;
    }
    get(l + 1, r, tot + cand[l], v, cand); 
    get(l + 1, r, tot, v, cand);
}
void solve(){
    cin >> n;
    ll tot = 0;
    for(int i = 1; i <= n; i++){
        cin >> a[i]; tot += (a[i] * (a[i] - 1));
    }
    for(int i = 2; i <= n; i++){
        cin >> p[i]; adj[p[i]].pb(i); adj[i].pb(p[i]);
    }
    dfs(1, -1);
    dfs2(1, -1);
    ll ans = 0;
    for(int i = 1; i <= n; i++){
        //cout << dp[i] << " ";
        vt<ll>cand;
        for(auto j: adj[i]){
            if(j == p[i]){
                cand.pb(siz[1] - siz[i]);
            }else{
                cand.pb(siz[j]);
            }
        }
        ll mx = *max_element(ALL(cand));
        ll all = siz[1] - a[i];
        if(mx >= all / 2){
            ans = max(ans, mx * (all - mx) + dp[i]);
        }else{
            // centroid
            vt<ll>fir, sec; 
            get(0, sz(cand) / 2, 0, fir, cand); 
            get(sz(cand) / 2, sz(cand), 0, sec, cand);
            sort(ALL(fir)); sort(ALL(sec));
            int rp = sz(sec) - 1;
            for(int j = 0; j < sz(fir); j++){
                while(rp >= 0 && fir[j] + sec[rp] >= all / 2){
                    ll tot = fir[j] + sec[rp];
                    ans = max(ans, tot * (all - tot) + dp[i]);
                    rp--;
                }
                if(rp != -1){
                    ll tot = fir[j] + sec[rp];
                    ans = max(ans, tot * (all - tot) + dp[i]);
                }
            }
        }
    }
    
    cout << ans + tot << "\n";
}



signed main(){
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    //freopen("THREE.inp", "r", stdin);
    //freopen("THREE.out", "w", stdout);
    int tt; tt = 1;
    while(tt--){
        solve();

    }
    return(0);
}
# Verdict Execution time Memory Grader output
1 Correct 10 ms 31320 KB Output is correct
2 Correct 6 ms 31324 KB Output is correct
3 Correct 8 ms 31320 KB Output is correct
4 Correct 7 ms 31324 KB Output is correct
5 Correct 9 ms 31324 KB Output is correct
6 Correct 7 ms 31324 KB Output is correct
7 Correct 7 ms 31324 KB Output is correct
8 Correct 7 ms 31344 KB Output is correct
9 Correct 7 ms 31324 KB Output is correct
10 Correct 7 ms 31324 KB Output is correct
11 Correct 7 ms 31408 KB Output is correct
12 Correct 7 ms 31324 KB Output is correct
13 Correct 7 ms 31324 KB Output is correct
14 Correct 6 ms 31324 KB Output is correct
15 Correct 7 ms 31324 KB Output is correct
16 Correct 6 ms 31324 KB Output is correct
17 Correct 7 ms 31320 KB Output is correct
18 Correct 7 ms 31324 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 7 ms 31324 KB Output is correct
2 Correct 8 ms 31324 KB Output is correct
3 Correct 7 ms 31324 KB Output is correct
4 Correct 7 ms 31460 KB Output is correct
5 Correct 6 ms 31320 KB Output is correct
6 Correct 7 ms 31320 KB Output is correct
7 Correct 7 ms 31320 KB Output is correct
8 Correct 8 ms 31576 KB Output is correct
9 Correct 7 ms 31320 KB Output is correct
10 Correct 7 ms 31324 KB Output is correct
11 Correct 7 ms 31352 KB Output is correct
12 Correct 7 ms 31324 KB Output is correct
13 Correct 7 ms 31324 KB Output is correct
14 Correct 7 ms 31320 KB Output is correct
15 Correct 7 ms 31464 KB Output is correct
16 Correct 7 ms 31464 KB Output is correct
17 Correct 7 ms 31324 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 7 ms 31324 KB Output is correct
2 Correct 8 ms 31324 KB Output is correct
3 Correct 7 ms 31324 KB Output is correct
4 Correct 7 ms 31460 KB Output is correct
5 Correct 6 ms 31320 KB Output is correct
6 Correct 7 ms 31320 KB Output is correct
7 Correct 7 ms 31320 KB Output is correct
8 Correct 8 ms 31576 KB Output is correct
9 Correct 7 ms 31320 KB Output is correct
10 Correct 7 ms 31324 KB Output is correct
11 Correct 7 ms 31352 KB Output is correct
12 Correct 7 ms 31324 KB Output is correct
13 Correct 7 ms 31324 KB Output is correct
14 Correct 7 ms 31320 KB Output is correct
15 Correct 7 ms 31464 KB Output is correct
16 Correct 7 ms 31464 KB Output is correct
17 Correct 7 ms 31324 KB Output is correct
18 Correct 84 ms 41248 KB Output is correct
19 Correct 69 ms 40672 KB Output is correct
20 Correct 83 ms 41376 KB Output is correct
21 Correct 104 ms 41300 KB Output is correct
22 Correct 102 ms 41456 KB Output is correct
23 Correct 100 ms 41300 KB Output is correct
24 Correct 62 ms 40020 KB Output is correct
25 Correct 62 ms 40020 KB Output is correct
26 Correct 91 ms 40224 KB Output is correct
27 Correct 66 ms 40016 KB Output is correct
28 Correct 75 ms 41060 KB Output is correct
29 Correct 83 ms 44156 KB Output is correct
30 Correct 47 ms 44880 KB Output is correct
31 Correct 55 ms 44940 KB Output is correct
32 Correct 51 ms 45012 KB Output is correct
33 Correct 47 ms 44828 KB Output is correct
34 Correct 47 ms 44884 KB Output is correct
35 Correct 46 ms 44832 KB Output is correct
36 Correct 52 ms 44884 KB Output is correct
37 Correct 47 ms 45000 KB Output is correct
38 Correct 70 ms 59732 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 45 ms 36560 KB Output is correct
2 Correct 46 ms 36556 KB Output is correct
3 Correct 45 ms 36560 KB Output is correct
4 Correct 45 ms 36560 KB Output is correct
5 Correct 45 ms 36560 KB Output is correct
6 Correct 45 ms 36556 KB Output is correct
7 Correct 52 ms 36560 KB Output is correct
8 Correct 45 ms 36560 KB Output is correct
9 Correct 49 ms 36560 KB Output is correct
10 Correct 50 ms 36336 KB Output is correct
11 Correct 46 ms 36556 KB Output is correct
12 Correct 45 ms 36560 KB Output is correct
13 Correct 45 ms 36560 KB Output is correct
14 Correct 46 ms 36556 KB Output is correct
15 Correct 45 ms 36760 KB Output is correct
16 Correct 45 ms 36556 KB Output is correct
17 Correct 44 ms 36556 KB Output is correct
18 Correct 45 ms 36556 KB Output is correct
19 Correct 46 ms 36560 KB Output is correct
20 Correct 44 ms 36556 KB Output is correct
21 Correct 44 ms 36556 KB Output is correct
22 Correct 46 ms 36560 KB Output is correct
23 Correct 44 ms 36560 KB Output is correct
24 Correct 52 ms 36556 KB Output is correct
25 Correct 45 ms 36560 KB Output is correct
26 Correct 45 ms 36556 KB Output is correct
27 Correct 45 ms 36556 KB Output is correct
28 Correct 48 ms 36556 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 10 ms 31320 KB Output is correct
2 Correct 6 ms 31324 KB Output is correct
3 Correct 8 ms 31320 KB Output is correct
4 Correct 7 ms 31324 KB Output is correct
5 Correct 9 ms 31324 KB Output is correct
6 Correct 7 ms 31324 KB Output is correct
7 Correct 7 ms 31324 KB Output is correct
8 Correct 7 ms 31344 KB Output is correct
9 Correct 7 ms 31324 KB Output is correct
10 Correct 7 ms 31324 KB Output is correct
11 Correct 7 ms 31408 KB Output is correct
12 Correct 7 ms 31324 KB Output is correct
13 Correct 7 ms 31324 KB Output is correct
14 Correct 6 ms 31324 KB Output is correct
15 Correct 7 ms 31324 KB Output is correct
16 Correct 6 ms 31324 KB Output is correct
17 Correct 7 ms 31320 KB Output is correct
18 Correct 7 ms 31324 KB Output is correct
19 Correct 7 ms 31324 KB Output is correct
20 Correct 8 ms 31324 KB Output is correct
21 Correct 7 ms 31324 KB Output is correct
22 Correct 7 ms 31460 KB Output is correct
23 Correct 6 ms 31320 KB Output is correct
24 Correct 7 ms 31320 KB Output is correct
25 Correct 7 ms 31320 KB Output is correct
26 Correct 8 ms 31576 KB Output is correct
27 Correct 7 ms 31320 KB Output is correct
28 Correct 7 ms 31324 KB Output is correct
29 Correct 7 ms 31352 KB Output is correct
30 Correct 7 ms 31324 KB Output is correct
31 Correct 7 ms 31324 KB Output is correct
32 Correct 7 ms 31320 KB Output is correct
33 Correct 7 ms 31464 KB Output is correct
34 Correct 7 ms 31464 KB Output is correct
35 Correct 7 ms 31324 KB Output is correct
36 Correct 84 ms 41248 KB Output is correct
37 Correct 69 ms 40672 KB Output is correct
38 Correct 83 ms 41376 KB Output is correct
39 Correct 104 ms 41300 KB Output is correct
40 Correct 102 ms 41456 KB Output is correct
41 Correct 100 ms 41300 KB Output is correct
42 Correct 62 ms 40020 KB Output is correct
43 Correct 62 ms 40020 KB Output is correct
44 Correct 91 ms 40224 KB Output is correct
45 Correct 66 ms 40016 KB Output is correct
46 Correct 75 ms 41060 KB Output is correct
47 Correct 83 ms 44156 KB Output is correct
48 Correct 47 ms 44880 KB Output is correct
49 Correct 55 ms 44940 KB Output is correct
50 Correct 51 ms 45012 KB Output is correct
51 Correct 47 ms 44828 KB Output is correct
52 Correct 47 ms 44884 KB Output is correct
53 Correct 46 ms 44832 KB Output is correct
54 Correct 52 ms 44884 KB Output is correct
55 Correct 47 ms 45000 KB Output is correct
56 Correct 70 ms 59732 KB Output is correct
57 Correct 45 ms 36560 KB Output is correct
58 Correct 46 ms 36556 KB Output is correct
59 Correct 45 ms 36560 KB Output is correct
60 Correct 45 ms 36560 KB Output is correct
61 Correct 45 ms 36560 KB Output is correct
62 Correct 45 ms 36556 KB Output is correct
63 Correct 52 ms 36560 KB Output is correct
64 Correct 45 ms 36560 KB Output is correct
65 Correct 49 ms 36560 KB Output is correct
66 Correct 50 ms 36336 KB Output is correct
67 Correct 46 ms 36556 KB Output is correct
68 Correct 45 ms 36560 KB Output is correct
69 Correct 45 ms 36560 KB Output is correct
70 Correct 46 ms 36556 KB Output is correct
71 Correct 45 ms 36760 KB Output is correct
72 Correct 45 ms 36556 KB Output is correct
73 Correct 44 ms 36556 KB Output is correct
74 Correct 45 ms 36556 KB Output is correct
75 Correct 46 ms 36560 KB Output is correct
76 Correct 44 ms 36556 KB Output is correct
77 Correct 44 ms 36556 KB Output is correct
78 Correct 46 ms 36560 KB Output is correct
79 Correct 44 ms 36560 KB Output is correct
80 Correct 52 ms 36556 KB Output is correct
81 Correct 45 ms 36560 KB Output is correct
82 Correct 45 ms 36556 KB Output is correct
83 Correct 45 ms 36556 KB Output is correct
84 Correct 48 ms 36556 KB Output is correct
85 Correct 76 ms 41300 KB Output is correct
86 Correct 72 ms 40528 KB Output is correct
87 Correct 84 ms 41556 KB Output is correct
88 Correct 85 ms 41296 KB Output is correct
89 Correct 76 ms 41296 KB Output is correct
90 Correct 81 ms 41504 KB Output is correct
91 Correct 65 ms 40028 KB Output is correct
92 Correct 71 ms 40020 KB Output is correct
93 Correct 69 ms 40148 KB Output is correct
94 Correct 61 ms 40164 KB Output is correct
95 Correct 74 ms 41056 KB Output is correct
96 Correct 94 ms 43856 KB Output is correct
97 Correct 82 ms 50028 KB Output is correct
98 Correct 82 ms 50120 KB Output is correct
99 Correct 83 ms 50116 KB Output is correct
100 Correct 80 ms 50116 KB Output is correct
101 Correct 95 ms 50136 KB Output is correct
102 Correct 82 ms 50136 KB Output is correct
103 Correct 84 ms 50112 KB Output is correct
104 Correct 88 ms 50096 KB Output is correct
105 Correct 84 ms 49948 KB Output is correct
106 Correct 81 ms 50112 KB Output is correct
107 Correct 82 ms 50164 KB Output is correct
108 Correct 83 ms 50124 KB Output is correct
109 Correct 81 ms 50116 KB Output is correct
110 Correct 82 ms 50116 KB Output is correct
111 Correct 84 ms 49980 KB Output is correct
112 Correct 81 ms 50116 KB Output is correct
113 Correct 74 ms 59860 KB Output is correct