二分法（整数和实数的二分）

是acwing的听课笔记

九.二分法（整数和实数的二分）

1.整数二分

本质：如果有单调性一定可以二分。但是能二分的题目不一定都要二分，所以二分本质并非单调性。

如果给定一个区间，可以找到一种区间使得范围能够被一分为二，则可以通过二分法寻找边界。假设左区间为红色，右区间为绿色。【五】给的两个二分模板即对应查找红或绿的情况。

红色情况：

设定中间值 mid=(l+r+1)/2。mid可以取到边界点

//表示检查mid是否满足给定的红色性质，true表示满足，否则false
if(check(mid)) true://则答案在区间[mid,r]，查找区间更新方式：l=mid
              false://则mid相当于取在了右边，不满足情况啊。区间更新为                          [l,mid-1]，更新方式为r=mid-1

绿色情况：

mid=(l+r)/2

//表示检查mid是否满足给定的绿色性质，true表示满足，否则false
if(check(mid)) true://则答案在区间[l,mid]，查找区间更新方式：r=mid
              false://则mid相当于取在了左边，不满足情况啊。区间更新为                          [mid+1,r]，更新方式为l=mid+1

如何去选择模板？

先写一个check模板，思考根据check的结果该如何更新区间。如果更新方式是l=mid，r=mid-1，则mid应等于(l+r+1)/2，选择模板一，否则二。

为什么要加上1？

因为相加/2的结果是向下取整。如果选择模板一，当l=r-1时，，mid=l+r/2=r-1=l。如果check成功，l=mid，则l实际上还是没变，最后就成了死循环。

例子：

给定一个已经按照升序排列的长度为n的整数数组，以及q个查询。对于每个查询，返回一个元素k的起始位置和终止位置，（位置从0开始计数），如果不存在则返回“-1，-1”

#include<iostream>
​
using namespace std;
​
const int N=100010;
​
int n,m,q[N];
​
int main(){
    scanf("%d%d",&n,&m);
    for(int i=0;i<n;i++) scanf("%d",&q[N]);
    
    while(m--){
        int x;
        scanf("%d",&x);
        
        int l=0,r=n-1;
        while(l<r){//二分查找
            int mid=l+r>>1;
            if(q[mid]>=x){//说明mid取到了右半边，更新为r=mid，模板二不加1
                r=mid;
            }
            else l=mid+1;
        }
        if(q[l]!=x) cout<<"-1,-1"<<endl;//说明数组中不存在x
        else{
            cout<<l<<' ';
            
            int l=0,r=n-1;
            while(l<r){//实际上两个模板都用到的话可以直接判断是否+1
                int mid=l+r+1>>1;
                if(q[mid]<=x) l=mid;//是模板一，补+1
                else r=mid-1;
            }
            cout<<l<<endl;
        }
    }
}

2.浮点数二分

例题：数的三次方根

给定一个浮点数，求他的三次方根。数据范围 -10000<=n<=10000，保留六位小数

首先确定区间[-10000,10000]。

mid=(l+r)/2,如果m³>=x，说明m取的大于最终答案了，则将区间更新为m左边即可。r=m，否则将l=m。

因为要保留六位小数，所以在写循环判定条件的时候，使r-l>1e-8（这是经验之谈，为了保证精确度，小数点要多往前挪两位）

#include <iostream>
​
using namespace std;
​
int main(){
    double x;
    cin>>x;
    
    double l=-10000,r=10000;
    while(r-l>1e-8){
        double mid=(l+r)/2;
        if(mid*mid*mid>=x)r=mid;
        else l=mid;
    }
    
    printf("%lf",l);
    
    return 0;
}