栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 3806 Incircle and Circumc...

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 3806 Incircle and Circumc...

#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>using namespace std ;#define REP( i , a , b ) for ( int i = a ; i < b ; ++ i )#define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )#define REV( i , a , b ) for ( int i = a ; i >= b ; -- i )#define CLR( a , x ) memset ( a , x , sizeof a )#define CPY( a , x ) memcpy ( a , x , sizeof a )const double eps = 1e-10 ;int sgn ( double x ) {return ( x > eps ) - ( x < -eps ) ;}int main () {double r , R ;while ( ~scanf ( "%lf%lf" , &r , &R ) ) {if ( sgn ( R - 2 * r ) < 0 ) printf ( "NO Solution!n" ) ;else {double low = 0 , high = sqrt ( 3.0 ) * R ;int cnt = 0 ;while ( high - low > eps ) {++ cnt ;double c = ( low + high ) / 2 ;double a = sqrt ( pow ( sqrt ( R * R - c * c / 4 ) + R , 2 ) + c * c / 4 ) ;if ( sgn ( pow ( ( sqrt ( r * r + pow ( a - c / 2 , 2 ) ) + r ) , 2 ) + c * c / 4 - a * a ) < 0 ) high = c ;else low = c ;}double c = low ;double a = sqrt ( pow ( sqrt ( R * R - c * c / 4 ) + R , 2 ) + c * c / 4 ) ;printf ( "%.15f %.15f %.15fn" , a , a , c ) ;}}return 0 ;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/374764.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号