C++实现二分法求方程近似解
二分法是一种求解方程近似根的方法。对于一个函数 f(x)f(x)，使用二分法求 f(x)f(x) 近似解的时候，我们先设定一个迭代区间（在这个题目上，我们之后给出了的两个初值决定的区间 [-20,20]），区间两端自变量 x 的值对应的 f(x) 值是异号的，之后我们会计算出两端 x 的中点位置 x' 所对应的 f(x') ，然后更新我们的迭代区间，确保对应的迭代区间的两端 x 的值对应的 f(x) 值还会是异号的。
重复这个过程直到我们某一次中点值 x' 对应的 f(x') < f(x′)<ϵ （题目中可以直接用EPSILON）就可以将这个 x′ 作为近似解返回给 main 函数了。
例如：
上面所示的一个迭代过程的第一次的迭代区间是 [a1,b1]，取中点 b2，然后第二次的迭代区间是 [a1,b2]，再取中点 a2，然后第三次的迭代区间是 [a2,b2]，然后取 a3，然后第四次的迭代区间是[a3,b2]，再取红色中点 c，我们得到发现 f(c) 的值已经小于ϵ，输出c 作为近似解。
在这里，我们将用它实现对形如 px+q=0 的一元一次方程的求解。
在这里，你完成的程序将被输入两个正整数 p 和 q（你可以认为测评机给出的 0<∣p∣≤1000 且0<∣q∣≤1000），程序需要用二分法求出 px+q=0 的近似解。
输入格式
测评机会反复运行你的程序。每次程序运行时，输入为一行，包括一组被空格分隔开的符合描述的正整数 p 和 q。你可以认为输入数据构成的方程 px+q=0 都是有解且解在[−20,20] 的区间内。
输出格式
输出为一行，包括一个数字。为方程 px+q=0 的近似解。请使用四舍五入的方式保留小数点后 4 位小数。
样例输入1
55 9
样例输出1
-0.1636
样例输入2
-22 4
样例输出2
0.1818
代码：
#include <cstdio>
#include <cmath>
#define EPSILON 1e-7
double bisection(int p, int q, double (*func)(int, int, double));
double f(int p, int q, double x);
int main() {
int p;
int q;
scanf("%d %d", &p, &q);
printf("%.4lf\n", bisection(p, q, f));
return 0;
}
double bisection(int p, int q, double (*func)(int, int, double)) {
double x1 = -20;
double x2 = 20;
double x = 0;
while(fabs((*func)(p,q,x))>EPSILON)
{
x = (x1+x2)/2.0;
double fx1 = (*func)(p,q,x1);
double fx = (*func)(p,q,x);
if(fx*fx1>0)
{x1 = x;}
else
{x2 = x;}
}
return x;
}
double f(int p, int q, double x) {
return p * x + q;
}
以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持猪先飞。
相关文章
- vector是表示可以改变大小的数组的序列容器，本文主要介绍了C++STL标准库std::vector的使用详解，文中通过示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下...2022-03-06
- 这篇文章主要介绍了C++中取余运算的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧...2021-02-23
- 这篇文章主要介绍了C++ string常用截取字符串方法，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- 本文通过例子，讲述了C++调用C#的DLL程序的方法，作出了以下总结，下面就让我们一起来学习吧。...2020-06-25
- 本篇文章主要介绍了C++中四种加密算法之AES源代码，小编觉得挺不错的，现在分享给大家，也给大家做个参考。...2020-04-25
- 整数拆分，指把一个整数分解成若干个整数的和。本文重点给大家介绍C++ 整数拆分方法详解，非常不错，感兴趣的朋友一起学习吧...2020-04-25
- 这篇文章主要介绍了C++中Sort函数详细解析，sort函数是algorithm库下的一个函数，sort函数是不稳定的，即大小相同的元素在排序后相对顺序可能发生改变...2022-08-18
- 这篇文章主要介绍了C++万能库头文件在vs中的安装步骤(图文)，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧...2021-02-23
- 这篇文章主要介绍了C++ bitset用法，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- 本篇文章小编并不是为大家讲解string类型的用法，而是讲解我个人比较好奇的问题，就是string 类型占几个字节...2020-04-25
- 这篇文章主要为大家详细介绍了C++ Eigen库计算矩阵特征值及特征向量，具有一定的参考价值，感兴趣的小伙伴们可以参考一下...2020-04-25
- 这篇文章主要介绍了C++ pair的用法实例详解的相关资料,需要的朋友可以参考下...2020-04-25
- 这篇文章主要介绍了VSCode C++多文件编译的简单使用方法,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下...2021-03-29
- 虽然C++11引入了智能指针的，但是开发人员在与内存的斗争问题上并没有解放，如果我门实用不当仍然有内存泄漏问题，其中智能指针的循环引用缺陷是最大的问题。下面通过实例代码给大家介绍c++中的循环引用，一起看看吧...2020-04-25
- 这篇文章主要给大家介绍了关于C++随机点名生成器的相关资料，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- map容器是C++ STL中的重要一员，删除map容器中value为指定元素的问题是我们经常与遇到的一个问题，下面这篇文章主要给大家介绍了关于利用C++如何删除map容器中指定值的元素的相关资料，需要的朋友可以参考借鉴，下面来一起看看吧。...2020-04-25
- 这篇文章主要介绍了C++ 约瑟夫环问题案例详解,本篇文章通过简要的案例,讲解了该项技术的了解与使用,以下就是详细内容,需要的朋友可以参考下...2021-08-15
- 这篇文章主要介绍了C++中cin的用法详细,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- 本篇文章是对C++中的常见编译错误进行了详细的分析介绍，需要的朋友参考下...2020-04-25
- 在本篇内容里小编给大家分享了关于C++实现递归函数的教学步骤，需要的朋友跟着参考下。...2020-04-25