uva 11731 - Ex-circles
题意:已知三角形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的更多相关文章
- UVa 247 - Calling Circles(Floyd求有向图的传递闭包)
链接: https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem& ...
- UVa 247 Calling Circles【传递闭包】
题意:给出n个人的m次电话,问最后构成多少个环,找出所有的环 自己想的是:用map来储存人名,每个人名映射成一个数字编号,再用并查集,求出有多少块连通块,输出 可是map不熟,写不出来,而且用并查集输 ...
- UVA - 247 Calling Circles Floyd判圈
思路:利用的Floyd判圈,如果i能到j,j也能到i说明i和j在同一个圈里.每个人的名字可用map编号.最后DFS打印答案即可. AC代码 #include <cstdio> #inclu ...
- UVa 247 Calling Circles (DFS+Floyd)
题意:如果两个人互通电话,那么他们就在一个电话圈里,现在给定 n 个人,并且给定 m 个通话记录,让你输出所有的电话圈. 析:刚开始没想到是Floyd算法,后来才知道是这个算法,利用这个算法进行连通性 ...
- UVA LIVE-4642 - Malfatti Circles
给出三角形三个顶点,求出三个互切的圆的半径 尽管大白鼠说能够推出公式,但是这个公式仅仅怕没那么easy推--我左看右看上看下看也推不出. 应该是要做辅助线什么的,那也-- 因为非常easy就推出了关于 ...
- UVA 247 - Calling Circles (Floyd)
互相可以打电话是一个传递关系,所以Floyd求传递封包,dfs找一个尽量大的圈. #include<bits/stdc++.h> using namespace std; ; map< ...
- UVA - 247 Calling Circles(Floyd求传递闭包)
题目: 思路: 利用Floyd求传递闭包(mp[i][j] = mp[i][j]||(mp[i][k]&&mp[k][j]);),当mp[i][j]=1&&mp[j][ ...
- UVA 247 Calling Circles —— (强连通分量模板题)
第一个强连通分量的题. 题意:有一堆人,a给b打电话表示a有一条向b的边,一个强连通分量代表一个电话圈,把每个电话圈里的人在一行内输出出来. 直接上模板即可,但是要注意把string用map映射一下的 ...
- UVA 247"Calling Circles"(floyd求传递闭包+SCC)
传送门 题意: 如果两个人相互打电话(直接或间接),则说他们在同一个电话圈里. (a,b) 表示 a 打给 b: 例如,(a,b),(b,c),(c,d),(d,a),则这四个人在同一个电话圈里: 输 ...
随机推荐
- The hacker's sanbox游戏
第一关:使用/usr/hashcat程序,对passwd中root的密码进行解密,得到gravity98 执行su,输入密码gravity98. 第二关:获取提供的工具,wget http://are ...
- MongoDB 3.0安全权限访问控制(Windows版)
MongoDB 3.0安全权限访问控制(Windows版) 1.首先,不使用 –auth 参数,启动 mongoDB: mongod --dbpath "d:\mongodb\data\db ...
- 在.Net中进行跨线程的控件操作(上篇:Control.Invoke)
本文的重点在于介绍如何在多线程编程中,从非UI线程上访问界面中的控件.有过多线程编程经验的人都知道,当我们在非UI线程上试图给一个界面中的控件赋值的时候,比如说label的Text属性,系统会抛出一个 ...
- 【基础】Oracle 表空间和数据文件
多个表空间的优势:1.能够将数据字典与用户数据分离出来,避免由于字典对象和用户对象保存在同一个数据文件中而产生的I/O冲突2.能够将回退数据与用户数据分离出来,避免由于硬盘损坏而导致永久性的数据丢失3 ...
- shiro认证
一.通过ini文件初始化一个用户 1.通过ini配置文件创建securityManager2.调用subject.login方法主体提交认证,提交的token3.securityManager进行认证 ...
- centos 彻底卸载mysql
yum remove mysql mysql-server mysql-libs compat-mysql51rm -rf /var/lib/mysqlrm /etc/my.cnf查看是否还有mysq ...
- UVA 10817 Headmaster's Headache(DP +状态压缩)
Headmaster's Headache he headmaster of Spring Field School is considering employing some new teacher ...
- windows注册表的基本使用——示例
上网找好多资料发现一个问题就是太老.例如只有RegCreateKey而没有RegCreateKeyEx用法详解,自己摸索了几个小时终于基本用到的几个函数都试用了一遍. 下面代码已经通过编译测试 // ...
- (一)跟我一起玩Linux网络服务:DNS服务——BIND(/etc/named.conf、/var/named)设置实现和解释
2015年3月24更新 添加了要加的配置域名解析器(否则会找不到域名) 一.创建该实验的的模型 配置完gate虚拟机的两张网卡后,就启动gate的转发 [root@localhost ro ...
- Linux下彻底卸载LibreOffice方法
Linux下彻底卸载LibreOffice方法 终端中输入命令: 对所有基于 Debian 的发行版(Debian.Ubuntu.Kubuntu.Xubuntu.*buntu.Sidux 等): su ...