# |
Submission time |
Handle |
Problem |
Language |
Result |
Execution time |
Memory |
917852 |
2024-01-29T02:05:06 Z |
CSQ31 |
New Home (APIO18_new_home) |
C++17 |
|
1371 ms |
291392 KB |
#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);
}
vector<pii>sweep;
void solve(){
vector<multiset<int>>store(k+1);
//coordinate compress
comp = {0};crd.clear();
for(auto z:sweep){
int t = z.fi;
int x = z.se;
if(t > k){
comp.pb(x);
continue;
}
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 : sweep){
int t = z.fi;
int x = z.se;
if(t > k){
t-=k+1;
if(open != k){
ans[t] = -1;
continue;
}
for(int i=1;i<=k;i++){
for(int x1 : store[i]){
ans[t] = max(ans[t],abs(x-x1));
}
}
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);
}
}*/
}
}
int main()
{
owo
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});
}
for(auto z:pt){
for(auto x:z.se){if(x.fi <= k)sweep.pb(x);}
for(auto x:z.se){if(x.fi > k)sweep.pb(x);}
}
solve();
for(pii &x:sweep)x.se = 100000000 - x.se;
//solve();
for(int i=0;i<q;i++)cout<<ans[i]<<'\n';
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
38 ms |
142428 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
38 ms |
142428 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
700 ms |
239756 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
1371 ms |
291392 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
38 ms |
142428 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
38 ms |
142428 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |