2016—9—24

var
        a:array[0..10000,1..50]of longint;
        f:array[0..10000,1..2]of longint;
        l,n,i,j,k:longint;
        c:char;
procedure dg(l,r:longint);//快排
var
        i,t,j,m,n:longint;
begin
        i:=l;
        j:=r;
        m:=f[i,1];
        n:=f[i,2];
        repeat
                while (f[i,1]<m) or (f[i,1]=m) and (f[i,2]<n) do
                        inc(i);
                while (f[j,1]>m) or (f[j,1]=m) and (f[j,2]>n) do
                        dec(j);
                if i<=j then
                begin
                        f[0]:=f[i];
                        f[i]:=f[j];
                        f[j]:=f[0];
                        inc(i);
                        dec(j);
                end;
        until i>j;
        if(j>l) then
                dg(l,j);
        if(i<r) then
                dg(i,r);
end;
begin
        assign(input,'dna.in');//文件输入输出
        assign(output,'dna.out');
        reset(input);
        rewrite(output);
        readln(l,n);
        for i:=1 to n do
        begin
                for j:=1 to l do
                begin
                        read(c);//把读入的换为字符
                        if c='A' then
                                a[i,j]:=1;
                        if c='C' then
                                a[i,j]:=2;
                        if c='G' then
                                a[i,j]:=3;
                        if c='T' then
                                a[i,j]:=4;
                end;
                readln;
        end;
        for i:=1 to n do
                f[i,2]:=i;//存原来的位置
        for i:=1 to n do
        begin
                for j:=2 to l do
                begin
                        for k:=1 to j-1 do
                        begin
                                if a[i,j]<a[i,k] then//看每行有多少逆序对
                                        inc(f[i,1]);
                        end;
                end;
        end;
        dg(1,n);//快排一遍逆序对
        for i:=1 to n do
        begin
                for j:=1 to l do
                begin
                        if a[f[i,2],j]=1 then//把原来转成数字的转回字符并输出
                                write('A');
                        if a[f[i,2],j]=2 then
                                write('C');
                        if a[f[i,2],j]=3 then
                                write('G');
                        if a[f[i,2],j]=4 then
                                write('T');
                end;
                writeln;
        end;
        close(input);
        close(output);
end.

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值