C语言DNA互补序列生成函数非0位置失效问题排查

现有代码的核心问题 #

• 目标数组未初始化：compliment_DNA_seq是栈上的局部数组，创建后没有初始化。当pos不为0时，数组前pos个位置是随机垃圾值，输出时会显示乱码，让你误以为函数没工作。
• 硬编码ASCII值可读性差：用65、84这类数字代表'A'、'T'，不仅难读，还容易出错，直接用字符常量更清晰。
• 无边界检查：调用函数时硬传length=22，如果输入的DNA序列长度小于pos+22，会越界访问原序列数组，引发未定义行为；同时也没判断pos是否在合法范围内（比如是否小于序列长度）。
• 变量拼写错误：lenght应为length，属于不良编码习惯。

修复后的代码 #

#include <cs50.h>
#include <stdio.h>
#include <string.h>
#include <strings.h>
#include <stdlib.h>
#include <ctype.h>

void compliment(string seq, string comp_seq, int pos, int length);

int main(void)
{
    string DNA_seq = get_string("Input: ");
    int seq_len = strlen(DNA_seq);
    char compliment_DNA_seq[seq_len + 1];
    
    // 先复制原序列到目标数组，未修改的位置保留原序列
    strcpy(compliment_DNA_seq, DNA_seq);
    
    // 示例：从位置1开始，修改2个字符的互补序列
    int target_pos = 1;
    int target_length = 2;
    compliment(DNA_seq, compliment_DNA_seq, target_pos, target_length);

    printf("\n%s\n%s\n\n", DNA_seq, compliment_DNA_seq);
}

void compliment(string seq, string comp_seq, int pos, int length)
{
    int seq_len = strlen(seq);
    // 边界检查：确保pos和length参数合法
    if (pos < 0 || pos >= seq_len || pos + length > seq_len) {
        printf("Invalid position or length!\n");
        return;
    }
    
    int j = length;
    for (int i = pos; j > 0; i++, j--)
    {
        // 处理大小写兼容，同时用switch提升可读性
        switch(toupper(seq[i])) {
            case 'A':
                comp_seq[i] = 'T';
                break;
            case 'C':
                comp_seq[i] = 'G';
                break;
            case 'G':
                comp_seq[i] = 'C';
                break;
            case 'T':
                comp_seq[i] = 'A';
                break;
            default:
                // 遇到非ATGC字符时保留原字符并提示
                comp_seq[i] = seq[i];
                printf("Warning: Invalid DNA character at position %d\n", i);
                break;
        }
    }
    // 确保字符串以'\0'结尾
    comp_seq[seq_len] = '\0';
}

修复说明 #

1. 先把原序列复制到目标数组，未指定修改的位置会保留原序列，避免乱码。
2. 添加了边界检查，防止越界访问数组。
3. 用字符常量替换ASCII数字，同时兼容大小写输入。
4. 修正了变量名拼写错误，提升代码可维护性。

内容的提问来源于stack exchange，提问作者Meet Patel