2019-08-09

环状序列(紫书例题3-6)

问题描述:
长度为n的环状串有n种表示法,分别为某个位置开始顺时针得到。CGAGTCAGCT,GAGTCAGCTC,AGTCAGCTCG等。在这些表示法中,字典序最小的称为“最小表示”。输入一个长度为n(n<=100)的环状DNA串(只包含A、C、G、T这4种字符)的一种表示法,你的任务是输出该环状串的最小表示。例如,CTCC的最小表示是CCCT,CGAGTCAGCT的最小表示为AGCTCGAGTC。


image.png

样例输入:
2
CGAGTCAGCT
CTCC

样例输出:
AGCTCGAGTC
CCCT

分析:

一开始理解这道题的时候很吃力,看了好久都不知道它想要干嘛,后面才知道原来是在一串成环的 DNA像字典排序一样找出其中最小的表示,也就是要按照字母排列寻找正确的起点。

难点:

1.当找到的起点并非是输入的DNA序列的第一个字母,怎么让“断掉“”的重新“续上”:
(起点s+增值i)%输入DNA的长度

2.当符合/不符合时,应怎么递增继续判断/切换起点:
符合时,记录当前i的值作为起点,用此起点继续跟该点以后做起点的序列作比较;不符合时,直接跳过;直到到达该DNA的结尾。

判断以不同起点的序列的字典序,从起点开始比较,当增加相同单位处的字母不同时,字母小的序列就小,相同时继续判断两序列的下一个位置字母的大小关系,直到回到原起点。

代码如下:
#include<iostream>
#include<string> 
using namespace std; 
int comp(string s,int s1,int s2){
    int num=s.size();
    for(int i=0;i<num;i++){
        if(s[(s1+i)%num]!=s[(s2+i)%num])
        return s[(s1+i)%num]>s[(s2+i)%num];
    }
    return 0;
}
int main(){
    int n;
    string dna;
    cin>>n;
    while(n--){
        cin>>dna;
        int start=0,num=dna.size();
        for(int i=1;i<num;i++){
            if(comp(dna,start,i))start=i;
        } 
        for(int i=0;i<num;i++){
            cout<<dna[(start+i)%num];
        }
        cout<<endl;
    } 
} 
运行截图

相似的,还有:

DNA序列(紫书习题3-7)

题目描述:

输入m个长度均为n的DNA序列,求一个DNA序列,到所有序列的总Hamming距离尽量小。两个等长字符串的Hamming距离等于字符不同的位置个数,例如,ACGT和GCGA的Hamming距离为2(左数第1,4字符不同)。
输入整数m和n(4<=m<=50,4<=n<=1000),以及m个长度为n的DNA序列(只是包含字母A,C,G,T),输出到m个序列的Hamming距离和最小的DNA序列和对应的距离。如有多解,要求为字典序最小的解。例如,对于下面5个DNA序列,最优解为TAAGATAC。


image.png
分析:

在上面那题对字典序有了了解以后,理解这道题就很简单了。所谓Hamming距离和最小意思就是跟所给出的所有序列的字符不同的位置个数总和最小,直接寻找在所有序列相同位置中出现次数最多的字符就可以找出这个序列了。又由于规定了多个解的时候,题目要求字典序最小,因此判断时可以用if-else结构按照字典顺序进行判断。

代码如下:

#include<iostream>
#include<string>
using namespace std;
int main(){
    int m,n,max,a,c,g,t;
    string dna;
    cin>>m>>n;
    char str[m][n];
    for(int i=0;i<m;++i){
        cin>>str[i];
    } 
    for(int i=0;i<n;++i){
        a=c=t=g=0;
        for(int j=0;j<m;++j){
            if(str[j][i]=='A')a++;
            else if(str[j][i]=='C')c++;
            else if(str[j][i]=='G')g++;
            else t++;
        }
        max=a;
        if(c>max)max=c;
        if(g>max)max=g;
        if(t>max)max=t;
        
        if(max==a)dna+='A';
        else if(max==c)dna+='C';
        else if(max==g)dna+='G';
        else if(max==t)dna+='T';
    }
    cout<<"\n"<<dna<<endl; 
}
运行截图
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,362评论 5 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,330评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,247评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,560评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,580评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,569评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,929评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,587评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,840评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,596评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,678评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,366评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,945评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,929评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,165评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 43,271评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,403评论 2 342

推荐阅读更多精彩内容

  • 1.dict字典 定义 字典 是一系列键—值对 。每个键 都与一个值相关联,你可以使用键来访问与之相关联的值。 与...
    棒棒糖117阅读 165评论 0 0
  • ElasticSearch 什么是ElasticSearch,根据其中文文档可以简单定义为以下三点,与传统的关系型...
    冰阔乐_b0a4阅读 107评论 0 0
  • 一、内置类属性 什么是内置类属性声明类的时候系统自动添加的属性(可能是字段也可能是对象属性) 定制当前类的对象的打...
    Wa_ngli阅读 235评论 0 0
  • 听过传染病会传染,听过自控力也会传染吗 根据大量数据跟踪发现如果你身边有某个人的体重增加了,那么你的体重增加的概率...
    baiying阅读 248评论 0 0
  • “美与悲剧向来不冲突,甚至悲剧中更能突出人性之美” 《边城》的故事本身是一个美丽的悲剧。沈从文通过自己优美的语言,...
    CHENYIFEI_阅读 349评论 1 3