#include "dna.h"
#include<bits/stdc++.h>
using namespace std;
long long s[100002][8], pref[100002][4], pref2[100002][4];
void init(std::string a, std::string b) {
long long i,j;
for(i=0;i<a.size();i++)
{
if(a[i]=='A' && b[i]=='T')s[i][1]++;
if(a[i]=='T' && b[i]=='A')s[i][2]++;
if(a[i]=='C' && b[i]=='T')s[i][3]++;
if(a[i]=='T' && b[i]=='C')s[i][4]++;
if(a[i]=='A' && b[i]=='C')s[i][5]++;
if(a[i]=='C' && b[i]=='A')s[i][6]++;
if(i!=0)for(j=1;j<=6;j++)s[i][j]+=s[i-1][j];
if(a[i]=='A')pref[i][1]++;
if(a[i]=='T')pref[i][2]++;
if(a[i]=='C')pref[i][3]++;
if(i!=0)for(j=1;j<=3;j++)pref[i][j]+=pref[i-1][j];
if(b[i]=='A')pref2[i][1]++;
if(b[i]=='T')pref2[i][2]++;
if(b[i]=='C')pref2[i][3]++;
if(i!=0)for(j=1;j<=3;j++)pref2[i][j]+=pref2[i-1][j];
}
}
int get_distance(int x, int y) {
long long A=pref[y][1], A2=pref2[y][1];
if(x>0){A-=pref[x-1][1]; A2-=pref2[x-1][1];}
long long B=pref[y][2], B2=pref2[y][2];
if(x>0){B-=pref[x-1][2]; B2-=pref2[x-1][2];}
long long C=pref[y][3], C2=pref2[y][3];
if(x>0){C-=pref[x-1][3]; C2-=pref2[x-1][3];}
if((A!=A2) || (B!=B2 )||(C!=C2))return -1;
long long AT = s[y][1]-(x>0?s[x-1][1]:0);
long long TA = s[y][2]-(x>0?s[x-1][2]:0);
long long CT = s[y][3]-(x>0?s[x-1][3]:0);
long long TC = s[y][4]-(x>0?s[x-1][4]:0);
long long AC = s[y][5]-(x>0?s[x-1][5]:0);
long long CA = s[y][6]-(x>0?s[x-1][6]:0);
long long ans=0;
ans+=min(AT,TA)+min(CT,TC)+min(AC,CA);
AT-=TA;
if(AT<0)AT=-AT;
ans+=2*AT;
return ans;
}
Compilation message
dna.cpp: In function 'void init(std::string, std::string)':
dna.cpp:9:14: warning: comparison of integer expressions of different signedness: 'long long int' and 'std::__cxx11::basic_string<char>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
9 | for(i=0;i<a.size();i++)
| ~^~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
37 ms |
15884 KB |
Output is correct |
2 |
Correct |
37 ms |
16028 KB |
Output is correct |
3 |
Correct |
37 ms |
14860 KB |
Output is correct |
4 |
Correct |
41 ms |
15992 KB |
Output is correct |
5 |
Correct |
0 ms |
308 KB |
Output is correct |
6 |
Correct |
1 ms |
212 KB |
Output is correct |
7 |
Correct |
1 ms |
212 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
300 KB |
Output is correct |
2 |
Correct |
0 ms |
308 KB |
Output is correct |
3 |
Correct |
1 ms |
212 KB |
Output is correct |
4 |
Correct |
8 ms |
13396 KB |
Output is correct |
5 |
Correct |
8 ms |
13524 KB |
Output is correct |
6 |
Correct |
8 ms |
13524 KB |
Output is correct |
7 |
Correct |
9 ms |
12604 KB |
Output is correct |
8 |
Correct |
8 ms |
13524 KB |
Output is correct |
9 |
Correct |
10 ms |
13524 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
300 KB |
Output is correct |
2 |
Correct |
0 ms |
308 KB |
Output is correct |
3 |
Correct |
1 ms |
212 KB |
Output is correct |
4 |
Correct |
8 ms |
13396 KB |
Output is correct |
5 |
Correct |
8 ms |
13524 KB |
Output is correct |
6 |
Correct |
8 ms |
13524 KB |
Output is correct |
7 |
Correct |
9 ms |
12604 KB |
Output is correct |
8 |
Correct |
8 ms |
13524 KB |
Output is correct |
9 |
Correct |
10 ms |
13524 KB |
Output is correct |
10 |
Correct |
42 ms |
15852 KB |
Output is correct |
11 |
Correct |
39 ms |
16032 KB |
Output is correct |
12 |
Correct |
38 ms |
15384 KB |
Output is correct |
13 |
Correct |
43 ms |
15824 KB |
Output is correct |
14 |
Correct |
40 ms |
16368 KB |
Output is correct |
15 |
Correct |
37 ms |
16264 KB |
Output is correct |
16 |
Correct |
38 ms |
15392 KB |
Output is correct |
17 |
Correct |
38 ms |
15704 KB |
Output is correct |
18 |
Correct |
36 ms |
16272 KB |
Output is correct |
19 |
Correct |
32 ms |
15340 KB |
Output is correct |
20 |
Correct |
32 ms |
15764 KB |
Output is correct |
21 |
Correct |
32 ms |
16264 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
300 KB |
Output is correct |
2 |
Correct |
0 ms |
308 KB |
Output is correct |
3 |
Correct |
1 ms |
212 KB |
Output is correct |
4 |
Correct |
8 ms |
13396 KB |
Output is correct |
5 |
Correct |
8 ms |
13524 KB |
Output is correct |
6 |
Correct |
8 ms |
13524 KB |
Output is correct |
7 |
Correct |
9 ms |
12604 KB |
Output is correct |
8 |
Correct |
8 ms |
13524 KB |
Output is correct |
9 |
Correct |
10 ms |
13524 KB |
Output is correct |
10 |
Correct |
8 ms |
12372 KB |
Output is correct |
11 |
Correct |
8 ms |
13572 KB |
Output is correct |
12 |
Correct |
8 ms |
12628 KB |
Output is correct |
13 |
Correct |
9 ms |
13496 KB |
Output is correct |
14 |
Correct |
9 ms |
13524 KB |
Output is correct |
15 |
Correct |
8 ms |
13504 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
37 ms |
15884 KB |
Output is correct |
2 |
Correct |
37 ms |
16028 KB |
Output is correct |
3 |
Correct |
37 ms |
14860 KB |
Output is correct |
4 |
Correct |
41 ms |
15992 KB |
Output is correct |
5 |
Correct |
0 ms |
308 KB |
Output is correct |
6 |
Correct |
1 ms |
212 KB |
Output is correct |
7 |
Correct |
1 ms |
212 KB |
Output is correct |
8 |
Correct |
1 ms |
300 KB |
Output is correct |
9 |
Correct |
0 ms |
308 KB |
Output is correct |
10 |
Correct |
1 ms |
212 KB |
Output is correct |
11 |
Correct |
8 ms |
13396 KB |
Output is correct |
12 |
Correct |
8 ms |
13524 KB |
Output is correct |
13 |
Correct |
8 ms |
13524 KB |
Output is correct |
14 |
Correct |
9 ms |
12604 KB |
Output is correct |
15 |
Correct |
8 ms |
13524 KB |
Output is correct |
16 |
Correct |
10 ms |
13524 KB |
Output is correct |
17 |
Correct |
42 ms |
15852 KB |
Output is correct |
18 |
Correct |
39 ms |
16032 KB |
Output is correct |
19 |
Correct |
38 ms |
15384 KB |
Output is correct |
20 |
Correct |
43 ms |
15824 KB |
Output is correct |
21 |
Correct |
40 ms |
16368 KB |
Output is correct |
22 |
Correct |
37 ms |
16264 KB |
Output is correct |
23 |
Correct |
38 ms |
15392 KB |
Output is correct |
24 |
Correct |
38 ms |
15704 KB |
Output is correct |
25 |
Correct |
36 ms |
16272 KB |
Output is correct |
26 |
Correct |
32 ms |
15340 KB |
Output is correct |
27 |
Correct |
32 ms |
15764 KB |
Output is correct |
28 |
Correct |
32 ms |
16264 KB |
Output is correct |
29 |
Correct |
8 ms |
12372 KB |
Output is correct |
30 |
Correct |
8 ms |
13572 KB |
Output is correct |
31 |
Correct |
8 ms |
12628 KB |
Output is correct |
32 |
Correct |
9 ms |
13496 KB |
Output is correct |
33 |
Correct |
9 ms |
13524 KB |
Output is correct |
34 |
Correct |
8 ms |
13504 KB |
Output is correct |
35 |
Correct |
0 ms |
312 KB |
Output is correct |
36 |
Correct |
40 ms |
14948 KB |
Output is correct |
37 |
Correct |
40 ms |
15948 KB |
Output is correct |
38 |
Correct |
38 ms |
15636 KB |
Output is correct |
39 |
Correct |
40 ms |
16372 KB |
Output is correct |
40 |
Correct |
41 ms |
16368 KB |
Output is correct |
41 |
Correct |
8 ms |
13524 KB |
Output is correct |
42 |
Correct |
36 ms |
15468 KB |
Output is correct |
43 |
Correct |
37 ms |
16360 KB |
Output is correct |
44 |
Correct |
36 ms |
16248 KB |
Output is correct |
45 |
Correct |
31 ms |
15512 KB |
Output is correct |
46 |
Correct |
35 ms |
16352 KB |
Output is correct |
47 |
Correct |
32 ms |
16300 KB |
Output is correct |