P1024【一元三次方程求解】

一元三次方程求解

题目

P1024

---

解析

这道题其实应该输出12位精度，开longdouble，结果只有2位精度，一堆直接枚举的，都不需要题目性质了
既然已保证两个实根之间距离不会小于1，那么我们枚举每一个i~i+1的区间(i是整数)，若f(i)==0则输出i，若区间两端f相乘<0则在区间内二分，若f(l)*f(mid)<0则r=mid，反之l=mid，最后直接输出右端点即可

code:

#include<cstdio>
#define db double
using namespace std;
db a,b,c,d,l,r,mid,eps=0.001;
inline db F(db x){return (((a*x+b)*x+c)*x)+d;}
int main()
{
	scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
	for(db i=-100;i<=100;++i)
	{
		if(F(i)==0)printf("%.2lf ",i);
		else if(F(i)*F(i+1)<0)
		{
			l=i,r=i+1;
			while(l+eps<r)
			{
				mid=(l+r)/2;
				if(F(mid)*F(l)<0)r=mid;
				else l=mid;
			}
			printf("%.2lf ",r);
		}
	}
	return 0;
}