【Noip2017】D2T1奶酪（cheese）

我终于把这道题给干掉了，介于这道题实在太皮，我这里详细讲解一下。

首先我对我的方法感到有些神奇，我一来就选择了并查集方法，我也不知道为什么，反正我现在想起来我也觉得当时非常的瓜，因为并查集这个东西听他们说很容易被卡，于是我就与这道题苦战了将近一个多月，今天早上终于把他A了。

下面我们来提一下本题的注意事项：

注意事项：

1、一定要sort一下，否则要被卡。

2、一定要做好强制转化（long long → double）

3、并查集好像还不是传统型（也就是因为这一点所以我的老并查集这道题一直A不了）

4、还有memset，不要忘！

好了差不多了，下面上代码了：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1010;
double dis[maxn][maxn];
struct sd{
	long long x,y,z;
}point[maxn];
long long n,h,r;
bool mainflag=false;
bool flagup=false,flagdown=false;
bool x[1010][3];
bool cmp(sd a,sd b)
{
	return a.z < b.z;
}
double dist(double x,double y,double z,double x1,double y1,double z1)
{
    return sqrt((x-x1)*(x-x1)+(y-y1)*(y-y1)+(z-z1)*(z-z1));
}
void init();
void tryit();
int main()
{
	int num;
	scanf("%d",&num);
	for(int i=1;i<=num;++i)
	{
		init();
		if(!mainflag)
		{
			printf("No\n");
		}
		mainflag=false;
		flagup=0;flagdown=0;
		memset(x,0,sizeof(x));
		memset(dis,0,sizeof(dis));
		memset(point,0,sizeof(point));
	}
}
void init()
{
	scanf("%lld%lld%lld",&n,&h,&r);
	for(int i=1;i<=n;++i)
	{
		scanf("%d%d%d",&point[i].x,&point[i].y,&point[i].z);
	}
	sort(point+1,point+1+n,cmp);
	for(int i=1;i<=n;++i)
	{
		if(point[i].z+r>=h)
		{
			x[i][2]=true; flagup=true;
		}
		if(point[i].z-r<=0)
		{
			x[i][1]=true; flagdown=true;
		}
	}
	if(!flagup&&!flagdown)
	{
		mainflag=true;
		printf("No\n");
		return;
	}
	for(int i=1;i<=n;++i)
	{
		for(int j=i;j<=n;++j)
		{
			dis[i][j]=dis[j][i]=dist(point[i].x,point[i].y,point[i].z,point[j].x,point[j].y,point[j].z);
		}
	}
	tryit();
	return;
}
void tryit()
{
	double maxlen=2.0*(double)r;
	for(int i=1;i<=n;++i)
	{
		for(int j=1;j<=n;++j)
		{
			if(i==j) continue;
			if(dis[i][j]<=maxlen&&x[i][1])
			{
				x[j][1]=true;
			}
			if(dis[i][j]<=maxlen&&x[i][2])
			{
				x[j][2]=true;
			}
		}
	}
	for(int i=1;i<=n;++i)
	{
		if(x[i][1]&&x[i][2])
		{
			mainflag=true;
			printf("Yes\n");
			return;
		}
	}
}

一直A不了的并查集：

#include<bits/stdc++.h>
using namespace std;
int num,n,h,r; 
long long x[1005];
int head,tail,up[1005],under[1005];
struct sd{
    long long x,y,z;
}loc[1005];
void init();
bool UFS();
double dis(double x,double y,double z,double x1,double y1,double z1)
{
    return sqrt((x-x1)*(x-x1)+(y-y1)*(y-y1)+(z-z1)*(z-z1));
}
int father(int v)
{
    if(x[v]!=v)
    {
        x[v]=father(x[v]);
    }
    return x[v];		
}
bool cmp(sd a,sd b)
{
    return a.z<b.z;
}
int main()
{
    scanf("%d",&num);	
    for(int i=1;i<=num;++i)
    {
        init();
        head=0; tail=0;
    }
}
void init()
{
    head=0;
    for(int i=1;i<=1005;++i)
    {
        x[i]=i;
    }
    scanf("%lld%lld%lld",&n,&h,&r);
    if(n==0)
    {
        printf("No\n");
        return ;
    }
    for(int i=1;i<=n;++i)
    {
        scanf("%lld%lld%lld",&loc[i].x,&loc[i].y,&loc[i].z);
    }
    sort(loc+1,loc+1+n,cmp);
    for(int i=1;i<=n;++i)
    {
        if(loc[i].z+r>=h)
        {
            head++;
            up[head]=i;
        }
        if(loc[i].z-r<=0)
        {
            tail++;
            under[tail]=i;
        }
    }
    for(int i=1;i<=n;++i)
    {
        for(int j=i-1;j>=1;--j)
        {
            if(dis(loc[i].x,loc[i].y,loc[i].z,loc[j].x,loc[j].y,loc[j].z)<=2*r)
            {
                x[i]=father(j);
            }
        }
    }
    if(!UFS())
    {
        printf("No\n");
    }
    else
    {
        printf("Yes\n");
    }
}
bool UFS()
{
    for(int i=1;i<=head;++i)
    {
        for(int j=1;j<=tail;++j)
        {
            if(x[up[i]]==x[under[j]])
            {
                return true;
            }
        }
    }
    return false;
}

2018.7.25

我还是回来了，因为一种不甘心，虽然我的并查集技术非常的菜，所以说这次我写的是BFS，然后居然一下子就过了，我还是非常的高兴的：

代码如下：

#include<bits/stdc++.h>
using namespace std;
struct sd{
    int x,y,z,num;
}loc[1005];
int flag[1005];
int qu[1000005];
int p1,p2;
bool vis[1005];
vector<int> edge[1005];
int t,n,h,r;
double dis(double x,double y,double z,double x1,double y1,double z1)
{
    return sqrt((x-x1)*(x-x1)+(y-y1)*(y-y1)+(z-z1)*(z-z1));
}
bool bfs()
{
    while(p1<=p2)
    {
        for(int i=edge[qu[p1]].size()-1;i>=0;--i)
        {
            
            if(flag[edge[qu[p1]][i]]) return true;
            if(!vis[edge[qu[p1]][i]])vis[edge[qu[p1]][i]]=true,p2++,qu[p2]=edge[qu[p1]][i];
        }
        p1++;
    }
    return false;
}
void init()
{
    p1=0;p2=0;
    for(int i=0;i<1005;++i) vector<int>().swap(edge[i]);//或者vec.swap(vector<int>());
    memset(qu,0,sizeof(qu));
    memset(flag,false,sizeof(flag));
    memset(vis,false,sizeof(vis));
    scanf("%d%d%d",&n,&h,&r);
    int x,y,z;
    bool c=false;
    for(int i=1;i<=n;++i)
    {
        scanf("%d%d%d",&x,&y,&z);
        loc[i].x=x;loc[i].y=y;loc[i].z=z;loc[i].num=i;
        for(int j=1;j<=i-1;++j)
        {
            if(dis(loc[i].x,loc[i].y,loc[i].z,loc[j].x,loc[j].y,loc[j].z)<=(double)r*2)
            edge[i].push_back(j),edge[j].push_back(i); 
        }
        if(loc[i].z-r<=0&&loc[i].z+r>=0&&loc[i].z+r>=h&&loc[i].z-r<=h){c=true;}
        if(loc[i].z-r<=0&&loc[i].z+r>=0)  vis[i]=true,p2++,qu[p2]=i;
        if(loc[i].z+r>=h&&loc[i].z-r<=h) flag[i]=true;
    }
    if(c){ printf("Yes\n");return;} 
    p1=1;
    if(bfs()) printf("Yes\n");
    else printf("No\n");
}
int main()
{
    scanf("%d",&t);
    for(int i=1;i<=t;++i)
    {
        init();
    } 
    return 0;
} 
/*
1
1 100 75
80 39 34
*/