码迷,mamicode.com
首页 > 其他好文 > 详细

UVA 10341 Solve It

时间:2014-04-28 19:00:17      阅读:1878      评论:0      收藏:0      [点我收藏+]

标签:com   http   class   blog   style   img   div   code   java   javascript   string   

Problem F
Solve It
Input: standard input
Output: standard output
Time Limit: 1 second
Memory Limit: 32 MB
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

首先对函数求导,发现是一个递减函数,这样就可以用二分求答案了,不过这个题精度要求好高啊,必须到10^-8次啊,一开始弄了一个10^-6,结果WA了,看来很久没有感觉有错误,就改了一下8就通过了,晕啊!

mamicode.com,码迷
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

const double esp=1e-8;

double p,q,r,s,t,u;

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

double find(double x,double y)
{
    while (y-x>esp)
    {
        double m=x+(y-x)/2;
        if (get_value(m)-esp>0) x=m;
        else y=m;
    }
    return x;
}

int main()
{
    while (scanf("%lf%lf%lf%lf%lf%lf",&p,&q,&r,&s,&t,&u)!=EOF)
    {
        if (get_value(0.0)<0 || get_value(1.0)>0)
        printf("No solution\n");
        else
        printf("%.4lf\n",find(0.0,1.0));
    }
    return 0;
}
mamicode.com,码迷

 

UVA 10341 Solve It,码迷,mamicode.com

UVA 10341 Solve It

标签:com   http   class   blog   style   img   div   code   java   javascript   string   

原文地址:http://www.cnblogs.com/chensunrise/p/3694254.html

(0)
(1)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!