P1024 一元三次方程求解

最新推荐文章于 2024-10-27 21:39:27 发布

MidJuly~

最新推荐文章于 2024-10-27 21:39:27 发布

阅读量114

点赞数

分类专栏：题解

本文链接：https://blog.csdn.net/plutozxc/article/details/108976507

版权

题解专栏收录该内容

23 篇文章 1 订阅

订阅专栏

今天来写一篇入门级的二分
题目
P1024 一元三次方程求解

在高中的时候就有学过这个方法，二分法。
题目也给了提示，相乘小于0的区间有解，或者直接等于0。
相乘小于0就可以用二分不断的缩小区间，求到最精确的一个数。从而得到答案。

#include<bits/stdc++.h>
using namespace std;
const int N = 1001;
typedef long long LL;
#define db double
int n, ans;
double a, b, c, d;
db f(db x) {
  return a*x*x*x+b*x*x+c*x+d;
}
db dfs(db l, db r) {
  if (r-l < 0.001) return l;
  db mid = (l+r)/2;
  if (f(l)*f(r) <= 0) return dfs(l, mid);
  else return dfs(mid, r);
} 
int main() {
  scanf("%lf%lf%lf%lf", &a, &b, &c, &d);
  double ans[5];
  int cnt = 0;
  for (int i = -100; i < 100; i++) {
    db l = i, r = i+1;
    if (f(l) == 0) ans[++cnt] = l;
    else if (f(l)*f(r) < 0) ans[++cnt] = dfs(l, r);
    if (cnt >= 3) break;
  }
  printf("%0.2lf %0.2lf %0.2lf\n",ans[1],ans[2],ans[3]);
  return 0;
  return 0;
}