求三角形外接圆 hdu4720

最新推荐文章于 2023-08-13 17:15:27 发布

tabo0

最新推荐文章于 2023-08-13 17:15:27 发布

阅读量280

点赞数

CC 4.0 BY-SA版权

分类专栏：计算几何文章标签：外接圆圆心

本文链接：https://blog.youkuaiyun.com/qq_25973789/article/details/80011821

计算几何专栏收录该内容

6 篇文章

订阅专栏

本文提供HDU4720题目的详细解答，针对钝角三角形特殊情况，介绍如何判断并处理以获得最小面积。代码实现采用C++，通过计算几何的方法找出三角形的外接圆心，进而判断指定点的安全性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

hdu4720

详细解释

坑点：如果为钝角三角形，不求外接圆，圆心直接在最大边中点，面积才为最小

#include<bits/stdc++.h>
#define inf  0x7fffffff
#define linf 0x7fffffffffffffff
#define fil(x,y) memset(x,y,sizeof(x))
#define fup(i,x,y) for(int i=(x);i<=(y);i++)
#define fdn(i,x,y) for(int i=(x);i>=(y);i--)
#define sp(x) setprecision(x)
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d%d",&n,&m)
#define sddd(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define sld(n) scanf("%lld",&n)
#define sldd(n,m) scanf("%lld%lld",&n,&m)
#define slddd(n,m,k) scanf("%lld%lld%lld",&n,&m,&k)
#define sf(n) scanf("%lf",&n)
#define sff(n,m) scanf("%lf%lf",&n,&m)
#define sfff(n,m,k) scanf("%lf%lf%lf",&n,&m,&k)
#define sc(n) scanf("%s",n)
#define lowbit(x) (x&(-x))
#define pf(x) printf("%d\n",x)
#define pfl(x) printf("%lld\n",x)
#define pff(x) printf("%lf\n",x)
#define debug printf("!!\n");
#define N 1000005
#define M 4000009
#define pi acos(-1)
#define eps 1e-2
//cout.setf(ios::fixed);
//freopen("out.txt","w",stdout);// freopen("in.txt","r",stdin);
using namespace std;
typedef long long  ll;
typedef double db;
using namespace std;
typedef long long ll;
using namespace std;
db dis(db a,db b,db c,db d)
{
    return pow(a-c,2)+pow(b-d,2);
}
db arccos(db a, db b,db c)
{
    return (a*a+b*b-c*c)/(2*a*b);
}
int main()
{
    int n,t;
    sd(t);
    int k=0;
    while(t--)
    {
        int x1,x2,x3,y1,y2,y3;
        db x4,y4;
        sdd(x1,y1);
        sdd(x2,y2);
        sdd(x3,y3);
        sff(x4,y4);
        int a1,a2,b1,b2,c1,c2;
        a1=2*(x2-x1);
        b1=2*(y2-y1);
        c1=x2*x2+y2*y2-x1*x1-y1*y1;
        a2=2*(x3-x2);
        b2=2*(y3-y2);
        c2=x3*x3+y3*y3-x2*x2-y2*y2;
        db x=((c1*b2)-(c2*b1))*1.0/((a1*b2)-(a2*b1));
        db y=((a1*c2)-(a2*c1))*1.0/((a1*b2)-(a2*b1));
        printf("Case #%d: ",++k);
        db dis1=dis(x1,y1,x2,y2);
        db dis2=dis(x1,y1,x3,y3);
        db dis3=dis(x2,y2,x3,y3);
        if(arccos(dis1,dis2,dis3)<0)
        {
            x=(x2+x3)/2;
            y=(y2+y3)/2;
        }
        else if(arccos(dis1,dis3,dis2)<0)
        {
            x=(x1+x3)/2;
            y=(y1+y3)/2;
        }
        else if(arccos(dis3,dis2,dis1)<0)
        {
            x=(x2+x1)/2;
            y=(y2+y1)/2;
        }
        if(dis(x1,y1,x,y)>=dis(x4,y4,x,y)) puts("Danger");
        else puts("Safe");
    }
    return 0;
}