哈夫曼樹的數組實現

（本篇博客是本人第一篇數據結構的博客，有什么不足還望各位看官指出！！）

---

題目來源：SOJ 1000. Huffman Coding V1,V3

題目描述

---

V3：

• Description
  對輸入的英文大寫字母序列進行統計概率，然后構建Huffman樹，得出每個字母的Huffman編碼，輸出字母序列的總編碼長度。

• Input
  第一行是大寫字母個數n(0<n<=100)
  第二行為n個字母，中間以一個空格分隔。

• Output
  輸出字母序列的總編碼長度。

• Sample Input
  10
  S S U U U S U L U U

• Sample Output
  14

---

V1:

• Description
  對輸入的英文大寫字母序列進行統計概率，然后構建Huffman樹，輸出按照概率降序排序輸出Huffman編碼。

• Input
  第一行是大寫字母個數n(0<n<=100)
  第二行為n個字母，中間以一個空格分隔。
  建樹過程把權值較小的子樹作為左孩子，數據保證建樹過程不會出現左右子樹權值一樣的情況。

• Output
  假設輸入中共有m個不同的字母，按照出現概率降序輸出，每個字母單獨一行輸出。格式如下：

　　字母1 出現次數 Huffman編碼
　　字母2 出現次數 Huffman編碼
　　字母3 出現次數 Huffman編碼
　　…
　　字母m 出現次數 Huffman編碼

• Sample Input
  Copy sample input to clipboard
  10
  S S U U U S U L U U
• Sample Output
  U 6 1
  S 3 01
  L 1 00

---

算法描述

其實哈夫曼樹的建樹規則的話網上都有不少資料可以看，此處不予贅述。講講個人的一些收獲和想法：數組這種實現方法也是我在網上學習來的，簡單講就是先計算輸入數據對應的字符的權重並進行記錄。主要的結構體是哈夫曼樹的節點，存儲的是每個字符的權重以及左右子樹的權重，還有就是很有用的一個數據：父節點的權重。這樣就可以以權重代替指針域進行上下的尋址，可以減少由於指針使用不當帶來的內存問題。然后寫代碼的過程中遇到的一個纏最久的bug就是在建立非字符節點時候查找無父節點的節點的函數select()中，遇到的很棘手的一個問題是已經標記為有父節點的節點仍未被識別，后來才發現問題是出現在查找權重最小的節點的過程中，for循環的邊界寫錯了= =
不多說了，直接上代碼吧

---

個人代碼實現

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;

#define N 26
#define M (2*N-1)

struct huffmanTreeNode
{
	int weight;
	int left, right, parent;
};

struct huffmanCode
{
	char data;
	int weight;
	char code[N];
};

int initialize(huffmanCode hfmCodeSet[], int n)
{
	int set[N + 1];
	char inputStr[5];
	memset(set, 0, sizeof(set));
	memset(inputStr, 0, sizeof(inputStr));
	int i, j;
	
	huffmanCode cd;
	cd.data = 0;
	cd.weight = 0;
	memset(cd.code, 0, sizeof(cd.code));
	for (i = 0; i < N + 1; i ++)
		hfmCodeSet[i] = cd;

	for (i = 0; i < n; i ++) {
		scanf("%s", inputStr);
		set[inputStr[0] - 'A']++;
	}

	j = 1;
	for (i = 0; i < N + 1; i ++) {
		if (set[i] > 0) {
			hfmCodeSet[j].data = 'A' + i;
			hfmCodeSet[j].weight = set[i];
			j++;
		}
	}

	return (j - 1);
}

void select(huffmanTreeNode hfmTree[],int idx, int* i1, int* i2)
{
	int i, j;
	for (i = 1; i <= idx; i ++) {
		if (hfmTree[i].parent == 0) {
			*i1 = i;
			break;
		}
	}

	for (i = 1; i <= idx; i ++) {
		if (hfmTree[i].parent == 0 &&
			hfmTree[i].weight < hfmTree[*i1].weight) {
			*i1 = i;
		}
	}

	for (i = 1; i <= idx; i ++) {
		

		if (i != *i1 && hfmTree[i].parent == 0) {
			*i2 = i;
			break;
		}
	}

	for (i = 1; i <= idx; i ++) {
		

		if (i != *i1) {
			if (hfmTree[i].parent == 0 &&
				hfmTree[i].weight < hfmTree[*i2].weight) {
				*i2 = i;
			} 
		}
	}
}

void hfmCoding(huffmanCode hfmCodeSet[], huffmanTreeNode hfmTree[], int n)
{
	int i;
	char tempCode[N + 1];
	for (i = 1; i <= 2 * n - 1; i ++) {
		hfmTree[i].weight = (i <= n ? hfmCodeSet[i].weight : 0);
		hfmTree[i].parent = hfmTree[i].left = hfmTree[i].right = 0;
	}

	int minIdx1, minIdx2;
	for (i = n + 1; i <= 2 * n - 1; i ++) {
		select(hfmTree, i - 1, &minIdx1, &minIdx2);
		hfmTree[i].weight = hfmTree[minIdx1].weight + hfmTree[minIdx2].weight;
		hfmTree[i].left = minIdx1;
		hfmTree[i].right = minIdx2;
		hfmTree[minIdx1].parent = i;
		hfmTree[minIdx2].parent = i;
	}

	int start, childIdx, parentIdx;
	for (i = 1; i <= n; i ++) {
		start = n - 1;
		tempCode[n] = '\0';		

		childIdx = i;
		parentIdx = hfmTree[childIdx].parent;
		while (parentIdx) {
			if (hfmTree[parentIdx].left == childIdx) {
				tempCode[--start] = '0';
			} else if (hfmTree[parentIdx].right == childIdx) {
				tempCode[--start] = '1';
			}
			childIdx = parentIdx;
			parentIdx = hfmTree[childIdx].parent;
		}

		strcpy(hfmCodeSet[i].code, &tempCode[start]);

	}
}

int totalLength(huffmanCode hfmCodeSet[])
{
	int i, sum = 0;
	for (i = 1; i <= N; i ++) {
		if (hfmCodeSet[i].weight > 0) {
			sum += ((hfmCodeSet[i].weight) * strlen(hfmCodeSet[i].code));
		}
	}
	return sum;
}

void sortDisplayCode(huffmanCode hfmCodeSet[], int n)
{
	int i, j;
	huffmanCode tempCode;
	for (i = 1; i <= n - 1; i ++) {
		for (j = 1; j <= n - 1; j ++) {
			if (hfmCodeSet[j].weight < hfmCodeSet[j + 1].weight) {
				tempCode = hfmCodeSet[j];
				hfmCodeSet[j] = hfmCodeSet[j + 1];
				hfmCodeSet[j + 1] = tempCode;
			}
		}
	}

	for (i = 1; i <= n; i ++) {
		printf("%c %d %s\n", hfmCodeSet[i].data, hfmCodeSet[i].weight, hfmCodeSet[i].code);
	}
}

int main(int argc, char const *argv[])
{
	huffmanCode hfmCodeSet[N + 1];
	huffmanTreeNode hfmTree[M + 1];
	int n;

	// input the number of letters
	scanf("%d", &n);
	
	// input the letters
	n = initialize(hfmCodeSet, n);
	
	// build a huffman tree using arrays and calculate the huffman codes
	hfmCoding(hfmCodeSet, hfmTree, n);

	if (n == 1) {
		printf("%d\n", strlen(hfmCodeSet[1].code));
		return 0;
	}
	// output the total length of huffman Code
	printf("%d\n", totalLength(hfmCodeSet));
	// output the weight and huffman Code of corresponding char
	sortDisplayCode(hfmCodeSet, n);
	return 0;
}