1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81
|
typedef struct
{
char nom[4];
int nombre;
} Codon;
void ExtractFirstCodon(char codon[], char chaine[])
{
strncpy(codon, chaine, 3);
codon[3] = '\0';
}
void InsertCodon(Codon listCodon[], char codon[], int position)
{
listCodon[position].nombre=1;
strcpy(listCodon[position].nom, codon);
}
int CompteCodons(char chaine[], Codon listCodon[])
{
int nb_codon = strlen(chaine)-2;
char codon[4];
int indiceCodon = 1;
for(int i=0; i<nb_codon; i++)
{
ExtractFirstCodon(codon, chaine);
if(i==0)
{
InsertCodon(listCodon, codon, 0);
chaine++;
continue;
}
int codonAbsent = 1;
for(int j=0; j<i; j++)
{
if(strcmp(listCodon[j].nom, codon)==0)
{
listCodon[j].nombre++;
codonAbsent=0;
}
}
if(codonAbsent)
{
InsertCodon(listCodon, codon, indiceCodon);
indiceCodon++;
}
chaine++;
}
return indiceCodon;
}
int Tri_occ(const void *cd1, const void *cd2)
{
int compar= ((Codon*)cd1)->nombre - ((Codon*)cd2)->nombre;
if (compar==0)
{
return strcmp(((Codon*)cd1)->nom , ((Codon*)cd2)->nom) ;
}
else
return -compar;
}
int main(void)
{
char chaine [100] = "TACGACGCGACTACGTCGTACG";
Codon listCodon[64];
int nbCodon = CompteCodons(chaine, listCodon);
qsort(listCodon, nbCodon, sizeof(Codon), Tri_occ);
for(int i=0; i<nbCodon; i++)
printf("%s --> %d\n", listCodon[i].nom, listCodon[i].nombre);
return 0;
} |
Partager