第三周练习——二分法2 方程求解

Problem B

Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 32768/32768K (Java/Other)

Total Submission(s) : 95   Accepted Submission(s) : 37

Problem Description

Now,given the equation 8*x^4 + 7*x^3 + 2*x^2 + 3*x + 6 == Y,can you find its solution between 0 and 100; Now please try your lucky.

 

Input

The first line of the input contains an integer T(1<=T<=100) which means the number of test cases. Then T lines follow, each line has a real number Y (fabs(Y) <= 1e10);

 

Output

For each test case, you should just output one real number(accurate up to 4 decimal places),which is the solution of the equation,or “No solution!”,if there is no solution for the equation between 0 and 100.

 

Sample Input

  
  

   
   2
100
-4

  
  

 

Sample Output

  
  

   
   1.6152
No solution!

  
  
  
  


  
  
  
  

   
   题目分析：
  
  
  
  

   
       1）这是简单的二分求解方程的题目，求导可知，方程的图像为一个增函数，因而可以用二分法
  
  
  
  

   
       2）题目总结：  
  
  
  
  

   
              1.实数输出时，用f替代lf（GCC的编译环境下）
  
  
  
  

   
              2.错误原因： break的使用出现错误 
  
  
  
  

   
                     原先时，在 NO solution 语句后面添加Break 语句，
  
  
  
  

   
                     将break语句删掉以后，用if  else 语句替代，出现了正确结果
  
  
  
  

   
           稍后分析原因
  
  
  
  


  
  
  
  

   
   代码如下
  
  
  
  


  
  
  
  

   
   #include<stdio.h>
#include<math.h>
#define eps 1e-9
double f(double x,double y)
{

   
   	return 8*pow(x,4)+7*pow(x,3)+2*pow(x,2)+3*x+6-y;
}
int main()
{

   
   	int t;

   
   	scanf("%d",&t);

   
   	double y,l_val,r_val,mid_val,mid;

   
   	while(t--)

   
   	{

   
   		scanf("%lf",&y);

   
   		double l=0.0,r=100.0;

   
   		l_val=f(l,y);

   
   		r_val=f(r,y);

   
   		if(l_val* r_val>0)

   
   		{

   
   			printf("No solution!\n");

   
   		

   
   		}
        else

   
   		{

   
   		while(r-l>eps)

   
   		{

   
   		  mid=(l+r)/2;

   
   		  mid_val=f(mid,y);
          l_val=f(l,y);

   
   		  r_val=f(r,y);

   
   		  if(mid_val==0)

   
   			  break;

   
   		  else

   
   		  {

   
   			  if(l_val*mid_val<0)

   
   				  r=mid;

   
   			  else

   
   				  l=mid;

   
   		  }

   
   		  

   
   		}



   
   		printf("%.4f\n",mid);

   
   		}







   
   	}

   
   	return 0;
}