给定方程 8*x^4 + 7*x^3 + 2*x^2 + 3*x + 6 == Y,请计算x在[0,100]范围内的解。
输入数据首先是一个正整数T(1<=T<=100),表示有T组测试数据。
接下来T行,每行包含一个实数Y ( fabs(Y) <= 1e10 )。
请计算并输出方程在范围[0,100]内的解,结果精确到小数点后4位。
如果无解,则请输出“No solution!”
2 100 -4
1.6152 No solution!
#include
#include
#include
#include
using namespace std;
double fun(double x) {
return 8 * pow(x, 4.0) + 7 * pow(x, 3.0) + 2 * pow(x, 2.0) + 3 * x + 6;
}
int main()
{
int T; cin >> T;
while (T--) {
double Y; cin >> Y;
double left, right, mid;
if (fun(0) > Y || fun(100) < Y) {
cout << "No solution!" << endl;
}
else {
left = 0; right = 100;
while (right - left > 1e-8) {
mid = (right + left) / 2;
if (fun(mid) == Y) { break; }
else if (fun(mid) < Y) left = mid;
else right = mid;
}
printf("%.4f\n", mid);
}
}
return 0;
}