Cody

# Problem 969. Genome Sequence 003: DNA Sequence with random positioned segments

Solution 320514

Submitted on 18 Sep 2013
This solution is locked. To view this solution, you need to provide a solution of the same size or smaller.

### Test Suite

Test Status Code Input and Output
1   Pass
tic L=1024; while 1 % Loop until create a Gtruth with no repeated 16 segments Gtruth=randi(4,1,L,'uint8')-1; Gseg=zeros(L-15,16,'uint8'); for j=1:L-15 Gseg(j,:)=Gtruth(j:j+15); end Gsegv=uint32(Gseg(:,1)); for j=2:16 Gsegv=4*Gsegv+uint32(Gseg(:,j)); end if size(unique(Gsegv),1)==L-15 break end end % while create valid Gtruth segs=Gtruth(1:48); % First segment is a required segs=[segs;Gtruth(L-47:end)]; % Last Segment is required ptr=0; while 1 ptr=ptr+randi(32); if ptr+48>L,break;end segs=[segs;Gtruth(ptr:ptr+47)]; end v=randperm(size(segs,1)); % Segs Scrambled Over Easy segs(:,:)=segs(v,:); Gout = Genome_003(segs); toc %assert(isequal(Gout,Gtruth)||isequal(Gout,fliplr(Gtruth))) assert( isequal(Gout,Gtruth) )

Elapsed time is 0.515770 seconds.

2   Fail
%% tic L=2000+randi(200); while 1 % Loop until create a Gtruth with no repeated 16 segments Gtruth=randi(4,1,L,'uint8')-1; Gseg=zeros(L-15,16,'uint8'); for j=1:L-15 Gseg(j,:)=Gtruth(j:j+15); end Gsegv=uint32(Gseg(:,1)); for j=2:16 Gsegv=4*Gsegv+uint32(Gseg(:,j)); end if size(unique(Gsegv),1)==L-15 break end end % while create valid Gtruth segs=Gtruth(1:48); % First segment is a required segs=[segs;Gtruth(L-47:end)]; % Last Segment is required ptr=0; while 1 ptr=ptr+randi(32); if ptr+48>L,break;end segs=[segs;Gtruth(ptr:ptr+47)]; end v=randperm(size(segs,1)); % Segs Scrambled Over Easy segs(:,:)=segs(v,:); Gout = Genome_003(segs); toc %assert(isequal(Gout,Gtruth)||isequal(Gout,fliplr(Gtruth))) assert( isequal(Gout,Gtruth) )

Error: Assertion failed.

3   Pass
%% tic L=4000+randi(200); while 1 % Loop until create a Gtruth with no repeated 16 segments Gtruth=randi(4,1,L,'uint8')-1; Gseg=zeros(L-15,16,'uint8'); for j=1:L-15 Gseg(j,:)=Gtruth(j:j+15); end Gsegv=uint32(Gseg(:,1)); for j=2:16 Gsegv=4*Gsegv+uint32(Gseg(:,j)); end if size(unique(Gsegv),1)==L-15 break end end % while create valid Gtruth segs=Gtruth(1:48); % First segment is a required segs=[segs;Gtruth(L-47:end)]; % Last Segment is required ptr=0; while 1 ptr=ptr+randi(32); if ptr+48>L,break;end segs=[segs;Gtruth(ptr:ptr+47)]; end v=randperm(size(segs,1)); % Segs Scrambled Over Easy segs(:,:)=segs(v,:); Gout = Genome_003(segs); toc %assert(isequal(Gout,Gtruth)||isequal(Gout,fliplr(Gtruth))) assert( isequal(Gout,Gtruth) )

Elapsed time is 1.530525 seconds.

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!