浮点数二分
用二分的思想可以找出这个分界点。
假设我们想二分出来红色(不满足条件的范围)这边的分界点:
- 找一个中间值mid :l+r >> 1(因为是浮点数,所以可以严格的二分,即mid就是中点)
- 写一个函数check()判断时mid时在满足条件(绿色)范围内,还是在不满足条件(红色)范围内。
- 当区间长度很小的时候,可以近似的认为找到了答案,例:r-l ≤ 10 − 6 ^{-6} −6 (如果题目要求保留4位小数,则应该写到-6次方,比有效小数的位数多2)
例题:用二分法做一个sqrt函数 (只能计算≥1的开方)
#include<iostream>
#include<cstdio>
#define IOS ios::sync_with_stdio(false); cin.tie(0); cout.tie(0)
using namespace std;
int main()
{
double x;
cin >> x;
double l = 0 , r = x;
//如果题目要求保留4位小数,则应该写到-6次方,比有效小数的位数多2
while (r - l > 1e-6)
//还可以写成 for(int i = 0; i < 100; ++i) 循环100次相当于整个区间/ 2的100次方,最后的结果也是很精确的。
{
double mid = (l + r) / 2;
if (mid * mid >= x)
r = mid;
else
l = mid;
}
printf("%lf\n",l);
return 0;
}