Submission #423782

#	Time	Username	Problem	Language	Result	Execution time	Memory
423782		jtnydv25	Shortcut (IOI16_shortcut)	C++17	97 / 100	2084 ms	55912 KiB

shortcut

#include "shortcut.h"
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int, int>
#define all(c) ((c).begin()), ((c).end())
#define sz(x) ((int)(x).size())
 
#ifdef LOCAL
#include <print.h>
#else
#define trace(...)
#define endl '\n'
#endif
const int N = 1 << 20;
const ll INF = 1e18;
struct Data{
    ll a, b;
    Data() : a(-INF), b(INF){}
    Data(ll a, ll b) : a(a), b(b){}
};
 
inline void merge(Data & X, const Data & Y){
    X.a = max(X.a, Y.a);
    X.b = min(X.b,Y.b);
}
 
Data prefix[N];
ll p[N];

int perm1[N], perm2[N], where[N], position[N];
long long find_shortcut(int n, std::vector<int> l, std::vector<int> d, int c){
    for(int i = 1; i < n; i++) p[i] = p[i - 1] + l[i - 1];
    iota(perm1, perm1 + n, 0);
    iota(perm2, perm2 + n, 0);
    sort(perm1, perm1 + n, [&](int i, int j){return d[i] - p[i] > d[j] - p[j];});
    sort(perm2, perm2 + n, [&](int i, int j){return d[i] + p[i] < d[j] + p[j];});
 
    for(int i = 0; i < n; i++) where[perm1[i]] = i;
    int mx = 0, mx2 = 0, ind_mx = 0;
    for(int i = 1; i <= n; i++){
        int ind = perm1[i - 1];
        prefix[i] = prefix[i - 1];
        merge(prefix[i], Data(p[ind] + d[ind], p[ind] - d[ind]));
        if(d[ind] > mx){
            mx2 = mx;
            mx = d[ind];
            ind_mx = ind;
        } else if(d[ind] > mx2) mx2 = d[ind];
    }
    ll lo = mx + mx2, hi = 1000000000LL * 10000000LL;
    while(lo < hi){
        ll D = (lo + hi) / 2;
        Data d1, d2;
        int cur = 0;
        for(int ind = 0; ind < n; ind++){
            int i = perm2[ind];
            while(cur < n && d[perm1[cur]] + d[i] + p[i] - p[perm1[cur]] > D) cur++;
            position[i] = cur;
        }
        
        for(int j = 0; j < n; j++){
            Data temp = prefix[position[j]];
            if(j == ind_mx && D < 2 * mx){
                temp = Data();
                int k = position[j];
                for(int r = 0; r < k; r++){
                    int _r = perm1[r];
                    if(_r == j) continue;
                    merge(temp, Data(p[_r] + d[_r], p[_r] - d[_r]));
                }
            }
            d1.a = max(d1.a, temp.a + d[j] + c - p[j] - D);
            d1.b = min(d1.b, temp.b - d[j] - c - p[j] + D);
            d2.a = max(d2.a, temp.a + p[j] - D + d[j] + c);
            d2.b = min(d2.b, temp.b + p[j] + D - d[j] - c);
        }
        ll L1 = d1.a, R1 = d1.b;
        ll L2 = d2.a, R2 = d2.b;
        if(L1 > R1 || L2 > R2){
            lo = D + 1;
            continue;
        } 
        bool ok = false;
        int ptr1 = 0;
        int ptr2 = n;
        for(int v = 1; v < n; v++){
            ll R = min(p[v] + R1, R2 - p[v]);
            while(ptr1 < n && p[ptr1] < L1 + p[v]) ptr1++;
            while(ptr2 > 0 && p[ptr2 - 1] >= L2 - p[v]) ptr2--;
            int u = max(ptr1, ptr2);
 
            if(u >= v) continue;
            if(p[u] > R) continue;
            ok = true;
            break;
        }
        if(ok) hi = D;
        else lo = D + 1;
    }
    return lo;
}

• Subtask 1: 9 / 9
• Subtask 2: 14 / 14
• Subtask 3: 8 / 8
• Subtask 4: 7 / 7
• Subtask 5: 33 / 33
• Subtask 6: 22 / 22
• Subtask 7: 4 / 4
• Subtask 8: 0 / 3