一:起因
(1)STL中关于二分查找的函数有三个:lower_bound 、upper_bound 、binary_search —— 这三个函数都运用于有序区间(当然这也是运用二分查找的前提),下面记录一下这两个函数;
(2)ForwardIter lower_bound(ForwardIter first, ForwardIter last,const _Tp& val)算法返回一个非递减序列[first, last)中的第一个 大于等于值val的位置;(3)ForwardIter upper_bound(ForwardIter first, ForwardIter last, const _Tp& val)算法返回一个非递减序列[first, last)中的第一个大于值val的位置。
二:lower_bound和upper_bound如下图所示:
(1)lower_bound函数源代码://这个算法中,first是最终要返回的位置分析:既然lower_bound()函数返回的是first的指针,那就看first指针变化的位置——first = middle + 1;的条件是 严格的小于(<),即大于等于(>=)时first的位置是不变的。或者可以这样理解,当middle指向的值 < key时,first 指向middle的下一个位置,而下一个位置可能等于或者大于key;因此 该函数是返回一个非递减序列[first, last)中的第一个大于等于值val的位置。
int lower_bound(int *array, int size, int key)
{
int first = 0, middle;
int half, len;
len = size;
while(len > 0) {
half = len >> 1;
middle = first + half;
if(array[middle] < key) {
first = middle + 1;
len = len-half-1; //在右边子序列中查找
}
else
len = half; //在左边子序列(包含middle)中查找
}
return first;
}
(2)upper_bound函数源代码:
int upper_bound(int *array, int size, int key)
{
int len = size-1;
int half, middle;
while(len > 0){
half = len >> 1;
middle = first + half;
if(array[middle] > key) //中位数大于key,在包含last的左半边序列中查找。
len = half;
else{
first = middle + 1; //中位数小于等于key,在右半边序列中查找。
len = len - half - 1;
}
}
return first;
}
分析:既然upper_bound()函数返回的是first的指针,那就看first指针变化的位置——first = middle + 1;的条件是 小于等于(<=),即只有严格大于(>)时first的位置是不变的。或者可以这样理解,当middle指向的值 <= key时,first 指向middle的下一个位置,而下一个位置可能等于或者大于key,但是要是仍是等于的话,first会继续指向middle的下一位置,直到大于为止;因此 该函数是返回一个非递减序列[first, last)中的第一个大于值val的位置。
(3)binary_search函数源代码
//int BinSearch(SeqList*R,int n,KeyType K)
int bin_search(int arr[],const int &n,const int &k)
{
//在有序表R[0..n-1]中进行二分查找,成功时返回结点的位置,失败时返回-1
int low=0,high=n-1,mid;//置当前查找区间上、下界的初值
while(low<=high)
{
mid=low+((high-low)/2);
//使用(low+high)/2会有整数溢出的问题
//(问题会出现在当low+high的结果大于表达式结果类型所能表示的最大值时,
//这样,产生溢出后再/2是不会产生正确结果的,而low+((high-low)/2)不存在这个问题
if(arr[mid] == k)
return mid;//查找成功返回
if(arr[mid] > k)
high=mid-1;//继续在R[low..mid-1]中查找
else
low=mid+1;//继续在R[mid+1..high]中查找
}
if(low>high)
return low;//当low>high时表示所查找区间内没有结果,查找失败
}
int bin_search2(int arr[],const int &n,const int &k)
{
if(arr && n>0)// 对数组和数组长度进行了判断,非常好的改进,上面也应该由此个改进的
{
int low,mid,high;
float rate;
low = 0;
high = n-1;
while(low<=high)
{
rate = (k-arr[low])*1.0/(arr[high]-arr[low]);
if(rate>1 || rate <0)// 对应数组的边界,不在数组中
return -1;
mid = low+(high-low)*rate;
if(arr[mid] == k)
return mid;
else if(arr[mid] > k)
high = mid - 1;
else
low = mid + 1;
}
return -1;// 查找失败
}
return -1;
}
分析:既然binary_search()函数返回的是middle 或 未找到时的指针,那就看middle指针变化的位置——mid = low + (high - low)/2;的条件是最外层while的条件(没有条件限制)第一次相等时就返回middle或者查找失败;因此 该函数是返回一个非递减序列[first, last)中的第一次找到的等于值val的位置,有可能是第一个中间的最后一个(如有序序列 ……3 3 3 ……)。
三:比较说明
(1)分析:既然lower_bound()函数返回的是first的指针,那就看first指针变化的位置——first = middle + 1;的条件是 严格的小于(<),即大于等于(>=)时first的位置是不变的。或者可以这样理解,当middle指向的值 < key时,first 指向middle的下一个位置,而下一个位置可能等于或者大于key;因此 该函数是返回一个非递减序列[first, last)中的第一个大于等于值val的位置。
(2)分析:既然upper_bound()函数返回的是first的指针,那就看first指针变化的位置——first = middle + 1;的条件是 小于等于(<=),即只有严格大于(>)时first的位置是不变的。或者可以这样理解,当middle指向的值 <= key时,first 指向middle的下一个位置,而下一个位置可能等于或者大于key,但是要是仍是等于的话,first会继续指向middle的下一位置,直到大于为止;因此 该函数是返回一个非递减序列[first, last)中的第一个大于值val的位置。
(3)分析:既然binary_search()函数返回的是middle 或 未找到时的指针,那就看middle指针变化的位置——mid = low + (high - low)/2;的条件是最外层while的条件(没有条件限制)第一次相等时就返回middle,或者查找失败;因此 该函数是返回一个非递减序列[first, last)中的第一次找到的等于值val的位置,有可能是第一个中间的最后一个(如有序序列 ……3 3 3 ……)。
(4)更正说明(优化说明)形式统一,格式统一,这样比较好的,就是不要太脱离你比较熟悉的binary_search()函数的形式,或者核心思想
#include <iostream>(5)运行结果(当时没有upper()函数的)
using namespace std;
//这个算法中,first是最终要返回的位置
int lower_bound(int *array, int size, int key)
{
int first = 0, middle;
int half, len;
len = size;
while(len > 0) {
half = len >> 1;
middle = first + half;
if(array[middle] < key) {
first = middle + 1;
len = len-half-1; //在右边子序列中查找
}
else
len = half; //在左边子序列(包含middle)中查找
}
return first;
}
int upper_bound(int *array, int size, int key)
{
int first = 0;
int len = size-1;
int half, middle;
while(len > 0){
half = len >> 1;
middle = first + half;
if(array[middle] > key) //中位数大于key,在包含last的左半边序列中查找。
len = half;
else{
first = middle + 1; //中位数小于等于key,在右半边序列中查找。
len = len - half - 1;
}
}
return first;
}
int lower(int *,int,int);
int upper(int*,int,int);
int main()
{
int key;
cout << "Hello world!" << endl;
int arr[] = {1,2,2,3,4,4,4,4,5,6,7,9,9,10};
while(cin >> key)
{
cout << "1ziji: " << lower(arr,sizeof(arr)/sizeof(arr[0]),key) << endl;
cout << "2: " << lower_bound(arr,sizeof(arr)/sizeof(arr[0]),key) << endl;
cout << "3: " << upper_bound(arr,sizeof(arr)/sizeof(arr[0]),key) << endl;
cout << "4ziji: " << upper_bound(arr,sizeof(arr)/sizeof(arr[0]),key) << endl;
}
return 0;
}
int lower(int arr[],int size,int key)
{
if(arr && size > 0)
{
int l = 0;
int r = size -1;
int mid;
while(l < r)
{
mid = l + (r-l)/2;
if(*(arr+mid)>=key)
r = mid;
else
l = mid + 1;
}
return l;
}
return -1;
}
int upper(int*arr,int size,int key)
{
if(arr && size > 0)
{
int l = 0;
int r = size -1;
int mid;
while(l < r)
{
mid = l + (r-l)/2;
if(arr[mid]>key)
r = mid;
else
l = mid+1;
}
return l;
}
return -1;
}
加上upper()函数之后的效果如下: