题意:解方程
p*e-x + q*sin(x) + r*cos(x) + s*tan(x) + t*x2 + u = 0
where 0 <= x <= 1.(0 <= p,r <= 20 and -20 <= q,s,t <= 0)
——>>仔细观察可发现每一项在定义的域上都是减函数,所以,直接二分吧……
第一条数值题,学了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;
}