uva10341Solve It（二分＋精度问题）

最新推荐文章于 2021-09-07 16:36:53 发布

zoe_ed

最新推荐文章于 2021-09-07 16:36:53 发布

阅读量403

点赞数

分类专栏：二分-精度文章标签： uva

本文链接：https://blog.csdn.net/zyd8888102/article/details/49258265

版权

二分-精度专栏收录该内容

1 篇文章 0 订阅

订阅专栏

链接：
http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=113&page=show_problem&problem=1282

原题：
Solve the equation:
p*e-x + q*sin(x) + r*cos(x) + s*tan(x) + t*x2 + u = 0
where 0 <= x <= 1.
Input
Input consists of multiple test cases and terminated by an EOF. Each test case consists of 6 integers in a single line: p, q, r, s, t and u(where 0 <= p,r <= 20 and -20 <= q,s,t <= 0). There will be maximum 2100 lines in the input file.
Output
For each set of input, there should be a line containing the value of x, correct upto 4 decimal places, or the string “No solution”, whichever is applicable.
Sample Input
0 0 0 0 -2 1
1 0 0 0 -1 2
1 -1 1 -1 -1 1
Sample Output
0.7071
No solution
0.7554
感悟：
这个题窝第一次做的时候一看只保留4位小数也就是说只有0.0000-1.0000 这一万零一种情况，于是就想，如果有解那可能是其中之一。题目也没说如果有多组解应该怎么输出那肯定只有一个解了。。。。想到这里，如果能枚举，肯定就可以二分。。于是这个函数肯定是个单调函数。。。只是枚举的话，将枚举的数与零比较的时候，无法准确判断正确答案，所以不能用枚举。
－－－－－－－－－－－－－－－－－－－－－－－－－－
首先，泥得发现这是一个单调的函数，然后就可以用二分。。。但是发现wa的原因是精确度的问题。理论上 0.1234561 0.1234569 这俩数差值就小于1e－6，无论窝return哪一个答案都是一样的可是，，，事实上窝错了。（r－l）<1e-6这个判断要改成（r-l）<1e-7 甚至更小才能过。。。目前还在找原因，，先暂时搁置一下(:з」∠)

#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=1e4;
const double minc=1e-7;
// const double e=2.718281929;
double p,q,r,s,t,u;
double js(double x)
{
    return (p*exp(-x)+q*sin(x)+r*cos(x)+s*tan(x)+t*x*x+u); 
}
double erfen(double l,double r)
{
    if (r-l<minc) return l; 
    double mid=(l+r) /2.0;
    // double ansm=js(mid);ansl=js(l);
    if (js(mid)*js(l)>0) return erfen(mid,r);
    else return erfen(l,mid);
}
int main()
{

    freopen("in.txt","r",stdin);
    while(~scanf("%lf%lf%lf%lf%lf%lf",&p,&q,&r,&s,&t,&u))
    {
        // cout<<js(0)<<' '<<js(1)<<' ';
        if(js(1)*js(0)>0) 
        {
            printf("No solution\n");
        }
        else 
        {
            printf("%.4f\n",erfen(0,1));
        }

    }
    return 0;

}