霍夫曼编码

霍夫曼编码_第1张图片 霍夫曼编码(Huffman Coding)是一种编码方式,是一种用于无损数据压缩的熵编码(权编码)算法。1952年,David A. Huffman在麻省理工攻读博士时所发明的,并发表于《一种构建极小多余编码的方法》(A Method for the Construction of Minimum-Redundancy Codes)一文。

在计算机数据处理中,霍夫曼编码使用变长编码表对源符号(如文件中的一个字母)进行编码,其中变长编码表是通过一种评估来源符号出现机率的方法得到的,出现机率高的字母使用较短的编码,反之出现机率低的则使用较长的编码,这便使编码之后的字符串的平均长度、期望值降低,从而达到无损压缩数据的目的。
例如,在英文中,e的出现机率最高,而z的出现概率则最低。当利用霍夫曼编码对一篇英文进行压缩时,e极有可能用一个位来表示,而z则可能花去25个位(不是26)。用普通的表示方法时,每个英文字母均占用一个字节(byte),即8个位。二者相比,e使用了一般编码的1/8的长度,z则使用了3倍多。倘若我们能实现对于英文中各个字母出现概率的较准确的估算,就可以大幅度提高无损压缩的比例。
霍夫曼树又称最优二叉树,是一种带权路径长度最短的二叉树。所谓树的带权路径长度,就是树中所有的叶结点的权值乘上其到根结点的路径长度(若根结点为0层,叶结点到根结点的路径长度为叶结点的层数)。树的路径长度是从树根到每一结点的路径长度之和,记为WPL=(W1*L1+W2*L2+W3*L3+...+Wn*Ln),N个权值Wi(i=1,2,...n)构成一棵有N个叶结点的二叉树,相应的叶结点的路径长度为Li(i=1,2,...n)。可以证明霍夫曼树的WPL是最小的。
代码:
//本来参照维基百科里的示例写的,可以向量容器始终搞不定,所以参照CSDN使用动态数组来完成
#include <iostream>    
 
#include <cstdlib>  
#include <string> 
 
using namespace std; 
 
typedef struct  
{  
unsigned int weight;//结点权值  
unsigned int parent,lchild,rchild;//结点的父指针,左右孩子指针  
}HTNode,*HuffmanTree;//动态分配数组存储哈夫曼树  
typedef char **HuffmanCode;//动态分配数组存储哈夫曼编码表 
HuffmanTree HT;//哈夫曼树HT  
HuffmanCode HC;//哈夫曼编码表HC  
unsigned int *w;//w存放叶子结点权值      
 
void CreateHuffmanTree(HuffmanTree &,unsigned int*,int);//生成一棵哈夫曼树  
void HuffmanCoding(HuffmanTree,HuffmanCode&,int);//对哈夫曼树进行编码   
void PrintHuffmanCode(HuffmanCode,unsigned int*,int);//显示哈夫曼编码  
void Select(HuffmanTree,int,int&,int&);//在数组中寻找权值最小的两个结点 
 
int main()  
{  
int   n,i;                  //n是哈夫曼树叶子结点数   
char   choose='y';  //用于选择程序是否退出     
 
while(choose!='N' && choose!='n')  
{  
   cout << "请输入叶子结点数目:";
  
   do  
   {
    cin >> n;//输入叶子结点数
    cin.clear();
    cin.ignore();
   } while( n<=1 );
 
   w=(unsigned int*)malloc(n*sizeof(unsigned int));//开辟空间存放权值  
    for(i=0;i<n;i++)    
   {   
    cout << "W" << i+1 << "的权值:";  
    cin >> w[i];//输入各叶子结点权值  
   }          
   CreateHuffmanTree(HT,w,n);//生成哈夫曼树  
   HuffmanCoding(HT,HC,n);//进行哈夫曼编码  
   cout << endl;
   PrintHuffmanCode(HC,w,n);//显示哈夫曼编码  
   cout << "还要继续吗?(Y/N)";  
  cin >> choose;  
}
 
system("pause");
return 0;
 
void CreateHuffmanTree(HuffmanTree &HT,unsigned int *w,int n)  
{//w存放n个结点的权值,将构造一棵哈夫曼树HT  
int i,m;  
int s1,s2;  
HuffmanTree p;  
if(n<=1)
   return;  
   
m=2*n-1;//n个叶子结点的哈夫曼树,有2*n-1个结点  
HT=(HuffmanTree)malloc((m+1)*sizeof(HTNode));//开辟2*n各结点空间,0号单元不用  
   
for(p=HT+1,i=1;i<=n;++i,++p,++w)//进行初始化  
{  
   p->weight=*w;  
   p->parent=0;  
   p->lchild=0;  
   p->rchild=0;  
}
 
for(;i<=m;++i,++p)  
{  
   p->weight=0;  
   p->parent=0;  
   p->lchild=0;  
   p->rchild=0;  
}
 
for(i=n+1;i<=m;++i)//建哈夫曼  
{  
  //从HT[1...i-1]中选择parent为0且weight最小的两个结点,其序号分别为s1和s2  
   Select(HT,i-1,s1,s2);      
   HT[s1].parent=i;   HT[s2].parent=i;//修改s1和s2结点的父指针parent  
   HT[i].lchild=s1;   HT[i].rchild=s2;//修改i结点的左右孩子指针  
   HT[i].weight=HT[s1].weight+HT[s2].weight;//修改权值  
 
 
}  
 
void   HuffmanCoding(HuffmanTree HT,HuffmanCode &HC,int n)  
{  
//将有n个叶子结点的哈夫曼树HT进行编码,所编的码存放在HC中  
//方法是从叶子到根逆向求每个叶子结点的哈夫曼编码   
int i,c,f,start;  
char *cd;  
HC=(HuffmanCode)malloc((n+1)*sizeof(char*));//分配n个编码的头指针向量  
cd=(char *)malloc(n*sizeof(char));//开辟一个求编码的工作空间  
cd[n-1]='/0';//编码结束符  
for(i=1;i<=n;++i)//逐个地求哈夫曼编码   
{  
   start=n-1;//编码结束位置  
   for(c=i,f=HT[i].parent;f!=0;c=f,f=HT[f].parent)//从叶子到根逆向求编码   
    if(HT[f].lchild==c)      
     cd[--start]='0';//若是左孩子编为'0'  
    else    
     cd[--start]='1';//若是右孩子编为'1'    
            HC[i]=(char *)malloc((n-start)*sizeof(char));//为第i个编码分配空间   
            strcpy(HC[i],&cd[start]);//将编码从cd复制到HC中  
}  
free(cd);//释放工作空间  
 
void PrintHuffmanCode(HuffmanCode HC,unsigned int *w,int n)  
{  
//显示有n个叶子结点的哈夫曼树的编码表  
int i;  
cout << "哈夫曼编码如下:/n";  
for(i=1;i<=n;i++)  
{  
   cout << HC[i]<<"/t" <<"W"<<i<< "ȨֵΪ£º" <<w[i-1]<<endl;     
}  
cout<<endl;  
 
void Select(HuffmanTree   HT,int   t,int&s1,int&s2)  
{  
//在HT[1...t]中选择parent不为0且权值最小的两个结点,其序号分别为s1和s2,s1存放最小的,s2存放次小的    
int i=0;  
int j=0;  
int k=0;  
int least=0;  
int second=0;  
 
for(i=1;i<=t;i++)  
{  
   if(HT[i].parent==0)
    break;  
}          
 
for(j=i+1;j<=t;j++)  
{  
   if(HT[j].parent==0)
    break;  
}      
 
if(HT[i].weight<HT[j].weight)  
{  
   least=i;  
   second=j;  
}              
 
else    
{  
   least=j;  
   second=i;  
}          
 
for(k=j+1;k<=t;k++)  
{  
   if(HT[k].parent==0)  
   {  
    if(HT[k].weight<HT[least].weight)  
    {  
     second=least;  
     least=k;  
    }          
  
    else
     if(HT[k].weight>=HT[least].weight&&HT[k].weight<HT[second].weight)  
     second=k;
   }          
}     
 
s1=least;  
s2=second;
 
 
运行示例:

 

你可能感兴趣的:(霍夫曼编码)