#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 |