二分算法

二分算法是解决形如，查找1-1000之间的某个数之类的查找问题。

如果从1到1000依次询问，则平均要用500次。

而如果采用二分的方法，每次取中间的数，每一次都能使查找范围缩小一半，则仅仅只需要10次！

 

二分查找函数模板

int binarysearch(int *a,int len,int p)
{
    int lo=0;//查找的左端点
    int hi=len-1; //查找的右端点
    while(lo<=hi){
        int mid=lo+(hi-lo)/2;
        if(p==a[mid])
            return mid;
        else if(p>a[mid])
            lo=mid+1;
        else
            hi=mid-1;
    }
    return -1;
}

时间复杂度

 

二次查找函数(小于p下标最大的元素)模板

数组a已从小到大排序

int lowersearch(int *a,int len,int p)
{
    int lo=0;
    int hi=len-1;
    int lastpost=-1;
    while(lo<=hi){
        int mid=lo+(hi-lo)/2;
        if(a[mid]>=p)
            hi=mid-1;
        else{
            lastpost=mid;
            lo=mid+1;
        }
    }
    return lastpost;
}

时间复杂度

 

注意

不要写成mid=(lo+hi)/2，要写成mid=lo+(hi-lo)/2。防止数过大而溢出。

 

它的一些应用

例题一 求下面方程的 一个根： 根：f(x) = x3 -5x2 +10x-80 = 0
若求出的根是 若求出的根是a ，则要求 |f(a)| <= 10 -6

因为求导之后f'(x)=3x 2 -10x+10，其无解所以它为单调递增函数，通过判断发现f(0)<0,f(100)>0，所以函数在[0,100]区间内存在一个解。此时二分就可以大显身手了！

#include<iostream>
#include<cstdio>
#include<cmath>
double eps=1e-6;
double f(double x)
{
    return x*x*x-5*x*x+10*x-80;
}
int main()
{
    double root,x1=0,x2=100,y;
    root=x1+(x2-x1)/2;
    y=f(root);
    int triedtime=1;
    while(fabs(y)>eps){
        if(y>0)
            x2=root;
        else
            x1=root;
        triedtime++;
        root=x1+(x2-x1)/2;
        y=f(root);
    }
    printf("%.8f\n",root);
    printf("%d\n",triedtime);
    return 0;
}

注 1e-6为。因为root为浮点数，而且题目分析出在[0,100]之间肯定存在一个答案，所以这里面没有用root-1，root+1而可以直接用x2=root,x1=root。怎么判断浮点数等于0 。fabs()表示给浮点数取绝对值，abs()表示给整型取绝对值。

 

例题二 输入n ( n<= 100000)个整数， 找出其中的两个数，它们之和等于整数m(假定肯定有解)。 题中所有整数都能用int 表示。

如果用二层for循环那么复杂度为，高达会超时。

那么可以用二分的方法。

解法1，先将存储这些数的数组排序，时间复杂度为。

之后对数组中的每个元素a[i],在数组中二分查找m-a[i]，看能否找到。复杂度log(n)，，最坏要查找n-2次，所以查找这部分的复杂度也是O(n×log(n))。

所以总的时间复杂度为。

解法2，先将存储这些数的数组排序，时间复杂度为。

之后设置两个变量，lo和hi，lo=0，hi=n-1。看a[lo]+a[hi]如果它的值大于m，那么hi--，如果值小于m，那会lo++，知道a[lo]+a[hi]=m。

时间复杂度为

 

例题三 百练 2456:Aggressive cows