Submission #799354

#	Time	Username	Problem	Language	Result	Execution time	Memory
799354		dreamoon	Meteors (POI11_met)	C++14	74 / 100	1269 ms	60492 KiB

met

#include <bits/stdc++.h>
#define ll long long
#define SZ(X) (int)(X.size())
#define MAX 300005
#define pii pair <int,ll>
#define ind first
#define val second
#define INF (ll)(1e18)
using namespace std;
 
struct Query{
    int le, ri, amo;
};
 
int typeNum, len, queNum;
int arr[MAX];
int ans[MAX];
vector <int> sta[MAX]; // stations of every type
vector <Query> allQ;
vector <pii> allT;
int a, b, c;
 
void solve(int le, int ri, vector <pii> & type){ // type: station types and how much is needed, pos = le
    //cout << le << " " << ri << " " << SZ(type) << "\n";
    if (type.empty()) return;
    vector <pii> loc; // considered locations and dif
    vector <pii> tCopy; // copy of type
    for (pii t : type){
        for (int s : sta[t.ind])
            loc.push_back({s, 0});
        tCopy.push_back(t);
    }
    sort(loc.begin(), loc.end());
    if (loc.empty()) return;
    
    int mid = le + (ri-le+2)/2;
    for (int i = le; i < mid; i++){ // set up prefix difference
        Query q = allQ[i];
        int add = q.le, sub = q.ri;
        if (q.le > q.ri) { // ring-shaped array
            loc[0].val += q.amo;
        }
        auto it = lower_bound(loc.begin(), loc.end(), make_pair(add, -INF));
        if (it != loc.end()) it->val += q.amo;
        it = lower_bound(loc.begin(), loc.end(), make_pair(sub+1, -INF));
        if (it != loc.end()) it->val -= q.amo;
    }
    ll sum = 0;
    for (pii l : loc){ // use prefix difference to calculate how much more we need for each type
        sum += l.val;
        auto it = lower_bound(tCopy.begin(), tCopy.end(), make_pair(arr[l.ind], -INF));
        it->val -= sum;
        //cout << sum << " ";
    }
    
    if (le == ri){ // this is the final step
        for (pii t : tCopy){
            if (t.val <= 0) ans[t.ind] = le+1;
        }
        return;
    }
    
    vector <pii> typeA, typeB;
    for (int i = 0; i < SZ(type); i++){
        if (tCopy[i].val <= 0) typeA.push_back(type[i]);
        else typeB.push_back(tCopy[i]);
    }
    {
        vector <pii> tmp;
        loc.swap(tmp);
    }
    {
        vector <pii> tmp;
        tCopy.swap(tmp);
    }
    solve(le, mid-1, typeA);
    solve(mid, ri, typeB);
}
 
int main() {
    ios::sync_with_stdio(0); cin.tie(0);
    cin >> typeNum >> len;
    for (int i = 1; i <= len; i++) {
        cin >> arr[i];
        sta[arr[i]].push_back(i);
    }
    for (int i = 1; i <= typeNum; i++){
        cin >> a;
        allT.push_back({i, a});
    }
    cin >> queNum;
    for (int i = 0; i < queNum; i++){
        cin >> a >> b >> c;
        allQ.push_back({a, b, c});
    }
    solve(0, queNum-1, allT);
    for (int i = 1; i <= typeNum; i++){
        if (ans[i] == 0) cout << "NIE\n";
        else cout << ans[i] << "\n";
    }
    
    return 0;
}
/*
 10 11
 1 2 3 4 5 6 7 8 9 10 11
 1 2 3 4 5 6 7 8 9 10
 10
 1 3 1
 3 4 1
 7 1 1
 1 5 1
 1 5 1
 6 6 1
 7 7 1
 8 8 1
 9 9 1
 10 10 1
 */

• Subtask 1: 12 / 12
• Subtask 2: 12 / 12
• Subtask 3: 12 / 12
• Subtask 4: 12 / 12
• Subtask 5: 13 / 13
• Subtask 6: 13 / 13
• Subtask 7: 0 / 13
• Subtask 8: 0 / 13