我们现在来看一下斐波拉契查找。
1.斐波拉契数列介绍
2.斐波拉契查找原理
3.代码实现
4.与二分查找的渊源
5.在向量模板中实现
斐波拉契数列(Fibonacci sequence)是指这样一个数列 fib={0,1,1,2,3,5,8,13,21,34·······}。其在数学上的定义为Fib[0]=0,Fib[1]=1,Fib[n]=Fib[n-1]+Fib[n-2]。
设V为被查找的数列
假设函数fib(k)是第k个斐波拉契数,令n=fib(k)-1,若要被查找的有序数列的元素个数正好等于n,e为被查找的数。则可以将查找区间分成如下三部分:
V[0]<=e<V[fib(k-1)-1]
e=V[fib(k-1)-1]
V[fib(k-1)-1)<e<V(fib(k)-1)
可用下图形象表示:
可知前一段被分割的长度为fib(k-1)-1,后一段长度为fib(k-2)-1。如果e属于左半段,我们可以深入左半段将fib(k-1)-1再次按如此方式分割,e属于右半段同理。若e等于V[fib(k-1)-1],则命中。
#include<iostream> using namespace std; class Fib{ private: int fibNumber;//获取斐波拉契数 int Rank;//斐波拉契数下标 int FibSquence[100];//用于存储斐波拉契数列 public: Fib(int n) { fibNumber=n+1; FibSquence[1]=1; FibSquence[2]=1; for(int i=3;i<100;i++)//生成斐波拉契数列 { FibSquence[i]=FibSquence[i-1]+FibSquence[i-2]; } for(int i=1;i<100;i++) { if(FibSquence[i]==fibNumber) { Rank=i; break; } } } ~Fib(){ } int get()//获取当前斐波拉契数 { return FibSquence[Rank]; } void prev()//前一项斐波拉契数 { Rank--; } }; int fibSearch(int *A,int e,int lo,int hi) { Fib fib(hi-lo); while(lo<hi) { while(lo<hi) { while(hi-lo<fib.get())//通过向前顺序查找,确定形如Fib(k)-1的轴点 { fib.prev(); } int mi=lo+fib.get()-1;//按黄金比例切分 if(e<A[mi]) { hi=mi;//深入前半段[lo,mi)查找 } else if(A[mi]<e) { lo=mi+1;//深入后半段(mi,hi)查找 } else { return mi;//在mi处命中 } } } return -1; //查找失败 } int main() { /***测试***/ int a[7]; for(int i=0;i<7;i++) { a[i]=i;//生成数列{0,1,2,3,4,5,6}。 } cout<<"Search(3)="<<fibSearch(a,3,0,7)<<endl; /********/ return 0; }我们不难理解,上面我写了一个Fib类,它的作用是获取数列的长度n,并获取n+1这个斐波拉契数,它的prev()函数可以让它指向前一个斐波拉契数,get()可返回此时的斐波拉契数。
由上述介绍我们不难看出,其实斐波拉契查找就是一个特殊的二分查找。只是分的不是1/2而已。它按斐波拉契数列分割,前段长后段短。这样就弥补了一点二分查找向前和向后查找的代价不一样的问题(二分查找向前查找需判断一次,向后查找需判断两次)。而又有数学式子可以证明,用按斐波拉契数列分割时,查找的效率最高,因此它又被称为黄金分割。但是由于它又有条件的限制(被查找的数列必须满足其长度n为一个斐波拉契数fib(k)减1,当n不足时需补齐),比较麻烦,所以大家还是少用为好,其实它也快不了多少,特别是在如今硬件计算速度如此快的条件下,可以忽略不计。
由于本页为向量这一页的分支,所以最终这个查找还是会与向量对接。
代码实现:
template<typename T> Rank myVector<T>:: fibSrarch(T*A,T const &e,Rank lo,Rank hi) { Fib fib(hi-lo); while(lo<hi) { while(hi-lo<fib.get()) { fib.prev();//返回前一项斐波拉契数值 } Rank mi=lo+fib.get()-1;//按黄金比例划分 if(e<A[mi]) { hi=mi;//深入前半段[lo,mi) } else if(A[mi]<e) { lo=mi+1;//深入后半段(mi,hi) } else { return mi;//在mi处命中 } } return -1;//查找失败 }
想要学习更多关于向量的知识吗?请点击!!