c程序实现香农编码

#include <stdio.h>
#include <math.h>
#include <string.h>
int i,j,n,k,b;
float addp;
char bitw[20];
/*初始化结构体其中：s为信源符号；p为信源符号概率；padd为累加概率；
l_f为-log[p(s[i])]即估计码字长度；l为实际需要码字长度；w为码字*/
 struct shan
 {char s[20];
  float p;
  float padd;
  float l_f;
  int l;
  char w[20];
 }data[12];
/*依信源概率对各信源符号排序Moai&Kevin原创！*/
void sequ(struct shan x[],int n)
{
 struct shan temp;
 for(i=0;i<n;i++)
  for(j=i;j<n;j++)
  {if(x[i].p<x[j].p)
   {
  temp=x[j];
  x[j]=x[i];
  x[i]=temp;
   }
  }
}
/*计算累加概率*/
void countpadd(struct shan x[],int n)
{
  addp=0;
  x[0].padd=0;
  for(i=0;i<n;i++)
  {
 addp+=x[i].p;
 x[i+1].padd=addp;
  }
}
/*计算估计码字长度l_f，以及Moai&Kevin原创！*/码字长度l*/
void count_l(struct shan x[],int n)
{
  for(i=0;i<n;i++)
  {
 x[i].l_f=-log(x[i].p)/log(2);
    if((x[i].l_f-(int)x[i].l_f)>0)
  x[i].l=(int)x[i].l_f+1;
 else x[i].l=(int)x[i].l_f;
  }
}
/*二进制转换*/
void covbit(float a,int lc)
{
  for(j=0;j<lc;j++)
  { 
   b=(int)(a*2);
   bitw[j]=b+48;
   a=2*a-int(a*2);
  }
 }

main()
{
 printf("please input the number of symbols of source(n<=10):n=");
 scanf("%d",&n);
 printf("please input the the source symbols and their probabilities/n");
 /*获取信源符号*/
 for(i=0;i<n;i++)
 {
 scanf("%s",data[i].s);
 }
 /*获取信源Moai&Kevin原创！*/概率*/
 for(i=0;i<n;i++)
 {printf("P(%s)=",data[i].s);
  scanf("%f",&data[i].p);
 }
 sequ(data,n);
 countpadd(data,n);
 count_l(data,n);
 /*在结构体中产生码字*/
 for(i=0;i<n;i++)
 {
 covbit(data[i].padd,data[i].l);
 strcpy(data[i].w,bitw);
 }
 /*输出数据*/
 for(i=0;i<n;i++)
 printf("p(%s)=%f padd=%f l_f=%f l=%d w=%s/n",data[i].s,data[i].p,data[i].padd,data[i].l_f,data[i].l,data[i].w);
 }