hdu 2199 Can you solve this equation?

最新推荐文章于 2020-11-11 19:50:07 发布

永远热血沸腾

最新推荐文章于 2020-11-11 19:50:07 发布

阅读量329

点赞数

CC 4.0 BY-SA版权

分类专栏：数学

本文链接：https://blog.youkuaiyun.com/xinag578/article/details/45437299

数学专栏收录该内容

56 篇文章

订阅专栏

本文介绍了一种在一维区间内寻找方程根的数值计算方法。通过输入方程的目标值，采用二分法逐步缩小可能包含根的区间范围，最终达到较高的精度找到满足条件的根。该算法适用于形如8x^4 + 7x^3 + 2x^2 + 3x + 6 = y类型的方程求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<queue>
#include<algorithm>
#include<string.h>
#include<iostream>
using namespace std;

double y;
double fcn(double x)
{
    double t=8*x*x*x*x+7*x*x*x+2*x*x+3*x+6-y;
    return t;
}

int main()
{
    int t;
    double x1,x2,x;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf",&y);
        x1=0,x2=100;

        if(fcn(x1)>0&&fcn(x2)>0)
        {
            printf("No solution!\n");
            continue;
        }
        if(fcn(x1)<0&&fcn(x2)<0)
        {
            printf("No solution!\n");
            continue;
        }

        if(fcn(x1)==0)
        {
            printf("%.4f\n",x1);
            continue;
        }
        if(fcn(x2)==0)
        {
            printf("%.4f\n",x2);
            //printf("%.2lf\n",fcn(x2));
            continue;
        }
        int f=0;
        while(x2-x1>1e-10)
        {
            x=(x1+x2)/2;
            double t1,t2,t3;
            t1=fcn(x1);
            t2=fcn(x);
            t3=fcn(x2);
            if(fcn(x1)==0)
            {
                f=1;
                break;
            }
            if(fcn(x2)==0)
            {
                f=2;
                break;
            }
            if(t1>0&&t2<0) x2=x;
            else if(t1<0&&t2>0) x2=x;
            else if(t2>0&&t3<0) x1=x;
            else if(t2<0&&t3>0) x1=x;
            //printf("%.4f\n",fcn(x1));
        }
        if(f==1) printf("%.4f\n",x1);
        else printf("%.4f\n",x2);
    }
    return 0;
}