Noip2017 Day2 T1 奶酪

题目描述

现有一块大奶酪，它的高度为 h，它的长度和宽度我们可以认为是无限大的，奶酪中间有许多半径相同的球形空洞。我们可以在这块奶酪中建立空间坐标系，在坐标系中，奶酪的下表面为z =0，奶酪的上表面为z=h。

现在，奶酪的下表面有一只小老鼠 Jerry，它知道奶酪中所有空洞的球心所在的坐 标。如果两个空洞相切或是相交，则 Jerry 可以从其中一个空洞跑到另一个空洞，特别地，如果一个空洞与下表面相切或是相交，Jerry 则可以从奶酪下表面跑进空洞；如果 一个空洞与上表面相切或是相交，Jerry 则可以从空洞跑到奶酪上表面。

位于奶酪下表面的 Jerry 想知道，在不破坏奶酪的情况下，能否利用已有的空洞跑 到奶酪的上表面去?

空间内两点P1(x1,y1,z1)、P2(x2​,y2​,z2​)的距离公式如下：

　　　DIS(P1,P2)=√[(x1-x2)²+(y1-y2)²+(z1-z2)²]
输入格式：

每个输入文件包含多组数据。

输入文件的第一行，包含一个正整数T，代表该输入文件中所含的数据组数。

接下来是T组数据，每组数据的格式如下： 第一行包含三个正整数 n,h 和r，两个数之间以一个空格分开，分别代表奶酪中空洞的数量，奶酪的高度和空洞的半径。

接下来的n行,每行包含三个整数 x,y,z，两个数之间以一个空格分开，表示空洞球心坐标为(x,y,z)。

输出格式：

输出文件包含 T 行，分别对应 T 组数据的答案，如果在第 i 组数据中,如果Jerry 能从下表面跑到上表面，则输出 Yes，如果不能，则输出 No 。

本题是用并查集做的（好像很多人都是用搜索做的），而且这道题好像也并不是一道小水题了（还需要判断一下距离==），但是这并不能难道蒟蒻的我←_←

思想：除了并查集外，我们就只需要判断距离了，首先，我们现将能与x-y平面连通的球与 0点 并到一起，高度能达到h的球与 n+1 点并到一起，然后如果这两个球之间的距离 l<=2*r ，那么我们就把他们并到一起，最后只去要查询一下 0点 与 n+1点 是否连同就能得出小老鼠是不是可以跑到奶酪上面去了。

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long LL;
 const int maxn=;
 struct cp
 {
 　　LL x,y,z;
 }e[maxn];
 int f[maxn];
 bool cmp(cp a,cp b)
 {
 　　return a.z<b.z;
 }
 LL square(LL x)
 {
 　　return x*x;
 }
 LL DIS(cp a,cp b)
 {
 　　return square(a.z-b.z)+square(a.x-b.x)+square(a.y-b.y);
 }
 int find(int x)
 {
 　　if(f[x]==x) return x;
 　　return f[x]=find(f[x]);
 }
 void Union(int a,int b)
 {
 　　int x=find(a),y=find(b);
 　　f[x]=y;
 }
 int main()
 {
 　　LL n,m;
 　　LL t,h;
 　　LL r;
 　　scanf("%lld",&t);
 　　while(t--){
 　　　　scanf("%lld%lld%lld",&n,&h,&r);
 　　　　for(int i=;i<=n+;i++) f[i]=i;
 　　　　for(int i=;i<=n;i++){
 　　　　　　scanf("%lld%lld%lld",&e[i].x,&e[i].y,&e[i].z);
 　　　　}
 　　　　sort(e+,e+n+,cmp);
 　　　　for(int i=;i<=n;i++)
 　　　　　for(int j=i+;j<=n;j++)
 　　　　　　if(DIS(e[i],e[j])<=*r*r) Union(i,j);
 　　　　for(int i=;i<=n;i++){
 　　　　　　if(e[i].z-r<=) Union(i,);
 　　　　　　if(e[i].z+r>=h) Union(i,n+);
 　　　　}
 　　　　if(find()==find(n+)) printf("Yes\n");
 　　　　else printf("No\n");
 　　}
 　　return ;
 }