Submission #1019120

# Submission time Handle Problem Language Result Execution time Memory
1019120 2024-07-10T13:44:43 Z GrindMachine Homecoming (BOI18_homecoming) C++17
100 / 100
909 ms 193352 KB
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2")

#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

/*

refs:
edi

*/

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

#include "homecoming.h"

template<typename T>
struct segtree {
    // https://codeforces.com/blog/entry/18051

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

    struct data {
        ll a[2][2];
        data(){
            rep(i,2) rep(j,2) a[i][j] = -inf2;
        }
    };

    data neutral = data();

    data merge(data &left, data &right) {
        data curr;
        
        rep(a,2){
            rep(b,2){
                amax(curr.a[a][b],left.a[a][b]);
                amax(curr.a[a][b],right.a[a][b]);

                rep(c,2){
                    amax(curr.a[a][c],left.a[a][b]+right.a[b][c]);
                }
            }
        }

        return curr;
    }

    void create(int i, T v) {
        auto [x,y] = v;
        tr[i].a[0][0] = max(x+y,0ll);
        tr[i].a[0][1] = x;
        tr[i].a[1][0] = y;
        tr[i].a[1][1] = -inf2;
    }

    void modify(int i, T v) {
        auto [x,y] = v;
        tr[i].a[0][0] = max(x+y,0ll);
        tr[i].a[0][1] = x;
        tr[i].a[1][0] = y;
        tr[i].a[1][1] = -inf2;
    }

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

    int n;
    vector<data> tr;

    segtree() {

    }

    segtree(int siz) {
        init(siz);
    }

    void init(int siz) {
        n = siz;
        tr.assign(2 * n, neutral);
    }

    void build(vector<T> &a, int siz) {
        rep(i, siz) create(i + n, a[i]);
        rev(i, n - 1, 1) tr[i] = merge(tr[i << 1], tr[i << 1 | 1]);
    }

    void pupd(int i, T v) {
        modify(i + n, v);
        for (i = (i + n) >> 1; i; i >>= 1) tr[i] = merge(tr[i << 1], tr[i << 1 | 1]);
    }

    data query(int l, int r) {
        data resl = neutral, resr = neutral;

        for (l += n, r += n; l <= r; l >>= 1, r >>= 1) {
            if (l & 1) resl = merge(resl, tr[l++]);
            if (!(r & 1)) resr = merge(tr[r--], resr);
        }

        return merge(resl, resr);
    }
};

ll solve(int n, int k, int *A, int *B){
    // all subjects win (so all books taken)
    ll all_win = 0;
    rep(i,n) all_win += A[i]-B[i];

    ll ans = max(all_win,0ll);

    // at least 1 subject doesnt win (so at least 1 book not taken)
    segtree<pll> st(n);
    ll s = 0, cost2 = 0;
    rep(i,k) s += B[i];
    vector<pll> init(n);

    rep(i,2*n){
        if(i+k-1 < 2*n){
            if(i){
                s -= B[(i-1)%n];
                s += B[(i+k-1)%n];
            }
     
            ll cost1 = A[i%n]-s;
            cost2 += A[i%n]-B[(i+k-1)%n];
            cost1 -= cost2;
            if(i < n){
                init[i] = {cost1,cost2};
            }
            else{
                st.pupd(i%n,{cost1,cost2});
            }
        }

        if(i == n-1){
            st.build(init,n);
        }

        if(i >= n){
            int f = i%n;
            int l = f+1, r = f+n-k;
            if(l > r) conts;
            l %= n, r %= n;
            ll val = 0;
     
            if(l <= r){
                val = st.query(l,r).a[0][0];
            }
            else{
                auto d1 = st.query(l,n-1);
                auto d2 = st.query(0,r);
                val = st.merge(d1,d2).a[0][0];
            }

            amax(ans,val);
        }
    }

    return ans;
}
# Verdict Execution time Memory Grader output
1 Correct 0 ms 344 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 1 ms 348 KB Output is correct
4 Correct 0 ms 348 KB Output is correct
5 Correct 0 ms 348 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 344 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 1 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 856 KB Output is correct
7 Correct 1 ms 860 KB Output is correct
8 Correct 1 ms 568 KB Output is correct
9 Correct 2 ms 856 KB Output is correct
10 Correct 1 ms 348 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 208 ms 48912 KB Output is correct
2 Correct 7 ms 860 KB Output is correct
3 Correct 902 ms 193352 KB Output is correct
4 Correct 8 ms 2396 KB Output is correct
5 Correct 34 ms 3436 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 344 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 1 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 856 KB Output is correct
7 Correct 1 ms 860 KB Output is correct
8 Correct 1 ms 568 KB Output is correct
9 Correct 2 ms 856 KB Output is correct
10 Correct 1 ms 348 KB Output is correct
11 Correct 208 ms 48912 KB Output is correct
12 Correct 7 ms 860 KB Output is correct
13 Correct 902 ms 193352 KB Output is correct
14 Correct 8 ms 2396 KB Output is correct
15 Correct 34 ms 3436 KB Output is correct
16 Correct 777 ms 192764 KB Output is correct
17 Correct 274 ms 31788 KB Output is correct
18 Correct 476 ms 33272 KB Output is correct
19 Correct 833 ms 100596 KB Output is correct
20 Correct 909 ms 184092 KB Output is correct
21 Correct 886 ms 101236 KB Output is correct