Submission #1043007

# Submission time Handle Problem Language Result Execution time Memory
1043007 2024-08-03T17:42:42 Z vjudge1 Meteors (POI11_met) C++17
100 / 100
5255 ms 65536 KB
#include<bits/stdc++.h>
#include <unordered_map>

#pragma GCC target("avx2")
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

using namespace std;

#define fast ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define file(s) freopen(s".in", "r", stdin);freopen(s".out", "w", stdout);
#define all(v) v.begin(),v.end()
#define rall(v) v.rbegin(),v.rend()
#define pb push_back
#define pf push_front
#define ppb pop_back
#define ft first
#define sd second
#define ll long long
#define ull unsigned ll
#define pii pair<int,int>
#define pll pair<ll,ll>
const int N = 3e5 + 10;
const int M = 1e6 + 5;
const int B = 316;
const ll msize = 2;
const ll mod1 = 1e9 + 7;
const ll mod2 = 998244353;
const long double Phi = acos(-1);
const long long inf = 2e9;
const vector <pair<int, int>> dir = {{-1, 0}, {1, 0}, {0, -1}, {0, 1}};

ll binmul(ll x, ll ti, ll m);
ll binpow(ll x, ll ti, ll m);

struct meteor {
    ll l, r, x;
} a[N];

struct Government {
    ll l, r, ans, id, need;
};

ll n, m, k;
ll t[N * 4], fl[N * 4];
vector <ll> gov[N];

void push(ll v, ll tl, ll tr) {
    t[v] += fl[v] * (tr - tl + 1);
    t[v] = min(t[v], inf);
    if (tl < tr) {
        fl[v + v] += fl[v];
        fl[v + v] = min(fl[v + v], inf);
        fl[v + v + 1] += fl[v];
        fl[v + v + 1] = min(fl[v + v + 1], inf);
    }
    fl[v] = 0;
}

void upd(ll l, ll r, ll x, ll v = 1, ll tl = 1, ll tr = m) {
    push(v, tl, tr);
    if (r < tl || tr < l) return;
    if (l <= tl && tr <= r) {
        fl[v] += x;
        fl[v] = min(fl[v], inf);
        push(v, tl, tr);
        return;
    }
    ll tm = (tl + tr) / 2;
    upd(l, r, x, v + v, tl, tm);
    upd(l, r, x, v + v + 1, tm + 1, tr);
}

ll get(ll pos, ll v = 1, ll tl = 1, ll tr = m) {
    push(v, tl, tr);
    if (tl == tr) return t[v];
    int tm = (tl + tr) / 2;
    if (pos <= tm) return get(pos, v + v, tl, tm);
    else return get(pos, v + v + 1, tm + 1, tr);
}

const void solve() {
    cin >> n >> m;
    for (ll i = 1; i <= m; ++ i) {
        int id; cin >> id;
        gov[id].pb(i);
    }

    vector <Government> p;
    p.pb({0, 0, 0, 0});
    for (ll i = 1; i <= n; ++ i) {
        ll x; cin >> x;
        p.pb({0, 0, -1, i, x});
    }

    cin >> k;
    for (ll i = 1; i <= k; ++ i) {
        cin >> a[i].l >> a[i].r >> a[i].x;
    }

    for (ll i = 1; i <= n; ++ i) p[i].l = 1, p[i].r = k;

    for (ll i = 1; i <= 20; ++ i) {
        vector <ll> g[k + 1];
        for (ll j = 1; j <= m * 4; ++ j) t[j] = fl[j] = 0;

        for (ll j = 1; j <= n; ++ j) {
            int mid = (p[j].l + p[j].r) / 2;
            g[mid].pb(p[j].id);
        }

        for (ll j = 1; j <= k; ++ j) {
            if (a[j].l <= a[j].r) {
                upd(a[j].l, a[j].r, a[j].x);
            } else {
                upd(a[j].l, m, a[j].x);
                upd(1, a[j].r, a[j].x);
            }

            for (auto id : g[j]) {
                ll sum = 0;
                for (auto location : gov[id]) {
                    sum += get(location);
                    sum = min(sum, inf);
                }

                ll mid = (p[id].l + p[id].r) / 2;
                if (sum >= p[id].need) {
                    p[id].r = mid - 1;
                    p[id].ans = mid;
                } else {
                    p[id].l = mid + 1;
                }
            }

        }
    }

    for (ll i = 1; i <= n; ++ i) {
        if (p[i].ans == -1) cout << "NIE\n";
        else cout << p[i].ans << "\n";
    }
}

int main() {
    ios_base::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);
    srand(time(NULL));

    // file("promote");

    // freopen("input.txt", "r", stdin);
    // freopen("output.txt", "w", stdout);

    int tt = 1;
    // cin >> tt;
    for (int i = 1; i <= tt; ++ i) {
        // cout << "Caso #" << i << "\n";
        solve();
    }

    #ifndef ONLINE_JUDGE
    cerr << "\n" << 1.0 * clock() / CLOCKS_PER_SEC << " s.\n";
    #endif

    return 0;
}


// Template functions
ll binmul(ll x, ll ti, ll m) { ll res = 0;while (ti){if(ti & 1)res += x;x += x;ti >>= 1; x %= m; res %= m;} return res;}
ll binpow(ll x, ll ti, ll m) { ll res = 1;while (ti){if(ti & 1)(res *= x)%=m;(x*=x)%=m;ti >>= 1; x %= m; res %= m;} return res;}
# Verdict Execution time Memory Grader output
1 Correct 7 ms 7516 KB Output is correct
2 Correct 9 ms 7716 KB Output is correct
3 Correct 8 ms 7692 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 7 ms 7512 KB Output is correct
2 Correct 7 ms 7708 KB Output is correct
3 Correct 8 ms 7528 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 395 ms 14808 KB Output is correct
2 Correct 467 ms 16584 KB Output is correct
3 Correct 466 ms 15692 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 442 ms 15328 KB Output is correct
2 Correct 437 ms 15436 KB Output is correct
3 Correct 498 ms 16552 KB Output is correct
4 Correct 119 ms 12888 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 228 ms 14740 KB Output is correct
2 Correct 353 ms 16336 KB Output is correct
3 Correct 191 ms 10332 KB Output is correct
4 Correct 446 ms 15996 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 441 ms 14292 KB Output is correct
2 Correct 397 ms 15344 KB Output is correct
3 Correct 405 ms 14796 KB Output is correct
4 Correct 479 ms 16836 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3783 ms 59712 KB Output is correct
2 Correct 2218 ms 50508 KB Output is correct
3 Correct 959 ms 21840 KB Output is correct
4 Correct 4685 ms 65536 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 4379 ms 58740 KB Output is correct
2 Correct 1655 ms 49192 KB Output is correct
3 Correct 762 ms 20380 KB Output is correct
4 Correct 5255 ms 65536 KB Output is correct