Uva - 10341 - Solve It（二分）

来源：互联网发布：数据挖掘原理与实践编辑：程序博客网时间：2024/05/24 04:02

题意：解方程

p*e^-x+ q*sin(x) + r*cos(x) + s*tan(x) + t*x² + u = 0
where 0 <= x <= 1.（0 <= p,r <= 20 and -20 <= q,s,t <= 0）

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

——>>仔细观察可发现每一项在定义的域上都是减函数，所以，直接二分吧……

第一条数值题，学了3样东西：

1、二分迭代100次；

2、浮点数比较用个误差eps = 1e-14，判正数用eps，判负数用-eps；

3、函数式可以用#define来定义喔……

#include <cstdio>#include <cmath>using namespace std;#define F(x) (p*exp(-x) + q*sin(x) + r*cos(x) + s*tan(x) + t*x*x + u)const double eps = 1e-14;int main(){    int p, q, r, s, t, u, i;    while(~scanf("%d%d%d%d%d%d", &p, &q, &r, &s, &t, &u))    {        double f0 = F(0), f1 = F(1);        if(f0 < -eps || f1 > eps) printf("No solution\n");        else        {            double m, L = 0, R = 1;            for(i = 0; i < 100; i++)            {                m = L + (R-L) / 2;                if(F(m) >= 0) L = m;                else R = m;            }            printf("%.4lf\n", L);        }    }    return 0;}