Submission #1078294

# Submission time Handle Problem Language Result Execution time Memory
1078294 2024-08-27T14:58:11 Z GrindMachine Bitaro's travel (JOI23_travel) C++17
5 / 100
3000 ms 84184 KB
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

using namespace std;
using namespace __gnu_pbds;

template<typename T> using Tree = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
typedef long long int ll;
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

#define fastio ios_base::sync_with_stdio(false); cin.tie(NULL)
#define pb push_back
#define endl '\n'
#define sz(a) (int)a.size()
#define setbits(x) __builtin_popcountll(x)
#define ff first
#define ss second
#define conts continue
#define ceil2(x,y) ((x+y-1)/(y))
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
#define yes cout << "Yes" << endl
#define no cout << "No" << endl

#define rep(i,n) for(int i = 0; i < n; ++i)
#define rep1(i,n) for(int i = 1; i <= n; ++i)
#define rev(i,s,e) for(int i = s; i >= e; --i)
#define trav(i,a) for(auto &i : a)

template<typename T>
void amin(T &a, T b) {
    a = min(a,b);
}

template<typename T>
void amax(T &a, T b) {
    a = max(a,b);
}

#ifdef LOCAL
#include "debug.h"
#else
#define debug(...) 42
#endif

/*

read the solution a long time ago, remember the key ideas

*/

const int MOD = 1e9 + 7;
const int N = 1e5 + 5;
const int inf1 = int(1e9) + 5;
const ll inf2 = ll(1e18) + 5;

template<typename T>
struct sparse_table {
    /*============================*/

    T merge(T a, T b) {
        return max(a,b);
    }

    /*============================*/

    vector<vector<T>> table;
    vector<int> bin_log;
    int LOG = 0;

    sparse_table() {

    }

    sparse_table(vector<T> &a, int n) {
        while ((1 << LOG) <= n) LOG++;

        table = vector<vector<T>>(n, vector<T>(LOG));
        bin_log = vector<int>(n + 1);

        rep(i, n) table[i][0] = a[i];

        rep1(j, LOG - 1) {
            rep(i, n) {
                int jump = 1 << (j - 1);
                if (i + jump >= n) {
                    break;
                }

                table[i][j] = merge(table[i][j - 1], table[i + jump][j - 1]);
            }
        }

        bin_log[1] = 0;
        for (int i = 2; i <= n; ++i) {
            bin_log[i] = bin_log[i / 2] + 1;
        }
    }

    T query(int l, int r) {
        int len = r - l + 1;
        int k = bin_log[len];

        T val1 = table[l][k];
        T val2 = table[r - (1 << k) + 1][k];

        return merge(val1, val2);
    }
};

void solve(int test_case)
{
    ll n; cin >> n;
    vector<ll> a(n);
    rep(i,n) cin >> a[i];
    
    vector<ll> val1(n);
    rep(i,n-1) val1[i] = 2*a[i+1]-a[i];

    vector<ll> val2(n);
    rep1(i,n-1) val2[i] = a[i]-2*a[i-1];

    sparse_table<ll> st1(val1,n), st2(val2,n);

    vector<ll> b(n);
    rep(i,n){
        ll l = i, r = i, pos = i;
        ll cost = 0;

        while(!(l == 0 and r == n-1)){
            ll left_cost = inf2, right_cost = inf2;
            if(l-1 >= 0){
                left_cost = a[pos]-a[l-1];
            }
            if(r+1 < n){
                right_cost = a[r+1]-a[pos];
            }

            if(left_cost <= right_cost){
                cost += left_cost;
                ll lo = 0, hi = l-2;
                ll best = l-1;
                ll x = inf2;
                if(r+1 < n) x = a[r+1];

                while(lo <= hi){
                    ll mid = (lo+hi) >> 1;
                    if(st1.query(mid,l-2) <= x){
                        best = mid;
                        hi = mid-1;
                    }
                    else{
                        lo = mid+1;
                    }
                }

                cost += a[l-1]-a[best];
                pos = best;
                l = best;
            }
            else{
                cost += right_cost;
                ll lo = r+2, hi = n-1;
                ll best = r+1;
                ll x = -inf2;
                if(l-1 >= 0) x = a[l-1];

                while(lo <= hi){
                    ll mid = (lo+hi) >> 1;
                    if(st1.query(r+2,mid) < -x){
                        best = mid;
                        lo = mid+1;
                    }
                    else{
                        hi = mid-1;
                    }
                }

                cost += a[best]-a[r+1];
                pos = best;
                r = pos;
            }
        }

        b[i] = cost;
    }

    ll q; cin >> q;
    while(q--){
        ll x; cin >> x;
        ll p = lower_bound(all(a),x)-a.begin();
        ll right_cost = inf2, left_cost = inf2;
        if(p < n){
            right_cost = abs(a[p]-x);
        }
        if(p-1 >= 0){
            left_cost = abs(a[p-1]-x);
        }

        ll ans = min(left_cost,right_cost);

        if(left_cost <= right_cost){
            ans += b[p-1];
        }
        else{
            ans += b[p];
        }

        cout << ans << endl;
    }
}

int main()
{
    fastio;

    int t = 1;
    // cin >> t;

    rep1(i, t) {
        solve(i);
    }

    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 22 ms 860 KB Output is correct
3 Correct 13 ms 856 KB Output is correct
4 Correct 13 ms 1024 KB Output is correct
5 Correct 13 ms 864 KB Output is correct
6 Correct 38 ms 860 KB Output is correct
7 Correct 37 ms 860 KB Output is correct
8 Correct 1 ms 344 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 344 KB Output is correct
11 Correct 0 ms 348 KB Output is correct
12 Correct 0 ms 348 KB Output is correct
13 Correct 0 ms 348 KB Output is correct
14 Correct 19 ms 860 KB Output is correct
15 Correct 24 ms 860 KB Output is correct
16 Correct 20 ms 860 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 22 ms 860 KB Output is correct
3 Correct 13 ms 856 KB Output is correct
4 Correct 13 ms 1024 KB Output is correct
5 Correct 13 ms 864 KB Output is correct
6 Correct 38 ms 860 KB Output is correct
7 Correct 37 ms 860 KB Output is correct
8 Correct 1 ms 344 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 344 KB Output is correct
11 Correct 0 ms 348 KB Output is correct
12 Correct 0 ms 348 KB Output is correct
13 Correct 0 ms 348 KB Output is correct
14 Correct 19 ms 860 KB Output is correct
15 Correct 24 ms 860 KB Output is correct
16 Correct 20 ms 860 KB Output is correct
17 Execution timed out 3035 ms 81996 KB Time limit exceeded
18 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 0 ms 348 KB Output is correct
5 Correct 0 ms 348 KB Output is correct
6 Correct 1 ms 348 KB Output is correct
7 Correct 1457 ms 83992 KB Output is correct
8 Correct 1494 ms 84000 KB Output is correct
9 Correct 1502 ms 84052 KB Output is correct
10 Correct 1448 ms 84000 KB Output is correct
11 Correct 1583 ms 84184 KB Output is correct
12 Correct 1581 ms 84004 KB Output is correct
13 Correct 31 ms 4180 KB Output is correct
14 Correct 19 ms 1560 KB Output is correct
15 Execution timed out 3071 ms 81888 KB Time limit exceeded
16 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 22 ms 860 KB Output is correct
3 Correct 13 ms 856 KB Output is correct
4 Correct 13 ms 1024 KB Output is correct
5 Correct 13 ms 864 KB Output is correct
6 Correct 38 ms 860 KB Output is correct
7 Correct 37 ms 860 KB Output is correct
8 Correct 1 ms 344 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 344 KB Output is correct
11 Correct 0 ms 348 KB Output is correct
12 Correct 0 ms 348 KB Output is correct
13 Correct 0 ms 348 KB Output is correct
14 Correct 19 ms 860 KB Output is correct
15 Correct 24 ms 860 KB Output is correct
16 Correct 20 ms 860 KB Output is correct
17 Execution timed out 3035 ms 81996 KB Time limit exceeded
18 Halted 0 ms 0 KB -