【每日刷题】 PAT 基础习题集 7-18 二分法求多项式单根 （20 分)

题目描述：


代码如下：

//7-18
#include <stdio.h>
double ret( double m, double a0, double a1, double a2, double a3 );
int main ( void )
{
	double a0, a1, a2, a3;
	double a, b, root, mid;
	scanf( "%lf%lf%lf%lf%lf%lf", &a3, &a2, &a1, &a0, &a, &b );
	while( b - a > 0.001 ){
		mid = ( b + a ) / 2;
		double fa = ret( a, a0, a1, a2, a3 );
		double fb = ret( b, a0, a1, a2, a3 );
		double fm = ret( mid, a0, a1, a2, a3 );
		if( fa * fb < 0  ){
			if( fm == 0 ){
				printf( "%.2f", mid );
				return 0;
			}
			else if( fm * fa > 0 )
				a = mid;
			else if( fm * fb > 0 )
				b = mid;
		}
		else if( fa * fb == 0 ){
			if( fa == 0 ){
				printf( "%.2f", a );
			return 0;
			}
			else if( fb == 0 ){
				printf( "%.2f", b );
				return 0;
			}
		}	
	} 
	printf( "%.2f", ( b + a ) / 2 );	
	return 0;
} 

double ret( double m, double a0, double a1, double a2, double a3 )
{
	return a3 * m * m * m + a2 * m * m + a1 * m + a0;
}

题目本身比较简单，但有两点需要注意：
（1）关于精度，由于题目要求保留小数点后两位，所以实际应该计算到后三位，所以区间长度应该大于0.001；
（2）在计算多项式时，不要使用pow函数，我之前就是使用pow函数，结果输出运行超时，虽然具体原因还不是很清楚（以后来补），但需要注意。