C与2D char数组的连接

时间:2012-12-09 01:24:59

标签: c char concatenation dna-sequence

我正在逐行读取文本文件到2D数组中。我想连接char数组,所以我有一个长char数组。我遇到了麻烦,我可以使用两个char数组,但是当我尝试做很多时我就错了。

目前char数组看起来像这样:

AGCTTTTCATTC

我想得到这样的东西:

AGCTTTTCATTCAGCTTTTCATTC

我已经包含了一些代码。

int counter = 0; 
fid = fopen("dna.fna","r");
while(fgets(line, sizeof(line), fid) != NULL && counter!=66283 ) {
    if (strlen(line)==70) {
        strcpy(dna[counter], line);        
    counter++;
    }
}
int dnaSize = 6628; 
//Concatenating the DNA into a single char array.
int i;
char DNA[dnaSize];
for(i = 0; i<66283;i++){
   strcpy(DNA[i],dna[i]);
   strcat(DNA[i+1],dna[i+1]);
}

1 个答案:

答案 0 :(得分:1)

您只需要循环到< counter 那么,你是在复制还是连接?你只需要做一个或另一个。

我建议在循环中使用strcat,但初始化DNA。

char DNA[dnaSize] = ""; //initalise so safe to pass to strcat
for(i = 0; i<counter;i++)
{
   strcat(DNA,dna[i]); //no need for indexer to DNA
}

此外,您需要考虑两个阵列的大小。我相信(希望)dna是一个char数组的数组。如果是的话,我猜它只有66283长的第一个维度。因此它不适合DNA6628 long),即使每行的长度为1个字符。

以下是关于如何准确分配适量内存的想法:

#define MAXLINELENGTH (70)
#define MAXDNALINES (66283)

//don't copy this line, it will not work because of the sizes involved (over 4MB)
//it will likely stack overflow
//just do what you are currently doing as long as it's a 2-d array.
char dna[MAXDNALINES][MAXLINELENGTH + 1];

int counter = 0; 
int totalSize = 0;
fid = fopen("dna.fna","r");
while(fgets(line, sizeof(line), fid) != NULL && counter!=MAXDNALINES ) {
    const int lineLength = strlen(line);
    if (lineLength==MAXLINELENGTH) {
        strcpy(dna[counter], line);        
        counter++;
        totalSize += lineLength;
    }
}

//Concatenating the DNA into a single char array (of exactly the right length)
int i;
char *DNA = malloc(totalSize+1); // the + 1 is for the final null, and putting on heap so don't SO
DNA[0] = '\0'; //the initial null is so that the first strcat works
for(i = 0; i<counter;i++){
   strcat(DNA,dna[i]);
}

//do work with DNA here

//finally free it  
free(DNA);