202009-4 星际旅行【CCF-CSP认证真题】

最新推荐文章于 2023-10-01 10:54:04 发布

WinnieXiao

最新推荐文章于 2023-10-01 10:54:04 发布

阅读量255

点赞数

分类专栏： CCF-CSP 文章标签： c语言 c++ 算法

本文链接：https://blog.csdn.net/qq_46084754/article/details/120178153

版权

CCF-CSP 专栏收录该内容

4 篇文章 0 订阅

订阅专栏

注意：

（1）能一遍算完存下来的就要及时存储，节约时间。

（2）尽可能简化判断条件。

满分代码

#include<bits/stdc++.h>
using namespace std;
int n,m;
double r;
struct point
{
    double d[110];
}p[2010];
double ans[2010];
point o;
double x[2010],diso[2010];

bool intersect(int x,int y,double c)
{
    double b=diso[x],a=diso[y];
    double p=(a+b+c)/2;
    double s=sqrt(p*(p-a)*(p-b)*(p-c));
    double h=2*s/c;
    if(h<r)
    {
        if((b*b+c*c-a*a)>0&&(a*a+c*c-b*b)>0)
        return 1;
    }
    return 0;
}

int main()
{
    cin>>n>>m>>r;
    for(int i=1;i<=n;i++) cin>>o.d[i];
    for(int i=1;i<=m;i++)
    {
        double s=0;
        for(int j=1;j<=n;j++)
        {
            cin>>p[i].d[j];
            p[i].d[j]-=o.d[j];
            s+=p[i].d[j]*p[i].d[j];
        } 
        diso[i]=sqrt(s);
        x[i]=sqrt(s-r*r);
    }
    for(int i=1;i<=m;i++)
    {
        for(int j=i+1;j<=m;j++)
        {
            double c=0,ang;
            double d;
            for(int k=1;k<=n;k++) c+=(p[i].d[k]-p[j].d[k])*(p[i].d[k]-p[j].d[k]);
            c=sqrt(c);
            if(intersect(i,j,c))
            {
                double angj=acos(r/diso[j]),angi=acos(r/diso[i]);
                ang=acos((diso[i]*diso[i]+diso[j]*diso[j]-c*c)/(2*diso[i]*diso[j]));
                double anga=ang-angj-angi;
                d=x[i]+x[j]+anga*r;
            }
            else d=c;
            //printf("%d %d %.16f\n",i,j,d);
            ans[i]+=d,ans[j]+=d;
        }
    }
    for(int i=1;i<=m;i++)
        printf("%.15lf\n",ans[i]);
    return 0;
}