Submission #240718

# Submission time Handle Problem Language Result Execution time Memory
240718 2020-06-21T00:36:17 Z WhaleVomit Lightning Conductor (POI11_pio) C++14
45 / 100
1000 ms 13048 KB
#include <bits/stdc++.h>
 
using namespace std;
 
#define f first
#define s second
#define pb push_back
#define mp make_pair
#define all(v) v.begin(), v.end()
#define sz(v) (int)v.size()
 
#define MOO(i, a, b) for(int i=a; i<b; i++)
#define M00(i, a) for(int i=0; i<a; i++)
#define MOOd(i,a,b) for(int i = (b)-1; i >= a; i--)
#define M00d(i,a) for(int i = (a)-1; i>=0; i--)
 
#define FAST ios::sync_with_stdio(0); cin.tie(0);
#define finish(x) return cout << x << '\n', 0;
#define dbg(x) cerr << ">>> " << #x << " = " << x << "\n";
#define _ << " _ " <<
 
typedef long long ll;
typedef long double ld;
typedef vector<int> vi;
typedef pair<int,int> pi;
typedef pair<ld,ld> pd;
typedef complex<ld> cd;

const int MAX_N = 500500;
int n;

ld eval(ld x, pair<ll,ll> curve) {
    return curve.s + pow(curve.f - x, 0.5);
}

ll evalLong(ll x, pair<ll,ll> curve) { // ceil(eval(x, curve))
    ld y = eval(x,curve);
    ll y1 = y;
    if((y1-curve.s)*(y1-curve.s) >= curve.f-x) return y1;
    return y1+1;
}

ld intersect(pair<ll,ll> c1, pair<ll,ll> c2) { // c2 comes after c1
    if(eval(c1.f, c2) <= c1.s) return c1.f;
    if(eval(0, c2) >= eval(0, c1)) return -1;
    ld lo = 0;
    ld hi = c1.f;
    while(hi-lo > 1e-8) {
        ld mid = (lo + hi)/2;
        if(eval(mid, c1) < eval(mid, c2)) {
            hi = mid;
        } else {
            lo = mid;
        }
    }
    return lo;
}

vector<ll> solve(vector<ll> v) {
    deque<pair<ll,ll>> dq;
    MOO(i, 1, n) {
        pair<ll,ll> curve = mp(i, v[i]);
        while(1) {
            if(dq.empty()) break;
            else if(intersect(dq[sz(dq)-1], curve) < 0) dq.pop_back();
            else {
                if(sz(dq) >= 2) {
                    if(intersect(dq[sz(dq)-2], curve) < intersect(dq[sz(dq)-2], dq[sz(dq)-1])) {
                        dq.pop_back();
                    } else {
                        break;
                    }
                } else {
                    break;
                }
            }
        }
        dq.pb(curve);
    }
    vector<ll> ans;
    M00(i, n-1) {
        while(dq[0].f <= i) dq.pop_front();
        while(sz(dq) >= 2 && eval(i, dq[1]) > eval(i, dq[0])) dq.pop_front();
        ans.pb(max(0LL,evalLong(i, dq[0])-v[i]));
    }
    ans.pb(-1);
    return ans;
}
 
int main() { FAST
    cin >> n;
    vector<ll> arr(n);
    M00(i, n) cin >> arr[i];
    vector<ll> res1 = solve(arr);
    reverse(all(arr));
    vector<ll> res2 = solve(arr);
    reverse(all(res2));
    M00(i, n) {
        cout << max(res1[i], res2[i]) << "\n";
    }
}

# Verdict Execution time Memory Grader output
1 Correct 6 ms 384 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 6 ms 384 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 32 ms 384 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 261 ms 1860 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 432 ms 2920 KB Output is correct
# Verdict Execution time Memory Grader output
1 Execution timed out 1085 ms 2808 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1096 ms 4056 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1090 ms 6136 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1092 ms 10448 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1094 ms 13048 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1090 ms 13048 KB Time limit exceeded
2 Halted 0 ms 0 KB -