Vijos P1116 一元三次方程求解【多解,暴力,二分】

简介: 一元三次方程求解 描述 有形如:ax^3+bx^2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间),且根与根之差的绝对值>=1。

一元三次方程求解

描述

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

格式

输入格式

输入该方程中各项的系数(a,b,c,d 均为实数),

输出格式

由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后2位。

样例1

样例输入1

1 -5 -4 20

样例输出1

-2.00 2.00 5.00

限制

每个测试点1s

提示

提示:记函数f(x),若存在2个实数x1和x2,且x1<x2,使f(x1)*(x2)<0,则在(x1, x2)之间一定存在实数x0使得f(x0)=0。

来源

NOIP2001第一题

题目链接:https://vijos.org/p/1116

分析:又来一道暴力题,听说这是一道很经典的二分?贴个二分的代码啊,暴力肯定可以,但是但是,一定要注意精度QAQ

暴力代码:【一个三元函数的性质,不懂的自己翻高中课本】

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 int main()
 4 {
 5     double a,b,c,d;
 6     cin>>a>>b>>c>>d;
 7     double j=-100;
 8     double i=j;
 9     for(j;j<=100;j+=0.01)
10     {
11         if(a*i*i*i+b*i*i+c*i+d<0)
12         {
13             if(a*j*j*j+b*j*j+c*j+d>0)
14             {
15                 printf("%.2lf ",j);
16                 i=j;
17             }
18         }
19         if(a*i*i*i+b*i*i+c*i+d>0)
20         {
21             if(a*j*j*j+b*j*j+c*j+d<0)
22             {
23                 printf("%.2lf ",j);
24                 i=j;
25             }
26         }
27     }
28     printf("\n");
29     return 0;
30 }

直接贴了二分代码:

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <algorithm>
 4 #include <cstring>
 5 #include <iomanip>
 6 #include <cstdlib>
 7 using namespace std;
 8 
 9 float ans[5];
10 float a,b,c,d;
11 int n=0;
12 
13 float f(float x)
14 {
15     return ((a*x+b)*x+c)*x+d;
16 }
17 
18 void solve(float l,float r)
19 {
20     if(f(l)*f(r)>0&&(((r-l)<1)||n>=2))
21         return;
22     float mid=(l+r)/2;
23 
24     if(f(mid)<=1e-4 && f(mid)>=-1e-4)
25     {
26         ans[n++]=mid;
27         return;
28     }
29 
30     solve(l,mid),solve(mid,r);
31 }
32 
33 int main()
34 {
35     cin>>a>>b>>c>>d;
36     solve(-100,100);
37     printf("%.2lf %.2lf %.2lf",ans[0],ans[1],ans[2]);
38     return 0;
39 }

 

目录
相关文章
|
5月前
|
算法 测试技术 C++
【数学归纳法 组合数学】容斥原理
【数学归纳法 组合数学】容斥原理
|
5月前
|
算法 测试技术 C#
【状态压缩 动态规划 数论】1799. N 次操作后的最大分数和
【状态压缩 动态规划 数论】1799. N 次操作后的最大分数和
|
机器学习/深度学习 存储 算法
回溯法求解N皇后问题
回溯法求解N皇后问题
130 0
蓝桥杯:暴力求解四平方和
蓝桥杯:暴力求解四平方和
55 0
|
存储 缓存 移动开发
动态规划法
动态规划是在20世纪50年代由美国数学家贝尔曼为研究最优控制问题而提出的,当该方法在应用数学中的价值被大家认同以后,在计算机学界,动态规划法成为一种通用的算法设计技术用来求解多阶段决策最优化问题。
|
人工智能 算法
【动态规划】矩阵连乘
完全加括号的矩阵连乘积可递归地定义为: • 单个矩阵是完全加括号的 • 矩阵连乘积A是完全加括号的,则A可表示为2个完全加括号的矩阵连乘积B和C 的乘积并加括号,即A=(BC) 设有四个矩阵A, B, C, D ,它们的维数分别是: A = 50*10 B = 10*40 C = 40*30 D = 30*5 总共有五种完全加括号的方式:
304 0
|
算法
数学:组合数算法模板
数学:组合数算法模板
185 0
553. 最优除法 : 数学类贪心运用题
553. 最优除法 : 数学类贪心运用题
|
机器学习/深度学习
【组合数学】递推方程 ( 常系数线性非齐次递推方程求解 | 递推方程标准型及通解 | 递推方程通解证明 )
【组合数学】递推方程 ( 常系数线性非齐次递推方程求解 | 递推方程标准型及通解 | 递推方程通解证明 )
188 0
|
算法 Java Python
【算法导论】最大子数组——暴力求解
1. 暴力方法求解最大子数组问题:求出所有子数组的和并比较; 2. 伪代码 FIND-MAXIMUM-SUBARRAY(A) max = -∞ for i = 1 to A.
1431 0