哈夫曼编码程序设计

算法与数据结构课程设计哈夫曼编码/译码器设计学生姓名：学号：专业：（计算机科学与技术）年级：（大二）指导教师：（汪洋）2009年6月19日哈夫曼编码/译码器问题描述：利用哈夫曼编码进行通信可以大大提高信道利用率，缩短信息传输时间，降低传输成本，但是，这要求在发送端通过一个编码系统对待传数据预先编码，在接收端将传来的数据进行译码（复原）。

对于双工信道（即可以双向传输信息的信道）每端都需要一个完整的编/译码系统。

试为这样的信息收发站写一哈夫曼编/译码系统。

基本要求：I：初始化（Initialization）。

从终端读入字符集大小n，以及n个字符和n个权值，建立哈夫曼树，并将它存于文件hfmTree中。

E：编码（Encoding）。

利用已建好的哈夫曼树（如不在内存，则从文件hfmTree中读入），对文件ToBeTran中的正文进行编码，然后将结果存入文件CodeFile中。

D：译码（Decoding）。

利用已建好的哈夫曼树将文件CodeFile中的代码进行译码，结果存入文件TextFile中。

P：打印代码文件（Print）。

将文件CodeFile以紧凑格式显示在终端上，每行50个代码。

同时将此字符形式的编码文件写入文件CodePrin中。

T：打印哈夫曼树（Tree printing）。

将已在中的哈夫曼树以直观的方式（树或凹入表形式）显示在终端上，同时将此字符形式的哈夫曼树写入文件TreePrint中。

大体解题思路:（1）对输入的一段欲编码的字符串进行统计各个字符出现的次数，并它们转化为权值{w1,w2,……，wN}构成n棵二叉树的集合F={T1,T2,……，Tn}把它们保存到结构体数组HT[n]中,其中{Ti是按它们的ASCⅡ码值先后排序。

其中每棵二叉树Ti中只有一个带权为Wi的根结点的权值为其左、右子树上根结点的权值之和。

（2）在HT[1..i]中选取两棵根结点的权值最小且没有被选过的树作为左右子树构造一棵新的二叉树，且置新的二叉树的根结点的权值为左、右子树上根结点的权值之和。

（3）哈夫曼树已经建立后，从叶子到根逆向求每一个字符的哈夫曼编码。

概要设计：实现的功能：1.查看原文（showpassage()），2.字符统计(showdetail())，3.哈夫曼编码并输出内容(HuffmanCoding(HT,HC,30,w))，4.输出整篇文章的码字(printcode())，5.求最小权值(minweight())，6.最码字进行解码(decode())，7.测试解码(testdecode())，8.推出(break)。

（1）定义结构体HTNOde,*HuffmanTree;typedefchar**HuffmanCode;定义堆结构体RedType;定义全局变量；（2）编程序：测试编码(void testdecode())；解码(void decode)；求最小权值(void minweight())；打印码字(voidprintcode())；堆排序(void HeapAdjust(SqList&L,ints,int m))；哈夫曼编码(voidHuffmanCoding(HuffmanTree &HT,HuffmanCode&HC,float *w,int n ))；输出个字符的次数比例(voidshowdetail())；输出原英文文章并做统计(voidshowpassage())；（3）主函数详细设计：（1）定义结构体1.哈夫曼结构体typedef struct{float weight;unsigned int parent, lchild, rchild;}HTNode, *HuffmanTree;typedef char **HuffmanCode ;2.堆结构体typedef struct{float key; //关键字项int otherinfo; //其他数据项（此题目中用不到）}RedType;typedef struct{RedType r[MAXSIZE+1]; //r[0]闲置用作哨兵int length; //顺序表长度}SqList;3.定义全局变量HuffmanTree HT; //赫夫曼树HuffmanCode HC; //码值FILE *fp, *fp1, *fp2;int a[30] = {0};float b[30];float *w; //权（2）程序代码1.测试解码（可以输入一个不正确的二进制码串）void testdecode(){char str[200]; //存放自己输入的码子int p, p1, i; //解码的线索char ch;printf("\n请根据以上各个字符的编码输入一串二进制码字(200个以内):\n");gets(str);printf("以上码子解码后为:\n");p = 59; //最初令p为树根整数，p由根顺着树枝一直到树叶i = 0;ch = str[i++];while (ch!='\0'){for ( ; ch!='\0' && (HT[p].lchild!=0||HT[p].rchild!=0) ; ){if (ch == '0'){p = HT[p].lchild;}else if(ch == '1'){p = HT[p].rchild;}else{printf("\n你输入了'0','1'之外的字符，无法正确译码，请检查!\n");return ; //直接结束}ch = str[i++]; //下一个码字不断的取下一个}if(p <= 30) //小于等于30的时候才正确，有可能最后一位p还没有在1-30范围内的时候就没有二进制码了，也就是说二进制码最后不完整。

switch (p){case 27 : printf(",");break;case 28 : printf(".");break;case 29 : printf(" ");break;case 30 : printf("\n");break;default : printf("%c", p+96);}p1 = p; //让p1记住pp = 59; //这里p要重置为59，因为经过上面的程序p已经变化了，不重置为1则HT[p].lchild!=0||HT[p].rchild!=0所以for语句无法进行} //while循环if (p1 > 30)printf("\n以上正确译出了前面的字符，由于你输入的二进制码不完整，最后一位字符无法译出，请检查!\n");}2.下面是解码void decode(){int p;char ch;printf("\n\n对码子解码后的如下:\n");fp1 = fopen("bianma.txt","r");if(!fp1){printf("can not open this file!\n");exit(0);}p = 59; //最初令p为任意一个非零整数，p由根顺着树枝一直到树叶ch = fgetc(fp1);fp2 = fopen("jiema.txt","w");if (!fp2){printf("can not open this file!\n");exit(0);}while (ch!=EOF){for ( ; HT[p].lchild!=0||HT[p].rchild!=0 ; ){if (ch == '0'){p = HT[p].lchild;}else{p = HT[p].rchild;}ch = fgetc(fp1); //下一个码字不断的取下一个}switch (p){case 27 : printf(",");fprintf(fp2,",");break;case 28 : printf(".");fprintf(fp2, ".");break;case 29 : printf(" ");fprintf(fp2, " ");break;case 30 : printf("\n");fprintf(fp2, "\n");break;default : printf("%c", p+96);fprintf(fp2, "%c", p+96);} p = 59; //这里p要重置为59，因为经过上面的程序p已经变化了，不重置为1则HT[p].lchild!=0||HT[p].rchild!=0所以for语句无法进行} //while循环printf("\n");fprintf(fp2, "\n");fclose(fp1);fclose(fp2);}3.求最小权值void minweight(){float Weight = 0;int i;for (i = 0 ; i < 30 ; i++)Weight = Weight + strlen(HC[i+1])*b[i];printf("最小权值是:%f\n",Weight);}4.打印码子void printcode(){char ch;fp = fopen("Huffman.txt","r");if (!fp){printf("can not open this file!\n");exit(0);}fp1 = fopen("bianma.txt","w");if (!fp1){printf("can not open this file!\n");exit(0);}printf("\n原英文文章经编码后如下:\n");ch = fgetc(fp);while (ch!=EOF){if (ch > 96 && ch < 123) //小写字母{printf("%s", HC[ch-96]);fputs(HC[ch-96], fp1); //输出到文件里面ch = fgetc(fp);continue;}if (ch > 64 && ch < 91) //大小字母{ printf("%s", HC[ch-64]);fputs(HC[ch-64], fp1); //输出到文件里面ch = fgetc(fp);continue;}if (ch == ','){printf("%s", HC[27]);fputs(HC[27], fp1);ch = fgetc(fp);continue;}if (ch == '.'){printf("%s", HC[28]);fputs(HC[28], fp1);ch = fgetc(fp);continue;}if (ch == ' '){printf("%s", HC[29]);fputs(HC[29], fp1);ch = fgetc(fp);continue;}if (ch == '\n'){printf("%s", HC[30]);fputs(HC[30], fp1);ch = fgetc(fp);continue;}}printf("\n\n");fclose(fp);fclose(fp1);}5.堆排序void HeapAdjust(SqList &L, int s, int m) {RedType rc;int j;rc = L.r[s];for (j = 2 * s ; j <= m ; j*=2){if (j < m && L.r[j].key > L.r[j+1].key) //即使等于也不要动，不用加1 ++j;if (rc.key <= L.r[j].key) //即使等于也不要动，直接跳出来break;L.r[s] = L.r[j];s = j;}L.r[s] = rc;}void select(HuffmanTree t, int i, int &s1, int &s2){ //此函数被调用一次则就用堆排序选择两个最小的赋给s1和s2SqList L;RedType rc;int j, k, n = 1;L.length = 0;for (j = 1 ; j <= i ; j++){ if (t[j].parent!=0)continue;L.r[n].key = t[j].weight; //赋值好,用堆排序L.r[n].otherinfo = j; //存放序号,j是一直在加一的，循环一次加1，但是n不是的只有在符合条件的情况下才加1n++;L.length++; //这样写很巧妙的}for (k = L.length/2 ; k > 0 ; --k)HeapAdjust(L,k,L.length);s1 = L.r[1].otherinfo; //最小的选出来了！/***把最小的换到最下面***/rc = L.r[1];L.r[1] = L.r[L.length]; //此次的select函数，进行堆排序的只有L.length 个（parent！=0的没有在里面），所以这里是L.length 而不是iL.r[L.length] = rc;HeapAdjust(L,1,L.length-1);s2 = L.r[1].otherinfo; //次小的选出来了(这里当输入1 4 2 8四个数时，第一次选出的s1=1,s2=3是对的，但第二次选出的s1=5，但s2是随机数)}6.赫夫曼编码void HuffmanCoding(HuffmanTree &HT, HuffmanCode &HC, float *w, int n) // 算法6.12{ // w存放n个字符的权值(均>0),构造赫夫曼树HT,并求出n个字符的赫夫曼编码HCint m, i, s1, s2, start, k;unsigned c, f;HuffmanTree p;char *cd;if (n <= 1)return;m = 2 * n - 1;w = (float *)malloc(30*sizeof(float));for (i = 0 ; i < 30 ; i++)*(w+i) = b[i];HT = (HuffmanTree)malloc((m+1)*sizeof(HTNode)); // 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 = n + 1 ; i <= m ; ++i, ++p)(*p).parent=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; //最小的和次小的双亲已经不为0了，下次就不在它两中间找最小的和次小的了。