哈夫曼树(Huffman Tree)
一种特殊的二叉树,用于编码和解码数据。它是一种最优二叉树,也称为最优前缀编码树或最优二进制编码树。
哈夫曼树的构建过程是基于贪心算法的,即每次选择权值最小的两个节点合并成一个新节点,直到所有节点都合并成一个根节点为止。在构建过程中,每个节点的权值是其子树中所有叶子节点的权值之和。
哈夫曼树的主要应用是数据压缩,通过将频率高的字符用较短的编码表示,可以大大减小数据的存储空间。在压缩过程中,需要先构建哈夫曼树,然后根据每个字符在哈夫曼树中的路径生成对应的编码表,最后将原始数据按照编码表进行编码。解压缩时,根据编码表将编码后的数据还原成原始数据。
先上代码:
#define _CRT_SECURE_NO_WARNINGS 1
#include <iostream>
#include <fstream>
#include <string.h>
using namespace std;
#define MaxSize 1024 // 读入文件的上限
#define OK 1
#define ERROR 0
typedef int Status;
typedef struct wordcnt { // 统计字符和对应的次数
char ch;
int cnt = 0;
}Count;
typedef struct NumCount { // 统计次数的外部封装
Count count[MaxSize];
int length = 0;
}NumCount;
typedef struct HTree { // 哈夫曼树结构
char data;
int weight;
int parent, lchild, rchild;
}HTNode, * HuffmanTree;
typedef struct HCode { // 编码结构
char data;
char* str;
}*HuffmanCode;
Status ReadData(char* source); // 读入文件
Status WordCount(char* data, NumCount* paraCnt); // 统计次数
Status Show(NumCount* paraCnt); // 展示次数
Status CreateHuffmanTree(HuffmanTree& HT, int length, NumCount cntarray); // 创建哈夫曼树
Status select(HuffmanTree HT, int top, int* s1, int* s2); // 选择权重最小的两个节点
Status CreateHuffmanCode(HuffmanTree HT, HuffmanCode& HC, int length); // 创建哈夫曼编码
Status Encode(char* data, HuffmanCode HC, int length); // 将读入的文件编码,写到txt文件
Status Decode(HuffmanTree HT, int length); //读入编码文件,解码
int main(int argc, char** argv) {
char data[MaxSize];
NumCount Cntarray;
ReadData(data); // 读入数据
WordCount(data, &Cntarray); // 统计次数
HuffmanTree tree;
CreateHuffmanTree(tree, Cntarray.length, Cntarray); // 建树
HuffmanCode code;
CreateHuffmanCode(tree, code, Cntarray.length); // 创建编码
Encode(data, code, Cntarray.length); // 生成编码文件
Decode(tree, Cntarray.length); // 解码
cout << "Please view the generated TXT file to check the result" << endl;
return 0;
}
Status ReadData(char* source)
{
//打开文件读入数据
ifstream infile;
infile.open("in.txt");
cout << "Reading..." << endl;
cout << "the input file is:" << endl;
infile.getline(source, MaxSize);
cout << source << endl;
infile.close();
cout << endl;
return OK;
}
Status WordCount(char* data, NumCount* paraCnt)
{
int flag;// 标识是否已经记录
int len = strlen(data);
for (int i = 0; i < len; ++i)
{
flag = 0;
for (int j = 0; j < paraCnt->length; ++j)
{
if (paraCnt->count[j].ch == data[i]) // 若已有记录,直接++
{
++paraCnt->count[j].cnt;
flag = 1;
break;
}
}
if (!flag) // 没有记录,则新增
{
paraCnt->count[paraCnt->length].ch = data[i];
++paraCnt->count[paraCnt->length].cnt;
++paraCnt->length;
}
}
return OK;
}
Status Show(NumCount* paraCnt)
{
cout << "the length is " << paraCnt->length << endl;
for (int i = 0; i < paraCnt->length; ++i)
{
cout << "The character " << paraCnt->count[i].ch << " appears " << paraCnt->count[i].cnt << endl;
}
cout << endl;
return OK;
}
Status CreateHuffmanTree(HuffmanTree& HT, int length, NumCount cntarray)
{
if (length <= 1) return ERROR;
int s1, s2;
int m = length * 2 - 1; // 没有度为1的节点,则总结点是2*叶子节点数-1个
HT = new HTNode[m + 1];
for (int i = 1; i <= m; ++i) // 初始化
{
HT[i].parent = 0;
HT[i].lchild = 0;
HT[i].rchild = 0;
}
for (int i = 1; i <= length; ++i)
{
HT[i].data = cntarray.count[i - 1].ch;
HT[i].weight = cntarray.count[i - 1].cnt;
}
for (int i = length + 1; i <= m; ++i)
{
select(HT, i - 1, &s1, &s2); // 从前面的范围里选择权重最小的两个节点
HT[s1].parent = i;
HT[s2].parent = i;
HT[i].lchild = s1;
HT[i].rchild = s2;
HT[i].weight = HT[s1].weight + HT[s2].weight; // 得到一个新节点
}
return OK;
}
Status select(HuffmanTree HT, int top, int* s1, int* s2)
{
int min = INT_MAX;
for (int i = 1; i <= top; ++i) // 选择没有双亲的节点中,权重最小的节点
{
if (HT[i].weight < min && HT[i].parent == 0)
{
min = HT[i].weight;
*s1 = i;
}
}
min = INT_MAX;
for (int i = 1; i <= top; ++i) // 选择没有双亲的节点中,权重次小的节点
{
if (HT[i].weight < min && i != *s1 && HT[i].parent == 0)
{
min = HT[i].weight;
*s2 = i;
}
}
return OK;
}
Status CreateHuffmanCode(HuffmanTree HT, HuffmanCode& HC, int length)
{
HC = new HCode[length + 1];
char* cd = new char[length]; // 存储编码的临时空间
cd[length - 1] = '\0'; // 方便之后调用strcpy函数
int c, f, start;
for (int i = 1; i <= length; ++i)
{
start = length - 1; // start表示编码在临时空间内的起始下标,由于是从叶子节点回溯,所以是从最后开始
c = i;
f = HT[c].parent;
while (f != 0)
{
--start; // 由于是回溯,所以从临时空间的最后往回计
if (HT[f].lchild == c)
cd[start] = '0';
else
cd[start] = '1';
c = f;
f = HT[c].parent;
}
HC[i].str = new char[length - start]; // 最后,实际使用的编码空间大小是length-start
HC[i].data = HT[i].data;
strcpy(HC[i].str, &cd[start]); // 从实际起始地址开始,拷贝到编码结构中
}
delete cd;
return 1;
}
Status Encode(char* data, HuffmanCode HC, int length)
{
ofstream outfile;
outfile.open("code.txt");
for (int i = 0; i < strlen(data); ++i) // 依次读入数据,查找对应的编码,写入编码文件
{
for (int j = 1; j <= length; ++j)
{
if (data[i] == HC[j].data)
{
outfile << HC[j].str;
}
}
}
outfile.close();
cout << "the code txt has been written" << endl;
cout << endl;
return OK;
}
Status Decode(HuffmanTree HT, int length)
{
char codetxt[100001];
ifstream infile;
infile.open("code.txt");
infile.getline(codetxt, MaxSize * length);
infile.close();
ofstream outfile;
outfile.open("out.txt");
int root = 2 * length - 1; // 从根节点开始遍历
for (int i = 0; i < strlen(codetxt); ++i)
{
if (codetxt[i] == '0') root = HT[root].lchild; //为0表示向左遍历
else if (codetxt[i] == '1') root = HT[root].rchild; //为1表示向右遍历
if (HT[root].lchild == 0 && HT[root].rchild == 0) // 如果已经是叶子节点,输出到输出文件中,然后重新回到根节点
{
outfile << HT[root].data;
root = 2 * length - 1;
}
}
outfile.close();
cout << "the output txt has been written" << endl;
cout << endl;
return OK;
}
运行结果
不知道为什么代码运行有误,在VS上运行错误,读不到文本内容
心得:
哈夫曼树是一种特殊的二叉树,它的每个叶子节点都有一个权值,而且权值越大的叶子节点越靠近根节点。哈夫曼树的构造方法是:
首先将所有的叶子节点按照权值从小到大排序,然后每次取出最小的两个节点,合并成一个新的节点,新节点的权值为两个子节点的权值之和,然后将新节点插入到原来的序列中,重复这个过程直到只剩下一个节点,这个节点就是哈夫曼树的根节点。
其次,从根节点开始遍历这棵树,并且给每条边赋予一个bit值:左边为0,右边为1。最后我们可以得到每个叶子节点对应的编码:从根节点到叶子节点经过的所有边上的bit值组成。