A - Can you solve this equation?
Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u
Submit
Status
Description
现在,给出等式8* X^4+ 7* X^3+ 2* X^2+ 3 * X +6= Y,请找出他在0和100之间的解(包含0和100)。
现在,请你试试运气。。。。
Input
输入的第一行包含一个整数T(1 <= T <=100),表示测试用例的数目。接下来T个数字,每一行都有一个实数Y(abs(Y)<=10^10);
Output
对于每个测试用例,如果有解,你应该输出一个实数(精确到小数点后4位,四舍五入),如果在0到100之间无解,就输出“No solution!”。
Sample Input
2
100
-4
Sample Output
1.6152
No solution!
- 题意:求上述方程在0-100方位内的解;
- 思路:求导可知,上述方程在给定范围内单调递增,符合二分的条件,对区间内的数进行枚举,找解(解方程的问题一般都用二分,还有在一个确定区间内的函数,一般都用枚举,只不过有序了就可以排除许多不必要枚举的元素了,了,这时枚举就成了二分);
- 失误:注意控制精度,浮点数不能直接用==;还有二分的方向,大了该向哪移动要想清楚;
- 代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
#define esp 1e-8
double judge(double mid)
{
return (8*pow(mid,4.0)+7*pow(mid,3)+2*pow(mid,2)+3*mid+6);
}
int main()
{
int t;
double y,mid,l,r,cnt;
cin>>t;
while(t--)
{
cin>>y;
if(y<6||y>judge(100.0))
cout<<"No solution!"<<endl;
else
{
r=100.0,l=0.0;
cnt=50;//每次少一半那就是2^-50倍 10^15级的,精确度足够
while(cnt--) //r-l>esp 精度控制x的误差小于10^-4
{
mid=(l+r)/2.0;
if(y>judge(mid))
{
l=mid;
}
else
{
r=mid;
}
}
printf("%.4lf\n",mid);
}
}
return 0;
}