首先給出一個對字符串比較好的散列函數,在有些地方把該算法稱為“均勻哈希算法”。
//提供一個對string進行散列的函數
int hashString(const string &str){
string s;
if(str.size()>1024) //如果str太長,則只取前1024個字符
s=str.substr(0,1024);
else
s=str;
int rect=0;
for(int i=0;i<s.size();++i)
rect=rect*37+s[i]; //Horner法則
return rect;
}
這個散列函數可能會溢出,導致返回值為負數。注意這里選擇的數字是37,好在哪里我也無法證明,至少可以看到素數在散列函數中十分有用。
再好的散列函數也會產生沖突(collision),而且沖突還時有發生,所以一個HashTable必須考慮沖突如何解決。方法有:
- 分離鏈接法(separate chaining),后面會詳細介紹。
- 線性探測法。就用采用散列函數h0發現有沖突時就采用h1,如果還有沖突就采用h2……其中hi(x)=[ h(x) + f(i) ] mod tableSize。 f(i)為i的線性函數,通常取f(i)=i。
- 平方探測法。與線性探測法法類似, f(i)為i的二次函數,通常取f(i)=i^2。
- 雙散列。f(i)=i·h2(x)。
分離鏈接法就是把沖突的元素存放在一個鏈表里。當然除鏈表外也可采用其他數據結構,如二叉查找樹、甚至是另外一個散列表,它們的查找速度都比鏈表要快。但是我們期望的是散列函數足夠地好、槽足夠地多,所以對應的鏈表都應該很短,不值得去嘗試更復雜的結構。
槽數為素數時,HashTable的性能會比較好,這就引用一個問題:如何確定一個數是素數?一種快速的判別方法起源於費馬小定理。
如果對於任意滿足1 < b < p的b下式都成立:
//模c下求a的b次方
int powermod(int a,int b,int c){
if(b==0)
return 1;
if(b==1)
return a%c;
int t=powermod(a,b/2,c);
if(b&1){ //如果b是奇數
return t*t*a%c;
}
else{
return t*t%c;
}
}
//判斷數p是否為素數
bool isPrime(int p){
for(int i=0;i<100;i++){;
if(powermod(rand()%(p-1)+1,p-1,p)!=1)
return false;
}
return true;
}
假設一個散列表能容納n個元素、具有m個槽,定義其裝載因子(load factor)α為n/m。假定散列函數足夠地好,任何一個元素散列到m個槽位的可能性是相同的,且與其他元素已被散列到什么位置上獨立無關(這個假設稱為簡單一致散列simple uniform hashing)。則平均情況下查找一個元素是否在散列表中的時間復雜度是O(1+α)。所以裝載因子α成為我們關注的焦點。
當一個HashTable太滿后,發生沖突的概率就會大大增加。我們的策略是:當達到事先設定的裝載因子時,就把槽位擴展成原先的2倍以上(取最小的素數)。這叫再散列。原先的HashTable完全釋放,申請新的更大的空間,然后把已有的元素重新散列到新的HashTable。再散列開銷很大,但我們期望的是發生再散列的次數很少。
template<typename HashObj>
void HashTable<HashObj>::rehash(){
int oldSlot=getSlots();
int newSlot=oldSlot*2;
while(!isPrime((unsigned long)newSlot)){
newSlot++;
}
vector<list<HashObj> > oldVector = table;
for(int i=0;i<table.size();++i)
table[i].clear();
table.resize(newSlot);
setSlots(newSlot);
capacity=0;
for(int i=0;i<oldVector.size();++i){
typename list<HashObj>::iterator itr=oldVector[i].begin();
while(itr!=oldVector[i].end())
insert(*itr++);
}
}
有一個語法點:注意17行,當使用含有模板類型的迭代器時,要在前面加一個typename。
再散列需要把之前散列過的key重新做一次散列,這稱為非一致性哈希。一致性哈希要求提供一個hashtable,它能在節點加入離開時不會導致映射關系的重大變化。在某些實際應用中需要使用一致性哈希。
下面給出完整代碼:
hash.h
#ifndef _HASH_H
#define _HASH_H
#include<vector>
#include<list>
#include<string>
using namespace std;
//HashObj類型的數據必須提供函數:int hash();operator ==
template<typename HashObj>
class HashTable{
private:
int capacity; //已容納的元素的個數
int slots; //槽的個數(取素數散列性比較好)
double alpha; //裝載因子
vector<list<HashObj> > table; //哈希表
int myhash(const HashObj& ele); //內部哈希函數,負責把數據映射到[0,slots-1]
void rehash(); //當達到裝載因子時,槽數擴大為原來的2倍以上(取最小的素數),重新進行再散列
public:
//判斷一個數據是否在HashTable中
bool contain(const HashObj &ele);
//插入一個元素到HashTable。插入成功則返回true,如果元素原先就存在於HashTable中,則返回false
bool insert(const HashObj& ele);
//從HashTable中刪除一個元素。如果指定元素不在HashTable中,則刪除失敗,返回false
bool remove(const HashObj &ele);
HashTable(int slots=10007,double alpha=0.7):slots(slots),alpha(alpha){
capacity=0;
table.resize(slots);
}
int getSlots(){
return slots;
}
void setSlots(int num){
slots=num;
}
};
//提供一個對string進行散列的函數
int hashString(const string &str){
string s;
if(str.size()>1024) //如果str太長,則只取前1024個字符
s=str.substr(0,1024);
else
s=str;
int rect=0;
for(int i=0;i<s.size();++i)
rect=rect*37+s[i]; //Horner法則
return rect;
}
//由於散列的槽數最好是素數,所以提供一個判別素數的函數
int fmod(int a, int b, int c)//快速模取冪
{
if(b == 1) return a;
int t = fmod(a, b / 2, c);
t = (t * t) % c;
if(b & 1) t = (t * a) % c;
return t;
}
bool isPrime(int n)//米勒-拉賓算法
{
for(int i = 0; i < 100; ++ i)
{
if(fmod(rand() % (n - 1) + 1, n - 1, n) != 1)//a的取值為[1,n-1],a的值需要變化,所以用到隨機函數
return false;
}
return true;
}
//提供一個滿足HashObj要求的類
class Employee{
private:
string name;
int id;
public:
Employee(string name="",int id=-1):name(name),id(id){}
bool operator == (const Employee & obj) const{
return this->name==obj.name;
}
int hash()const{
return hashString(name);
}
void setName(string str){
name=str;
}
void setId(int i){
id=i;
}
};
#endif
hash.cpp
#include<iostream>
#include<cassert>
#include<algorithm>
#include"hash.h"
template<typename HashObj>
int HashTable<HashObj>::myhash(const HashObj & ele){
int rect=ele.hash();
rect%=slots;
if(rect<0)
rect+=slots;
return rect;
}
template<typename HashObj>
bool HashTable<HashObj>::contain(const HashObj & ele){
int index=myhash(ele);
const list<HashObj> & whichList=table[index];
return find(whichList.begin(),whichList.end(),ele)!=whichList.end();
}
template<typename HashObj>
bool HashTable<HashObj>::insert(const HashObj & ele){
int index=myhash(ele);
list<HashObj> & whichList=table[index];
if(find(whichList.begin(),whichList.end(),ele)!=whichList.end())
return false;
whichList.push_back(ele);
capacity++;
if(capacity*1.0/slots>alpha)
rehash();
return true;
}
template<typename HashObj>
bool HashTable<HashObj>::remove(const HashObj & ele){
int index=myhash(ele);
list<HashObj> & whichList=table[index];
typename list<HashObj>::iterator itr=find(whichList.begin(),whichList.end(),ele);
if(itr==whichList.end())
return false;
whichList.erase(itr);
capacity--;
return true;
}
template<typename HashObj>
void HashTable<HashObj>::rehash(){
int oldSlot=getSlots();
int newSlot=oldSlot*2;
while(!isPrime((unsigned long)newSlot)){
newSlot++;
}
vector<list<HashObj> > oldVector = table;
for(int i=0;i<table.size();++i)
table[i].clear();
table.resize(newSlot);
setSlots(newSlot);
capacity=0;
for(int i=0;i<oldVector.size();++i){
typename list<HashObj>::iterator itr=oldVector[i].begin();
while(itr!=oldVector[i].end())
insert(*itr++);
}
}
int main(){
const int arrSize=9;
HashTable<Employee> hashTable(7,1.0); //剛開始設槽數為7
string names[arrSize]={"hujintao","jiangzeming","heizeming",
"chaogai","jingchengwu","liangchaowei",
"weijiabao","zhamgsanbao","zengxiaoxian"};
Employee employee[arrSize];
for(int i=0;i<arrSize;++i){
employee[i].setName(names[i]);
employee[i].setId(i+1);
}
for(int i=0;i<arrSize;++i){
hashTable.insert(employee[i]);
}
assert(hashTable.getSlots()==17); //擴容后槽數應該為17
for(int i=0;i<arrSize;++i){
assert(hashTable.contain(employee[i]));
hashTable.remove(employee[i]);
assert(!hashTable.contain(employee[i]));
}
return 0;
}
