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);
}