范式Huffman编码的实现

传统的Huffman 编码依赖于二叉树，一方面Huffman编码生成的速度不够快，另一方面在没有指针的程序设计语言中要实现二叉树结构又非常复杂。所以懂得构造范式Huffman编码去解决一些问题在某些方面就很有必要。范式Huffman编码相比与传统Huffman编码，不需要建造Huffman树，只需求出每个符号在Huffman树中从树根到到该符号节点的路径长度即可，此路径长度即为该符号的编码长度。在求出各符号的编码长度之后，按如下规则进行编码：

1）求出编码长度相同的符号的数目，如

　　Length　　　　Num
　　　5　　　　　　2
　　　4　　　　　　3
　　　3　　　　　　2
　　　2　　　　　　2
　　　1　　　　　　0

2）生成编码，char1,char2,char3,char4,char5,char6,char7,char8,char9的深度为5,5,4,4,4,3,3,2,2所以char1 的编码为00000,长度为5；char2的编码为 char1的编码加1,再右移(char1长度减char2长度)位,即00001；同理char3,char4,char5的编码为0001,0010,0011；到char6时为0100右移4(char4的长度)减去3(char5的长度),即010；char7则为011;char8,char9的编码为10，11

 

以上为算法，下面给出代码：

#include<stdio.h> #include<math.h> #define n 26 /****************************************************/ int HuffmanCode(int *w);//计算符号的范式Huffman编码 void select(int a[][51],int &s1,int &s2,int boundary);//从a数组中选出权值最小的两个数 void printb(int x,int m);//将x以二进制形式输出，位数为m位 /****************************************************/ int HuffmanCode(char *w) { int i,s1,s2,f,m; int a[2][2*n-1]={0},b[n]={0};/*这里使用二维数组a[][]模拟Huffman树的建造过程，其中第一行存放符号的出现次数，第二行存放它的父节点*/ m=2*n-1; for(i=0;i<n;i+=1)/*将符号的出现次数copy到a[][]数组的第一行*/ a[0][i]=*(w+i); for(i=0;i<n;i++) if(a[0][i]>0)printf("%c ,%d/n",i+97,a[0][i]); for(;i<m;i+=1) { select(a,s1,s2,i-1);/*选择a[0][0……i-1]中父节点为0且权值最小的两个节点，其位置为s1、s2*/ a[0][i]=a[0][s1]+a[0][s2]; a[1][s1]=i;a[1][s2]=i; } for(i=0;i<n;i+=1) {/*求每个符号在Huffman树中深度，即从树根到表示该符号的叶子的路径长度，也就是每个符号的编码长度*/ int dis=0; for(f=a[1][i];f;f=a[1][f]) dis+=1; a[1][i]=dis; } for(i=0;i<n;i++) if(a[1][i]<=3)printf("%c ,%d/n",i+97,a[1][i]); /*下面进行编码*/ int e=0,p,p0,flag=0,c=0; while(c<n) { p=0; for(i=0;i<n;i+=1) if(a[1][i]>p)p=a[1][i];//计算a[1][i]中最大的数 for(i=0;i<n;i+=1)//对编码长度为p的符号进行编码 if(a[1][i]==p) { e+=1; if(flag){e=e/pow(2,p0-p);flag=0;}/*将e进行右移(p0-p)位的操作，这里的右移是对e的二进制的数值而言的*/ printf("%c:",i+97); printb(e,p); printf("/n"); c+=1; a[1][i]=-1;//对已编码的符号进行标记，本程序里标记为-1 }//if、for p0=p;flag=1; }//while return 1; } void select(int a[][51],int &s1,int &s2,int boundary) {//求出a[0][0]~a[0][i-1]中权值最小的两个元素，返回其位置s1、s2 int k=0,s[2],t; for(int i=0;i<2;i+=1) { for(int l=0;l<=boundary;l+=1) if(!a[1][l]){t=a[0][l];k=l;break;} for(int j=0;j<=boundary;j+=1) if(!a[1][j]&&a[0][j]<t){t=a[0][j];k=j;} s[i]=k;a[1][k]=1; } s1=s[0];s2=s[1]; //a[1][s1]=0;a[1][s2]=0; } void printb(int x,int m) /*输入整数x的二进制数的后n位*/ { if(m>0){ putchar('0'+((unsigned)(x&(1<<(m-1)))>>(m-1))); printb(x,m-1); /*递归调用，输出后n-1位*/ } } int main() {//主函数 char ch[n]="",c; printf("Please input a string:/n"); while((c=getchar())!=10)ch[c-97]++;//统计输入的各种字符的次数，即统计各种符号出现的次数 HuffmanCode(ch); return 1; }