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

最新推荐文章于 2023-09-21 14:12:42 发布

BlackJack_

最新推荐文章于 2023-09-21 14:12:42 发布

阅读量526

点赞数 1

分类专栏：切比雪夫距离 —————————计算几何平衡树 —————————中级数据结构

本文链接：https://blog.csdn.net/BlackJack_/article/details/78270419

版权

—————————中级数据结构同时被 3 个专栏收录

70 篇文章 0 订阅

订阅专栏

—————————计算几何

22 篇文章 0 订阅

订阅专栏

平衡树

18 篇文章 0 订阅

订阅专栏

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行输入N和C，之后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 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<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
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.x<y.x;}}p[N];

int fa[N],sonsize[N];

inline int find(int x)
{
	int t=x;
	while(t^fa[t])t=fa[t];
	while(x^fa[x])x=fa[x],fa[x]=t;
	return x;
}

int n,root,sz;
struct balanced_tree{int ls,rs,sz,val,pos,rnd;}tr[N];

inline void lturn(int &k)
{
	int t=tr[k].rs;tr[k].rs=tr[t].ls;tr[t].ls=k;
	tr[t].sz=tr[k].sz;tr[k].sz=tr[tr[k].ls].sz+tr[tr[k].rs].sz+1;k=t;
}

inline void rturn(int &k)
{
	int t=tr[k].ls;tr[k].ls=tr[t].rs;tr[t].rs=k;
	tr[t].sz=tr[k].sz;tr[k].sz=tr[tr[k].ls].sz+tr[tr[k].rs].sz+1;k=t;
}

void insert(int &k,int val,int pos)
{
	if(!k)
	{
		k=++sz;tr[k].sz=1;tr[k].val=val;
		tr[k].pos=pos;tr[k].rnd=rand();
		return ;
	}
	tr[k].sz++;
	if(val<=tr[k].val)
	{
		insert(tr[k].ls,val,pos);
		if(tr[k].rnd>tr[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].rnd<tr[tr[k].rs].rnd)rturn(k),del(k,val,pos);
		else lturn(k),del(k,val,pos);
		return ;
	}
	val<tr[k].val?del(tr[k].ls,val,pos):del(tr[k].rs,val,pos);
}

int X,Y;

void query_pre(int k,int val)
{
	if(!k)return ;
	if(tr[k].val<=val)
	{
		X=k;
		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=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(head<tail&&p[i].x-p[q[head]].x>c)
		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<<num<<" "<<ans<<endl;return 0;
}
/*
4 2
1 1
3 3
2 2
10 10

2 3
*/

最开始写的

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

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

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
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.x<y.x;}}p[N];

int fa[N],sonsize[N];

inline int find(int x)
{
	int t=x;
	while(t^fa[t])t=fa[t];
	while(x^fa[x])x=fa[x],fa[x]=t;
	return x;
}

int n,root,sz;
struct balanced_tree{int ls,rs,sz,val,pos,rnd;}tr[N];

inline void lturn(int &k)
{
	int t=tr[k].rs;tr[k].rs=tr[t].ls;tr[t].ls=k;
	tr[t].sz=tr[k].sz;tr[k].sz=tr[tr[k].ls].sz+tr[tr[k].rs].sz+1;k=t;
}

inline void rturn(int &k)
{
	int t=tr[k].ls;tr[k].ls=tr[t].rs;tr[t].rs=k;
	tr[t].sz=tr[k].sz;tr[k].sz=tr[tr[k].ls].sz+tr[tr[k].rs].sz+1;k=t;
}

int insert(int &k,int val,int pos)
{
	if(!k)
	{k=++sz;tr[k].sz=1;tr[k].val=val;tr[k].pos=pos;tr[k].rnd=rand();return 1;}
	if(tr[k].val==val){tr[k].pos=pos;return 0;}
	if(val<tr[k].val)
	{
		if(insert(tr[k].ls,val,pos));tr[k].sz++;
		if(tr[k].rnd>tr[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].rnd<tr[tr[k].rs].rnd)rturn(k),del(tr[k].rs,val,pos);
		else lturn(k),del(tr[k].ls,val,pos);
		return 1;
	}
	if(val<tr[k].val)if(del(tr[k].ls,val,pos))tr[k].sz--;
	if(val>tr[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(head<tail&&(p[i].x-p[q[head]].x)>c)
		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<<num<<" "<<ans<<endl;return 0;
}
/*
4 2
1 1
3 3
2 2
10 10

2 3
*/

BlackJack_

关注

1
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
BZOJ 1604: [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居曼哈顿距离转切比雪夫距离 Treap

1604: [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1131 Solved: 463[Submit][Status][Discuss]Description了解奶牛们的人都知道，奶牛喜欢成群结队．观察约翰的N(1≤N≤100000)只奶牛，你会
复制链接

扫一扫