题意:链接

方法: Treap

解析:

前几道资格赛的题水的不行,这道Gold的题就够分量辣。

首先这个曼哈顿距离啥的肯定能做文章,怎么转化是个问题,自己玩了一会没玩出来,就查了查曼哈顿距离的转化,发现这个玩意转化之后就变得有思路多了,所以这数学本领还是非常重要啊=-=

先看曼哈顿距离的定义

|x1−x2|+|y1−y2|

拆绝对值

x1−x2+y1−y2或x1−x2+y2−y1

x2−x1+y1−y2或x2−x1+y2−y1

即|x1+y1−(x2+y2)|或|x1−y1−(x2−y2)|

设x1+y1为x′,x1−y1为y′

则|x1′−x2′|或|y1′−y2′|

所以原要求1转化为

max(|x1′−x2′|,|y1′−y2′|)<=c

这样的二维的东西显然排序一下降一维。

按x’排序后。维护一个x’的队列,再对y’维护一个平衡树即可了。

至于要求2,即是并查集,也就是说平衡树每一次拿出来前驱后继维护下并查集即可。

y’显然可能反复,又维护并查集我们须要拿出来标号。所以平衡树须要多维护一个no,所以再删除的时候我们要找到v与no都跟要删除的目标节点同样的节点删除。

(前驱写挫WA一次= =!

代码:

#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 100010
using namespace std;
int fa[N];
int q[N];
int tot;
int root;
int cnt[N];
struct node
{
int l,r,w,v,no,rnd,siz;
}tr[N];
void pushup(int rt)
{
tr[rt].siz=tr[tr[rt].l].siz+tr[tr[rt].r].siz+tr[rt].w;
}
void rturn(int &rt)
{
int t=tr[rt].l;
tr[rt].l=tr[t].r;
tr[t].r=rt;
tr[t].siz=tr[rt].siz;
pushup(rt);
rt=t;
}
void lturn(int &rt)
{
int t=tr[rt].r;
tr[rt].r=tr[t].l;
tr[t].l=rt;
tr[t].siz=tr[rt].siz;
pushup(rt);
rt=t;
}
void insert(int &rt,int v,int no)
{
if(!rt)
{
rt=++tot;
tr[rt].siz=1,tr[rt].no=no,tr[rt].rnd=rand();
tr[rt].v=v,tr[rt].w=1;
return;
}
tr[rt].siz++;
if(v<=tr[rt].v)
{
insert(tr[rt].l,v,no);
if(tr[tr[rt].l].rnd<tr[rt].rnd)rturn(rt);
}else
{
insert(tr[rt].r,v,no);
if(tr[tr[rt].r].rnd<tr[rt].rnd)lturn(rt);
}
}
void del(int &rt,int v,int no)
{
if(!rt)return;
tr[rt].siz--;
if(tr[rt].v==v&&tr[rt].no==no)
{
if(tr[rt].l*tr[rt].r==0){rt=tr[rt].l+tr[rt].r;return;}
else
{
if(tr[tr[rt].l].rnd<tr[tr[rt].r].rnd)
{
rturn(rt);
del(rt,v,no);
}else
{
lturn(rt);
del(rt,v,no);
}
}
}else if(v<tr[rt].v)
{
del(tr[rt].l,v,no);
}else del(tr[rt].r,v,no);
}
int ans;
void q_pre(int rt,int v)
{
if(!rt)return;
if(v>=tr[rt].v)
{
ans=rt;
q_pre(tr[rt].r,v);
}else q_pre(tr[rt].l,v);
}
void q_sub(int rt,int v)
{
if(!rt)return;
if(v<tr[rt].v)
{
ans=rt;
q_sub(tr[rt].l,v);
}else q_sub(tr[rt].r,v);
}
struct point
{
int x,y;
}pt[N];
int n,c;
int find(int x)
{
if(x!=fa[x])return fa[x]=find(fa[x]);
return x;
}
int cmp(point a,point b)
{
if(a.x==b.x)return a.y<b.y;
return a.x<b.x;
}
int main()
{
scanf("%d%d",&n,&c);
for(int i=1;i<=n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
pt[i].x=x+y,pt[i].y=x-y;
fa[i]=i;
}
sort(pt+1,pt+n+1,cmp);
int head=1,tail=0;
for(int i=1;i<=n;i++)
{
while(head<=tail&&pt[i].x-pt[q[head]].x>c)
{
del(root,pt[q[head]].y,q[head]);
head++;
}
ans=0;
q_pre(root,pt[i].y);
int tmp=ans;
if(tmp!=0)
{
if(pt[i].y-tr[tmp].v<=c)
{
if(find(i)!=find(tr[tmp].no))
{
fa[find(i)]=find(tr[tmp].no);
}
}
}
ans=0;
q_sub(root,pt[i].y);
tmp=ans;
if(tmp!=0)
{
if(tr[tmp].v-pt[i].y<=c)
{
if(find(i)!=find(tr[tmp].no))
{
fa[find(i)]=find(tr[tmp].no);
}
}
}
insert(root,pt[i].y,i);
q[++tail]=i;
}
int ma=0,print=0;
for(int i=1;i<=n;i++)
{
int fx=find(i);
if(!cnt[fx])
{
print++;
}
cnt[fx]++;
if(cnt[fx]>ma)ma=cnt[fx];
}
printf("%d %d\n",print,ma);
}

最新文章

  1. javascript 获取滚动条高度+常用js页面宽度与高度
  2. sip协议音视频性能测试
  3. 《C专家编程》第三章——分析C语言的声明
  4. python---dnspython
  5. Expression Tree Basics 表达式树原理
  6. C++宽窄字符串转换
  7. javascript数组浅谈2
  8. 制作nginx和php的rpm包
  9. Android(java)学习笔记266:Android线程形态之 IntentService
  10. 239. Sliding Window Maximum
  11. html5 做游戏 Quintus Sublime Text牛逼的神器
  12. delphi7开发webservice部属在apache服务器中 转
  13. Qt中OpenGL的初步使用
  14. 【C#】Smtp发送邮件
  15. Python内置函数(47)——vars
  16. Java8学习笔记(四)--接口增强
  17. Linux指令之netstat
  18. js把预定义的html字符串转换为 HTML 实体 htmlspecialchars 输出html实体内容包括标签,而不自动转义标签,只显示内容,类似php的htmlspecialchars
  19. 智能家居中的物联网网关的可信计算平台模块(TPM)设计
  20. android入门问题--R文件丢失

热门文章

  1. [NOIP2015提高组]子串
  2. jumpserver 安装python 报错
  3. 01-JS起步
  4. 分别改动Cube每一个面的贴图UV(Unity3D开发之十八)
  5. PHP Apache shutdown unexpectedly启动错误解释及解决的方法
  6. java实习生的成长之路&lt;转&gt;
  7. jquery常规选择器再学习_1123
  8. Linux 串口终端调试工具minicom
  9. selenium 窗口句柄之间的切换
  10. Yeslab华为安全HCIE七门之--防火墙高级技术(17篇)