欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

洛谷P1024一元三次方程求解题解--zhengjun

程序员文章站 2022-05-09 09:57:36
...

题目描述

有形如:ax3+bx2+cx1+dx0=0ax^3+bx^2+cx^1+dx^0=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,da,b,c,d均为实数),并约定该方程存在三个不同实根(根的范围在100-100100100之间),且根与根之差的绝对值 1\ge 1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后22位。

提示:记方程f(x)=0f(x)=0,若存在22个数x1x_1x2x_2,且x1<x2x_1<x_2f(x1)×f(x2)<0f(x_1) \times f(x_2)<0,则在(x1,x2)(x_1,x_2)之间一定有一个根。

输入格式

一行,44个实数A,B,C,DA,B,C,D

输出格式

一行,33个实根,并精确到小数点后22位。

输入输出样例

输入 #1 复制
1 -5 -4 20
输出 #1 复制
-2.00 2.00 5.00

思路

难得的一道水题,因为题中说根和根之差的绝对值 1\ge 1,也就是说,在一段长度为 11 的区间中,不会有两个根。那么枚举这段区间在按照题中的方法做就是了

代码

#include<bits/stdc++.h>
using namespace std;
double a,b,c,d;
double f(double x){
    return a*x*x*x+b*x*x+c*x+d;
}
int tot;//用于计数
int main()
{
	cin>>a>>b>>c>>d;
	for(int i=-100;i<=100;i++){//枚举区间[i,i+1)
		if(tot>=3)//已经找到三个根
		    return 0;
		if(f(i)==0){//如果i本身就是一个根
			printf("%0.2lf ",double(i));
			tot++;
			continue;
		}
		if(f(i)*f(i+1)<0){//这段区间中有一个根
		    double l=i,r=i+1,m;
		    while(l+0.001<r){//因为要精确到2位小数,所以要算到3位小数
		    	m=(l+r)/2;
		    	if(f(l)*f(m)<0)//根在左边区间
		    	    r=m;
		    	else//根在右边区间
		    	    l=m;
			}
			printf("%0.2lf ",l);
			//因为r始终在(x1,x2)当中,所以根不在r
			tot++;
		}
	}
	//前面已经结束了,这里就不用了(我可能有点懒)
}

谢谢–zhengjun