【PA2014】【BZOJ3719】Plemiona

Description

远古时代，在吉丽王国的版图上分布着n个部落。建立平面直角坐标系后，每个部落都是一个边平行于坐标轴的矩形。有些地盘可能同时属于多个部落。随着时间推移，部落之间会发生融合。具体来说，若两个部落的公共面积严格大于零，它们会合并成一个新的部落，新部落的形状是包含原来两个部落的最小矩形（边平行于坐标轴）。
数百万年后，部落之间终于达到了稳定状态（任两个部落都不能再合并了），然而吉丽也已经老了。他想知道最终还剩下几个部落，以及各个部落的位置。你能替他完成遗业吗？

Input

第一行一个整数n(1<=n<=100000)，表示远古时代的部落数量。
接下来n行，每行四个整数x1,x2,y1,y2(0<=x1< x2<=1000000,0<=y1< y2<=1000000)，表示部落的坐标。

Output

第一行输出一个整数m，表示稳定后还剩下的部落数量。
接下来m行，每行四个整数x1,x2,y1,y2，表示部落的坐标。请按照字典序（先比较x1，若x1相等则比较x2，以此类推）从小到大输出。

Sample Input

5

7 8 1 4

1 5 2 3

4 5 2 7

2 3 5 9

4 6 8 9
Sample Output

2

1 6 2 9

7 8 1 4
HINT

Source

鸣谢Jcvb

因为坐标范围很小,所以考虑对x坐标建线段树
线段树每个节点维护两个vector,配合并查集来维护矩形的合并信息
代码又长..跑的又慢…
光荣垫底了..

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define MAXN 100010
#define SIZE 1000010
#define GET (ch>='0'&&ch<='9')
#define lchild rt<<1,l,mid
#define rchild rt<<1|1,mid+1,r
#define ln rt<<1
#define rn rt<<1|1
using namespace std;
const int size=1000000;
int n,ans;
int st[MAXN],sta[MAXN],top;
struct Matrix
{
    int x1,y1,x2,y2;
    inline bool operator <(const Matrix& a)const    {   return x1==a.x1?(x2==a.x2?(y1==a.y1?y2<a.y2:y1<a.y1):x2<a.x2):x1<a.x1;  }
}s[MAXN],Ans[MAXN];
int f[MAXN];
int find(int x) {   return f[x]==x?x:f[x]=find(f[x]);   }
void Union(int a,int b)
{
    int x=find(a),y=find(b);
    if (x==y)   return;
    s[a].x1=min(s[a].x1,s[b].x1);s[a].x2=max(s[a].x2,s[b].x2);
    s[a].y1=min(s[a].y1,s[b].y1);s[a].y2=max(s[a].y2,s[b].y2);
    f[find(b)]=a;
}
inline bool cmp(int a,int b)    {   return s[a].y2==s[b].y2?a<b:s[a].y2<s[b].y2;    }
struct seg  {   vector<int> q1,q2;  }tree[SIZE<<2];
void insert(int rt,int l,int r,int L,int R)
{
    if (l>R||r<L)   return;
    if (l>=L&&r<=R) {   sta[++top]=rt;return;   }
    int mid=(l+r)>>1;insert(lchild,L,R);insert(rchild,L,R);
}
void update(int rt,int x)   {   for (int i=rt;i;i>>=1)  tree[i].q1.push_back(x);    }
int find_rt(int rt,int l,int r,int x)
{
    if (l==r)   return rt;int mid=(l+r)>>1;
    return x<=mid?find_rt(lchild,x):find_rt(rchild,x);
}
inline void in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
}
int main()
{
    in(n);
    for (int i=1;i<=n;i++)
        in(s[i].x1),in(s[i].x2),in(s[i].y1),in(s[i].y2),
        s[i].x1++,s[i].y1++,f[i]=st[i]=i;
    sort(st+1,st+n+1,cmp);
    for (int i=1;i<=n;i++)
    {
        bool flag=1;
        while (flag)
        {
            flag=0;top=0;insert(1,1,size,s[st[i]].x1,s[st[i]].x2);
            for (int j=1;j<=top;++j)    for (int k=tree[sta[j]].q1.size()-1;~k;--k)
            {
                int y=tree[sta[j]].q1[k];
                if (y!=f[y]||find(y)==find(st[i]))  {   tree[sta[j]].q1.pop_back();continue;    }
                if (s[y].y2<s[st[i]].y1)    break;
                Union(st[i],y);flag=1;tree[sta[j]].q1.pop_back();
            }
            for (int now=find_rt(1,1,size,s[st[i]].x1);now;now>>=1) for (int k=tree[now].q2.size()-1;~k;--k)
            {
                int y=tree[now].q2[k];
                if (y!=f[y]||find(y)==find(st[i]))  {   tree[now].q2.pop_back();continue;   }
                if (s[y].y2<s[st[i]].y1)    break;
                Union(st[i],y);flag=1;tree[now].q2.pop_back();
            }
            for (int now=find_rt(1,1,size,s[st[i]].x2);now;now>>=1) for (int k=tree[now].q2.size()-1;~k;--k)
            {
                int y=tree[now].q2[k];
                if (y!=f[y]||find(y)==find(st[i]))  {   tree[now].q2.pop_back();break;  }
                if (s[y].y2<s[st[i]].y1)    break;
                Union(st[i],y);flag=1;tree[now].q2.pop_back();
            }
        }
        update(find_rt(1,1,size,s[st[i]].x1),st[i]);
        if (s[st[i]].x1!=s[st[i]].x2)   update(find_rt(1,1,size,s[st[i]].x2),st[i]);
        insert(1,1,size,s[st[i]].x1,s[st[i]].x2);
        for (int j=1;j<=top;++j)    tree[sta[j]].q2.push_back(st[i]);
    }
    for (int i=1;i<=n;++i)  if (f[i]==i)    Ans[++ans]=s[i],Ans[ans].x1--,Ans[ans].y1--;
    printf("%d\n",ans);sort(Ans+1,Ans+ans+1);
    for (int i=1;i<=ans;++i)    printf("%d %d %d %d\n",Ans[i].x1,Ans[i].x2,Ans[i].y1,Ans[i].y2);
}