解题思考F. 解方程(二分法解方程)
题目描述
Problem Description
给定方程 8x^4 + 7x^3 + 2x^2 + 3x + 6 == Y,请计算x在[0,100]范围内的解。
Input
输入数据首先是一个正整数T(1<=T<=100),表示有T组测试数据。
接下来T行,每行包含一个实数Y ( fabs(Y) <= 1e10 )。
Output
请计算并输出方程在范围[0,100]内的解,结果精确到小数点后4位。
如果无解,则请输出“No solution!”
输入样例
2
100
-4
输出样例
1.6152
No solution!
思考:二分控制解方程,设置一头一尾为范围,然后逐步取中间值比较来控制精度。
代码如下:
#include<bits/stdc++.h>
using namespace std;
double equation(double x);
int main()
{int n;cin>>n;double y;double left,right,mid;while(n--){cin>>y;if(equation(0)<=y&&y<=equation(100))//判断方程是否有解 {left=0;right=100;while(right-left>1e-6)//精度控制 {mid=(left+right)/2;double ans=equation(mid);if(ans>y)right=mid-1e-7;elseleft=mid+1e-7;}printf("%.4lf\n",(left+right)/2);}elsecout<<"No solution!"<<endl;}
}
double equation(double x)
{return 8*pow(x,4.0)+7*pow(x,3.0)+2*pow(x,2.0)+3*x+6;
}
解题思考F. 解方程(二分法解方程)相关推荐
- 二分法解具有单调性的方程
解方程的手段有很多,如二分法.牛顿迭代法等等,本次介绍的是采用二分法,在使用二分法解方程时,有一点应该特别注意,就是那种"具有单调性"的函数才可以,否则是会有问题的 ** HDU ...
- 5.C语言二分法求方程解
用二分法求方程解 #include <stdio.h> #include <math.h> double f(double x){return x*x-2*x-1; } dou ...
- python二分法求方程的根_Python查找函数f(x)=0根的解决方法
线性代数分享方程f(x)=0的根 函数F(x)=0的重根与F'(x)=0的根有什么关系?有些人一旦错过了,就是一辈子不再主动联系,不愿打扰你的生活,连偶尔的寒暄都没有,成长就是这样的,不断的告别,不断 ...
- 如何用matlab解异或方程,Matlab-6:解非线性方程组newton迭代法
函数文件: function x=newton_Iterative_method(f,n,Initial) x0=Initial; tol=1e-11; x1=x0-Jacobian(f,n,x0)\ ...
- C语言 牛顿法 解方程,牛顿法解方程
数值法解方程广泛应用于各种工程科学领域,它可以求出解的近似值,而对于高次方程没有求根公式,而且对于各种困难的方程,数值法就显得尤为重要. 我们的老师向我们介绍了"二分法",这种方法 ...
- scipy.optimize.fsolve:用Python求解方程的解
例1: 求解方程组的一个解: x0*cos(x1) = 4, x1*x0 - x1 = 5. 需要注意两点: 1.定义方程组,方程组要写出f(x)=0的形式(=0不需要写出来),所以原方程右边4和5都 ...
- NOIP模拟测试7「方程的解·visit」
visit 由于一些不可预知的错误导致我一直WA 错误最后说 思路 方案一 假设终点在出发点右上方(这样假设只是为了方便) 假设向左走了a步,向右下了b布,那么相应的我们要向右走m+a,向上n+b步 ...
- 1308 方程的解(组合计数--隔板法)
1. 问题描述: 佳佳碰到了一个难题,请你来帮忙解决.对于不定方程 a1+a2+⋯+ak−1+ak=g(x),其中 k ≥ 1 且 k∈N∗,x 是正整数,g(x)=x ^ x mod 1000(即 ...
- 4.1 简单方程的解
前面章节介绍了向量组和矩阵理论,利用这些理论可以解决线性方程 Amnx=bA_{mn}\mathbf{x}=\mathbf{b}Amnx=b 解的存在性和唯一性问题.向量组理论如下:向量 b\mat ...
最新文章
- Eigen向量化内存对齐/Eigen的SSE兼容,内存分配/EIGEN_MAKE_ALIGNED_OPERATOR_NEW
- Fiddler之弱网测试(Web)
- resoult在python啥意思,python - __lt__而不是__cmp__
- 阿里巴巴发送短信接口的使用
- 输出1/n(是循环小数的,只输出第一个循环节)
- JVM(2)——JVM类加载机制
- 使用Google Test的一个简单例子
- jQuery水印插件 - Watermark 和 FormWatermark
- ndarray是什么_python数据分析用什么软件?
- Ubuntu 左边栏和顶栏都不见了,ctrl+alt+t 也调用不出terminal
- bzoj2146 Construct
- WineQQ2012 最新下载
- 闪客工具:好玩的GifCam录制gif软件
- LeetCode - 1419 数青蛙
- FPGA之IP核的使用、搞懂功能仿真和时序仿真
- [创业之路-50] :动态股权机制 -5- 创业公司股权分配常见的坑
- 如何实现外网访问内网ip?公网端口映射或内网映射来解决
- 优化器:torch.optimizer
- Matlab学习(台大郭彦甫)第5节-初阶绘图
- 打印机显示服务器磁盘已满,打印机无法打印显示内存已满是怎么回事,怎么解决?...