Алгоритм Шеннона-Фано

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
 
int main(int argc, char *argv[])
{
    unsigned int c[256], mal2; // массив c[256] - количество повторений байта
    char filename[50]; // массив для имени файла
    FILE* f; 
    unsigned char b; 
    int i=0,j=0; 
    int n[256],mal1; // массив n[256] - массив под индексы байтов
    double p[256], mal; // массив p[256] - массив вероятностей
    double s=0; 
    double L=1/log(2); // необходимая переменная для подсчёта ср.кол-ва инф.
    double I=0; // обнуляем значение среднего количества информации
    printf("Input a filename:\n"); 
    scanf ("%s",&filename); // ввод имени файла
    
    for(i=0;i<256;i++) // обнуляем элементы массива количества повторений байтов по циклу
    {c[i]=0;}
    f = fopen(filename,"rb"); // открываем файл в бинарном режиме чтения
 
 /*Считаем количество повторений*/   
    do
{
    fread (&b,1,1, f); 
    c[b]++; 
}
    while (!feof(f));
    c[b]--;
    /*Считаем вероятности повторений байтов*/
    for(i=0; i<256;i++)
    {s=s+c[i];}
    for(i=0; i<256;i++)
    {p[i]=c[i]/s;}
   for(i=0; i<256;i++)
    {
    if (p[i]!=0) 
    {I = I+p[i]*log(1/(p[i]))*L;} // считаем среднее кол-во информации по формуле
    }  
    fclose(f); 
    
    /*Сортируем данные по невозрастанию пузырьковой сортировкой*/    
    for (i=0;i<256;i++) //Данный цикл нужен для того, чтобы у каждого значения байта оставались свои значения вероятности и повторений после сортировки
          { n[i]=i; }  // иначе значения повторения и вероятности у байтов поменяются при выводе результатов, так как индексы остались на своих местах
    for (i=255;i>=0;i--) // сама сортировка
    {
        for (j=i;j<255;j++)
        {
            if (p[j]<p[j+1])
            {
                            mal=p[j];
                            p[j]=p[j+1];
                            p[j+1]=mal;
                            mal1=n[j];
                            n[j]=n[j+1];
                            n[j+1]=mal1;
                            mal2=c[j];
                            c[j]=c[j+1];
                            c[j+1]=mal2;                         
                      }}}
 for (i=0; i<256; i++)   // вывод индексов байтов, кол-ва повторений и вероятности повторений по циклу
{ printf("indeks=%d, kolichestvo=%d, veroyatnost=%f\n", n[i],c[i],p[i]); }
printf("Kol-vo informacii= %lf\n",I); //вывод среднего количества информации    
    system("PAUSE");
    return EXIT_SUCCESS;
}