Submission #917836

#	Time	Username	Problem	Language	Result	Execution time	Memory
917836		CSQ31	New Home (APIO18_new_home)	C++17	0 / 100	5065 ms	579644 KiB

new_home

#include <bits/stdc++.h>
using namespace std;
#define pb push_back
#define fi first
#define se second
#define sz(a) (int)(a.size())
#define all(a) a.begin(),a.end()
#define lb lower_bound
#define ub upper_bound
#define owo ios_base::sync_with_stdio(0);cin.tie(0);
#define debug(...) fprintf(stderr, __VA_ARGS__),fflush(stderr)
#define time__(d) for(long blockTime = 0; (blockTime == 0 ? (blockTime=clock()) != 0 : false);\
debug("%s time : %.4fs\n", d, (double)(clock() - blockTime) / CLOCKS_PER_SEC))
typedef long long int ll;
typedef long double ld;
typedef pair<ll,ll> PII;
typedef pair<int,int> pii;
typedef vector<vector<int>> vii;
typedef vector<vector<ll>> VII;
ll gcd(ll a,ll b){if(!b)return a;else return gcd(b,a%b);}
//solving only for negative slopes
//two adjacent stores at a and b causes a line at x = (a+b+1)/2 with y = b-x
const int MAXT = 3e6,MAXN = 3e5+1;
int n,q,k;
vector<int>ans(MAXN,-1),T;
multiset<int,greater<int>>lines[MAXT];
map<int,vector<pii>>pt;
map<int,int>crd;
int N;
void upd(int v,int l,int r,int pos){
    if(l==r){
        if(lines[l].empty())T[v] = -1e9;
        else T[v] = *lines[l].begin();
        return;
    }
    int tm = (l+r)/2;
    if(pos<=tm)upd(2*v,l,tm,pos);
    else upd(2*v+1,tm+1,r,pos);
    T[v] = max(T[2*v],T[2*v+1]);
}
int query(int v,int l,int r,int tl,int tr){
    if(l>r)return -1e9;
    if(l==tl && r==tr)return T[v];
    int tm = (tl+tr)/2;
    return max(query(2*v,l,min(r,tm),tl,tm),
               query(2*v+1,max(l,tm+1),r,tm+1,tr));
}
void add(int l,int r){
    int x = (l+r+1)/2;
    int y = r - x;
    lines[crd[x]].insert(x+y);
    //upd(1,0,N,crd[x]);
}
void del(int l,int r){
    int x = (l+r+1)/2;
    int y = r - x;
    lines[crd[x]].erase(lines[crd[x]].find((x+y)));
    //upd(1,0,N,crd[x]);
}
bool debug = 0;
vector<int>comp;
void addcomp(int l,int r){
    int x = (l+r+1)/2;
    comp.pb(x);
}
void solve(){
    vector<multiset<int>>store(k+1);
    //coordinate compress
    comp = {0};crd.clear();
    for(auto z:pt){
        vector<pii>v = z.se;
        for(pii st : v){
            int t = st.fi;
            int x = st.se;
            if(t > k){
                comp.pb(x);
                continue;
            }
            //cout<<t<<" "<<x<<'\n';
            if(t > 0){ 
                //add this guy
                if(store[t].find(x) != store[t].end()){
                    store[t].insert(x);
                    continue;
                }
                auto r = store[t].lb(x);
                bool lf = (r!=store[t].begin());
                bool rg = (r!=store[t].end());
                if(lf){
                    auto l = r;
                    l--;
                    addcomp(*l,x);
                }
                if(rg)addcomp(x,*r);
                store[t].insert(x);
            }else{ 
                t*=-1;
                //remove this guy
                store[t].erase(store[t].find(x));
                if(store[t].find(x) != store[t].end())continue;
                
                auto r = store[t].lb(x);
                bool lf = (r!=store[t].begin());
                bool rg = (r!=store[t].end());
                if(lf && rg){
                    auto l = r;
                    l--;
                    addcomp(*l,*r);
                }
            }
        }
    }
    sort(all(comp));
    comp.resize(unique(all(comp)) - comp.begin());
    for(int i=0;i<sz(comp);i++)crd[comp[i]] = i;
    N = sz(comp)-1;
    T.assign(4*(N+1),-1e9);

    for(int i=1;i<=k;i++)store[i].clear();

    int open = 0;
    for(auto z:pt){
        vector<pii>v = z.se;
        for(pii st : v){
            int t = st.fi;
            int x = st.se;
            if(t > k)continue;
            //cout<<"store "<<t<<" "<<x<<'\n';
            if(t > 0){ 
                if(store[t].empty())open++;
                if(store[t].find(x) != store[t].end()){
                    store[t].insert(x);
                    continue;
                }
                auto r = store[t].lb(x);
                bool lf = (r!=store[t].begin());
                bool rg = (r!=store[t].end());
                if(!lf){
                    add(-x,x); //line at 0
                    if(rg)add(x,*r);
                }else{
                    auto l = r;
                    add(*l,x);
                    if(rg){
                        add(x,*r);
                        del(*l,*r);
                    }
                }
                store[t].insert(x);
            }else{ 
                t*=-1;
                if(sz(store[t]) == 1)open--;
                store[t].erase(store[t].find(x));
                if(store[t].find(x) != store[t].end())continue;

                auto r = store[t].lb(x);
                bool lf = (r!=store[t].begin());
                bool rg = (r!=store[t].end());
                if(lf){
                    auto l = r;
                    l--;
                    if(rg){
                        add(*l,*r);
                        del(*l,x);
                        del(x,*r);
                    }else{
                        del(*l,x);
                    }
                }else{
                    del(-x,x);
                    if(rg)add(-*r,*r);
                }
            }
            
        }
        for(pii st : v){
            int id = st.fi - k - 1;
            int x = st.se;
            if(id < 0)continue;
            if(open != k){
                ans[id] = -1;
                continue;
            }
            for(int i=0;i<=crd[x];i++){
                if(lines[i].empty())continue;
                ans[id] = max(ans[id],*lines[i].begin() - x);
            }
        }
    }


}
int main()
{
    cin>>n>>k>>q;
    for(int i=0;i<n;i++){
        int x,t,a,b;
        cin>>x>>t>>a>>b;
        pt[a].pb({t,x}); //add
        pt[b+1].pb({-t,x}); //remove
    }
    for(int i=0;i<q;i++){
        int l,y;
        cin>>l>>y;
        pt[y].pb({i+k+1,l});
    }
    solve();
    for(auto &z:pt){
        for(pii &x : z.se)x.se = 100000000 - x.se;
    }
    solve();
    for(int i=0;i<q;i++)cout<<ans[i]<<'\n';


}

• Subtask 1: 0 / 5
• Subtask 2: 0 / 7
• Subtask 3: 0 / 10
• Subtask 4: 0 / 23
• Subtask 5: 0 / 35
• Subtask 6: 0 / 20