问题 D: 解方程

题目描述

求方程f(x)=2^x+3^x-4^x=0在[1，2]内的根，精确到10位小数。

输入

输入m(0<=m<=8)，输出m位小数

输出

输出方程f(x)=0的根，精确到10位小数，输出m位小数。

样例输入 

3

样例输出 

1.507

问题分析

核心思想：使用二分搜索找根

情况1：如果mid代入后函数值为0，即f(l)=0，则说明是根，循环结束

情况2：如果f(l)*f(mid)<0，那么说明l和mid是异侧的，也就是根在左区间，r=mid。

情况3：如果f(l)*f(mid)>0，那么说明l和mid是同侧的，也就是根在右区间，l=mid。

#include <bits/stdc++.h>
using namespace std;

// 定义函数 f，它接受一个双精度浮点数 x 并返回 2^x + 3^x - 4^x 的值。
double f(double x) {
    return pow(2, x) + pow(3, x) - pow(4, x);
}

int main() {
    // 声明整数变量 m 用于存储小数位数
    int m;
    cin >> m; // 从标准输入读取 m 的值

    // 初始化二分搜索的边界，l 为左边界，r 为右边界，mid 为中间值
    double l = 1, r = 2, mid;

    // 当右边界和左边界的差大于 10^-10 时继续循环
    // 这保证了我们找到的解的精度至少为 10 位小数
    while (r - l > pow(10, -10)) {
        mid = (l + r) / 2; // 计算中点

        // 如果在中点 f(mid) 的值为 0，那么我们找到了方程的根
        if (f(mid) == 0) {
            l = mid;
            break;
        } 
        // 如果 f(l) 和 f(mid) 乘积小于 0，说明根在左半区间
        else if (f(l) * f(mid) < 0) {
            r = mid;
        } 
        // 否则，根在右半区间
        else {
            l = mid;
        }
    }

    // 输出找到的根，精度为 m 位小数
    // fixed 和 setprecision 用于控制输出格式
    cout << fixed << setprecision(m) << l;
}