hdu-1140(求距离，精度判断)

卫星风暴消除算法解析

最新推荐文章于 2025-09-05 16:58:34 发布

转载最新推荐文章于 2025-09-05 16:58:34 发布 · 130 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/2018zxy/p/9905051.html

文章标签：

#php

本文详细解析了一种算法，该算法用于确定卫星能否消除地面风暴，通过计算卫星与风暴的距离来判断是否在卫星作用范围内。文章包括完整的代码实现，考虑了PI的精度问题，并提供了清晰的思路说明。

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1140

思路：卫星只能消灭地面上一部分的风暴，即风暴与卫星的距离最大是卫星到地球的切线的距离，大于这个距离就不行了。

（注意PI的精度）

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
const double PI = 3.1415926;
const double R = 20000.0/PI;
struct Node{
    double x,y,z,dis;
};
Node a[120],tp;
int main(void)
{
    int n,m,i,j,ans;
    double len;
    while(~scanf("%d%d",&n,&m)&&(n+m))
    {
        for(i=0;i<n;i++) 
        {
            scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].z);
            a[i].dis=a[i].x*a[i].x+a[i].y*a[i].y+a[i].z*a[i].z-R*R;
        }
        ans=0;
        for(j=0;j<m;j++)
        {
            scanf("%lf%lf%lf",&tp.x,&tp.y,&tp.z);
            
            for(i=0;i<n;i++)
            {
                len=(a[i].x-tp.x)*(a[i].x-tp.x)+(a[i].y-tp.y)*(a[i].y-tp.y)+(a[i].z-tp.z)*(a[i].z-tp.z);
                if(sqrt(len)<=sqrt(a[i].dis))
                {
                    ans++;break;
                }
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}