Uva - 10341 - Solve It (二分)

/**
 *  二分法。。。
 *  一开始还纠结单调性,看了看个参数的bound, 对原式求导下,可以求出是单调递减的。
 *  这题找错找N久!! 因为答案一直不对,可是算法都没错,网上搜代码,算法也一样
 *  答案就是不对!!  最后才发现,输入的参数中有个r!! 而我用二分的时候的变量也是r!!!
 *  粗心的哟。。。。。
 */
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#define INF 0x7fffffff)
#define PRE 1e-14

using namespace std;
double p, q, r, s, t, u;

double solve(double ans) {
    return p * exp(-ans) + q * sin(ans) + r * cos(ans) + s * tan(ans) + t * ans * ans + u;
}

int main()
{
    double L, R, cur, ans;
    while(scanf("%lf%lf%lf%lf%lf%lf", &p, &q, &r, &s, &t, &u) != EOF) {
        L = 0.0, R = 1.0;
        if(solve(L) * solve(R) > 0) {
            printf("No solution\n");
            continue;
        }
        while(R - L > PRE) {
            ans = L + (R - L) / 2.0;
            cur = solve(ans);
            if(cur > 0) L = ans;
            else        R = ans;
        }
        printf("%.4lf\n", ans);
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值