#include "dna.h"
#include<bits/stdc++.h>
using namespace std;
string s,ss;
int dna[2][3][100001],jem[100001],res[100001][6];
map<char,int>mp;
void init(string s,string ss) {
mp['A']=0;
mp['T']=1;
mp['C']=2;
dna[0][mp[s[0]]][0]=1;
dna[1][mp[ss[0]]][0]=1;
if(s[0]==ss[0])jem[0]=1;
for(int i=1;i<s.size();i++){
jem[i]=jem[i-1];
if(s[i]==ss[i])jem[i]++;
for(int j=0;j<3;j++){
dna[0][j][i]=dna[0][j][i-1];
dna[1][j][i]=dna[1][j][i-1];
}
dna[0][mp[s[i]]][i]++;
dna[1][mp[ss[i]]][i]++;
}
for(int i=0;i<s.size();i++){
if(i)for(int j=0;j<6;j++)res[i][j]=res[i-1][j];
if(s[i]=='A'){
if(ss[i]=='T'){
res[i][0]++;
}
else if(ss[i]=='C'){
res[i][2]++;
}
}else if(s[i]=='T'){
if(ss[i]=='A'){
res[i][1]++;
}else if(ss[i]=='C'){
res[i][5]++;
}
}else if(s[i]=='C'){
if(ss[i]=='A'){
res[i][3]++;
}else if(ss[i]=='T'){
res[i][4]++;
}
}
}
}
int get_distance(int x,int y) {
int cnt=jem[y]-(x-1<0?0:jem[x-1]);
int a1=dna[0][0][y]-(x-1<0?0:dna[0][0][x-1]);
int t1=dna[0][1][y]-(x-1<0?0:dna[0][1][x-1]);
int c1=dna[0][2][y]-(x-1<0?0:dna[0][2][x-1]);
int a2=dna[1][0][y]-(x-1<0?0:dna[1][0][x-1]);
int t2=dna[1][1][y]-(x-1<0?0:dna[1][1][x-1]);
int c2=dna[1][2][y]-(x-1<0?0:dna[1][2][x-1]);
int len=y-x+1;
if(a1!=a2 || t1!=t2 || c1!=c2)return -1;
int at=min(res[y][0]-(x-1<0?0:res[x-1][0]),res[y][1]-(x-1<0?0:res[x-1][1]));
int ac=min(res[y][2]-(x-1<0?0:res[x-1][2]),res[y][3]-(x-1<0?0:res[x-1][3]));
int tc=min(res[y][4]-(x-1<0?0:res[x-1][4]),res[y][5]-(x-1<0?0:res[x-1][5]));
int atc=len-(cnt+(at*2)+(ac*2)+(tc*2));
return (at+ac+tc+((atc/3)*2));
}
Compilation message
dna.cpp: In function 'void init(std::string, std::string)':
dna.cpp:14:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::__cxx11::basic_string<char>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
14 | for(int i=1;i<s.size();i++){
| ~^~~~~~~~~
dna.cpp:24:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::__cxx11::basic_string<char>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
24 | for(int i=0;i<s.size();i++){
| ~^~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
38 ms |
7308 KB |
Output is correct |
2 |
Correct |
36 ms |
7508 KB |
Output is correct |
3 |
Correct |
37 ms |
6900 KB |
Output is correct |
4 |
Correct |
39 ms |
7524 KB |
Output is correct |
5 |
Correct |
0 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
308 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
0 ms |
212 KB |
Output is correct |
4 |
Correct |
8 ms |
6012 KB |
Output is correct |
5 |
Correct |
7 ms |
6100 KB |
Output is correct |
6 |
Correct |
7 ms |
6072 KB |
Output is correct |
7 |
Correct |
7 ms |
5716 KB |
Output is correct |
8 |
Correct |
8 ms |
6096 KB |
Output is correct |
9 |
Correct |
6 ms |
6100 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
308 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
0 ms |
212 KB |
Output is correct |
4 |
Correct |
8 ms |
6012 KB |
Output is correct |
5 |
Correct |
7 ms |
6100 KB |
Output is correct |
6 |
Correct |
7 ms |
6072 KB |
Output is correct |
7 |
Correct |
7 ms |
5716 KB |
Output is correct |
8 |
Correct |
8 ms |
6096 KB |
Output is correct |
9 |
Correct |
6 ms |
6100 KB |
Output is correct |
10 |
Correct |
41 ms |
7312 KB |
Output is correct |
11 |
Correct |
35 ms |
7432 KB |
Output is correct |
12 |
Correct |
36 ms |
7320 KB |
Output is correct |
13 |
Correct |
37 ms |
7576 KB |
Output is correct |
14 |
Correct |
38 ms |
7704 KB |
Output is correct |
15 |
Correct |
33 ms |
7620 KB |
Output is correct |
16 |
Correct |
33 ms |
7324 KB |
Output is correct |
17 |
Correct |
36 ms |
7496 KB |
Output is correct |
18 |
Correct |
33 ms |
7728 KB |
Output is correct |
19 |
Correct |
34 ms |
7464 KB |
Output is correct |
20 |
Correct |
32 ms |
7592 KB |
Output is correct |
21 |
Correct |
32 ms |
7816 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
308 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
0 ms |
212 KB |
Output is correct |
4 |
Correct |
8 ms |
6012 KB |
Output is correct |
5 |
Correct |
7 ms |
6100 KB |
Output is correct |
6 |
Correct |
7 ms |
6072 KB |
Output is correct |
7 |
Correct |
7 ms |
5716 KB |
Output is correct |
8 |
Correct |
8 ms |
6096 KB |
Output is correct |
9 |
Correct |
6 ms |
6100 KB |
Output is correct |
10 |
Correct |
7 ms |
5588 KB |
Output is correct |
11 |
Correct |
8 ms |
6100 KB |
Output is correct |
12 |
Correct |
8 ms |
5852 KB |
Output is correct |
13 |
Correct |
8 ms |
6076 KB |
Output is correct |
14 |
Correct |
8 ms |
6120 KB |
Output is correct |
15 |
Correct |
6 ms |
6080 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
38 ms |
7308 KB |
Output is correct |
2 |
Correct |
36 ms |
7508 KB |
Output is correct |
3 |
Correct |
37 ms |
6900 KB |
Output is correct |
4 |
Correct |
39 ms |
7524 KB |
Output is correct |
5 |
Correct |
0 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
1 ms |
308 KB |
Output is correct |
9 |
Correct |
1 ms |
340 KB |
Output is correct |
10 |
Correct |
0 ms |
212 KB |
Output is correct |
11 |
Correct |
8 ms |
6012 KB |
Output is correct |
12 |
Correct |
7 ms |
6100 KB |
Output is correct |
13 |
Correct |
7 ms |
6072 KB |
Output is correct |
14 |
Correct |
7 ms |
5716 KB |
Output is correct |
15 |
Correct |
8 ms |
6096 KB |
Output is correct |
16 |
Correct |
6 ms |
6100 KB |
Output is correct |
17 |
Correct |
41 ms |
7312 KB |
Output is correct |
18 |
Correct |
35 ms |
7432 KB |
Output is correct |
19 |
Correct |
36 ms |
7320 KB |
Output is correct |
20 |
Correct |
37 ms |
7576 KB |
Output is correct |
21 |
Correct |
38 ms |
7704 KB |
Output is correct |
22 |
Correct |
33 ms |
7620 KB |
Output is correct |
23 |
Correct |
33 ms |
7324 KB |
Output is correct |
24 |
Correct |
36 ms |
7496 KB |
Output is correct |
25 |
Correct |
33 ms |
7728 KB |
Output is correct |
26 |
Correct |
34 ms |
7464 KB |
Output is correct |
27 |
Correct |
32 ms |
7592 KB |
Output is correct |
28 |
Correct |
32 ms |
7816 KB |
Output is correct |
29 |
Correct |
7 ms |
5588 KB |
Output is correct |
30 |
Correct |
8 ms |
6100 KB |
Output is correct |
31 |
Correct |
8 ms |
5852 KB |
Output is correct |
32 |
Correct |
8 ms |
6076 KB |
Output is correct |
33 |
Correct |
8 ms |
6120 KB |
Output is correct |
34 |
Correct |
6 ms |
6080 KB |
Output is correct |
35 |
Correct |
0 ms |
340 KB |
Output is correct |
36 |
Correct |
35 ms |
8036 KB |
Output is correct |
37 |
Correct |
38 ms |
8556 KB |
Output is correct |
38 |
Correct |
40 ms |
8852 KB |
Output is correct |
39 |
Correct |
36 ms |
8968 KB |
Output is correct |
40 |
Correct |
41 ms |
8908 KB |
Output is correct |
41 |
Correct |
6 ms |
6084 KB |
Output is correct |
42 |
Correct |
33 ms |
8564 KB |
Output is correct |
43 |
Correct |
34 ms |
8864 KB |
Output is correct |
44 |
Correct |
35 ms |
8840 KB |
Output is correct |
45 |
Correct |
31 ms |
8592 KB |
Output is correct |
46 |
Correct |
33 ms |
8840 KB |
Output is correct |
47 |
Correct |
32 ms |
8840 KB |
Output is correct |