Skip to content

Latest commit

ย 

History

History
110 lines (76 loc) ยท 4.15 KB

2019-08-13-AT-GenomicRangeQuery.md

File metadata and controls

110 lines (76 loc) ยท 4.15 KB

2019๋…„ 8์›” 13์ผ

Codility - GenomicRangeQuery (Prefix Sums) {docsify-ignore-all}

๏ฟฝ

์ถœ์ฒ˜: https://app.codility.com/programmers/lessons/5-prefix_sums/genomic_range_query/

๋ฌธ์ œ

A DNA sequence can be represented as a string consisting of the letters A, C, G and T, which correspond to the types of successive nucleotides in the sequence. Each nucleotide has an impact factor, which is an integer. Nucleotides of types A, C, G and T have impact factors of 1, 2, 3 and 4, respectively. You are going to answer several queries of the form: What is the minimal impact factor of nucleotides contained in a particular part of the given DNA sequence?

The DNA sequence is given as a non-empty string S = S[0]S[1]...S[N-1] consisting of N characters. There are M queries, which are given in non-empty arrays P and Q, each consisting of M integers. The K-th query (0 โ‰ค K < M) requires you to find the minimal impact factor of nucleotides contained in the DNA sequence between positions P[K] and Q[K] (inclusive).

For example, consider string S = CAGCCTA and arrays P, Q such that:

P[0] = 2    Q[0] = 4
P[1] = 5    Q[1] = 5
P[2] = 0    Q[2] = 6

The answers to these M = 3 queries are as follows:

  • The part of the DNA between positions 2 and 4 contains nucleotides G and C (twice), whose impact factors are 3 and 2 respectively, so the answer is 2.
  • The part between positions 5 and 5 contains a single nucleotide T, whose impact factor is 4, so the answer is 4.
  • The part between positions 0 and 6 (the whole string) contains all nucleotides, in particular nucleotide A whose impact factor is 1, so the answer is 1. Write a function:

def solution(S, P, Q)

that, given a non-empty string S consisting of N characters and two non-empty arrays P and Q consisting of M integers, returns an array consisting of M integers specifying the consecutive answers to all queries.

Result array should be returned as an array of integers.

For example, given the string S = CAGCCTA and arrays P, Q such that:

P[0] = 2    Q[0] = 4
P[1] = 5    Q[1] = 5
P[2] = 0    Q[2] = 6

the function should return the values [2, 4, 1], as explained above.

Write an efficient algorithm for the following assumptions:

  • N is an integer within the range [1..100,000];
  • M is an integer within the range [1..50,000];
  • each element of arrays P, Q is an integer within the range [0..N โˆ’ 1];
  • P[K] โ‰ค Q[K], where 0 โ‰ค K < M;
  • string S consists only of upper-case English letters A, C, G, T.

๋‚ด ํ’€์ด1 (์‹œ๊ฐ„ ์ดˆ๊ณผ)

์‹œ๊ฐ„๋ณต์žก๋„ O(N^2)

  • ์ •ํ™•๋„๋Š” 100%์ง€๋งŒ ์‹œ๊ฐ„๋ณต์žก๋„๊ฐ€ ๋†’์•„ ์„ฑ๋Šฅ ๊ฒ€์‚ฌ๋ฅผ ํ†ต๊ณผํ•˜์ง€ ๋ชป ํ•œ๋‹ค.
def genomicRangeQuery(S,P,Q):
    nucleotide = [0]*len(S)

    for i in range(len(S)):
        if S[i] =='A': nucleotide[i] = 1
        elif S[i] =='C': nucleotide[i] = 2
        elif S[i] =='G': nucleotide[i] = 3
        else: nucleotide[i] = 4

    z = zip(P,Q)
    answer =[]
    for s,e in z:
        answer.append(min(nucleotide[s:e+1]))
    return answer

์ ‘๊ทผ ๋ฐฉ๋ฒ•

  • ๋ฌธ์ž์—ด S ์˜ ํ•œ ๋ฌธ์ž๋ฅผ ์ฝ๋Š” ๊ณผ์ •์„ ํ•œ ๋‹จ๊ณ„๋ผ๊ณ  ํ•  ๋•Œ, ๊ฐ ๋‹จ๊ณ„๋งˆ๋‹ค A,C,G,T์˜ ๊ฐœ์ˆ˜๋ฅผ ๋ฐฐ์—ด์— ์ €์žฅํ•œ๋‹ค.

  • P,Q์˜ ๋ฒ”์œ„์— ์žˆ๋Š” ๋‹จ๊ณ„์—์„œ A,C,G,T์˜ ๊ฐœ์ˆ˜๋ฅผ ๋น„๊ตํ•œ ๋’ค, ๊ฐœ์ˆ˜์— ๋ณ€ํ™”๊ฐ€ ์žˆ๋Š” A,C,G,T ์ค‘ minimal impact ๊ฐ’์„ ๊ฒฐ๊ณผ(result)๋ฐฐ์—ด์— ๋‹ด๋Š”๋‹ค.

๋‚ด ํ’€์ด2

def genomicRangeQuery(S,P,Q):
    nucleotide = {'A': 1, 'C': 2, 'G': 3, 'T': 4}

    curr_counts = [0, 0, 0, 0]
    counts = [curr_counts[:]]

    for s in S:
        curr_counts[nucleotide[s] - 1] += 1
        counts.append(curr_counts[:])

    result = []
    for i in range(len(P)):
        count_p = counts[P[i]]
        count_q = counts[Q[i]+1]

        if P[i] == Q[i]: result.append(nucleotide[S[P[i]]])
        elif count_p[0] < count_q[0]:
            result.append(1)
        elif count_p[1] < count_q[1]:
            result.append(2)
        elif count_p[2] < count_q[2]:
            result.append(3)
        elif count_p[3] < count_q[3]:
            result.append(4)
    return result

๋Š๋‚€์ 

  • ์–ด๋ ค์šด ๋ฌธ์ œ์˜€๋‹ค. ์‹œ๊ฐ„๋ณต์žก๋„๋ฅผ ์ค„์ด๋Š” ๊ณผ์ •์ด ์–ด๋ ค์›Œ์„œ stackoverflow ๋“ฑ์„ ์ฐธ๊ณ ํ•ด์„œ ํ’€์—ˆ๋‹ค.