Can you solve this equation?(简单二分)

没有分析函数的单调性,但是直觉告诉我这是单调函数,至于单增还是单减就没有判断了,直接二分。


#include<iostream>
#include<cstdlib>
#include<cstdio>
#define long long LL
const double dev = 1E-9;
//const double dev2= 1E-5;
using namespace std;

double dabs(double a)
{
    return a>0? a : -a;
}

double fun(double x,double y)<span style="white-space:pre">//构造方程为一个函数,如果函数值很接近于0,则x就是方程的解
{
    return 8*x*x*x*x + 7*x*x*x + 2*x*x + 3*x + 6 - y;
}

int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int flag = 0;
        double y;
        cin>>y;
        double low=0,high=100,mid=50;
        if(fun(0,y)/fun(100,y)>0)          //高中学过的二分法,怕溢出,用的除法
            cout<<"No solution!"<<endl,flag = 1;
        else
            while(dabs(low-high)>dev)        //控制精度
            {
                if(fun(low,y)/fun(mid,y) < 0)          //如果解在这个范围内
                {
                    high = mid;           //更改区间
                    mid = (low+high)/2;
                }
                else          //如果不在这个范围内,换一种方式更改解所在的区间
                {
                    low = mid;
                    mid = (low+high)/2;
                }
            }
        if(flag == 0)printf("%.4lf\n",mid);
    }
    return EXIT_SUCCESS;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值