1,、問題描述
哈夫曼編碼是廣泛地用于數(shù)據(jù)文件壓縮的十分有效的編碼方法。其壓縮率通常在20%~90%之間,。哈夫曼編碼算法用字符在文件中出現(xiàn)的頻率表來建立一個用0,,1串表示各字符的最優(yōu)表示方式。一個包含100,000個字符的文件,,各字符出現(xiàn)頻率不同,,如下表所示。
有多種方式表示文件中的信息,,若用0,1碼表示字符的方法,,即每個字符用唯一的一個0,1串表示。若采用定長編碼表示,,則需要3位表示一個字符,,整個文件編碼需要300,000位;若采用變長編碼表示,,給頻率高的字符較短的編碼,;頻率低的字符較長的編碼,,達到整體編碼減少的目的,,則整個文件編碼需要(45×1+13×3+12×3+16×3+9×4+5×4)×1000=224,000位,由此可見,,變長碼比定長碼方案好,,總碼長減小約25%。
前綴碼:對每一個字符規(guī)定一個0,1串作為其代碼,,并要求任一字符的代碼都不是其他字符代碼的前綴,。這種編碼稱為前綴碼。編碼的前綴性質(zhì)可以使譯碼方法非常簡單,;例如001011101可以唯一的分解為0,0,101,1101,,因而其譯碼為aabe。
譯碼過程需要方便的取出編碼的前綴,,因此需要表示前綴碼的合適的數(shù)據(jù)結(jié)構(gòu),。為此,可以用二叉樹作為前綴碼的數(shù)據(jù)結(jié)構(gòu):樹葉表示給定字符,;從樹根到樹葉的路徑當(dāng)作該字符的前綴碼,;代碼中每一位的0或1分別作為指示某節(jié)點到左兒子或右兒子的“路標(biāo)”。
從上圖可以看出,,表示最優(yōu)前綴碼的二叉樹總是一棵完全二叉樹,,即樹中任意節(jié)點都有2個兒子。圖a表示定長編碼方案不是最優(yōu)的,,其編碼的二叉樹不是一棵完全二叉樹,。在一般情況下,若C是編碼字符集,表示其最優(yōu)前綴碼的二叉樹中恰有|C|個葉子,。每個葉子對應(yīng)于字符集中的一個字符,,該二叉樹有|C|-1個內(nèi)部節(jié)點。
給定編碼字符集C及頻率分布f,即C中任一字符c以頻率f(c)在數(shù)據(jù)文件中出現(xiàn),。C的一個前綴碼編碼方案對應(yīng)于一棵二叉樹T,。字符c在樹T中的深度記為dT(c)。dT(c)也是字符c的前綴碼長,。則平均碼長定義為:使平均碼長達到最小的前綴碼編碼方案稱為C的最優(yōu)前綴碼,。
2、構(gòu)造哈弗曼編碼
哈夫曼提出構(gòu)造最優(yōu)前綴碼的貪心算法,,由此產(chǎn)生的編碼方案稱為哈夫曼編碼,。其構(gòu)造步驟如下:
(1)哈夫曼算法以自底向上的方式構(gòu)造表示最優(yōu)前綴碼的二叉樹T。
(2)算法以|C|個葉結(jié)點開始,,執(zhí)行|C|-1次的“合并”運算后產(chǎn)生最終所要求的樹T,。
(3)假設(shè)編碼字符集中每一字符c的頻率是f(c)。以f為鍵值的優(yōu)先隊列Q用在貪心選擇時有效地確定算法當(dāng)前要合并的2棵具有最小頻率的樹,。一旦2棵具有最小頻率的樹合并后,,產(chǎn)生一棵新的樹,其頻率為合并的2棵樹的頻率之和,,并將新樹插入優(yōu)先隊列Q,。經(jīng)過n-1次的合并后,優(yōu)先隊列中只剩下一棵樹,,即所要求的樹T,。
構(gòu)造過程如圖所示:
具體代碼實現(xiàn)如下:
(1)4d4.cpp,程序主文件
- //4d4 貪心算法 哈夫曼算法
- #include "stdafx.h"
- #include "BinaryTree.h"
- #include "MinHeap.h"
- #include <iostream>
- using namespace std;
-
- const int N = 6;
-
- template<class Type> class Huffman;
-
- template<class Type>
- BinaryTree<int> HuffmanTree(Type f[],int n);
-
- template<class Type>
- class Huffman
- {
- friend BinaryTree<int> HuffmanTree(Type[],int);
- public:
- operator Type() const
- {
- return weight;
- }
- //private:
- BinaryTree<int> tree;
- Type weight;
- };
-
- int main()
- {
- char c[] = {'0','a','b','c','d','e','f'};
- int f[] = {0,45,13,12,16,9,5};//下標(biāo)從1開始
- BinaryTree<int> t = HuffmanTree(f,N);
-
- cout<<"各字符出現(xiàn)的對應(yīng)頻率分別為:"<<endl;
- for(int i=1; i<=N; i++)
- {
- cout<<c[i]<<":"<<f[i]<<" ";
- }
- cout<<endl;
-
- cout<<"生成二叉樹的前序遍歷結(jié)果為:"<<endl;
- t.Pre_Order();
- cout<<endl;
-
- cout<<"生成二叉樹的中序遍歷結(jié)果為:"<<endl;
- t.In_Order();
- cout<<endl;
-
- t.DestroyTree();
- return 0;
- }
-
- template<class Type>
- BinaryTree<int> HuffmanTree(Type f[],int n)
- {
- //生成單節(jié)點樹
- Huffman<Type> *w = new Huffman<Type>[n+1];
- BinaryTree<int> z,zero;
-
- for(int i=1; i<=n; i++)
- {
- z.MakeTree(i,zero,zero);
- w[i].weight = f[i];
- w[i].tree = z;
- }
-
- //建優(yōu)先隊列
- MinHeap<Huffman<Type>> Q(n);
- for(int i=1; i<=n; i++) Q.Insert(w[i]);
-
- //反復(fù)合并最小頻率樹
- Huffman<Type> x,y;
- for(int i=1; i<n; i++)
- {
- x = Q.RemoveMin();
- y = Q.RemoveMin();
- z.MakeTree(0,x.tree,y.tree);
- x.weight += y.weight;
- x.tree = z;
- Q.Insert(x);
- }
-
- x = Q.RemoveMin();
-
- delete[] w;
-
- return x.tree;
- }
(2)BinaryTree.h 二叉樹實現(xiàn)
- #include<iostream>
- using namespace std;
-
- template<class T>
- struct BTNode
- {
- T data;
- BTNode<T> *lChild,*rChild;
-
- BTNode()
- {
- lChild=rChild=NULL;
- }
-
- BTNode(const T &val,BTNode<T> *Childl=NULL,BTNode<T> *Childr=NULL)
- {
- data=val;
- lChild=Childl;
- rChild=Childr;
- }
-
- BTNode<T>* CopyTree()
- {
- BTNode<T> *nl,*nr,*nn;
-
- if(&data==NULL)
- return NULL;
-
- nl=lChild->CopyTree();
- nr=rChild->CopyTree();
-
- nn=new BTNode<T>(data,nl,nr);
- return nn;
- }
- };
-
-
- template<class T>
- class BinaryTree
- {
- public:
- BTNode<T> *root;
- BinaryTree();
- ~BinaryTree();
-
- void Pre_Order();
- void In_Order();
- void Post_Order();
-
- int TreeHeight()const;
- int TreeNodeCount()const;
-
- void DestroyTree();
- void MakeTree(T pData,BinaryTree<T> leftTree,BinaryTree<T> rightTree);
- void Change(BTNode<T> *r);
-
- private:
- void Destroy(BTNode<T> *&r);
- void PreOrder(BTNode<T> *r);
- void InOrder(BTNode<T> *r);
- void PostOrder(BTNode<T> *r);
-
- int Height(const BTNode<T> *r)const;
- int NodeCount(const BTNode<T> *r)const;
- };
-
- template<class T>
- BinaryTree<T>::BinaryTree()
- {
- root=NULL;
- }
-
- template<class T>
- BinaryTree<T>::~BinaryTree()
- {
-
- }
-
- template<class T>
- void BinaryTree<T>::Pre_Order()
- {
- PreOrder(root);
- }
-
- template<class T>
- void BinaryTree<T>::In_Order()
- {
- InOrder(root);
- }
-
- template<class T>
- void BinaryTree<T>::Post_Order()
- {
- PostOrder(root);
- }
-
- template<class T>
- int BinaryTree<T>::TreeHeight()const
- {
- return Height(root);
- }
-
- template<class T>
- int BinaryTree<T>::TreeNodeCount()const
- {
- return NodeCount(root);
- }
-
- template<class T>
- void BinaryTree<T>::DestroyTree()
- {
- Destroy(root);
- }
-
- template<class T>
- void BinaryTree<T>::PreOrder(BTNode<T> *r)
- {
- if(r!=NULL)
- {
- cout<<r->data<<' ';
- PreOrder(r->lChild);
- PreOrder(r->rChild);
- }
- }
-
- template<class T>
- void BinaryTree<T>::InOrder(BTNode<T> *r)
- {
- if(r!=NULL)
- {
- InOrder(r->lChild);
- cout<<r->data<<' ';
- InOrder(r->rChild);
- }
- }
-
- template<class T>
- void BinaryTree<T>::PostOrder(BTNode<T> *r)
- {
- if(r!=NULL)
- {
- PostOrder(r->lChild);
- PostOrder(r->rChild);
- cout<<r->data<<' ';
- }
- }
-
- template<class T>
- int BinaryTree<T>::NodeCount(const BTNode<T> *r)const
- {
- if(r==NULL)
- return 0;
- else
- return 1+NodeCount(r->lChild)+NodeCount(r->rChild);
- }
-
- template<class T>
- int BinaryTree<T>::Height(const BTNode<T> *r)const
- {
- if(r==NULL)
- return 0;
- else
- {
- int lh,rh;
- lh=Height(r->lChild);
- rh=Height(r->rChild);
- return 1+(lh>rh?lh:rh);
- }
- }
-
- template<class T>
- void BinaryTree<T>::Destroy(BTNode<T> *&r)
- {
- if(r!=NULL)
- {
- Destroy(r->lChild);
- Destroy(r->rChild);
- delete r;
- r=NULL;
- }
- }
-
- template<class T>
- void BinaryTree<T>::Change(BTNode<T> *r)//將二叉樹bt所有結(jié)點的左右子樹交換
- {
- BTNode<T> *p;
- if(r){
- p=r->lChild;
- r->lChild=r->rChild;
- r->rChild=p; //左右子女交換
- Change(r->lChild); //交換左子樹上所有結(jié)點的左右子樹
- Change(r->rChild); //交換右子樹上所有結(jié)點的左右子樹
- }
- }
-
- template<class T>
- void BinaryTree<T>::MakeTree(T pData,BinaryTree<T> leftTree,BinaryTree<T> rightTree)
- {
- root = new BTNode<T>();
- root->data = pData;
- root->lChild = leftTree.root;
- root->rChild = rightTree.root;
- }
(3)MinHeap.h 最小堆實現(xiàn)
3,、貪心選擇性質(zhì)
二叉樹T表示字符集C的一個最優(yōu)前綴碼,,證明可以對T作適當(dāng)修改后得到一棵新的二叉樹T”,在T”中x和y是最深葉子且為兄弟,,同時T”表示的前綴碼也是C的最優(yōu)前綴碼,。設(shè)b和c是二叉樹T的最深葉子,且為兄弟,。設(shè)f(b)<=f(c),,f(x)<=f(y)。由于x和y是C中具有最小頻率的兩個字符,,有f(x)<=f(b),,f(y)<=f(c),。首先,在樹T中交換葉子b和x的位置得到T',,然后再樹T'中交換葉子c和y的位置,,得到樹T''。如圖所示:
由此可知,,樹T和T'的前綴碼的平均碼長之差為:
因此,,T''表示的前綴碼也是最優(yōu)前綴碼,且x,y具有相同的碼長,,同時,,僅最優(yōu)一位編碼不同。
4,、最優(yōu)子結(jié)構(gòu)性質(zhì)
二叉樹T表示字符集C的一個最優(yōu)前綴碼,,x和y是樹T中的兩個葉子且為兄弟,z是它們的父親,。若將z當(dāng)作是具有頻率f(z)=f(x)+f(y)的字符,,則樹T’=T-{x,y}表示字符集C’=C-{x, y} ∪ { z}的一個最優(yōu)前綴碼。因此,,有:
如果T’不是C’的最優(yōu)前綴碼,,假定T”是C’的最優(yōu)前綴碼,那么有,,顯然T”’是比T更優(yōu)的前綴碼,,跟前提矛盾,!故T'所表示的C'的前綴碼是最優(yōu)的,。
由貪心選擇性質(zhì)和最優(yōu)子結(jié)構(gòu)性質(zhì)可以推出哈夫曼算法是正確的,即HuffmanTree產(chǎn)生的一棵最優(yōu)前綴編碼樹,。
程序運行結(jié)果如圖:
|