哈夫曼树要解决的问题是对不同频率元素查找的效率最高,哈夫曼编码代码如下:
首先要明白,哈夫曼树要查找的元素都在叶结点上,而且哈夫曼树的最终元素是确定的为2*n-1个。这样我们可以用结点数组来表示树,每一个结点的 child和parent指向自己的父亲和儿子的索引。首先实现select函数,从结点数组选出权重最小的两个合并,存储新结点,构造出哈夫曼树,然后从叶结点开始向上查找,左儿子为0,右儿子为1得到哈夫曼编码,最后再把编码反转即可。
//huffmanCoding.c
#include <stdio.h>
#include <limits.h>
#include <string.h>
#include <stdlib.h>
#define N 6
typedef struct huffNode
{
int weight; //权重
int lchild, rchild, parent; //左右子节点和父节点
}HTNode, *HuffTree;
typedef char **HuffCode;
//找出数组中无父节点且权值最小的两个节点下标,分别用s1和s2保存
void select(const HuffTree &HT, int n, int &s1, int &s2);
//HT:哈夫曼树,HC:哈夫曼编码,w:构造哈夫曼树节点的权值,n:构造哈夫曼树节点的个数
void HuffmanCode(HuffTree &HT, HuffCode &HC, int *w, int n);
int main()
{
int i;
char key[N] = { '0','A','B','C','D','E' };//第0个元素保留不用
int w[N] = { 0,1,2,4,5,6 }; //第0个元素保留不用
HuffTree HT;
HuffCode HC;
HuffmanCode(HT, HC, w, N - 1);
for (i = 1; i < N; i++)
printf("%c:%s\n", key[i], HC[i]);
printf("\n");
return 0;
}
//找出数组中权值最小的两个节点下标,分别用s1和s2保存
void select(const HuffTree &HT, int n, int &s1, int &s2)
{
int i;
s1 = s2 = 0;
int min1 = INT_MAX;//最小值,INT_MAX在<limits.h>中定义的
int min2 = INT_MAX;//次小值
for (i = 1; i <= n; ++i)
{
if (HT[i].parent == 0)
{//筛选没有父节点的最小和次小权值下标
if (HT[i].weight < min1)
{//如果比最小值小
min2 = min1;
s2 = s1;
min1 = HT[i].weight;
s1 = i;
}
else if ((HT[i].weight >= min1) && (HT[i].weight < min2))
{//如果大于等于最小值,且小于次小值
min2 = HT[i].weight;
s2 = i;
}
else
{//如果大于次小值,则什么都不做
;
}
}
}
}
//HT:哈夫曼树,HC:哈夫曼编码,w:构造哈夫曼树节点的权值,n:构造哈夫曼树节点的个数
void HuffmanCode(HuffTree &HT, HuffCode &HC, int *w, int n)
{
int s1;
int s2;
int m = 2 * n - 1; //容易知道n个节点构造的哈夫曼树是2n-1个节点
int i, c, f, j;
char *code; //暂存编码的
HT = (HuffTree)malloc((m + 1) * sizeof(HTNode)); //0单元未使用
for (i = 1; i <= n; i++)
HT[i] = { w[i],0,0,0 };//初始化前n个节点(构造哈夫曼树的原始节点)
for (i = n + 1; i <= m; i++)
HT[i] = { 0,0,0,0 }; //初始化后n-1个节点
//构建哈夫曼树
for (i = n + 1; i <= m; i++)
{
select(HT, i - 1, s1, s2);//找出前i-1个节点中权值最小的节点下标
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;
}
//哈夫曼编码
HC = (char **)malloc((n) * sizeof(char *));
//暂存编码
code = (char *)malloc(n * sizeof(char));//使用了第0单元
for (i = 1; i <= n; i++)
{
for (c = i, f = HT[c].parent, j = 0; f != 0; c = HT[c].parent, f = HT[c].parent, j++)
{//从叶子扫描到根
if (HT[f].lchild == c)
{
code[j] = '0';
}
else if (HT[f].rchild == c)
{
code[j] = '1';
}
else
{//否则什么也不做
;
}
}
code[j] = '\0';
for (int k = 0; k < j/2; k++)
{
char tmp = code[k];
code[k] = code[j - 1- k];
code[j - 1- k] = tmp;
}
HC[i] = (char *)malloc(strlen(code) * sizeof(char));
strcpy(HC[i], code);
}
}