【CF1028C】Rectangles（线段树）

题意：

n<=1e5,abs(x[i]),abs(y[i]<=1e9

思路：这是正解

离散后线段树强打，数据结构越学越傻

 #include<cstdio>
 #include<cstring>
 #include<string>
 #include<cmath>
 #include<iostream>
 #include<algorithm>
 #include<map>
 #include<set>
 #include<queue>
 #include<vector>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef pair<int,int> PII;
 typedef vector<int> VI;
 #define fi first
 #define se second
 #define MP make_pair
 #define N   1100000
 #define MOD 1000000007
 #define eps 1e-8
 #define pi acos(-1)
 struct arr1
 {
     int x1,y1,x2,y2;
 }a[N];

 struct arr2
 {
     int x,y1,y2,z;
 }c[N];

 struct tree
 {
     int a,s,x;
 }t[N<<];

 int b[N],d[N],q[N],m;

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 void swap(int &x,int &y)
 {
     int t=x;x=y;y=t;
 }

 int hash1(int x)
 {
     int l=;
     int r=m;
     while(l<=r)
     {
         int mid=(l+r)>>;
         if(b[mid]==x) return d[mid];
         if(b[mid]<x) l=mid+;
          else r=mid-;
     }
 }

 bool cmp(arr2 a,arr2 b)
 {
     return a.x<b.x;
 }

 void pushdown(int p)
 {
     int l=p<<;
     int r=l+;
     if(!t[p].a) return;
     t[l].s=t[l].s+t[p].a;
     t[l].a=t[l].a+t[p].a;
     t[r].s=t[r].s+t[p].a;
     t[r].a=t[r].a+t[p].a;
     t[p].a=;
 }

 void pushup(int p)
 {
     if(t[p<<].s<t[p+p+].s)
     {
         t[p].s=t[p+p+].s;
         t[p].x=t[p+p+].x;
     }
      else
      {
          t[p].s=t[p<<].s;
          t[p].x=t[p<<].x;
      }
 }

 void build(int l,int r,int p)
 {
     if(l>r) return;
      t[p].a=;
      t[p].x=l;
      if(l==r) return;
      int mid=(l+r)>>;
      build(l,mid,p<<);
     build(mid+,r,p<<|);
 }

 int query(int l,int r,int x,int y,int p)
 {
     if(x<=l&&r<=y) return t[p].s;
     pushdown(p);
     int mid=(l+r)>>;
     int ret=-;
     if(x<=mid) ret=query(l,mid,x,y,p<<);
     if(y>mid)
     {
         int t=query(mid+,r,x,y,p<<|);
         if(t>ret) ret=t;
     }
     return ret;
 }

 void update(int l,int r,int x,int y,int v,int p)
 {
     if(x<=l&&r<=y)
     {
         t[p].s+=v;
         t[p].a+=v;
         return;
     }
     pushdown(p);
     int mid=(l+r)>>;
     if(x<=mid) update(l,mid,x,y,v,p<<);
     if(y>mid) update(mid+,r,x,y,v,p<<|);
     pushup(p);
 }

 int main()
 {
    // freopen("cf1028C.in","r",stdin);
    // freopen("cf1028C.out","w",stdout);
     int x;
     scanf("%d",&x);
     int k1=x;
     m=;
     int n=;
     for(int i=;i<=x;i++)
     {
          scanf("%d%d%d%d",&a[i].x1,&a[i].y1,&a[i].x2,&a[i].y2);
          b[++m]=a[i].x1;
          b[++m]=a[i].y1;
          b[++m]=a[i].x2;
          b[++m]=a[i].y2;
     }
     sort(b+,b+m+);
     d[]=;
     for(int i=;i<=m;i++)
      if(b[i]==b[i-]) d[i]=d[i-];
       else d[i]=d[i-]+;
     q[]=b[];
     int len=;
     for(int i=;i<=m;i++)
      if(b[i]!=b[i-]) q[++len]=b[i];

     for(int i=;i<=x;i++)
     {
         c[++n].x=hash1(a[i].x1);
         c[n].y1=hash1(a[i].y1);
         c[n].y2=hash1(a[i].y2);
         c[n].z=;
         c[++n].x=hash1(a[i].x2)+;
         c[n].y1=hash1(a[i].y1);
         c[n].y2=hash1(a[i].y2);
         c[n].z=-;
     }
     sort(c+,c+n+,cmp);
     int line=;
     int row=;
     for(int i=;i<=n;i++)
     {
         line=max(line,c[i].x);
         row=max(row,c[i].y2);
     }
     build(,row,);

     int j=;
     for(int i=;i<=line;i++)
     {
         while(j<=n&&c[j].x==i)
         {
             update(,row,c[j].y1,c[j].y2,c[j].z,);
         //    printf("%d %d %d\n",c[j].y1,c[j].y2,c[j].z);
             j++;
         }
         int k=query(,row,,row,);
         //printf("%d %d\n",i,k);
         if(k>=k1-)
         {
         //    int p=find(1,row,1,row,1);
             int p=t[].x;
         //    printf("%d %d\n",i,p);
             int x=q[i];
             int y=q[p];
             printf("%d %d\n",x,y);
             return ;
         }
     }
     return ;
 }