HDU 2457 DNA repair AC 로봇 + 동적 계획
바로 현재 n(1<=n<=50)개의 DNA 서열이 질병이 있는데, 현재 하나의 DNA 서열(길이가 1000을 넘지 않음)이 있는데, 그 중의 문자를 교체해야 한다(A, C, G, T로만 바꿀 수 있다). 이 DNA 서열은 질병 서열을 포함하지 않는다. 최소한 몇 글자를 교체해야 하는지 묻는다. 이런 해가 존재하지 않으면 출력-1
대략적인 사고방식:
간단한 AC 자동기 + 동적 기획으로 어려운 점이 없어요...
상태 이동 방정식 설명은 코드를 보십시오.
Result : Accepted Memory : 4328 KB Time : 62 ms
/*
* Author: Gatevin
* Created Time: 2014/11/22 15:00:55
* File Name: Kagome.cpp
*/
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;
int n;
const int inf = 9999999;
int dp[1010][1010];
/*
* dp[i][j] i Trie j
* dp[i][j]
* dp[i + 1][next[j][t]] = min(dp[i + 1][next[j][t]], dp[i][j]) (M[s[i]] == t) M[s[i]] s[i]
* dp[i + 1][next[j][t]] = min(dp[i + 1][next[j][t]], dp[i][j] + 1) (M[s[i]] != t)
*/
map <char, int> M;
char s[1010];
struct Trie
{
int next[1010][4], fail[1010];
bool end[1010];
int L, root;
int newnode()
{
for(int i = 0; i < 4; i++)
next[L][i] = -1;
end[L++] = 0;
return L - 1;
}
void init()
{
L = 0;
root = newnode();
return;
}
void insert(char *s)
{
int now = root;
for(; *s; s++)
{
if(next[now][M[*s]] == -1)
next[now][M[*s]] = newnode();
now = next[now][M[*s]];
}
end[now] = 1;
return;
}
void build()
{
queue <int> Q;
fail[root] = root;
Q.push(root);
while(!Q.empty())
{
int now = Q.front();
Q.pop();
if(end[fail[now]]) end[now] = 1;
for(int i = 0; i < 4; i++)
if(next[now][i] == -1)
next[now][i] = now == root ? root : next[fail[now]][i];
else
{
fail[next[now][i]] = now == root ? root : next[fail[now]][i];
Q.push(next[now][i]);
}
}
return;
}
void solve(int cas)
{
int len = strlen(s);
for(int i = 0; i <= len; i++)
for(int j = 0; j < L; j++)
dp[i][j] = inf;
dp[0][0] = 0;
for(int i = 0; i < len; i++)
for(int j = 0; j < L; j++)
{
if(dp[i][j] == inf) continue;
for(int k = 0; k < 4; k++)
{
if(!end[next[j][k]])//
{
if(k == M[s[i]])//
dp[i + 1][next[j][k]] = min(dp[i + 1][next[j][k]], dp[i][j]);
else
dp[i + 1][next[j][k]] = min(dp[i + 1][next[j][k]], dp[i][j] + 1);
}
}
}
int ans = inf;
for(int i = 0; i < L; i++)
ans = min(ans, dp[len][i]);
printf("Case %d: %d
", cas, ans == inf ? -1 : ans);
return;
}
};
Trie AC;
int main()
{
char ts[22];
int cas = 0;
M['A'] = 0; M['C'] = 1; M['G'] = 2; M['T'] = 3;
while(scanf("%d", &n), n)
{
cas++;
AC.init();
while(n--)
{
scanf("%s", ts);
AC.insert(ts);
}
scanf("%s", s);
AC.build();
AC.solve(cas);
}
return 0;
}
이 내용에 흥미가 있습니까?
현재 기사가 여러분의 문제를 해결하지 못하는 경우 AI 엔진은 머신러닝 분석(스마트 모델이 방금 만들어져 부정확한 경우가 있을 수 있음)을 통해 가장 유사한 기사를 추천합니다:
01 가방, 완전 가방, 다중 가방 dp(동적 기획 입문 dp)01 가방은 2진법으로 직접 표시할 수 있지만 데이터 양이 너무 많으면 시간을 초과하는 것이 폭력이다.01 가방의 사상은 바로 이 물품에 대해 내가 넣은 가치가 큰지 안 넣은 가치가 큰지 비교하여 방정식 f[i][v...
텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
CC BY-SA 2.5, CC BY-SA 3.0 및 CC BY-SA 4.0에 따라 라이센스가 부여됩니다.