BZOJ 1604: [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居 曼哈顿距离转切比雪夫距离 Treap

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

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1131  Solved: 463
[Submit][Status][Discuss]

Description

了解奶牛们的人都知道,奶牛喜欢成群结队.观察约翰的N(1≤N≤100000)只奶牛,你会发现她们已经结成了几个“群”.每只奶牛在吃草的时候有一个独一无二的位置坐标Xi,Yi(l≤Xi,Yi≤[1..10^9];Xi,Yi∈整数.当满足下列两个条件之一,两只奶牛i和j是属于同一个群的:
  1.两只奶牛的曼哈顿距离不超过C(1≤C≤10^9),即lXi - xil+IYi - Yil≤C.
  2.两只奶牛有共同的邻居.即,存在一只奶牛k,使i与k,j与k均同属一个群.
    给出奶牛们的位置,请计算草原上有多少个牛群,以及最大的牛群里有多少奶牛

Input

   1行输入NC,之后N行每行输入一只奶牛的坐标.

Output

仅一行,先输出牛群数,再输出最大牛群里的牛数,用空格隔开.

Sample Input

4 2
1 1
3 3
2 2
10 10
* Line 1: A single line with a two space-separated integers: the
number of cow neighborhoods and the size of the largest cow
neighborhood.

Sample Output

2 3
OUTPUT DETAILS:
There are 2 neighborhoods, one formed by the first three cows and
the other being the last cow. The largest neighborhood therefore
has size 3.


首先觉得翻译有点问题 给个英文题面

BZOJ 1604: [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居 曼哈顿距离转切比雪夫距离 Treap_第1张图片

在昨天做过两个题之后 觉得这样的就是套路题了(这题也是昨天写的 没交 /捂脸熊)

对切比雪夫和曼哈顿的理解又加深了些

我在BZOJ 3382: [Usaco2004 Open]Cave Cows 3 洞穴里的牛之三中

说了一下曼哈顿距离的小小转化

是用去绝对值符号得到的结论

在BZOJ 3170: [Tjoi 2013]松鼠聚会 切比雪夫距离中

为了解决切比雪夫距离转成了曼哈顿距离(具体请点开链接)

再做完这道题 觉得已经可以【融会贯通】

我们得到结论,按照这个顺序做最好(一本正经的扯淡)


好了我们说题解:


题意就是只要两个点曼哈顿距离<=c那么就连边

问最后有多少联通块,最大的大小多少


看到曼哈顿距离,第一反应就是kdtree(也不管是啥题 /捂脸熊)

所以就略懵逼 然后突然想起是为了做切比雪夫才看这个题 之后就知道先转切比雪夫(mdzz)

按照上面BZOJ 3382: [Usaco2004 Open]Cave Cows 3 洞穴里的牛之三或BZOJ 3170: [Tjoi 2013]松鼠聚会 切比雪夫距离给出的方式

曼哈顿转切比雪夫都可以

就是转个坐标/开个绝对值 (x,y) 就变成了 (x+y,x-y)

然后用切比雪夫搞一搞

具体怎么做呢?

(0,0) (x,y) 的切比雪夫距离就是max(x,y)

所以 x<=c&&y<=c 时两点就连边

这时发现它二维的结构,直接排序开心降一维,单调队列扫一遍就好了

然后剩下一维的怎么搞呢?

考虑对于队列中插入一个点

因为已经保证内部连通性的正确,它只需要连接这里的前驱后继就可以了

所以写个平衡树就好了


写完就WA,一顿狂调,结果发现根本不是平衡树的锅 GG


在放弃调自己的平衡树后参考各路神犇敲的平衡树(这应该是一棵名次树??)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;

typedef double db;
typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=1001000;

struct point{int x,y;friend bool operator <(const point &x,const point &y){return x.xtr[tr[k].ls].rnd)rturn(k);
	}
	else
	{
		insert(tr[k].rs,val,pos);
		if(tr[k].rnd>tr[tr[k].rs].rnd)lturn(k);
	}
}

void del(int &k,int val,int pos)
{
	if(!k)return ;
	tr[k].sz--;
	if(tr[k].val==val&&tr[k].pos==pos)
	{
		if(tr[k].ls*tr[k].rs==0)k=tr[k].ls+tr[k].rs;
		else if(tr[tr[k].ls].rndval)
	{
		Y=k;
		query_sub(tr[k].ls,val);
	}
	else query_sub(tr[k].rs,val);
}

int q[N];

int main()
{
	n=read();int c=read();
	register int i,x,y;
	for(i=1;i<=n;++i)
	{
		x=read();y=read();
		p[i].x=x+y;p[i].y=x-y;
		fa[i]=i;
	}
	sort(p+1,p+1+n);
	register int head=1,tail=1;
	for(i=1;i<=n;++i)
	{
		while(headc)
		del(root,p[q[head]].y,q[head]),head++;
		X=Y=0;
		query_pre(root,p[i].y);
		if(X&&(p[i].y-tr[X].val<=c))fa[find(i)]=find(tr[X].pos);
		query_sub(root,p[i].y);
		if(Y&&(tr[Y].val-p[i].y<=c))fa[find(i)]=find(tr[Y].pos);
		insert(root,p[i].y,i);
		q[tail++]=i;
	}
	int ans=0,num=0;
	for(i=1;i<=n;++i)
	{
		if(!sonsize[find(i)])num++;
		sonsize[find(i)]++;ans=max(ans,sonsize[find(i)]);
	}
	cout<

最开始写的

然后最后知道是怎么错的呢?

带权并查集错了。。。单调队列错了(要记住,函数传参的时候不要用"head++"这种东东啊)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;

typedef double db;
typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=100100;

struct point{int x,y;friend bool operator <(const point &x,const point &y){return x.xtr[tr[k].ls].rnd)rturn(k);
	}
	else
	{
		if(insert(tr[k].rs,val,pos));tr[k].sz++;
		if(tr[k].rnd>tr[tr[k].rs].rnd)lturn(k);
	}
}

int del(int &k,int val,int pos)
{
	if(tr[k].val==val)
	{
		if(pos^tr[k].pos)return 0;
		if(tr[k].ls*tr[k].rs==0)k=tr[k].ls+tr[k].rs;
		else if(tr[tr[k].ls].rndtr[k].val)if(del(tr[k].rs,val,pos))tr[k].sz--;
}

int X,Y;

void query_pre(int k,int val)
{
	if(!k)return ;
	if(tr[k].val<=val)
	{
		X=tr[k].pos;
		query_pre(tr[k].rs,val);
	}
	else query_pre(tr[k].ls,val);
}

void query_sub(int k,int val)
{
	if(!k)return ;
	if(tr[k].val>=val)
	{
		Y=tr[k].pos;
		query_sub(tr[k].ls,val);
	}
	else query_sub(tr[k].rs,val);
}

int q[N];

int main()
{
	n=read();int c=read();
	register int i,x,y;
	for(i=1;i<=n;++i)
	{
		x=read();y=read();
		p[i].x=x+y;p[i].y=x-y;
		fa[i]=i;
	}
	sort(p+1,p+1+n);
	register int head=1,tail=1;
	for(i=1;i<=n;++i)
	{
		while(headc)
		del(root,p[q[head]].y,q[head]),head++;
		X=Y=0;
		query_pre(root,p[i].y);
		if(X&&(p[i].y-p[X].y<=c))fa[find(i)]=find(X);
		query_sub(root,p[i].y);
		if(Y&&(p[Y].y-p[i].y<=c))fa[find(i)]=find(Y);
		insert(root,p[i].y,i);
		q[tail++]=i;
	}
	int ans=0,num=0;
	for(i=1;i<=n;++i)
	{
		if(!sonsize[find(i)])num++;
		sonsize[find(i)]++;ans=max(ans,sonsize[find(i)]);
	}
	cout<

你可能感兴趣的:(切比雪夫距离,—————————计算几何,平衡树,—————————中级数据结构)