#include
 #include
 #include
 
void dna8_val2str( unsigned long val, char str[8] )
 {
     assert( val < (1ul<<(8*2)) );
 
    for( size_t i=0; i<8; ++i )
     {
         str[7-i] = "agct"[val%4];
         val >>= 2;
     }
 }
 
int main()
 {
     unsigned long words = 0;
     unsigned long numbers[65536] = { 0 };
 
    // 处理
     FILE* fin = fopen( "sequence.txt", "r" );
     if( !fin )
         return 1;
     int bav = 0;
     unsigned long val = 0;
     for( int c; c=fgetc(fin), c!=EOF; )
     {
         switch( c )
         {
         case '@': // 遇到@则结束
             break;
         case '\r':
         //case '\n':
             ++words;
             continue;
         case 'a':
         case 'g':
         case 'c':
         case 't':
             val = ((val<<2)&0xFFFF) | (c%36%5);
             if(bav<8) // 不足8个有效字符时先等等
                 ++bav;
             else
                 ++numbers[val];
             break;
         default: // 出现了agct之外的字符
             long pos = ftell(fin)-1;
             if( c>0x20 && c<0xFF ) // 可显示的字符，就显示其本身
                 printf( "--- 0x%08lX处出现非法字符\'%c\'\n", pos, (char)c );
             else // 不可显示的字符，就显示其对应的ASCII值
                 printf( "--- 0x%08lX处出现非法字符0x%02hhX\n", pos, (char)c );
         }
     }
     fclose( fin );
 
    // 输出
     FILE* fout = fopen( "result.txt", "w" );
     if( !fout )
         return 3;
     for( int i=0; i<65536; ++i )
     {
         // if( numbers[i] !=0 )
         {
             char str[8];
             dna8_val2str( i, str );
             fprintf( fout, "%.8s\t%ld\t%f\n", str, numbers[i], numbers[i]/65536.0 );
         }
     }
     fclose( fout );
     printf( "处理完毕\n" );
 
    return 0;
 }
程序没有问题，现在这个是统计8个字符的问题。我想改成统计四个字符的程序，自己改了半天也没成功。你可以先运行一下附件，看结果，然后就明白我说的意思了。

你看不懂原来的程序   自然就无法更改

此贴也可定义为求代码贴

那你帮我看一下吧，大兄弟！万分感谢！！！

#include
#include
#include 

头文件倒是可以修改修改...

犀利，残忍。。。

大兄弟，你帮我修改一下吧，我是学生物的很少用到编程，现在研三了，马上要毕业，却需要统计基因频数，我自己不太会编程，现学时间不够了，明年3月份论文都得写出来，现在却遇到了这个问题，出不来结果。万分感谢！！！

标头，我写的漏了，是这个。
#include <stdio.h>
#include <string.h>
 #include <assert.h>
 
void dna8_val2str( unsigned long val, char str[8] )
 {
     assert( val < (1ul<<(8*2)) );
 
    for( size_t i=0; i<8; ++i )
     {
         str[7-i] = "agct"[val%4];
         val >>= 2;
     }
 }
 
int main()
 {
     unsigned long words = 0;
     unsigned long numbers[65536] = { 0 };
 
    // 处理
     FILE* fin = fopen( "sequence.txt", "r" );
     if( !fin )
         return 1;
     int bav = 0;
     unsigned long val = 0;
     for( int c; c=fgetc(fin), c!=EOF; )
     {
         switch( c )
         {
         case '@': // 遇到@则结束
             break;
         case '\r':
         //case '\n':
             ++words;
             continue;
         case 'a':
         case 'g':
         case 'c':
         case 't':
             val = ((val<<2)&0xFFFF) | (c%36%5);
             if(bav<8) // 不足8个有效字符时先等等
                 ++bav;
             else
                 ++numbers[val];
             break;
         default: // 出现了agct之外的字符
             long pos = ftell(fin)-1;
             if( c>0x20 && c<0xFF ) // 可显示的字符，就显示其本身
                 printf( "--- 0x%08lX处出现非法字符\'%c\'\n", pos, (char)c );
             else // 不可显示的字符，就显示其对应的ASCII值
                 printf( "--- 0x%08lX处出现非法字符0x%02hhX\n", pos, (char)c );
         }
     }
     fclose( fin );
 
    // 输出
     FILE* fout = fopen( "result.txt", "w" );
     if( !fout )
         return 3;
     for( int i=0; i<65536; ++i )
     {
         // if( numbers[i] !=0 )
         {
             char str[8];
             dna8_val2str( i, str );
             fprintf( fout, "%.8s\t%ld\t%f\n", str, numbers[i], numbers[i]/65536.0 );
         }
     }
     fclose( fout );
     printf( "处理完毕\n" );
 
    return 0;
 }

我擦，这么严重啊，亚历山大啊...
好吧，那我就冒着被老师罚站的风险来改一改吧...
其实我本来就不想做作业...

这个程序来统计8个字符完全能用，运行一下我上传的附件就能看到结果，现在我要统计4个字符。统计8 个字符时，4的8次方=65536，是这么来的。当统计4个字符时，是4的4次方=256。同理4的6次方=4096。现在只要统计4个字符就行了。

这程序在我这儿有错误啊，你用的是哪个编辑器啊？