#include <bits/stdc++.h>
#include "dna.h"
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define MASK(i) (1ULL << (i))
#define GETBIT(mask, i) (((mask) >> (i)) & 1)
#define ALL(v) (v).begin(), (v).end()
ll max(ll a, ll b){return (a > b) ? a : b;}
ll min(ll a, ll b){return (a < b) ? a : b;}
ll LASTBIT(ll mask){return (mask) & (-mask);}
int pop_cnt(ll mask){return __builtin_popcountll(mask);}
int ctz(ull mask){return __builtin_ctzll(mask);}
int logOf(ull mask){return 63 - __builtin_clzll(mask);}
mt19937_64 rng(chrono::high_resolution_clock::now().time_since_epoch().count());
ll rngesus(ll l, ll r){return l + (ull) rng() % (r - l + 1);}
template <class T1, class T2>
bool maximize(T1 &a, T2 b){
if (a < b) {a = b; return true;}
return false;
}
template <class T1, class T2>
bool minimize(T1 &a, T2 b){
if (a > b) {a = b; return true;}
return false;
}
template <class T>
void printArr(T& container, string separator = " ", string finish = "\n", ostream &out = cout){
for(auto item: container) out << item << separator;
out << finish;
}
template <class T>
void remove_dup(vector<T> &a){
sort(ALL(a));
a.resize(unique(ALL(a)) - a.begin());
}
int convert(char c){
if (c == 'A') return 0;
if (c == 'T') return 1;
if (c == 'C') return 2;
return 3;
}
const int N = 1e5 + 69;
int n;
int cnt[N][3];
int cnt_pair[N][3][3];
void init(string a, string b){
n = a.size();
for(int i = 1; i<=n; ++i){
for(int j = 0; j<3; ++j) cnt[i][j] = cnt[i-1][j];
cnt[i][convert(a[i-1])]++;
cnt[i][convert(b[i-1])]--;
for(int x = 0; x < 3; ++x) for(int y = 0; y < 3; ++y) cnt_pair[i][x][y] = cnt_pair[i-1][x][y];
cnt_pair[i][convert(a[i-1])][convert(b[i-1])]++;
}
}
int get_distance(int l, int r){
l++; r++;
for(int x = 0; x< 3; ++x) if (cnt[r][x] - cnt[l-1][x] != 0) return -1;
int graph[3][3]; memset(graph, 0, sizeof graph);
for(int x = 0; x < 3; ++x) for(int y = 0; y < 3; ++y) graph[x][y] = cnt_pair[r][x][y] - cnt_pair[l-1][x][y];
int ans = 0;
for(int x = 0; x < 3; ++x) for(int y = x + 1; y < 3; ++y) {
int cur = min(graph[x][y], graph[y][x]);
ans += cur;
graph[x][y] -= cur; graph[y][x] -= cur;
}
for(int x = 0; x < 3; ++x) for(int y = 0; y < 3; ++y) if (y != x){
int z = 3 - x - y;
if (x == 0 || y == 1 || z == 2) continue;
int cur = min({graph[0][x], graph[1][y], graph[2][z]});
ans += cur * 2;
graph[0][x] -= cur; graph[1][y] -= cur; graph[2][z] -= cur;
}
return ans;
}
// int main(void){
// ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
// init("ATAC", "ACAT");
// cout << get_distance(1, 3) << "\n";
// return 0;
// }
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
32 ms |
6928 KB |
Output is correct |
2 |
Correct |
27 ms |
6928 KB |
Output is correct |
3 |
Correct |
27 ms |
6444 KB |
Output is correct |
4 |
Correct |
30 ms |
7516 KB |
Output is correct |
5 |
Correct |
0 ms |
344 KB |
Output is correct |
6 |
Correct |
1 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
4 ms |
5468 KB |
Output is correct |
5 |
Correct |
4 ms |
5724 KB |
Output is correct |
6 |
Correct |
3 ms |
5468 KB |
Output is correct |
7 |
Correct |
3 ms |
5212 KB |
Output is correct |
8 |
Correct |
4 ms |
5720 KB |
Output is correct |
9 |
Correct |
3 ms |
5724 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
4 ms |
5468 KB |
Output is correct |
5 |
Correct |
4 ms |
5724 KB |
Output is correct |
6 |
Correct |
3 ms |
5468 KB |
Output is correct |
7 |
Correct |
3 ms |
5212 KB |
Output is correct |
8 |
Correct |
4 ms |
5720 KB |
Output is correct |
9 |
Correct |
3 ms |
5724 KB |
Output is correct |
10 |
Correct |
26 ms |
6940 KB |
Output is correct |
11 |
Correct |
26 ms |
6932 KB |
Output is correct |
12 |
Correct |
27 ms |
6940 KB |
Output is correct |
13 |
Correct |
30 ms |
6936 KB |
Output is correct |
14 |
Correct |
34 ms |
7180 KB |
Output is correct |
15 |
Correct |
31 ms |
7096 KB |
Output is correct |
16 |
Correct |
27 ms |
6908 KB |
Output is correct |
17 |
Correct |
26 ms |
6948 KB |
Output is correct |
18 |
Correct |
27 ms |
7192 KB |
Output is correct |
19 |
Correct |
25 ms |
6940 KB |
Output is correct |
20 |
Correct |
25 ms |
7192 KB |
Output is correct |
21 |
Correct |
25 ms |
7192 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
4 ms |
5468 KB |
Output is correct |
5 |
Correct |
4 ms |
5724 KB |
Output is correct |
6 |
Correct |
3 ms |
5468 KB |
Output is correct |
7 |
Correct |
3 ms |
5212 KB |
Output is correct |
8 |
Correct |
4 ms |
5720 KB |
Output is correct |
9 |
Correct |
3 ms |
5724 KB |
Output is correct |
10 |
Correct |
4 ms |
5212 KB |
Output is correct |
11 |
Correct |
4 ms |
5724 KB |
Output is correct |
12 |
Correct |
4 ms |
5276 KB |
Output is correct |
13 |
Correct |
4 ms |
5624 KB |
Output is correct |
14 |
Correct |
4 ms |
5792 KB |
Output is correct |
15 |
Correct |
4 ms |
5724 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
32 ms |
6928 KB |
Output is correct |
2 |
Correct |
27 ms |
6928 KB |
Output is correct |
3 |
Correct |
27 ms |
6444 KB |
Output is correct |
4 |
Correct |
30 ms |
7516 KB |
Output is correct |
5 |
Correct |
0 ms |
344 KB |
Output is correct |
6 |
Correct |
1 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
344 KB |
Output is correct |
9 |
Correct |
1 ms |
344 KB |
Output is correct |
10 |
Correct |
0 ms |
348 KB |
Output is correct |
11 |
Correct |
4 ms |
5468 KB |
Output is correct |
12 |
Correct |
4 ms |
5724 KB |
Output is correct |
13 |
Correct |
3 ms |
5468 KB |
Output is correct |
14 |
Correct |
3 ms |
5212 KB |
Output is correct |
15 |
Correct |
4 ms |
5720 KB |
Output is correct |
16 |
Correct |
3 ms |
5724 KB |
Output is correct |
17 |
Correct |
26 ms |
6940 KB |
Output is correct |
18 |
Correct |
26 ms |
6932 KB |
Output is correct |
19 |
Correct |
27 ms |
6940 KB |
Output is correct |
20 |
Correct |
30 ms |
6936 KB |
Output is correct |
21 |
Correct |
34 ms |
7180 KB |
Output is correct |
22 |
Correct |
31 ms |
7096 KB |
Output is correct |
23 |
Correct |
27 ms |
6908 KB |
Output is correct |
24 |
Correct |
26 ms |
6948 KB |
Output is correct |
25 |
Correct |
27 ms |
7192 KB |
Output is correct |
26 |
Correct |
25 ms |
6940 KB |
Output is correct |
27 |
Correct |
25 ms |
7192 KB |
Output is correct |
28 |
Correct |
25 ms |
7192 KB |
Output is correct |
29 |
Correct |
4 ms |
5212 KB |
Output is correct |
30 |
Correct |
4 ms |
5724 KB |
Output is correct |
31 |
Correct |
4 ms |
5276 KB |
Output is correct |
32 |
Correct |
4 ms |
5624 KB |
Output is correct |
33 |
Correct |
4 ms |
5792 KB |
Output is correct |
34 |
Correct |
4 ms |
5724 KB |
Output is correct |
35 |
Correct |
0 ms |
348 KB |
Output is correct |
36 |
Correct |
27 ms |
7400 KB |
Output is correct |
37 |
Correct |
27 ms |
7916 KB |
Output is correct |
38 |
Correct |
28 ms |
7820 KB |
Output is correct |
39 |
Correct |
28 ms |
7948 KB |
Output is correct |
40 |
Correct |
28 ms |
8460 KB |
Output is correct |
41 |
Correct |
3 ms |
5720 KB |
Output is correct |
42 |
Correct |
27 ms |
7724 KB |
Output is correct |
43 |
Correct |
28 ms |
8468 KB |
Output is correct |
44 |
Correct |
28 ms |
8360 KB |
Output is correct |
45 |
Correct |
26 ms |
7864 KB |
Output is correct |
46 |
Correct |
26 ms |
8468 KB |
Output is correct |
47 |
Correct |
26 ms |
7944 KB |
Output is correct |