7-3 哈夫曼编码 (30 分)

				7-3 哈夫曼编码 (30 分)

1. 给定一段文字，如果我们统计出字母出现的频率，是可以根据哈夫曼算法给出一套编码，使得用此编码压缩原文可以得到最短的编码总长。然而哈夫曼编码并不是唯一的。例如对字符串"aaaxuaxz"，容易得到字母
   ‘a’、‘x’、‘u’、‘z’ 的出现频率对应为 4、2、1、1。我们可以设计编码 {‘a’=0, ‘x’=10, ‘u’=110,
   ‘z’=111}，也可以用另一套 {‘a’=1, ‘x’=01, ‘u’=001, ‘z’=000}，还可以用 {‘a’=0,
   ‘x’=11, ‘u’=100, ‘z’=101}，三套编码都可以把原文压缩到 14 个字节。但是 {‘a’=0, ‘x’=01,
   ‘u’=011, ‘z’=001} 就不是哈夫曼编码，因为用这套编码压缩得到 00001011001001
   后，解码的结果不唯一，“aaaxuaxz” 和 “aazuaxax” 都可以对应解码的结果。本题就请你判断任一套编码是否哈夫曼编码。

2. 输入格式： 首先第一行给出一个正整数 N（2≤N≤63），随后第二行给出 N 个不重复的字符及其出现频率，格式如下：

3. c[1] f[1] c[2] f[2] … c[N] f[N] 其中c[i]是集合{‘0’ - ‘9’, ‘a’ - ‘z’,
   ‘A’ - ‘Z’, ‘_’}中的字符；f[i]是c[i]的出现频率，为不超过 1000 的整数。再下一行给出一个正整数
   M（≤1000），随后是 M 套待检的编码。每套编码占 N 行，格式为：

   c[i] code[i] 其中c[i]是第i个字符；code[i]是不超过63个’0’和’1’的非空字符串。

   输出格式： 对每套待检编码，如果是正确的哈夫曼编码，就在一行中输出"Yes"，否则输出"No"。

   注意：最优编码并不一定通过哈夫曼算法得到。任何能压缩到最优长度的前缀编码都应被判为正确。

   7
   A 1 B 1 C 1 D 3 E 3 F 6 G 6
   4
   A 00000
   B 00001
   C 0001
   D 001
   E 01
   F 10
   G 11
   A 01010
   B 01011
   C 0100
   D 011
   E 10
   F 11
   G 00
   A 000
   B 001
   C 010
   D 011
   E 100
   F 101
   G 110
   A 00000
   B 00001
   C 0001
   D 001
   E 00
   F 10
   G 11
   输出样例：
   Yes
   Yes
   No
   No
   

//思路：先求出哈夫曼编码的的最短加权路径和然后再判断是不是前缀编码两者缺一不可
//注意：判断是不是前缀编码时一定不要忘记等长且相等的！！！血的教训！！！
#include
#include
#include
#include
typedef int Status;
typedef struct HTNode{
    unsigned int weight;
    unsigned int parent;
    unsigned int lchild,rchild;
}HTNode,*HuffmanTree;
typedef char ** HuffmanCode;
Status Select(HuffmanTree &HT,int n,int &s1,int &s2);
void HuffmanTree_Create(HuffmanTree &HT, int *w,int n)
{
    int i,m,s1,s2;
    m = 2*n-1;
    HT = (HTNode*)malloc((m+1)*sizeof(HTNode));
    HTNode *p;
    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 <= m;++i,++p,++w){
        p->weight = 0;
        p->parent = 0;
        p->lchild = 0;
        p->rchild = 0;
    }
    for(i = n+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;
    }
}

void Huffman_coding(HuffmanTree HT,HuffmanCode &HC,int n)
{
    int i,c,f,start;
    char *cd;
    HC = (HuffmanCode)malloc((n+1)*sizeof(char*));
    cd = (char*)malloc(n*sizeof(char));
    cd[n-1] = '\0';
    for(i = 1;i <= n;++i){
        start = n-1;
        for(c = i,f = HT[i].parent;f != 0;c = f,f = HT[f].parent){
            if(HT[f].lchild == c)
            cd[--start] = '0';
            else cd[--start] = '1';
			}
            HC[i] = (char*)malloc((n-start)*sizeof(char));
            strcpy(HC[i],&cd[start]);

    }
    free(cd);
}

Status Select(HuffmanTree &HT,int n,int &s1,int &s2)
{
	int i,min1 = -1,min2 = -1,j,temp;
	for(i = 1;i <= n;++i){
		if(HT[i].parent == 0){
			min1 = i;
			break;
		}
	}
	for(j = i+1;j <= n;++j){
		if(HT[j].parent == 0){
			min2 = j;
			break;
		}
	}
	for(i = j+1;i <= n;++i){
		if(HT[i].parent == 0&&HT[i].weight<HT[min1].weight)min1 = i;
		else if(HT[i].parent == 0&&HT[i].weight < HT[min2].weight&&HT[i].weight >= HT[min1].weight)min2 = i;
	}
	s1 = min1;
	s2 = min2;
	return 0;
}
int main()
{
    int a[2000],i,n,m,k,j,m1,i1,min,sum2,sum1,l;
    char s[2000],s2[2000];
    int s1[2000];
    scanf("%d",&n);
    for(i=0;i<n;++i){
        getchar();
        scanf("%s",s+i);
        scanf("%d",a+i);
    }
    scanf("%d",&k);
    HuffmanTree HT;
    HuffmanCode HC;
    HuffmanTree_Create(HT,a,n);
	Huffman_coding(HT,HC,n);
    HuffmanCode HC1;
    HC1=(HuffmanCode)malloc((n+1)*sizeof(char*));
    for(j=0;j<k;++j){
    for(i=1;i<=n;++i){
        HC1[i]=(char*)malloc(n*sizeof(char));
        getchar();
        scanf("%s",s2+i);
        scanf("%s",HC1[i]);

    }
    min=1;
    for(l = 1;l <= n;++l){
        if(min == 0)break;
    for(i=l+1;i<=n;++i){
        m = strlen(HC1[l]);
        m1 = strlen(HC1[i]);
        if(m >= m1)//一定不要忘记等于的情况
            if(strncmp(HC1[i],HC1[l],m1) == 0){min = 0;break;}
        else if(m1>m)
            if(strncmp(HC1[l],HC1[i],m) == 0){min = 0;break;}
    }
    }
    if(min == 1)s1[j] = 1;
    else s1[j] = 0;
        sum1 = 0;
        sum2 = 0;
    for(i = 1;i <= n;++i){
        m = strlen(HC1[i]);
        m1 = strlen(HC[i]);
        sum1+=m1*a[i-1];
        sum2+=m*a[i-1];
    }
    if(sum1 == sum2);
    else s1[j] = 0;
}
    for(i = 0;i<k;++i){
        if(s1[i] == 1)printf("Yes\n");
        else printf("No\n");
    }
    return 0;
}

`注：俺第一次写如有不对的的地方还请赐教哈！