【Biocode】产生三行的seq+01序列

代码说明：

sequence.txt与site.txt整合 如下图：

sequence.txt:

site.txt:

整理之后如下：

蛋白质序列中发生翻译后修饰的位置标记为“1”，其他的位置标记为“0”

代码如下：

#include<stdio.h>
#include<stdlib.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<time.h>
#include<math.h>
using namespace std;
char str[];
char s[];
int map[];
int main()
{
    //序列与位点结合 输出 一行序列 一行0 1 样本
    int i,j,k,l,m,n;
    //序列文件
    FILE *p=fopen("sequence.txt","r");
    //位置文件
    FILE *p1=fopen("site_after.txt","r");
    //输出文件
    freopen("3lines_sequence.txt","w",stdout);
    while(fgets(str,,p))
    {
        printf("%s",str);
        fgets(str,,p);
        l=strlen(str);
        printf("%s",str);
        fgets(s,,p1);
        char *pp=strtok(s," ");
        i=;
        memset(map,,sizeof(map));
        while(pp!=NULL)
        {
            k=atoi(pp);
        //    printf("**%d\n",k);
            map[k-]=;
            pp=strtok(NULL," ");
        }
        for(i=;i<l-;i++)
        {
            if(map[i]==)
            printf("");
            else
            printf("");
        }
        printf("\n");
    }
    return();
}

注：

①在初始化字符数组的时候，空间要大于所读取文件最长的长度，否则会出现篡位的问题。