#include <bits/stdc++.h>
#define pb push_back
#define f first
#define sc second
using namespace std;
typedef long long int ll;
typedef string str;
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
int next(int a, int b){
return uniform_int_distribution<int>(a, b)(rng);
}
int n, m, k;
bitset<4100> candidate;
char s[4100][4100];
int cnt[4100][4];
int dna[255];
int comp(int a, int b){
int c = 0;
for(int i = 0; i < m; i++) if(s[a][i] != s[b][i]) c++;
return c;
}
int ans = -1;
void dfs(int nd){
if(ans != -1) return;
if(!candidate[nd]) return;
int nxt = -1;
bool ok = 1;
for(int i = 0; i < n; i++) if(i != nd && candidate[i]){
int c = comp(i, nd);
if(c == k){
nxt = i;
if(!ok){
dfs(i);
break;
}
}
else{
ok = 0;
candidate[nd] = 0;
candidate[i] = 0;
if(nxt != -1){
dfs(nxt);
break;
}
}
}
if(!ok) return;
for(int i = 0; i < n; i++) if(i != nd && !candidate[i]){
int c = comp(i, nd);
if(c != k){
ok = 0;
break;
}
}
if(ok) ans = nd;
else candidate[nd] = 0, dfs(nxt);
}
int main(){
ios_base::sync_with_stdio(0);
cin.tie(0);
dna['A'] = 0, dna['T'] = 1, dna['C'] = 2, dna['G'] = 3;
cin >> n >> m >> k;
candidate.flip();
for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) cin >> s[i][j];
for(int T = 0; T <= 5; T++){
vector<int> cur;
vector<bool> in(n, 0);
for(int i = 0; i < n; i++){
if(T == 0 || next(0, 1) == 1) in[i] = 1, cur.pb(i);
}
if(cur.empty()){
cur = {0};
in[0] = 1;
}
for(int i: cur) for(int j = 0; j < m; j++){
if(i == cur[0]) cnt[j][0] = 0, cnt[j][1] = 0, cnt[j][2] = 0, cnt[j][3] = 0;
cnt[j][dna[s[i][j]]]++;
}
for(int i = 0; i < n; i++){
int c = 0;
for(int j = 0; j < m; j++) for(int x = 0; x < 4; x++) if(x != dna[s[i][j]]) c+=cnt[j][x];
if(c != k*cur.size()-k*in[i]) candidate[i] = 0;
}
}
for(int i = 0; i < n; i++) if(candidate[i]){
dfs(i);
break;
}
cout << ans+1 << "\n";
}
Compilation message
genetics.cpp: In function 'int main()':
genetics.cpp:88:30: warning: array subscript has type 'char' [-Wchar-subscripts]
88 | cnt[j][dna[s[i][j]]]++;
| ~~~~~~^
genetics.cpp:93:85: warning: array subscript has type 'char' [-Wchar-subscripts]
93 | for(int j = 0; j < m; j++) for(int x = 0; x < 4; x++) if(x != dna[s[i][j]]) c+=cnt[j][x];
| ~~~~~~^
genetics.cpp:94:18: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
94 | if(c != k*cur.size()-k*in[i]) candidate[i] = 0;
| ~~^~~~~~~~~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
596 KB |
Output is correct |
2 |
Correct |
2 ms |
724 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
596 KB |
Output is correct |
5 |
Correct |
1 ms |
596 KB |
Output is correct |
6 |
Correct |
1 ms |
660 KB |
Output is correct |
7 |
Correct |
1 ms |
468 KB |
Output is correct |
8 |
Correct |
0 ms |
340 KB |
Output is correct |
9 |
Correct |
2 ms |
468 KB |
Output is correct |
10 |
Correct |
2 ms |
596 KB |
Output is correct |
11 |
Correct |
2 ms |
596 KB |
Output is correct |
12 |
Correct |
2 ms |
724 KB |
Output is correct |
13 |
Correct |
1 ms |
724 KB |
Output is correct |
14 |
Correct |
1 ms |
724 KB |
Output is correct |
15 |
Correct |
2 ms |
724 KB |
Output is correct |
16 |
Correct |
1 ms |
724 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
214 ms |
6388 KB |
Output is correct |
2 |
Correct |
280 ms |
7576 KB |
Output is correct |
3 |
Correct |
281 ms |
7372 KB |
Output is correct |
4 |
Correct |
64 ms |
4948 KB |
Output is correct |
5 |
Correct |
268 ms |
7532 KB |
Output is correct |
6 |
Correct |
284 ms |
7580 KB |
Output is correct |
7 |
Correct |
115 ms |
3356 KB |
Output is correct |
8 |
Correct |
113 ms |
3412 KB |
Output is correct |
9 |
Correct |
252 ms |
7292 KB |
Output is correct |
10 |
Correct |
253 ms |
7244 KB |
Output is correct |
11 |
Correct |
221 ms |
6432 KB |
Output is correct |
12 |
Correct |
246 ms |
6384 KB |
Output is correct |
13 |
Correct |
225 ms |
6432 KB |
Output is correct |
14 |
Correct |
190 ms |
5588 KB |
Output is correct |
15 |
Correct |
184 ms |
5640 KB |
Output is correct |
16 |
Correct |
196 ms |
6292 KB |
Output is correct |
17 |
Correct |
262 ms |
7224 KB |
Output is correct |
18 |
Correct |
261 ms |
7220 KB |
Output is correct |
19 |
Correct |
266 ms |
7280 KB |
Output is correct |
20 |
Correct |
264 ms |
7184 KB |
Output is correct |
21 |
Correct |
260 ms |
7372 KB |
Output is correct |
22 |
Correct |
267 ms |
7120 KB |
Output is correct |
23 |
Correct |
253 ms |
7144 KB |
Output is correct |
24 |
Correct |
273 ms |
7232 KB |
Output is correct |
25 |
Correct |
272 ms |
7188 KB |
Output is correct |
26 |
Correct |
260 ms |
7244 KB |
Output is correct |
27 |
Correct |
272 ms |
7204 KB |
Output is correct |
28 |
Correct |
264 ms |
7192 KB |
Output is correct |
29 |
Correct |
256 ms |
7240 KB |
Output is correct |
30 |
Correct |
279 ms |
7624 KB |
Output is correct |
31 |
Correct |
261 ms |
7584 KB |
Output is correct |
32 |
Correct |
259 ms |
7632 KB |
Output is correct |
33 |
Correct |
0 ms |
340 KB |
Output is correct |
34 |
Correct |
1 ms |
724 KB |
Output is correct |
35 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
214 ms |
6388 KB |
Output is correct |
2 |
Correct |
280 ms |
7576 KB |
Output is correct |
3 |
Correct |
281 ms |
7372 KB |
Output is correct |
4 |
Correct |
64 ms |
4948 KB |
Output is correct |
5 |
Correct |
268 ms |
7532 KB |
Output is correct |
6 |
Correct |
284 ms |
7580 KB |
Output is correct |
7 |
Correct |
115 ms |
3356 KB |
Output is correct |
8 |
Correct |
113 ms |
3412 KB |
Output is correct |
9 |
Correct |
252 ms |
7292 KB |
Output is correct |
10 |
Correct |
253 ms |
7244 KB |
Output is correct |
11 |
Correct |
221 ms |
6432 KB |
Output is correct |
12 |
Correct |
246 ms |
6384 KB |
Output is correct |
13 |
Correct |
225 ms |
6432 KB |
Output is correct |
14 |
Correct |
190 ms |
5588 KB |
Output is correct |
15 |
Correct |
184 ms |
5640 KB |
Output is correct |
16 |
Correct |
196 ms |
6292 KB |
Output is correct |
17 |
Correct |
262 ms |
7224 KB |
Output is correct |
18 |
Correct |
261 ms |
7220 KB |
Output is correct |
19 |
Correct |
266 ms |
7280 KB |
Output is correct |
20 |
Correct |
264 ms |
7184 KB |
Output is correct |
21 |
Correct |
260 ms |
7372 KB |
Output is correct |
22 |
Correct |
267 ms |
7120 KB |
Output is correct |
23 |
Correct |
253 ms |
7144 KB |
Output is correct |
24 |
Correct |
273 ms |
7232 KB |
Output is correct |
25 |
Correct |
272 ms |
7188 KB |
Output is correct |
26 |
Correct |
260 ms |
7244 KB |
Output is correct |
27 |
Correct |
272 ms |
7204 KB |
Output is correct |
28 |
Correct |
264 ms |
7192 KB |
Output is correct |
29 |
Correct |
256 ms |
7240 KB |
Output is correct |
30 |
Correct |
279 ms |
7624 KB |
Output is correct |
31 |
Correct |
261 ms |
7584 KB |
Output is correct |
32 |
Correct |
259 ms |
7632 KB |
Output is correct |
33 |
Correct |
0 ms |
340 KB |
Output is correct |
34 |
Correct |
1 ms |
724 KB |
Output is correct |
35 |
Correct |
1 ms |
340 KB |
Output is correct |
36 |
Correct |
1175 ms |
14796 KB |
Output is correct |
37 |
Correct |
1338 ms |
33296 KB |
Output is correct |
38 |
Correct |
1319 ms |
32540 KB |
Output is correct |
39 |
Correct |
493 ms |
20172 KB |
Output is correct |
40 |
Correct |
1331 ms |
33240 KB |
Output is correct |
41 |
Correct |
679 ms |
16812 KB |
Output is correct |
42 |
Correct |
667 ms |
16904 KB |
Output is correct |
43 |
Correct |
949 ms |
25984 KB |
Output is correct |
44 |
Correct |
1354 ms |
33296 KB |
Output is correct |
45 |
Correct |
1341 ms |
33316 KB |
Output is correct |
46 |
Correct |
1360 ms |
33308 KB |
Output is correct |
47 |
Correct |
1184 ms |
29140 KB |
Output is correct |
48 |
Correct |
1178 ms |
29216 KB |
Output is correct |
49 |
Correct |
1022 ms |
25040 KB |
Output is correct |
50 |
Correct |
1004 ms |
25036 KB |
Output is correct |
51 |
Correct |
1132 ms |
28204 KB |
Output is correct |
52 |
Correct |
1296 ms |
32604 KB |
Output is correct |
53 |
Correct |
1318 ms |
32400 KB |
Output is correct |
54 |
Correct |
1270 ms |
32084 KB |
Output is correct |
55 |
Correct |
1259 ms |
32084 KB |
Output is correct |
56 |
Correct |
1287 ms |
32204 KB |
Output is correct |
57 |
Correct |
1318 ms |
32552 KB |
Output is correct |
58 |
Correct |
1309 ms |
32588 KB |
Output is correct |
59 |
Correct |
1306 ms |
32456 KB |
Output is correct |
60 |
Correct |
1300 ms |
32548 KB |
Output is correct |
61 |
Correct |
1307 ms |
32344 KB |
Output is correct |
62 |
Correct |
1299 ms |
32316 KB |
Output is correct |
63 |
Correct |
1303 ms |
32716 KB |
Output is correct |
64 |
Correct |
1313 ms |
32556 KB |
Output is correct |
65 |
Correct |
1309 ms |
32280 KB |
Output is correct |
66 |
Correct |
1306 ms |
32412 KB |
Output is correct |
67 |
Correct |
1267 ms |
32552 KB |
Output is correct |
68 |
Correct |
1305 ms |
32528 KB |
Output is correct |
69 |
Correct |
1309 ms |
32420 KB |
Output is correct |
70 |
Correct |
1394 ms |
32480 KB |
Output is correct |
71 |
Correct |
1382 ms |
32340 KB |
Output is correct |
72 |
Correct |
1389 ms |
32596 KB |
Output is correct |
73 |
Correct |
1391 ms |
32380 KB |
Output is correct |
74 |
Correct |
1 ms |
340 KB |
Output is correct |
75 |
Correct |
2 ms |
728 KB |
Output is correct |
76 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
596 KB |
Output is correct |
2 |
Correct |
2 ms |
724 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
596 KB |
Output is correct |
5 |
Correct |
1 ms |
596 KB |
Output is correct |
6 |
Correct |
1 ms |
660 KB |
Output is correct |
7 |
Correct |
1 ms |
468 KB |
Output is correct |
8 |
Correct |
0 ms |
340 KB |
Output is correct |
9 |
Correct |
2 ms |
468 KB |
Output is correct |
10 |
Correct |
2 ms |
596 KB |
Output is correct |
11 |
Correct |
2 ms |
596 KB |
Output is correct |
12 |
Correct |
2 ms |
724 KB |
Output is correct |
13 |
Correct |
1 ms |
724 KB |
Output is correct |
14 |
Correct |
1 ms |
724 KB |
Output is correct |
15 |
Correct |
2 ms |
724 KB |
Output is correct |
16 |
Correct |
1 ms |
724 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
214 ms |
6388 KB |
Output is correct |
19 |
Correct |
280 ms |
7576 KB |
Output is correct |
20 |
Correct |
281 ms |
7372 KB |
Output is correct |
21 |
Correct |
64 ms |
4948 KB |
Output is correct |
22 |
Correct |
268 ms |
7532 KB |
Output is correct |
23 |
Correct |
284 ms |
7580 KB |
Output is correct |
24 |
Correct |
115 ms |
3356 KB |
Output is correct |
25 |
Correct |
113 ms |
3412 KB |
Output is correct |
26 |
Correct |
252 ms |
7292 KB |
Output is correct |
27 |
Correct |
253 ms |
7244 KB |
Output is correct |
28 |
Correct |
221 ms |
6432 KB |
Output is correct |
29 |
Correct |
246 ms |
6384 KB |
Output is correct |
30 |
Correct |
225 ms |
6432 KB |
Output is correct |
31 |
Correct |
190 ms |
5588 KB |
Output is correct |
32 |
Correct |
184 ms |
5640 KB |
Output is correct |
33 |
Correct |
196 ms |
6292 KB |
Output is correct |
34 |
Correct |
262 ms |
7224 KB |
Output is correct |
35 |
Correct |
261 ms |
7220 KB |
Output is correct |
36 |
Correct |
266 ms |
7280 KB |
Output is correct |
37 |
Correct |
264 ms |
7184 KB |
Output is correct |
38 |
Correct |
260 ms |
7372 KB |
Output is correct |
39 |
Correct |
267 ms |
7120 KB |
Output is correct |
40 |
Correct |
253 ms |
7144 KB |
Output is correct |
41 |
Correct |
273 ms |
7232 KB |
Output is correct |
42 |
Correct |
272 ms |
7188 KB |
Output is correct |
43 |
Correct |
260 ms |
7244 KB |
Output is correct |
44 |
Correct |
272 ms |
7204 KB |
Output is correct |
45 |
Correct |
264 ms |
7192 KB |
Output is correct |
46 |
Correct |
256 ms |
7240 KB |
Output is correct |
47 |
Correct |
279 ms |
7624 KB |
Output is correct |
48 |
Correct |
261 ms |
7584 KB |
Output is correct |
49 |
Correct |
259 ms |
7632 KB |
Output is correct |
50 |
Correct |
0 ms |
340 KB |
Output is correct |
51 |
Correct |
1 ms |
724 KB |
Output is correct |
52 |
Correct |
1 ms |
340 KB |
Output is correct |
53 |
Correct |
1175 ms |
14796 KB |
Output is correct |
54 |
Correct |
1338 ms |
33296 KB |
Output is correct |
55 |
Correct |
1319 ms |
32540 KB |
Output is correct |
56 |
Correct |
493 ms |
20172 KB |
Output is correct |
57 |
Correct |
1331 ms |
33240 KB |
Output is correct |
58 |
Correct |
679 ms |
16812 KB |
Output is correct |
59 |
Correct |
667 ms |
16904 KB |
Output is correct |
60 |
Correct |
949 ms |
25984 KB |
Output is correct |
61 |
Correct |
1354 ms |
33296 KB |
Output is correct |
62 |
Correct |
1341 ms |
33316 KB |
Output is correct |
63 |
Correct |
1360 ms |
33308 KB |
Output is correct |
64 |
Correct |
1184 ms |
29140 KB |
Output is correct |
65 |
Correct |
1178 ms |
29216 KB |
Output is correct |
66 |
Correct |
1022 ms |
25040 KB |
Output is correct |
67 |
Correct |
1004 ms |
25036 KB |
Output is correct |
68 |
Correct |
1132 ms |
28204 KB |
Output is correct |
69 |
Correct |
1296 ms |
32604 KB |
Output is correct |
70 |
Correct |
1318 ms |
32400 KB |
Output is correct |
71 |
Correct |
1270 ms |
32084 KB |
Output is correct |
72 |
Correct |
1259 ms |
32084 KB |
Output is correct |
73 |
Correct |
1287 ms |
32204 KB |
Output is correct |
74 |
Correct |
1318 ms |
32552 KB |
Output is correct |
75 |
Correct |
1309 ms |
32588 KB |
Output is correct |
76 |
Correct |
1306 ms |
32456 KB |
Output is correct |
77 |
Correct |
1300 ms |
32548 KB |
Output is correct |
78 |
Correct |
1307 ms |
32344 KB |
Output is correct |
79 |
Correct |
1299 ms |
32316 KB |
Output is correct |
80 |
Correct |
1303 ms |
32716 KB |
Output is correct |
81 |
Correct |
1313 ms |
32556 KB |
Output is correct |
82 |
Correct |
1309 ms |
32280 KB |
Output is correct |
83 |
Correct |
1306 ms |
32412 KB |
Output is correct |
84 |
Correct |
1267 ms |
32552 KB |
Output is correct |
85 |
Correct |
1305 ms |
32528 KB |
Output is correct |
86 |
Correct |
1309 ms |
32420 KB |
Output is correct |
87 |
Correct |
1394 ms |
32480 KB |
Output is correct |
88 |
Correct |
1382 ms |
32340 KB |
Output is correct |
89 |
Correct |
1389 ms |
32596 KB |
Output is correct |
90 |
Correct |
1391 ms |
32380 KB |
Output is correct |
91 |
Correct |
1 ms |
340 KB |
Output is correct |
92 |
Correct |
2 ms |
728 KB |
Output is correct |
93 |
Correct |
1 ms |
340 KB |
Output is correct |
94 |
Correct |
1605 ms |
31484 KB |
Output is correct |
95 |
Correct |
1549 ms |
33292 KB |
Output is correct |
96 |
Correct |
1597 ms |
32672 KB |
Output is correct |
97 |
Correct |
861 ms |
17996 KB |
Output is correct |
98 |
Correct |
580 ms |
20076 KB |
Output is correct |
99 |
Correct |
1421 ms |
33256 KB |
Output is correct |
100 |
Correct |
835 ms |
16876 KB |
Output is correct |
101 |
Correct |
805 ms |
16712 KB |
Output is correct |
102 |
Correct |
1185 ms |
25940 KB |
Output is correct |
103 |
Correct |
1432 ms |
33228 KB |
Output is correct |
104 |
Correct |
1593 ms |
33236 KB |
Output is correct |
105 |
Correct |
1669 ms |
33232 KB |
Output is correct |
106 |
Correct |
1853 ms |
31576 KB |
Output is correct |
107 |
Correct |
1205 ms |
29136 KB |
Output is correct |
108 |
Correct |
1073 ms |
25036 KB |
Output is correct |
109 |
Correct |
1435 ms |
28712 KB |
Output is correct |
110 |
Correct |
1336 ms |
27800 KB |
Output is correct |
111 |
Correct |
1592 ms |
33260 KB |
Output is correct |
112 |
Correct |
1711 ms |
32808 KB |
Output is correct |
113 |
Correct |
1540 ms |
32088 KB |
Output is correct |
114 |
Correct |
1375 ms |
32088 KB |
Output is correct |
115 |
Correct |
1552 ms |
32088 KB |
Output is correct |
116 |
Correct |
1659 ms |
32780 KB |
Output is correct |
117 |
Correct |
1670 ms |
32808 KB |
Output is correct |
118 |
Correct |
1701 ms |
32824 KB |
Output is correct |
119 |
Correct |
1693 ms |
32724 KB |
Output is correct |
120 |
Correct |
1775 ms |
32840 KB |
Output is correct |
121 |
Correct |
1499 ms |
29136 KB |
Output is correct |
122 |
Correct |
1513 ms |
33296 KB |
Output is correct |
123 |
Correct |
1493 ms |
32492 KB |
Output is correct |
124 |
Correct |
1 ms |
336 KB |
Output is correct |
125 |
Correct |
2 ms |
728 KB |
Output is correct |
126 |
Correct |
2 ms |
336 KB |
Output is correct |