Mathematics:DNA Sorting(POJ 1007)

　　　　　　　　　　　　　　　　

　　　　　　　　　　　　　　　　DNA排序

　　题目大意：给定多个ACGT序列，按照字母顺序算出逆序数，按逆序数从小到大排列

　　这题其实很简单，我们只要用一个归并排序算逆序数，然后快排就可以了（插入排序也可以，数据量不大），但是要注意的是他要保持原来的顺序，fcmop那里要写当a.inverse=b.inverse时，返回0

　　

 #include <iostream>
 #include <functional>
 #include <algorithm>

 using namespace std;
 typedef int Position;
 typedef struct _set
 {
     char input[];
     int inverse;
 }Set;

 void Merge_Sort(Position, Position, int *const);
 void Merge(Position, Position, int *const);

 static Set input[];
 static char A[], Tmp_Array[];

 int fcomp(const void *a, const void *b)
 {
     if ((*(Set *)a).inverse - (*(Set *)b).inverse == )
         return ;
     return (*(Set *)a).inverse - (*(Set *)b).inverse;
 }

 int main(void)
 {
     int line_sum, length, inverse_num;

     while (~scanf("%d%d", &length, &line_sum))
     {
         getchar();
         for (int i = ; i < line_sum; i++)
         {
             scanf("%s", input[i].input);
             strcpy(A, input[i].input);
             inverse_num = ;
             Merge_Sort(, length - , &inverse_num);
             input[i].inverse = inverse_num;
         }
         qsort(input, line_sum, sizeof(Set), fcomp);
         for (int i = ; i < line_sum; i++)
             printf("%s\n", input[i].input);
     }
     return ;
 }

 void Merge_Sort(Position Left, Position Right, int *const Inverse_Num)
 {
     if (Left < Right)
     {
         Position Mid = (Left + Right) / ;
         Merge_Sort(Left, Mid, Inverse_Num);
         Merge_Sort(Mid + , Right, Inverse_Num);
         Merge(Left, Right, Inverse_Num);
     }
 }

 void Merge(Position Left, Position Right, int *const Inverse_Num)
 {
     Position Mid = (Left + Right) / , lpos = Left, rpos = Mid + , pos = Left;

     for (; lpos <= Mid && rpos <= Right;)
     {
         if (A[lpos] <= A[rpos])
             Tmp_Array[pos++] = A[lpos++];
         else
         {
             Tmp_Array[pos++] = A[rpos++];
             (*Inverse_Num) += (Mid +  - lpos);
         }
     }
     while (lpos <= Mid)
         Tmp_Array[pos++] = A[lpos++];
     while (rpos <= Right)
         Tmp_Array[pos++] = A[rpos++];
     for (int i = Left; i <= Right; i++)
         A[i] = Tmp_Array[i];
 }