编程之美-程序改错及扩展问题

原文的大意是这样的,要求写一个二分查找算法,并且当要查找的数出现不止一次时,返回最后那个数的下标。然后给你一段有错的代码,请你改。

这段有错的代码是这样的(我作了简化):

/* b和e是查找区间的两端,v是要找的数 */  
int bsearch(int a[], int b, int e, int v)  
{  
    int l = b;  
    int u = e;  
    while (l < u) {  
        int m = (l + u) / 2;  
        if (a[m] <= v)   l = m;  
        else        u = m – 1;  
    }  
    return a[u]==v? u:-1;  
} 

文中指出了代码中的错误,一是第7行(l+u)可能溢出,二是当程序运行到l=u-1时,计算出的m就与l相等,如果这时不幸地第8行的条件成立,那么就会进入死循环,因为这次迭代没有改变l和u的值,无论再迭代多少遍也不会了。

        这两点分析都很正确,但接下来作者给出的修改版本,不仅完全不是在上面代码的基础上改,而且给出的代码逻辑复杂、混乱,缺乏一致性。

        下面我将分析上面这段有错代码,给出自己的修改,然后再分析作者给出的版本。

        上面代码的思想是清晰的,维护一个循环不变式:如果要找的数存在的话,它下标一定在区间[ l , u ]中,然后通过迭代,不断缩小区间,直到退化成一个点,最后检测这个点是不是目标值。循环体中对不变式的维护都没有错。问题就是怎样避免第8行的把m直接赋给l。题目要求返回最后一个满足的数,所以在循环体中,当检测到相等时,一不能立刻就返回,二不能缩减区间的后半段。要维护这样一个不变式,这个赋值在所难免。我能想到的唯一的办法就是,当l=u-1时,就退出循环,然后先后检查此时的a[u]和a[l]。

        所以我的修改是这样的:

int bsearch(int a[], int b, int e, int v)  
{  
        int l = b;  
    int u = e;  
    while (l < u -1) {  
        int m = l + (u – l) * 0.5;//原来的/2并没错,但*0.5可以减少CPU时钟周期  
        if (a[m] <= v)   l = m;  
        else        u = m – 1;  
    }  
    if (a[u] == v)          return u;  
    else if (a[l] == v) return l;  
    else return -1;

接下来,我要分析一下作者的修改,他的代码大意是这样的:

int bsearch(int a[], int b, int e, int v)  
{  
    int l = b;  
    int u = e;  
    //循环结束有两种情况:   
    //若l为偶数则l==u   
    //否则,l==u-1   
    while (l < u -1) {  
        int m = l + (u – l)/ 2;  
        if (a[m] <= v)   l = m;  
        else {  
            //不需要m-1,防止l==u  为什么呢?不明白 
            u = m;  
        }  
    }  
    if (a[u] == v)  return u;   //先判断序号最大的值   
    else if(a[l] == v)  return l;  
    else    return -1;  
}  

我们先无视注释,我完全搞不懂注释,文中也不解释清楚。从循环体的if-else语句判断,他要维护的不变式是考察区间[l, u),注意是“左闭右开”,要找的数下标一定在这个区间内,如果存在的话。而看他对区间端点的赋初值,u=e,把a[e]给排除在外了。那要找的数就是a[e]怎么办呢?他有这条语句把关呢,if (a[u] == v) return u;后面还十分堂皇地加了注释。其实如果初值正确的话,只需要检测a[l]即可,因为a[u]这个值根本不在我们要考察的区间里。

        所以我对作者的版本的改法是:

int bsearch(int a[], int b, int e, int v)  
{  
    int l = b;  
    int u = e+1;  //有改动   
    while (l < u -1) {  
        int m = l + (u – l)* 0.5;  
        if (a[m] <= v)   l = m;  
        else    u = m;  
    }  
    return a[l]==v? l:-1;  //有改动   
}  

 

 

 

 

    原文作者:zhuhuiby
    原文地址: https://blog.csdn.net/zhuhuiby/article/details/6742909
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞