hdu4720 三角形的外接圆
题意:
给你四个点,问你第四个点是否在前三个点围成的三角形的外接圆上.
思路:
水题,就是练练用魔板罢了,当该三角形是锐角三角形的时候,圆心是任意两条边中垂线的交点,半径是圆心到任意一点的距离,否则圆心就是最长的那条边的中点位置,半径就是最长的那条边的一半..
- #include <cstdio>
- #include <cmath>
- #include <algorithm>
- #define maxn 60
- #define eps 1e-7
- using namespace std;
- int dcmp(double x) //控制精度
- {
- if(fabs(x)<eps) return 0;
- else return x<0?-1:1;
- }
- double toRad(double deg) //角度转弧度
- {
- return deg/180.0*acos(-1.0);
- }
- struct Point
- {
- double x,y;
- Point(){}
- Point(double x,double y):x(x),y(y) {}
- void input()
- {
- scanf("%lf %lf",&x,&y);
- }
- };
- typedef Point Vector;
- Vector operator+( Vector A, Vector B ) //向量加
- {
- return Vector( A.x + B.x, A.y + B.y );
- }
- Vector operator-(Vector A,Vector B) //向量减
- {
- return Vector( A.x - B.x, A.y - B.y );
- }
- Vector operator*( Vector A, double p ) //向量数乘
- {
- return Vector( A.x * p, A.y * p );
- }
- Vector operator/( Vector A, double p ) //向量数除
- {
- return Vector( A.x / p, A.y / p );
- }
- bool operator<(const Point& A, const Point& B ) //两点比较
- {
- return dcmp( A.x - B.x ) < 0 || ( dcmp( A.x - B.x ) == 0 && dcmp( A.y - B.y ) < 0 );
- }
- bool operator==( const Point& a, const Point& b ) //两点相等
- {
- return dcmp( a.x - b.x ) == 0 && dcmp( a.y - b.y ) == 0;
- }
- struct Line
- {
- Point s,e;
- Vector v;
- Line() {}
- Line(Point s,Point v,int type)://法向量式
- s(s),v(v){}
- Line(Point s,Point e):s(s),e(e)//两点式
- {v=e-s;}
- };
- double Dot(Vector A,Vector B)//向量点乘
- {
- return A.x*B.x+A.y*B.y;
- }
- double Length(Vector A)//向量模
- {
- return sqrt(Dot(A,A));
- }
- double Angle(Vector A,Vector B)//向量夹角
- {
- return acos(Dot(A,B)/Length(A)/Length(B));
- }
- double Cross(Vector A,Vector B)//向量叉积
- {
- return A.x*B.y-A.y*B.x;
- }
- double Area2(Point A,Point B,Point C )//向量有向面积
- {
- return Cross(B-A,C-A);
- }
- double Dist(Point A,Point B)
- {
- return Length(A-B);
- }
- Vector Rotate(Vector A, double rad)//向量逆时针旋转
- {
- return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));
- }
- Vector Normal(Vector A)//向量单位法向量
- {
- double L=Length(A);
- return Vector(-A.y/L,A.x/L);
- }
- Point GetLineIntersection(Line l1,Line l2)//两直线交点
- {
- Point P=l1.s;
- Vector v=l1.v;
- Point Q=l2.s;
- Vector w=l2.v;
- Vector u=P-Q;
- double t=Cross(w,u)/Cross(v,w);
- return P+v*t;
- }
- double DistanceToLine(Point P,Line L)//点到直线的距离
- {
- Point A,B;
- A=L.s,B=L.e;
- Vector v1=B-A,v2=P-A;
- return fabs(Cross(v1,v2))/Length(v1);
- }
- double DistanceToSegment(Point P, Line L)//点到线段的距离
- {
- Point A,B;
- A=L.s,B=L.e;
- if(A==B) return Length(P-A);
- Vector v1=B-A,v2=P-A,v3=P-B;
- if (dcmp(Dot(v1,v2))<0) return Length(v2);
- else if (dcmp(Dot(v1,v3))>0) return Length(v3);
- else return fabs(Cross(v1,v2)) / Length(v1);
- }
- Point GetLineProjection(Point P,Line L)// 点在直线上的投影
- {
- Point A,B;
- A=L.s,B=L.e;
- Vector v=B-A;
- return A+v*(Dot(v,P-A)/Dot(v,v));
- }
- bool OnSegment(Point p,Line l)//点在线段上包括端点
- {
- Point a1=l.s;
- Point a2=l.e;
- return dcmp(Cross(a1-p,a2-p))==0&&dcmp(Dist(p,a1)+Dist(p,a2)-Dist(a1,a2))==0;
- }
- bool Paralled(Line l1,Line l2)//直线平行
- {
- return dcmp(Cross(l1.e-l1.s,l2.e-l2.s))==0;
- }
- bool SegmentProperIntersection(Line l1,Line l2)//线段相交
- {
- if(Paralled(l1,l2))
- {
- return false;
- }
- Point t=GetLineIntersection(l1,l2);
- if(OnSegment(t,l1))
- {
- return true;
- }
- return false;
- }
- int ConvexHull(Point *p,int n,Point *ch) //求凸包
- {
- sort(p,p+n);
- int m=0;
- for ( int i = 0; i < n; ++i )
- {
- while ( m > 1 && Cross( ch[m - 1] - ch[m - 2], p[i] - ch[m - 2] ) <= 0 ) --m;
- ch[m++] = p[i];
- }
- int k = m;
- for ( int i = n - 2; i >= 0; --i )
- {
- while ( m > k && Cross( ch[m - 1] - ch[m - 2], p[i] - ch[m - 2] ) <= 0 ) --m;
- ch[m++] = p[i];
- }
- if ( n > 1 ) --m;
- return m;
- }
- double PolygonArea(Point *p,int n) //多边形有向面积
- {
- double area=0;
- for (int i=1;i<n-1;++i)
- area+=Cross(p[i]-p[0],p[i+1]-p[0]);
- return area/2.0;
- }
- double dis(Point A ,Point B)
- {
- double tmp = (A.x - B.x) * (A.x - B.x) + (A.y - B.y) * (A.y - B.y);
- return sqrt(tmp);
- }
- typedef struct
- {
- double dis;
- Point A ,B;
- }EDGE;
- EDGE edge[5];
- bool campp(EDGE a ,EDGE b)
- {
- return a.dis < b.dis;
- }
- int main ()
- {
- int t ,i ,cas = 1;
- Point p1 ,p2 ,p3 ,p;
- Point O;
- double R;
- scanf("%d" ,&t);
- while(t--)
- {
- scanf("%lf %lf" ,&p1.x ,&p1.y);
- scanf("%lf %lf" ,&p2.x ,&p2.y);
- scanf("%lf %lf" ,&p3.x ,&p3.y);
- scanf("%lf %lf" ,&p.x ,&p.y);
- edge[1].A = p1 ,edge[1].B = p2;
- edge[2].A = p1 ,edge[2].B = p3;
- edge[3].A = p2 ,edge[3].B = p3;
- edge[1].dis = dis(p1 ,p2);
- edge[2].dis = dis(p1 ,p3);
- edge[3].dis = dis(p2 ,p3);
- sort(edge + 1 ,edge + 3 + 1 ,campp);
- if(edge[1].dis * edge[1].dis + edge[2].dis * edge[2].dis <= edge[3].dis * edge[3].dis)
- {
- O.x = (edge[3].A.x + edge[3].B.x) / 2;
- O.y = (edge[3].A.y + edge[3].B.y) / 2;
- R = edge[3].dis / 2;
- }
- else
- {
- Line L1 = Line((p1 + p2)/2 ,Normal(p1 - p2),1);
- Line L2 = Line((p1 + p3)/2 ,Normal(p1 - p3),1);
- O = GetLineIntersection(L1 ,L2);
- R = dis(O ,p1);
- }
- double diss = dis(p ,O);
- if(diss <= R) printf("Case #%d: Danger\n" ,cas ++);
- else printf("Case #%d: Safe\n" ,cas ++);
- }
- return 0;
- }
hdu4720 三角形的外接圆的更多相关文章
- HDU4720+三角形外接圆
/* 几何 求给定三角形的外接圆圆心 方法:求解二元方程组 */ #include<stdio.h> #include<string.h> #include<math.h ...
- CAD:计算三角形的外接圆圆心
条件:三个定点不共线
- poj1266Cover an Arc(三角形外接圆)
链接 求出三角形的外接圆,通过圆心和半径可以知道这个圆的上下左右最远点,分别判断这个四个点跟弧的两端点A,B的关系,假如判断P点,弧内给出点为C,判断PC是否与AB相交即可判断出P是否在弧上. 精度问 ...
- POJ 1329
三角外接圆
Circle Through Three Points Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 3169 Acce ...
- poj 1090:The Circumference of the Circle(计算几何,求三角形外心)
The Circumference of the Circle Time Limit: 2 Seconds Memory Limit: 65536 KB To calculate the c ...
- Delaunay剖分与平面欧几里得距离最小生成树
这个东西代码我是对着Trinkle的写的,所以就不放代码了.. Delaunay剖分的定义: 一个三角剖分是Delaunay的当且仅当其中的每个三角形的外接圆内部(不包括边界)都没有点. 它的存在性是 ...
- XVI Open Cup named after E.V. Pankratiev. GP of Ekaterinburg
A. Avengers, The 留坑. B. Black Widow 将所有数的所有约数插入set,然后求mex. #include<bits/stdc++.h> using names ...
- Delaunay Triangulation in OpenCascade
Delaunay Triangulation in OpenCascade eryar@163.com 摘要:本文简要介绍了Delaunay三角剖分的基础理论,并使用OpenCascade的三角剖分算 ...
- 探索性思维——How to Solve It
我觉得这篇文章和什么都能扯上点关系,比如编程. 很多人已经讨论过数学与编程的关系了,这里不想过多探讨,只是简单提一下:有些人把数学贬低地一文不值,认为做一般的应用软件用不到数学:而有些人则把数学拔高到 ...
随机推荐
- HDOJ-1024(动态规划+滚动数组)
Max Sum Plus Plus HDOJ-1024 动态转移方程:dp[i][j]=max(dp[i][j-1]+a[j],max(dp[i-1][k])+a[j]) (0<k<j) ...
- 漏洞复现-CVE-2014-3120-ElasticSearch 命令执行漏洞
0x00 实验环境 攻击机:Win 10 靶机也可作为攻击机:Ubuntu18 (docker搭建的vulhub靶场) 0x01 影响版本 < ElasticSearch 1.2的版本 ...
- Java基础语法学习
Java基础语法学习 1. 注释 单行注释: //单行注释 多行注释: /*多行注释 多行注释 多行注释 多行注释 */ 2. 关键字与标识符 关键字: Java所有的组成部分都需要名字.类名.变量名 ...
- cpu缓存和volatile
目录 CPU缓存的由来 CPU缓存的概念 CPU缓存的意义 缓存一致性协议-MESI协议 Store Buffers Store Forwarding Memory Barriers Invalida ...
- Apache配置 2.用户认证
1.用户认证用来对某些目录中的网页进行访问控制,当用户访问这些页面的时候需要输入用户名和密码进行认证. 2. 配置: # vim /usr/local/apache2.4/conf/extra/htt ...
- android分析之Binder 02
分析Java层的ServiceManager,看看Binder在Java层是如何实现的. public final class ServiceManager { private static fina ...
- VS添加dll引用
直接添加(CADImport.dll) 手动添加 (sgcadexp.dll) 直接放到项目bin的目录下
- PTA 二叉树的三种遍历(先序、中序和后序)
6-5 二叉树的三种遍历(先序.中序和后序) (6 分) 本题要求实现给定的二叉树的三种遍历. 函数接口定义: void Preorder(BiTree T); void Inorder(BiTr ...
- Redis入门到放弃系列-redis数据类型
Redis数据类型? Redis 提供一些常用的数据类型:Strings.Lists.Sets.Sorted sets.Hashes.Arrays.Bitmap.Streams Strings(字符串 ...
- 用 Go + WebSocket 快速实现一个 chat 服务
前言 在 go-zero 开源之后,非常多的用户询问是否可以支持以及什么时候支持 websocket,终于在 v1.1.6 里面我们从框架层面让 websocket 的支持落地了,下面我们就以 cha ...