2017.4.23 一元三次方程求解 思考记录

        实在太菜了、、竟然不知道浮点数有误差、、、

          另外：修改的时候一次要改全！！！

码：

#include<iostream>
#include<cstdio>
using namespace std;
#include<cmath>
double a,ans[6],b,c,d,j,i;
int qushi,cnt;

int main()
{
	scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
	for(i=-100;i<=100;i++)
	{//	cout<<i<<" "<<i*(i*(a*i+b)+c)+d<<"pp";
		if(i==-100)
		{
			if(i*(i*(a*i+b)+c)+d>0)qushi=1;
		//	printf("%.5lf",fabs(i*(i*(a*i+b)+c)+d));
			if(fabs(i*(i*(a*i+b)+c)+d)<=0.0001){ans[++cnt]=i;}
			if(i*(i*(a*i+b)+c)+d<0)qushi=0;	
	           continue;
		}
		if(fabs(i*(i*(a*i+b)+c)+d)<=0.001)
		{ans[++cnt]=i;continue;}
		if((i*(i*(a*i+b)+c)+d>0&&qushi==0)||(i*(i*(a*i+b)+c)+d<0&&qushi==1))
		{  qushi^=1;
			j=double(i-1)+0.001;
		    while(j<i)
		    {   if(fabs(j*(j*(a*j+b)+c)+d)<=0.001)
		        {
		        	ans[++cnt]=j;break;
				}	    
		    	if(((j*(j*(a*j+b)+c)+d)>0&&((j+0.001)*((a*j+0.001)*((j+0.001)+b)+c)+d)<0)||((j*(j*(a*j+b)+c)+d)<0&&((j+0.001)*((a*j+0.001)*((j+0.001)+b)+c)+d)>0))
		    	{
                       ans[++cnt]=j;break;		    
				}
		    	j+=0.001;
			}
		}
	 } 
	printf("%.2lf %.2lf %.2lf",ans[1],ans[2],ans[3]);
 }