P1024 [NOIP2001 提高组] 一元三次方程求解
题目描述
有形如:a x^3 + b x^2 + c x + d = 0a**x3+b**x2+c**x+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,da,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在 -100−100 至 100100 之间),且根与根之差的绝对值 \ge 1≥1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后 22 位。
提示:记方程 f(x) = 0f(x)=0,若存在 22 个数 x_1x1 和 x_2x2,且 x_1 < x_2x1<x2,f(x_1) \times f(x_2) < 0f(x1)×f(x2)<0,则在 (x_1, x_2)(x1,x2)之间一定有一个根。
输入格式
一行,44 个实数 a, b, c, da,b,c,d。
输出格式
一行,33 个实根,从小到大输出,并精确到小数点后 22 位。
输入输出样例
输入 #1复制
1 -5 -4 20
输出 #1复制
-2.00 2.00 5.00
package com.kk.luogu;
import java.util.Scanner;
//解一元三次方程组
//这里使用枚举 暴力的方式
public class P1024 {
public static void main(String[] args) {
Scanner scanner = new Scanner(System.in);
double a=scanner.nextDouble();
double b=scanner.nextDouble();
double c=scanner.nextDouble();
double d=scanner.nextDouble();
for (double i = -100; i <=100 ; i+=0.001) {//枚举每个答案
if (Math.abs(i*i*i*a+i*i*b+i*c+d)<0.0001)//一个一个的测试
System.out.println(i);
}
}
}
#include<bits/stdc++.h>
using namespace std;
double a,b,c,d;
int main(){
cin>>a>>b>>c>>d; // 输入
for(double i=-100;i<=100;i+=0.001){//枚举每个答案
if(fabs(i*i*i*a+i*i*b+i*c+d)<0.0001)//避免double精度错误
printf("%.02lf ",i);//两位小数输出
}
return 0;
}