用C语言写一份判断一段DNA序列是否为CpG岛的程序，并且可以输入一段DNA序列，然后输出结果

在C语言中，我们可以编写一个简单的程序来判断一段DNA序列是否符合CpG岛的特性。CpG岛通常是指富含胞嘧啶(C)和鸟嘌呤(G)，并且C/G比值较高的区域。这里我们假设CpG岛的C/G比率大于某个阈值，例如0.6。

#include <stdio.h>
#include <string.h>

// 定义CpG岛的C/G比率阈值
#define CpG_THRESHOLD 0.6

// 检查CpG岛函数
int isCpGIsland(char *dna_sequence) {
    int c_count = 0;
    int g_count = 0;
    double c_g_ratio;

    // 遍历DNA序列
    for (int i = 0; i < strlen(dna_sequence); i += 2) { // DNA的碱基对步长通常是两个
        if (dna_sequence[i] == 'C' && dna_sequence[i + 1] == 'G') {
            c_count++;
            g_count++;
        }
    }

    // 计算C/G比率并检查是否超过阈值
    if (g_count != 0) {
        c_g_ratio = (double)c_count / g_count;
    } else {
        c_g_ratio = 0.0; // 如果没有G，视作C/G比为0
    }

    return (c_g_ratio > CpG_THRESHOLD) ? 1 : 0; // 1表示是CpG岛，0表示不是
}

int main() {
    char dna_input[100];
    printf("请输入DNA序列（只接受CG组合，用空格分隔）：");
    fgets(dna_input, sizeof(dna_input), stdin);

    // 移除末尾的换行符
    dna_input[strlen(dna_input) - 1] = '\0';

    // 调用函数并输出结果
    int result = isCpGIsland(dna_input);
    if (result) {
        printf("%s 是CpG岛\n", dna_input);
    } else {
        printf("%s 不是CpG岛\n", dna_input);
    }

    return 0;
}