标签：lf% 一元 方程 实根 double 分治 P1024 端点 include

题目描述:有形如：ax^3+bx^2+cx^1+dx^0=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d均为实数)，并约定该方程存在三个不同实根(根的范围在-100至100之间)，且根与根之差的绝对值≥1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位。

提示：记方程f(x)=0，若存在2个数x1​和x2​，且x1​<x2​，f(x1​)×f(x2​)<0，则在(x1​,x2​)之间一定有一个根。

输入格式:一行，4个实数A,B,C,D

输出格式:一行，3个实根，并精确到小数点后2位。

题目的数据范围不大，且根与根的绝对值之差大于等于1，因此可以考虑枚举宽度为1的区间。只要这个区间的左端点和右端点的函数值乘积小于零（即两个函数值异号），则其中必定有一个根。然后可以使用二分算法，即每次取区间中点，如果中点和左端点的函数值乘积小于零，则在(左端点，中点)间一定有一个根。再将中点赋值为右端点，继续上述操作。否则，在(中点，右端点)间有根，将中点赋值为左端点，继续上述操作。直到左端点和右端点的差小于0.001（题目要求精确到小数点后2位），此时输出任意一个端点。这里使用右端点。注意：如果在枚举区间时，左端点的函数值恰好等于0，则直接输出，不需再二分。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
double a,b,c,d;
double func(int x)
{
    return a*x*x*x+b*x*x+c*x+d;//求函数值
}
int main()
{
    double left,right,mid;
    int p=0;
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);//输入系数
    for(double i=-100.00;i<=100.00;i+=1.00)
    {
        double j=i+1.00;//枚举区间
        if(func(i)==0.00) //左端点就是一个解
        {
            printf("%.2lf ",(double)i);
            p++;
        }
        if(func((double)i)*func((double)j)<0)//发现区间(left,right)中有解
        {
            left=i;right=j;
            while(right-left>=0.001)//二分
            {
                mid=(left+right)/2;//取中间值
                if(func(mid)*func(left)<0)//判定解的区域为(left,mid)
                {
                    right=mid;
                }
                else//判定解的区域为(mid,right)
                {
                    left=mid;
                }
            }
            printf("%.2lf ",right);
            p++;
        }
        if(p==3) break;//一元三次方程至多有三个解
    }
    return 0;
}

标签：lf%,一元,方程,实根,double,分治,P1024,端点,include
来源： https://www.cnblogs.com/jiangyuechen/p/12457871.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。