UESTC1713 无关青云路,无关诗书,无你处,无江湖(计算几何)
2017-08-09 14:02
197 查看
无关青云路,无关诗书,无你处,无江湖
Time Limit: 3000/1000MS (Java/Others) Memory Limit: 65535/65535KB (Java/Others)
Submit
Status
每个圆的属性可以使用三个参数来描述,xx-圆心横坐标,yy-圆心纵坐标,rr-半径,现给出两个圆分别的参数,求两个圆相交的面积。
Input
第一行是一个正整数TT(T≤12T≤12),表示有TT组测试数据接下来有T行,每行有六个正整数描述两个圆的属性x1,y1,r1,x2,y2,r2x1,y1,r1,x2,y2,r2。 0≤x1,y1,r1,x2,y2,r2≤5000≤x1,y1,r1,x2,y2,r2≤500
Output
对于每组测试数据输出一行,表示两个圆相交的面积大小,结果保留6位小数
Sample input and output
Sample Input | Sample Output |
---|---|
3 0 0 10 15 0 10 -10 -10 5 0 -10 10 100 100 20 100 110 20 | 45.331175 35.076661 860.843690 |
思路:两圆情况分析一下,分为相切,相交,内含,直接套模板就行
#include <stdio.h> #include <stdlib.h> #include <math.h> #define PI 3.1415926535898 struct circle { double x,y; double r; }; double dis(circle a,circle b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } double solve(circle a, circle b) { double d = dis(a, b); if (d >= a.r + b.r) return 0; if (d <= fabs(a.r - b.r)) { double r = a.r < b.r ? a.r : b.r; return PI * r * r; } double ang1 = acos((a.r * a.r + d * d - b.r * b.r) / 2. / a.r / d); double ang2 = acos((b.r * b.r + d * d - a.r * a.r) / 2. / b.r / d); double ret = ang1 * a.r * a.r + ang2 * b.r * b.r - d * a.r * sin(ang1); return ret; } int main() { int t; while(scanf("%d",&t)!=EOF) { for(int i=0;i<t;i++) { circle a,b; scanf("%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&a.r,&b.x,&b.y,&b.r); printf("%lf\n",solve(a,b)); } } return 0; }
相关文章推荐
- 计算几何基础模板
- 计算几何模板
- POJ1269---Intersecting Lines (计算几何基础:直线相交求交点)
- poj 2318 TOYS (计算几何)
- 【计算几何】bzoj1043 [HAOI2008]下落的圆盘
- poj1039 Pipe(计算几何叉积求交点)
- ACM 计算几何向量
- 计算几何专项:UVa 10566
- POJ 1410 计算几何
- POJ 计算几何入门题目推荐
- HDU 4173 Party Location(计算几何,枚举)
- UVA-1331 Minimax Triangulation 区间dp 计算几何 三角剖分 最大三角形最小化
- hdu 1077 Catching Fish(计算几何)
- 计算几何基础知识
- POJ 1696 Space Ant <计算几何>
- 计算几何专项:UVa 10522
- 3716: [PA2014]Muzeum 计算几何+贪心+set
- poj 2318 TOYS (计算几何)
- 留着以后慢慢做的计算几何(题表)
- 【bzoj 1043】下落的圆盘(计算几何)