【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;
}

持续更新······