heap的分析

以前对堆都没什么了解，今天看了书之后才终于明白了，堆其实就是一个完全二叉树，对于查找和排序时间复杂度可以很好的降低：

heap分为两种：一种是max-heap一种是min-heap，按照书上的，我们假设为 max-heap，即根节点为最大值：

所以：

我们的内存分配方式为线性存储，用的是vector。

父节点如果位置是 i，那么子节点位置就是2*i+1 和2*i+2；

所以如果我们插入一个元素的话，通过跟它的父节点比较，然后依次互换位置。这样就能实现插入操作：

inline void push_heap(Iterator first,Iterator last,T value) { size_t distance=last-first-1; value=*(last-1); size_t parent=(distance-1)/2; while(distance>0&&*(first+parent)<value) { *(first+distance)=*(first+parent); distance=parent; parent=(distance-1)/2; } *(first+distance)=value; } 

如上面，我们要在first 和 last之间插入元素，当然这个时候要插入的元素已经在 last-1 这个位置上。那么我们要做的就是比较它跟父节点的大小，如果大于父节点 那么就互换位置，然后继续循环比较，如果小于父节点，那么就返回，因为这个时候我们已经完成了插入操作，符合完全二叉树，父节点值大于子节点。

完整代码：

template<typename Iterator> inline void push_heap(Iterator first,Iterator last) { typedef typename Iterator::value_type T; push_heap(first,last,T()); } template<typename Iterator,typename T> inline void push_heap(Iterator first,Iterator last,T value) { size_t distance=last-first-1; value=*(last-1); size_t parent=(distance-1)/2; while(distance>0&&*(first+parent)<value) { *(first+distance)=*(first+parent); distance=parent; parent=(distance-1)/2; } *(first+distance)=value; } 

对于删除操作：

那么我们首先把跟节点和最后一个元素值互换，然后从跟节点开始依次跟他的子节点比较，如果小于直接点那么跟子节点互换，如果大于子节点那么就保持不同返回：

template<typename Iterator ,typename T> inline void pop_heap(Iterator first,Iterator last,T value) { value=*(last-1);//得到尾部节点的数值 *(last-1)=*first;//把头节点的值赋给尾部节点 size_t parent=0;//得到父节点的位置 size_t distance=last-first-1;//得到最大距离 size_t r_child=2*parent+2; size_t child; while(r_child<distance&&value<max(*(first+r_child-1),*(first+r_child))) { child=*(first+r_child)>*(first+r_child-1)?r_child:(r_child-1); *(first+parent)=*(first+child); parent=child; r_child=2*parent+2; } if(distance==r_child) { if(value<*(first+r_child-1)) *(first+parent)=*(first+r_child-1); *(first+r_child-1)=value; return ; } *(first+parent)=value; return ; }