typedef struct
 {
   unsigned int weight;
   unsigned int parent,lchild,rchild;
 }HTNode,*HuffmanTree; /* 動態分配數組存儲赫夫曼樹 */
 typedef char **HuffmanCode; /* 動態分配數組存儲赫夫曼編碼表 */
#include<string.h>
 #include<ctype.h>
 #include<malloc.h> /* malloc()等 */
 #include<limits.h> /* INT_MAX等 */
 #include<stdio.h> /* EOF(=^Z或F6),NULL */
 #include<stdlib.h> /* atoi() */
 #include<io.h> /* eof() */
 #include<math.h> /* floor(),ceil(),abs() */
 #include<process.h> /* exit() */
 /* 函數結果狀態代碼 */
 #define TRUE 1
 #define FALSE 0
 #define OK 1
 #define ERROR 0
 #define INFEASIBLE -1
 /* #define OVERFLOW -2 因為在math.h中已定義OVERFLOW的值為3,故去掉此行 */
 typedef int Status; /* Status是函數的類型,其值是函數結果狀態代碼，如OK等 */
 typedef int Boolean; /* Boolean是布爾類型,其值是TRUE或FALSE */


int min1(HuffmanTree t,int i)
 { /* 函數void select()調用 */
   int j,flag;
   unsigned int k=UINT_MAX; /* 取k為不小于可能的值 */
   for(j=1;j<=i;j++)
     if(t[j].weight<k&&t[j].parent==0)
       k=t[j].weight,flag=j;
   t[flag].parent=1;
   return flag;
 }

 void select(HuffmanTree t,int i,int *s1,int *s2)
 { /* s1為最小的兩個值中序號小的那個 */
   int j;
   *s1=min1(t,i);
   *s2=min1(t,i);
   if(*s1>*s2)
   {
     j=*s1;
     *s1=*s2;
     *s2=j;
   }
 }

 void HuffmanCoding(HuffmanTree *HT,HuffmanCode *HC,int *w,int n) /* 算法6.12 */
 { /* w存放n個字符的權值(均>0),構造赫夫曼樹HT,并求出n個字符的赫夫曼編碼HC */
   int m,i,s1,s2,start;
   unsigned c,f;
   HuffmanTree p;
   char *cd;
   if(n<=1)
     return;
   m=2*n-1;
   *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<=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=(*HT)[s2].parent=i;
     (*HT)[i].lchild=s1;
     (*HT)[i].rchild=s2;
     (*HT)[i].weight=(*HT)[s1].weight+(*HT)[s2].weight;
   }
   /* 從葉子到根逆向求每個字符的赫夫曼編碼 */
   *HC=(HuffmanCode)malloc((n+1)*sizeof(char*));
   /* 分配n個字符編碼的頭指針向量([0]不用) */
   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));
     /* 為第i個字符編碼分配空間 */
     strcpy((*HC)[i],&cd[start]); /* 從cd復制編碼(串)到HC */
   }
   free(cd); /* 釋放工作空間 */
 }

 void main()
 {
   HuffmanTree HT;
   HuffmanCode HC;
   int *w,n,i;
   printf("請輸入權值的個數(>1)：");
   scanf("%d",&n);
   w=(int*)malloc(n*sizeof(int));
   printf("請依次輸入%d個權值(整型)：\n",n);
   for(i=0;i<=n-1;i++)
     scanf("%d",w+i);
   HuffmanCoding(&HT,&HC,w,n);
   for(i=1;i<=n;i++)
     puts(HC[i]);
 }