#include "dna.h"
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
string s,ss;
ll dna[2][3][100001],jem[100001],res[100001][6];
map<char,ll>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++){
if(s[i]==ss[i])jem[i]=jem[i-1]+1;
else jem[i]=jem[i-1];
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) {
ll cnt=jem[y]-(x-1<0?0:jem[x-1]);
ll a1=dna[0][0][y]-(x-1<0?0:dna[0][0][x-1]);
ll t1=dna[0][1][y]-(x-1<0?0:dna[0][1][x-1]);
ll c1=dna[0][2][y]-(x-1<0?0:dna[0][2][x-1]);
ll a2=dna[1][0][y]-(x-1<0?0:dna[1][0][x-1]);
ll t2=dna[1][1][y]-(x-1<0?0:dna[1][1][x-1]);
ll c2=dna[1][2][y]-(x-1<0?0:dna[1][2][x-1]);
ll len=y-x+1;
if(a1!=a2 || t1!=t2 || c1!=c2)return -1;
ll at=min(res[y][0]-(x-1<0?0:res[x-1][0]),res[y][1]-(x-1<0?0:res[x-1][1]));
ll ac=min(res[y][2]-(x-1<0?0:res[x-1][2]),res[y][3]-(x-1<0?0:res[x-1][3]));
ll tc=min(res[y][4]-(x-1<0?0:res[x-1][4]),res[y][5]-(x-1<0?0:res[x-1][5]));
ll atc=len-(cnt+(at*2)+(ac*2)+(tc*2));
return (at+ac+tc+((2*atc)/3));
}
Compilation message
dna.cpp: In function 'void init(std::string, std::string)':
dna.cpp:15:19: warning: comparison of integer expressions of different signedness: 'int' and 'std::__cxx11::basic_string<char>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
15 | for(int i=1;i<s.size();i++){
| ~^~~~~~~~~
dna.cpp:25:19: warning: comparison of integer expressions of different signedness: 'int' and 'std::__cxx11::basic_string<char>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
25 | for(int i=0;i<s.size();i++){
| ~^~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
46 ms |
13532 KB |
Output is correct |
2 |
Correct |
48 ms |
13640 KB |
Output is correct |
3 |
Correct |
46 ms |
12756 KB |
Output is correct |
4 |
Correct |
44 ms |
13656 KB |
Output is correct |
5 |
Correct |
1 ms |
256 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
212 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
212 KB |
Output is correct |
4 |
Correct |
10 ms |
11100 KB |
Output is correct |
5 |
Correct |
9 ms |
11240 KB |
Output is correct |
6 |
Correct |
11 ms |
11200 KB |
Output is correct |
7 |
Correct |
9 ms |
10528 KB |
Output is correct |
8 |
Correct |
9 ms |
11240 KB |
Output is correct |
9 |
Correct |
11 ms |
11228 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
212 KB |
Output is correct |
4 |
Correct |
10 ms |
11100 KB |
Output is correct |
5 |
Correct |
9 ms |
11240 KB |
Output is correct |
6 |
Correct |
11 ms |
11200 KB |
Output is correct |
7 |
Correct |
9 ms |
10528 KB |
Output is correct |
8 |
Correct |
9 ms |
11240 KB |
Output is correct |
9 |
Correct |
11 ms |
11228 KB |
Output is correct |
10 |
Correct |
43 ms |
13488 KB |
Output is correct |
11 |
Correct |
42 ms |
13684 KB |
Output is correct |
12 |
Correct |
41 ms |
13304 KB |
Output is correct |
13 |
Correct |
43 ms |
13488 KB |
Output is correct |
14 |
Correct |
42 ms |
13944 KB |
Output is correct |
15 |
Correct |
44 ms |
13960 KB |
Output is correct |
16 |
Correct |
55 ms |
13160 KB |
Output is correct |
17 |
Correct |
48 ms |
13460 KB |
Output is correct |
18 |
Correct |
41 ms |
13936 KB |
Output is correct |
19 |
Correct |
33 ms |
13208 KB |
Output is correct |
20 |
Correct |
43 ms |
13460 KB |
Output is correct |
21 |
Correct |
41 ms |
13972 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
212 KB |
Output is correct |
4 |
Correct |
10 ms |
11100 KB |
Output is correct |
5 |
Correct |
9 ms |
11240 KB |
Output is correct |
6 |
Correct |
11 ms |
11200 KB |
Output is correct |
7 |
Correct |
9 ms |
10528 KB |
Output is correct |
8 |
Correct |
9 ms |
11240 KB |
Output is correct |
9 |
Correct |
11 ms |
11228 KB |
Output is correct |
10 |
Correct |
10 ms |
10324 KB |
Output is correct |
11 |
Correct |
11 ms |
11224 KB |
Output is correct |
12 |
Correct |
11 ms |
10532 KB |
Output is correct |
13 |
Correct |
11 ms |
11100 KB |
Output is correct |
14 |
Correct |
10 ms |
11200 KB |
Output is correct |
15 |
Correct |
9 ms |
11164 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
46 ms |
13532 KB |
Output is correct |
2 |
Correct |
48 ms |
13640 KB |
Output is correct |
3 |
Correct |
46 ms |
12756 KB |
Output is correct |
4 |
Correct |
44 ms |
13656 KB |
Output is correct |
5 |
Correct |
1 ms |
256 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
212 KB |
Output is correct |
8 |
Correct |
1 ms |
340 KB |
Output is correct |
9 |
Correct |
1 ms |
340 KB |
Output is correct |
10 |
Correct |
1 ms |
212 KB |
Output is correct |
11 |
Correct |
10 ms |
11100 KB |
Output is correct |
12 |
Correct |
9 ms |
11240 KB |
Output is correct |
13 |
Correct |
11 ms |
11200 KB |
Output is correct |
14 |
Correct |
9 ms |
10528 KB |
Output is correct |
15 |
Correct |
9 ms |
11240 KB |
Output is correct |
16 |
Correct |
11 ms |
11228 KB |
Output is correct |
17 |
Correct |
43 ms |
13488 KB |
Output is correct |
18 |
Correct |
42 ms |
13684 KB |
Output is correct |
19 |
Correct |
41 ms |
13304 KB |
Output is correct |
20 |
Correct |
43 ms |
13488 KB |
Output is correct |
21 |
Correct |
42 ms |
13944 KB |
Output is correct |
22 |
Correct |
44 ms |
13960 KB |
Output is correct |
23 |
Correct |
55 ms |
13160 KB |
Output is correct |
24 |
Correct |
48 ms |
13460 KB |
Output is correct |
25 |
Correct |
41 ms |
13936 KB |
Output is correct |
26 |
Correct |
33 ms |
13208 KB |
Output is correct |
27 |
Correct |
43 ms |
13460 KB |
Output is correct |
28 |
Correct |
41 ms |
13972 KB |
Output is correct |
29 |
Correct |
10 ms |
10324 KB |
Output is correct |
30 |
Correct |
11 ms |
11224 KB |
Output is correct |
31 |
Correct |
11 ms |
10532 KB |
Output is correct |
32 |
Correct |
11 ms |
11100 KB |
Output is correct |
33 |
Correct |
10 ms |
11200 KB |
Output is correct |
34 |
Correct |
9 ms |
11164 KB |
Output is correct |
35 |
Correct |
1 ms |
340 KB |
Output is correct |
36 |
Correct |
51 ms |
12760 KB |
Output is correct |
37 |
Correct |
48 ms |
13672 KB |
Output is correct |
38 |
Correct |
46 ms |
13408 KB |
Output is correct |
39 |
Correct |
54 ms |
14004 KB |
Output is correct |
40 |
Correct |
48 ms |
14016 KB |
Output is correct |
41 |
Correct |
9 ms |
11240 KB |
Output is correct |
42 |
Correct |
47 ms |
13296 KB |
Output is correct |
43 |
Correct |
50 ms |
13944 KB |
Output is correct |
44 |
Correct |
47 ms |
13944 KB |
Output is correct |
45 |
Correct |
38 ms |
13332 KB |
Output is correct |
46 |
Correct |
43 ms |
13972 KB |
Output is correct |
47 |
Correct |
44 ms |
14004 KB |
Output is correct |