POJ 3737

第一道三分题,有模板

#define eps 10e-6
double cal(){}//计算题目所需要的值
while(l+eps<r)
{
    m1=l+(r-l)/3;
    m2=r-(r-l)/3;
    v1=cal(m1);
    v2=cal(m2);
    if(v1<v2)l=m1;
    else r=m2;
}

  

感觉三分像是数学题,因为它求的是单峰极值的问题。

#include <iostream>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;

const double PI=3.141592653;
const double eps=1e-8;

double S,RR,LL,m,mm;

double cal(double r){
	return sqrt(((r*S*S)-(2*PI*S*r*r))/9.0);
}

int main(){
	while(scanf("%lf",&S)!=EOF){
		RR=S/PI;LL=0;
		while(LL+eps<RR){
			m=LL+(RR-LL)/3;
			mm=RR-(RR-LL)/3;
			if(cal(m)<cal(mm))
			LL=m;
			else RR=mm;
		}
		double v=cal(LL);
		printf("%.2lf\n%.2lf\n%.2lf\n",v,v*3/(PI*LL),sqrt(LL));
	}
	return 0;
}

  

转载于:https://www.cnblogs.com/jie-dcai/p/4273316.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值