BZOJ 1604 [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居 Treap

题意:链接

方法: 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);
}