二分查找算法在C/C++法式中的運用示例。本站提示廣大學習愛好者:(二分查找算法在C/C++法式中的運用示例)文章只能為提供參考,不一定能成為您想要的結果。以下是二分查找算法在C/C++法式中的運用示例正文
二分查找算法的思惟很簡略,《編程珠玑》中的描寫: 在一個包括t的數組內,二分查找經由過程對規模的跟綜來處理成績。開端時,規模就是全部數組。經由過程將規模中央的元素與t比擬並拋棄一半規模,規模就被減少。這個進程一向連續,直到在t被發明,或許誰人可以或許包括t的規模已成為空。
Donald Knuth在他的《Sorting and Searching》一書中指出,雖然第一個二分查找算法早在1946年就被揭橥,但第一個沒有bug的二分查找算法倒是在12年後才被揭橥出來。個中罕見的一個bug是對中央值下標的盤算,假如寫成(low+high)/2,當low+high很年夜時能夠會溢出,從而招致數組拜訪失足。改良的辦法是將盤算方法寫成以下情勢:low+ ( (high-low) >>1)便可。上面給出修正後的算法代碼:
int binarysearch1(int a[],int n,int x) { int l,u,m; l=0;u=n; while(l<u) { m=l+((u-l)>>1); if(x<a[m]) u=m; else if(x==a[m]) return m; else l=m+1; } return -1; }
這裡留意一點,因為應用的是纰謬稱區間,所以下標的調劑看上去有點不規整。一個是u=m,另外一個是l=m+1。其實很好懂得,調劑前區間的情勢應當是 [ )的情勢,假如中央值比查找值小,那末調劑的是右邊界,也就是閉的部門,所以加1;不然,調劑是左邊界,是開的部門,所以不消減1。調劑後還是[ )的情勢。固然也能夠寫成對稱的情勢。代碼以下:
int binarysearch1(int a[],int n,int x) { int l,u,m; l=0;u=n-1; while(l<=u) { m=l+((u-l)>>1); if(x<a[m]) u=m-1; else if(x==a[m]) return m; else l=m+1; } return -1; }
如許也看上去比擬規整,然則有個缺乏。假如想把法式改成“純指針”的情勢,就會有費事。修正成純指針的代碼以下:
int binarysearch2(int *a,int n,int x) { int *l,*u,*m; l=a;u=a+n-1; while(l<=u) { m=l+((u-l)>>1); if(x<*m) u=m-1; else if(x==*m) return m-a; else l=m+1; } return -1; }
當n為0時,會援用有效地址。而用非對稱區間則不會有這個成績。代碼以下:
int binarysearch2(int *a,int n,int x) { int *l,*u,*m; l=a;u=a+n; while(l<u) { m=l+((u-l)>>1); if(x<*m) u=m; else if(x==*m) return m-a; else l=m+1; } return -1; }
下面給出的二分查找是迭代法完成,固然也能夠用遞歸的方法完成。代碼以下:
int binarysearch3(int a[],int l,int u,int x) int m=l+((u-l)>>1); if(l<=u) { if(x<a[m]) return binarysearch3(a,l,m-1,x); else if(x==a[m]) return m; else return binarysearch3(a,m+1,u,x); } return -1;
上述這些二分算法,若數組元素反復,前往的是反復元素的某一個元素。假如願望前往被查找元素第一次湧現的地位,則須要修正代碼。上面給出了一種解法:
int binarysearch4(int a[],int n,int x) { int l,u,m; int flag=-1; l=0;u=n; while(l<u) { m=l+((u-l)>>1); if(x<a[m]) u=m; else if(x==a[m]) flag=u=m; else l=m+1; } return flag; }
上面是《編程珠玑》上的解法:
int binarysearch4(int a[],int n,int x) { int l,u,m; l=-1;u=n; while(l+1<u) { m=l+((u-l)>>1); if(a[m]<x) l=m; else u=m; } return (u>=n||a[u]!=x)?-1:u; }
至此二分算法的代碼評論辯論停止,上面評論辯論一下法式的測試成績。《代碼之美》有一章專門引見二分查找算法的測試,異常英俊。這裡布鼓雷門,簡略給出幾個測試用例。針對binarysearch1。測試法式以下:
#include <iostream> #include <cassert> #include <algorithm> #include <ctime> using namespace std; int calmid(int l,int u) { return l+((u-l)>>1); } int binarysearch1(int a[],int n,int x); #define bs1 binarysearch1 int main() { long start,end; start=clock(); int a[9]={-2147483648,-13,-10,-5,-3,0,1,400,2147483647}; //中值下標盤算的測試 assert(calmid(0,1)==0); assert(calmid(0,2)==1); assert(calmid(1000000,2000000)==1500000); assert(calmid(2147483646,2147483647)==2147483646); assert(calmid(2147483645,2147483647)==2147483646); //冒煙測試 assert(bs1(a,9,0)==5); assert(bs1(a,9,1)==6); assert(bs1(a,9,2)==-1); //界限測試 assert(bs1(a,0,1)==-1); //0個元素 assert(bs1(a,1,-2147483648)==0); //1個元素 勝利 assert(bs1(a,1,-2147483647)==-1); //1個元素 掉敗 assert(bs1(a,9,-2147483648)==0); //首個元素 assert(bs1(a,9,-3)==4); //中央元素 assert(bs1(a,9,2147483647)==8); //末尾元素 //主動化測試 int b[10000]; int i,j; for(i=0;i<10000;i++) { b[i]=i*10; for(j=0;j<=i;j++) { assert(bs1(b,i+1,j*10)==j); assert(bs1(b,i+1,j*10-5)==-1); } } //主動化測試 引入隨機數 srand(time(0)); for(i=0;i<10000;i++) { b[i]=rand()%1000000; sort(&b[0],&b[i]); for(j=0;j<=i;j++) { int x=rand(); int k=bs1(b,i+1,x); if(k!=-1) assert(b[k]==x); } } end=clock(); cout<<(end-start)/1000.0<<'s'<<endl; return 0; }
留意到數組的元素有負數,正數,零,最年夜值,最小值。平日會忘失落正數的測試,引入最年夜值和最小值,重要是為了界限測試。
第一,測試了中值下標的盤算。別的寫了一個小函數,零丁測試。斟酌到內存能夠放不下這麼年夜的數組,是以只是模仿測試,並沒有真正請求這麼年夜的空間,然則關於中值下標的測試足夠了。
第二,冒煙測試。即做一些最根本的測試。測試經由過程落後行界限測試。
第三,界限測試。這裡有三品種型,一是針對數組元素個數,分離是0個,1個。二是針對元素地位,分離是首個元素,中央元素,末尾元素。三是針對元素值,有最年夜值,最小值,0等測試。
第四,主動化測試。這裡主動生成測試的數組,然後針對每一個元素停止勝利查找測試。
第五,主動化測試,只不外數組的元素是隨機值。
第五,機能測試。這裡相干代碼沒有列出。以上測試都經由過程時,可以修正查找算法,添加機能測試的代碼。其實可以簡略添加一個比擬的計數器。前往值從本來的查找成果改成比擬的計數器值便可。代碼比擬簡略,就不列了。
Note:二分查找輕易疏忽的一個bug
關於二分查找算法,信任年夜家確定不會生疏。算法從一個排好序的數組中找指定的元素,假如找到了前往該元素在數組中的索引,不然前往-1。上面給出懂得法。
//a為排好序的數組,n為數組的年夜小,x為指定元素 int binarySearch(int a[], int n, int x) { int left = 0, right = n-1, middle = 0; int tmp = 0; while(left <= right) { middle = (left + right)/2; tmp = a[middle]; if(x < tmp) right = middle - 1; else if(x > tmp) left = middle + 1; else return middle; } return -1; }
乍看沒有毛病,然則不幸的是,該法式存在一個bug。當數組極年夜時,(left+right)能夠為正數,則數組下標溢出,法式瓦解。
處理的計劃:將middle=(left+right)/2改成middle=left+(right-left)/2便可。即應用減法取代加法,從而清除上溢。
參考自《代碼之美》