| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 708906 | Essa2006 | DNA 돌연변이 (IOI21_dna) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include<bits/stdc++.h>
//#include "dna.h"
using namespace std;
#define ll long long
#define endl '\n'
#define FF first
#define SS second
#define all(a) a.begin(), a.end()
#define mod (ll)(1000000007)
vector<array<array<int, 3>, 3>>AA;
void init(string a, string b) {
int n=a.size();
AA.resize(n);
for(int i=0;i<n;i++){
if(a[i]=='A')
a[i]='1';
if(a[i]=='C')
a[i]='3';
if(a[i]=='T')
a[i]='2';
if(b[i]=='A')
b[i]='1';
if(b[i]=='C')
b[i]='3';
if(b[i]=='T')
b[i]='2';
}
for(int i=0;i<n;i++){
AA[i][a[i]-'1'][b[i]-'1']++;
}
for(int i=1;i<n;i++){
for(int j=0;j<3;j++){
for(int k=0;k<3;k++){
AA[i][j][k]+=AA[i-1][j][k];
}
}
}
}
int get_distance(int x, int y) {
array<array<int, 3>, 3>b;
for(int j=0;j<3;j++){
for(int k=0;k<3;k++){
b[j][k]=AA[y][j][k];
b[j][k]-=(x?AA[x-1][j][k]:0);
}
}
map<int, int>v, v2;
for(int i=0;i<3;i++){
for(int j=i+1;j<3;j++){
v[i]+=b[i][j], v[j]+=b[j][i], v2[i]+=b[j][i], v2[j]+=b[i][j];
}
}
for(int i=0;i<3;i++){
if(v[i]!=v2[i]){
return -1;
}
}
int ans=0, cur=0;
for(int i=0;i<3;i++){
for(int j=i+1;j<3;j++){
int mn=min(b[i][j], b[j][i]);
ans+=mn;
b[i][j]-=mn, b[j][i]-=mn;
cur+=b[i][j]+b[j][i];
}
}
cur/=3;
cur*=2;
ans+=cur;
return ans;
}
int main() {
int n, q;
assert(scanf("%d %d", &n, &q) == 2);
char A[n+1], B[n+1];
assert(scanf("%s", A) == 1);
assert(scanf("%s", B) == 1);
std::string a = std::string(A);
std::string b = std::string(B);
std::vector<int> x(q), y(q);
for (int i = 0; i < q; i++) {
assert(scanf("%d %d", &x[i], &y[i]) == 2);
}
fclose(stdin);
std::vector<int> results(q);
init(a, b);
for (int i = 0; i < q; i++) {
results[i] = get_distance(x[i], y[i]);
}
for (int i = 0; i < q; i++) {
printf("%d\n", results[i]);
}
fclose(stdout);
return 0;
}