题意:已知三角形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. JavaScript与html5写的贪吃蛇完整代码

    JavaScript与html5写的贪吃蛇完整代码 查看运行效果可访问http://www.codesocang.com/texiao/youxitexiao/2014/0402/7045.html# ...

  2. RedHat7上安装MariaDB

    编译安装MariaDB 下载MariaDB# wget http://mirrors.opencas.cn/mariadb//mariadb-10.1.8/source/mariadb-10.1.8. ...

  3. sublime 2如何进入vim模式

    点击菜单栏[Preferences]——[Settings - Defaults] 查找: "ignored_packages": ["Vintage"] 改为 ...

  4. Spring读书笔记-----Spring的Bean之设置Bean值

    [java] view plaincopyprint? Java实例的属性值可以有很多种数据类型.基本类型值.字符串类型.java实例甚至其他的Bean实例.java集合.数组等.所以Spring允许 ...

  5. Mono For Android的项目发布

    Mono for Android的发布步骤: 1.选择项目名称鼠标右键,选择Options选项. 2.弹出工程选项框,选择Android Build,勾选 Supported ABIS里面的所有项. ...

  6. python基础知识一

    数 python中有4种类型的数--整数.长整数.浮点数和复数. --2是一个整数 --长整数不过是大一些的整数 --3.23和52.3E-4是浮点数的例子.E标记表示10的幂.在这里52.3E-4表 ...

  7. jmeter压测app

    使用代理的方式,录制app端脚本,之后用jmeter压测就没啥好说的了 1.电脑端谷歌设置本地代理(端口号为8888) 2.jmeter设置HTTP代理服务器(端口号为8888) 3.手机端wifi设 ...

  8. Android - This Handler class should be static or leaks might occur.

    今天学习了使用 HTTP协议,从Android客户端往Tomcat服务器端以GET发送请求,途中无意中发现自己写的Handler类被Android提示:This Handler class shoul ...

  9. javascript - 图片的幻灯片效果

    javascript 代码: <script type="text/javascript"> function select_play() { var select_p ...

  10. .NET下的加密解密大全(3):非对称加密

    本博文列出了.NET下常用的非对称加密算法,并将它们制作成小DEMO,希望能对大家有所帮助. RSA[csharp]static string EnRSA(string data,string pub ...