#include <bits/stdc++.h>
using namespace std;
#define int long long
#define maxn 400050
struct Line {
int m, b;
int operator()(int x) { return m * x + b; }
} a[maxn*4];
void insert(int l, int r, Line seg, int o) {
if(l + 1 == r) {
if(seg(l) > a[o](l)) a[o] = seg;
return;
}
int mid= (l + r) >> 1, lson = o * 2 + 1, rson = o * 2 + 2;
if(a[o].m > seg.m) swap(a[o], seg);
if(a[o](mid) < seg(mid)) {
swap(a[o], seg);
insert(l, mid, seg, lson);
}
else insert(mid, r, seg, rson);
}
int query(int l, int r, int x, int o) {
if(l + 1 == r) return a[o](x);
int mid = (l + r) >> 1, lson = o * 2 + 1, rson = o * 2 + 2;
if(x < mid) return max(a[o](x), query(l, mid, x, lson));
else return max(a[o](x), query(mid, r, x, rson));
}
int32_t main() {
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
for(int i=0;i<maxn*4;i++){
a[i] = {(int)0,(int)(-10000000000000000)};
}
int no_of_input;
cin >> no_of_input;
int input;
int s = 0;
vector<int> vect1;
vector<int> vect2;
vector<int> dp(no_of_input+1,0);
for(int i=0;i<no_of_input;i++){
cin >> input;
vect1.push_back(input);
}
for(int i=0;i<no_of_input;i++){
cin >> input;
s += input;
vect2.push_back(input);
}
dp[0] = 0;
Line temp;
temp.m = 2*vect1[0];
temp.b = -(vect1[0]*vect1[0]-vect2[0]);
insert(0,maxn,temp,0);
for(int i=1;i<no_of_input;i++){
//cout << query(0,maxn,vect1[i],0) << endl;
assert(query(0,maxn,vect1[i],0)!=-10000000000000000);
dp[i] = -query(0,maxn,vect1[i],0)-vect2[i]+vect1[i]*vect1[i];
Line temp;
temp.m = 2*vect1[i];
temp.b = -(vect1[i]*vect1[i]+dp[i]);
//cout << temp.m << " " << temp.b << endl;
insert(0,maxn,temp,0);
}
/*for(int i=0;i<no_of_input;i++){
cout << dp[i] << ' ';
}
cout << endl;*/
cout << dp[no_of_input-1]+s;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
17 ms |
25216 KB |
Output is correct |
2 |
Correct |
17 ms |
25472 KB |
Output is correct |
3 |
Correct |
18 ms |
25472 KB |
Output is correct |
4 |
Correct |
19 ms |
25472 KB |
Output is correct |
5 |
Correct |
19 ms |
25472 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
66 ms |
29552 KB |
Output is correct |
2 |
Correct |
63 ms |
29552 KB |
Output is correct |
3 |
Correct |
65 ms |
29560 KB |
Output is correct |
4 |
Correct |
60 ms |
29424 KB |
Output is correct |
5 |
Correct |
59 ms |
29424 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
17 ms |
25216 KB |
Output is correct |
2 |
Correct |
17 ms |
25472 KB |
Output is correct |
3 |
Correct |
18 ms |
25472 KB |
Output is correct |
4 |
Correct |
19 ms |
25472 KB |
Output is correct |
5 |
Correct |
19 ms |
25472 KB |
Output is correct |
6 |
Correct |
66 ms |
29552 KB |
Output is correct |
7 |
Correct |
63 ms |
29552 KB |
Output is correct |
8 |
Correct |
65 ms |
29560 KB |
Output is correct |
9 |
Correct |
60 ms |
29424 KB |
Output is correct |
10 |
Correct |
59 ms |
29424 KB |
Output is correct |
11 |
Correct |
79 ms |
29680 KB |
Output is correct |
12 |
Correct |
72 ms |
29552 KB |
Output is correct |
13 |
Correct |
64 ms |
29552 KB |
Output is correct |
14 |
Correct |
75 ms |
29680 KB |
Output is correct |
15 |
Correct |
59 ms |
29468 KB |
Output is correct |
16 |
Correct |
55 ms |
29424 KB |
Output is correct |
17 |
Correct |
48 ms |
29560 KB |
Output is correct |
18 |
Correct |
52 ms |
29560 KB |
Output is correct |