题意:已知三角形ABC的3条边长,求三角形ABC 的面积,以及阴影部分的总面积。

#include<iostream>
#include<cstdio>
#include<cmath>
#define sqr(a) ((a)*(a))
#define pi 2.0*asin(1.0) using namespace std; typedef struct point
{
double x,y;
point(double xx=0,double yy=0):x(xx),y(yy){}
}vector; vector operator - (point a,point b)
{
return vector(a.x-b.x,a.y-b.y);
}
point operator + (point a,vector b)
{
return point(a.x+b.x,a.y+b.y);
}
vector operator * (vector a,double b)
{
return vector(a.x*b,a.y*b);
}
double dot(vector a,vector b)
{
return a.x*b.x+a.y*b.y;
}
double cross(vector a,vector b)
{
return a.x*b.y-a.y*b.x;
}
double len(vector a)
{
return sqrt(dot(a,a));
}
vector resiz(vector a,double l)
{
l/=len(a);
return vector(a.x*l,a.y*l);
}
vector rot(vector a,double rad)
{
double c=cos(rad),s=sin(rad);
return vector(a.x*c-a.y*s,a.x*s+a.y*c);
}
point inter(point p,vector v,point q,vector w)
{
vector u=p-q;
double t=cross(w,u)/cross(v,w);
return p+v*t;
}
point f(point o,vector v,point a,vector w,double r,double &rr,double rad)
{
vector u=resiz(rot(w,pi/2),2*r);
if(dot(v,u)>0) u=vector(-u.x,-u.y);
point b=a+u;
point p=inter(o,v,a,w),q=inter(o,v,b,w);
rr=len(p-o)/len(q-o)*r;
double d=rr/sin(rad);
v=resiz(v,d);
return o+v;
} int main()
{
int i=1;
double a,b,c,tha,thb,thc;
double r,r1,r2,r3;
point pa,pb,pc,p1,p2,p3;
vector va,vb,vc;
while(cin>>a>>b>>c && a+b+c)
{
tha=(sqr(b)+sqr(c)-sqr(a))/(2*b*c);
thb=(sqr(a)+sqr(c)-sqr(b))/(2*a*c);
thc=(sqr(b)+sqr(a)-sqr(c))/(2*b*a);
tha=acos(tha);
thb=acos(thb);
thc=acos(thc);
vc.x=pa.x+c;
pb=pa+vc;
vb=resiz(rot(vc,tha),b);
pc=pa+vb;
va=pc-pb;
r=cross(vc,vb)/(a+b+c);
p1=f(pa,rot(vc,tha/2),pb,va,r,r1,tha/2);
p2=f(pb,rot(va,thb/2),pc,vb,r,r2,thb/2);
p3=f(pc,rot(vb,thc/2-pi),pa,vc,r,r3,thc/2);
double s1=fabs(cross(p1-p2,p1-p3))/2;
a=len(p1-p2);
b=len(p2-p3);
c=len(p3-p1);
tha=(sqr(b)+sqr(c)-sqr(a))/(2*b*c);
thb=(sqr(a)+sqr(c)-sqr(b))/(2*a*c);
thc=(sqr(b)+sqr(a)-sqr(c))/(2*b*a);
tha=acos(tha);
thb=acos(thb);
thc=acos(thc);
double s2=tha/2*sqr(r3)+thb/2*sqr(r1)+thc/2*sqr(r2);
printf("Case %d: %.2lf %.2lf\n",i++,s1,s2);
}
return 0;
}

uva 11731 - Ex-circles的更多相关文章

  1. UVa 247 - Calling Circles(Floyd求有向图的传递闭包)

    链接: https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem& ...

  2. UVa 247 Calling Circles【传递闭包】

    题意:给出n个人的m次电话,问最后构成多少个环,找出所有的环 自己想的是:用map来储存人名,每个人名映射成一个数字编号,再用并查集,求出有多少块连通块,输出 可是map不熟,写不出来,而且用并查集输 ...

  3. UVA - 247 Calling Circles Floyd判圈

    思路:利用的Floyd判圈,如果i能到j,j也能到i说明i和j在同一个圈里.每个人的名字可用map编号.最后DFS打印答案即可. AC代码 #include <cstdio> #inclu ...

  4. UVa 247 Calling Circles (DFS+Floyd)

    题意:如果两个人互通电话,那么他们就在一个电话圈里,现在给定 n 个人,并且给定 m 个通话记录,让你输出所有的电话圈. 析:刚开始没想到是Floyd算法,后来才知道是这个算法,利用这个算法进行连通性 ...

  5. UVA LIVE-4642 - Malfatti Circles

    给出三角形三个顶点,求出三个互切的圆的半径 尽管大白鼠说能够推出公式,但是这个公式仅仅怕没那么easy推--我左看右看上看下看也推不出. 应该是要做辅助线什么的,那也-- 因为非常easy就推出了关于 ...

  6. UVA 247 - Calling Circles (Floyd)

    互相可以打电话是一个传递关系,所以Floyd求传递封包,dfs找一个尽量大的圈. #include<bits/stdc++.h> using namespace std; ; map< ...

  7. UVA - 247 Calling Circles(Floyd求传递闭包)

    题目: 思路: 利用Floyd求传递闭包(mp[i][j] = mp[i][j]||(mp[i][k]&&mp[k][j]);),当mp[i][j]=1&&mp[j][ ...

  8. UVA 247 Calling Circles —— (强连通分量模板题)

    第一个强连通分量的题. 题意:有一堆人,a给b打电话表示a有一条向b的边,一个强连通分量代表一个电话圈,把每个电话圈里的人在一行内输出出来. 直接上模板即可,但是要注意把string用map映射一下的 ...

  9. UVA 247"Calling Circles"(floyd求传递闭包+SCC)

    传送门 题意: 如果两个人相互打电话(直接或间接),则说他们在同一个电话圈里. (a,b) 表示 a 打给 b: 例如,(a,b),(b,c),(c,d),(d,a),则这四个人在同一个电话圈里: 输 ...

随机推荐

  1. springmvc使用@ResponseBody返回json乱码解决方法

    1.springmvc 3.2以上的版本解决乱码的方法: 第一步:在配置中加入: <mvc:annotation-driven> <mvc:message-converters re ...

  2. 【php】中【event】之实现方式

    这两天看了点事件机制,那么在php中,如何实现最简单的事件呢? 废话不多说,我们上代码. <?php class Event{ //事件名称 public $name; //存储hander p ...

  3. jquery插件下载地址

    以下是本人收集的jquery插件下载地址: .............版本自行选择. jquery官网:http://jquery.com/ jquery.validate.js 官网下载地址:htt ...

  4. dispatchkeyevent的调用机制

    http://blog.csdn.net/look85/article/details/23740761 dispatchKeyEvent和onKeyDown关系: 当键盘按下时 首先触发dispat ...

  5. 黑马程序员-集合(二)contains()方法的内部探索

    ------Java培训.Android培训.iOS培训..Net培训.期待与您交流! ------- 我们知道集合是用来存储对象的.在他们实现了众多的接口我们以Arraylist为列子 所有已实现的 ...

  6. 1.redis.3.2 下载,安装、配置、使用 - 1

    1.下载:   2.使用: 挤压之后,使用cmd执行,如下图 redis-server--service-installredis.windows.conf,执行安装 提示成功之后,剩下就好办了, 这 ...

  7. ubuntu出现有线已连接却无法上网

    或者直接追加到/etc/sysctl.conf 如果遇到“设备未托管”,一般是台式机默认移动ip后禁用网络. 那么修改/etc/NetworkManager/NetworkManager.conf,设 ...

  8. while循环的跳出

    今天在编码时突然产生一个疑问:程序中有一个while循环,循环体执行的是某个附带条件限制的操作.我现在想达到的目的是 => 条件成立,就执行操作,并跳出循环:条件不成立就跳出当次的while循环 ...

  9. web页面的优化

    众所周知,一个web页面通常会包括HTML(XHTML.XML).CSS.Javascript,而其中HTML(XHTML.XML)为结构化语言,用于构建页面结构和相关数据:CSS则负责页面的样式,即 ...

  10. Github上图文演示仓库创建

    1.注册账号 https://github.com/ 2.注册完账号 创建仓库 a.点击按钮 Start a project b.填写信息,点击按钮 Create repository c.仓库信息( ...