UVA 10341 Solve It (解方程 二分查找+精度)

来源:互联网 发布:stc51单片机的优点 编辑:程序博客网 时间:2024/05/16 07:56

Solve It

Input: standard input

Output: standard output

Time Limit: 1 second

Memory Limit: 32 MB

Solve the equation:
        p*e-xq*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: pqrst 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

题目大意:给你p,q,r,s,t,u六个变量的数值,要求你求出在区间(1,0)内,求x的值。

解析:对函数求导后导函数>0,可得函数在(1,0)上是单调递增的,所以直接用二分查找法求解。


#include <iostream>#include <string.h>#include <stdio.h>#include <math.h>using namespace std;double p,q,r,s,t,u;double f(double x) {double result = p*exp(-x) + q*sin(x) + r*cos(x) + s*tan(x) + t*x*x + u;return result;}int main() {double ans;double result;while(scanf("%lf %lf %lf %lf %lf %lf",&p,&q,&r,&s,&t,&u) != EOF) {double mid,left,right;left = 0, right = 1;if( f(left) * f(right) > 0) {printf("No solution\n");continue;}while(left < right) {mid = (left + right) / 2;result = f(mid);if(fabs(result) < 1e-6) {printf("%.4lf\n",mid);break;}if( result > 0) {left = mid;}else {right = mid;}}}return 0;}


0 0
原创粉丝点击