답안 #540509

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
540509 2022-03-20T15:40:24 Z skittles1412 Genetics (BOI18_genetics) C++17
100 / 100
1771 ms 61260 KB
#pragma GCC optimize("Ofast,unroll-loops")
#pragma GCC target("sse4.2,bmi2,avx2")

#include "bits/extc++.h"

using namespace std;

template <typename T>
void dbgh(const T& t) {
	cerr << t << endl;
}

template <typename T, typename... U>
void dbgh(const T& t, const U&... u) {
	cerr << t << " | ";
	dbgh(u...);
}

#ifdef DEBUG
#define dbg(...)                                           \
	cerr << "L" << __LINE__ << " [" << #__VA_ARGS__ << "]" \
		 << ": ";                                          \
	dbgh(__VA_ARGS__)
#else
#define cerr   \
	if (false) \
	cerr
#define dbg(...)
#endif

using u64 = uint64_t;

#define endl "\n"
#define long int64_t
#define sz(x) int((x).size())

const int maxn = 4105, maxb = maxn / 64 + 1;
const long mod = 1e9 + 7;
const string chars = "ACGT";

struct BS {
	u64 arr[maxb];

	BS() : arr() {}

	void set(int ind) {
		arr[ind >> 6] |= u64(1) << (ind & 63);
	}

	int operator^(const BS& y) const {
		int ans = 0;
		for (int i = 0; i < maxb; i++) {
			ans += __builtin_popcountll(arr[i] ^ y.arr[i]);
		}
		return ans;
	}
};

BS arrb[maxn][4];

void solve() {
	int n, m, k;
	cin >> n >> m >> k;
	string arr[n];
	for (auto& a : arr) {
		cin >> a;
	}
	auto hsh = [&](int ind) -> long {
		long h1 = 0, h2 = 0;
		for (int i = 0; i < m; i++) {
			h1 = (h1 * 7 + int(chars.find(arr[ind][i])) + 1) % mod;
			h2 = (h2 * 11 + int(chars.find(arr[ind][i])) + 1) % mod;
		}
		return h1 * mod + h2;
	};
	mt19937 cowng(chrono::steady_clock::now().time_since_epoch().count());
	int ord[n], conv[n], cind = 0;
	iota(ord, ord + n, 0);
	shuffle(ord, ord + n, cowng);
	bool bad[n] {}, checked[n][n] {};
	map<long, int> vis;
	for (auto& a : ord) {
		auto [it, inserted] = vis.emplace(hsh(a), cind);
		if (!inserted) {
			bad[it->second] = true;
			continue;
		}
		for (int c = 0; c < 4; c++) {
			for (int j = 0; j < m; j++) {
				if (arr[a][j] == chars[c]) {
					arrb[cind][c].set(j);
				}
			}
		}
		conv[cind] = a;
		cind++;
	}
	k *= 2;
	for (int i = 0; i < cind; i++) {
		if (bad[i]) {
			continue;
		}
		for (int j = 0; j < cind; j++) {
			if (i == j || checked[i][j]) {
				continue;
			}
			int cur = 0;
			for (int c = 0; c < 4; c++) {
				cur += arrb[i][c] ^ arrb[j][c];
			}
			checked[i][j] = checked[j][i] = true;
			if (cur != k) {
				bad[i] = bad[j] = true;
				goto loop;
			}
		}
		cout << conv[i] + 1 << endl;
		return;
	loop:;
	}
}

int main() {
	cin.tie(nullptr);
	ios_base::sync_with_stdio(false);
	cin.exceptions(ios::failbit);
	solve();
}
# 결과 실행 시간 메모리 Grader output
1 Correct 5 ms 8660 KB Output is correct
2 Correct 5 ms 8660 KB Output is correct
3 Correct 4 ms 8660 KB Output is correct
4 Correct 5 ms 8676 KB Output is correct
5 Correct 5 ms 8660 KB Output is correct
6 Correct 6 ms 8660 KB Output is correct
7 Correct 5 ms 8636 KB Output is correct
8 Correct 5 ms 8648 KB Output is correct
9 Correct 5 ms 8688 KB Output is correct
10 Correct 5 ms 8660 KB Output is correct
11 Correct 5 ms 8660 KB Output is correct
12 Correct 5 ms 8684 KB Output is correct
13 Correct 5 ms 8660 KB Output is correct
14 Correct 6 ms 8784 KB Output is correct
15 Correct 5 ms 8684 KB Output is correct
16 Correct 5 ms 8660 KB Output is correct
17 Correct 5 ms 8660 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 126 ms 16360 KB Output is correct
2 Correct 67 ms 18532 KB Output is correct
3 Correct 118 ms 17924 KB Output is correct
4 Correct 32 ms 11380 KB Output is correct
5 Correct 142 ms 18528 KB Output is correct
6 Correct 127 ms 18528 KB Output is correct
7 Correct 58 ms 12228 KB Output is correct
8 Correct 60 ms 12108 KB Output is correct
9 Correct 100 ms 17712 KB Output is correct
10 Correct 95 ms 17788 KB Output is correct
11 Correct 203 ms 16260 KB Output is correct
12 Correct 101 ms 16416 KB Output is correct
13 Correct 113 ms 16428 KB Output is correct
14 Correct 85 ms 14976 KB Output is correct
15 Correct 75 ms 15092 KB Output is correct
16 Correct 90 ms 15524 KB Output is correct
17 Correct 195 ms 18148 KB Output is correct
18 Correct 140 ms 17948 KB Output is correct
19 Correct 105 ms 18028 KB Output is correct
20 Correct 194 ms 17812 KB Output is correct
21 Correct 162 ms 17996 KB Output is correct
22 Correct 135 ms 17896 KB Output is correct
23 Correct 268 ms 17996 KB Output is correct
24 Correct 122 ms 17944 KB Output is correct
25 Correct 167 ms 17888 KB Output is correct
26 Correct 101 ms 17956 KB Output is correct
27 Correct 245 ms 17868 KB Output is correct
28 Correct 202 ms 17848 KB Output is correct
29 Correct 105 ms 17936 KB Output is correct
30 Correct 73 ms 18496 KB Output is correct
31 Correct 76 ms 18716 KB Output is correct
32 Correct 101 ms 18552 KB Output is correct
33 Correct 5 ms 8660 KB Output is correct
34 Correct 5 ms 8688 KB Output is correct
35 Correct 4 ms 8660 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 126 ms 16360 KB Output is correct
2 Correct 67 ms 18532 KB Output is correct
3 Correct 118 ms 17924 KB Output is correct
4 Correct 32 ms 11380 KB Output is correct
5 Correct 142 ms 18528 KB Output is correct
6 Correct 127 ms 18528 KB Output is correct
7 Correct 58 ms 12228 KB Output is correct
8 Correct 60 ms 12108 KB Output is correct
9 Correct 100 ms 17712 KB Output is correct
10 Correct 95 ms 17788 KB Output is correct
11 Correct 203 ms 16260 KB Output is correct
12 Correct 101 ms 16416 KB Output is correct
13 Correct 113 ms 16428 KB Output is correct
14 Correct 85 ms 14976 KB Output is correct
15 Correct 75 ms 15092 KB Output is correct
16 Correct 90 ms 15524 KB Output is correct
17 Correct 195 ms 18148 KB Output is correct
18 Correct 140 ms 17948 KB Output is correct
19 Correct 105 ms 18028 KB Output is correct
20 Correct 194 ms 17812 KB Output is correct
21 Correct 162 ms 17996 KB Output is correct
22 Correct 135 ms 17896 KB Output is correct
23 Correct 268 ms 17996 KB Output is correct
24 Correct 122 ms 17944 KB Output is correct
25 Correct 167 ms 17888 KB Output is correct
26 Correct 101 ms 17956 KB Output is correct
27 Correct 245 ms 17868 KB Output is correct
28 Correct 202 ms 17848 KB Output is correct
29 Correct 105 ms 17936 KB Output is correct
30 Correct 73 ms 18496 KB Output is correct
31 Correct 76 ms 18716 KB Output is correct
32 Correct 101 ms 18552 KB Output is correct
33 Correct 5 ms 8660 KB Output is correct
34 Correct 5 ms 8688 KB Output is correct
35 Correct 4 ms 8660 KB Output is correct
36 Correct 786 ms 52976 KB Output is correct
37 Correct 334 ms 60792 KB Output is correct
38 Correct 885 ms 59628 KB Output is correct
39 Correct 285 ms 32672 KB Output is correct
40 Correct 761 ms 60540 KB Output is correct
41 Correct 417 ms 30488 KB Output is correct
42 Correct 374 ms 30488 KB Output is correct
43 Correct 351 ms 45292 KB Output is correct
44 Correct 508 ms 61216 KB Output is correct
45 Correct 644 ms 61116 KB Output is correct
46 Correct 553 ms 61208 KB Output is correct
47 Correct 1109 ms 52868 KB Output is correct
48 Correct 1072 ms 52992 KB Output is correct
49 Correct 448 ms 45044 KB Output is correct
50 Correct 696 ms 44960 KB Output is correct
51 Correct 710 ms 50356 KB Output is correct
52 Correct 716 ms 59620 KB Output is correct
53 Correct 695 ms 59468 KB Output is correct
54 Correct 317 ms 57964 KB Output is correct
55 Correct 358 ms 57960 KB Output is correct
56 Correct 486 ms 57964 KB Output is correct
57 Correct 749 ms 59776 KB Output is correct
58 Correct 1047 ms 59716 KB Output is correct
59 Correct 685 ms 59400 KB Output is correct
60 Correct 688 ms 59628 KB Output is correct
61 Correct 569 ms 59592 KB Output is correct
62 Correct 506 ms 59184 KB Output is correct
63 Correct 663 ms 59828 KB Output is correct
64 Correct 1771 ms 59596 KB Output is correct
65 Correct 992 ms 59408 KB Output is correct
66 Correct 833 ms 59428 KB Output is correct
67 Correct 878 ms 59812 KB Output is correct
68 Correct 1283 ms 59800 KB Output is correct
69 Correct 495 ms 59520 KB Output is correct
70 Correct 769 ms 59728 KB Output is correct
71 Correct 1484 ms 59508 KB Output is correct
72 Correct 588 ms 59852 KB Output is correct
73 Correct 541 ms 59464 KB Output is correct
74 Correct 5 ms 8660 KB Output is correct
75 Correct 6 ms 8660 KB Output is correct
76 Correct 4 ms 8676 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 5 ms 8660 KB Output is correct
2 Correct 5 ms 8660 KB Output is correct
3 Correct 4 ms 8660 KB Output is correct
4 Correct 5 ms 8676 KB Output is correct
5 Correct 5 ms 8660 KB Output is correct
6 Correct 6 ms 8660 KB Output is correct
7 Correct 5 ms 8636 KB Output is correct
8 Correct 5 ms 8648 KB Output is correct
9 Correct 5 ms 8688 KB Output is correct
10 Correct 5 ms 8660 KB Output is correct
11 Correct 5 ms 8660 KB Output is correct
12 Correct 5 ms 8684 KB Output is correct
13 Correct 5 ms 8660 KB Output is correct
14 Correct 6 ms 8784 KB Output is correct
15 Correct 5 ms 8684 KB Output is correct
16 Correct 5 ms 8660 KB Output is correct
17 Correct 5 ms 8660 KB Output is correct
18 Correct 126 ms 16360 KB Output is correct
19 Correct 67 ms 18532 KB Output is correct
20 Correct 118 ms 17924 KB Output is correct
21 Correct 32 ms 11380 KB Output is correct
22 Correct 142 ms 18528 KB Output is correct
23 Correct 127 ms 18528 KB Output is correct
24 Correct 58 ms 12228 KB Output is correct
25 Correct 60 ms 12108 KB Output is correct
26 Correct 100 ms 17712 KB Output is correct
27 Correct 95 ms 17788 KB Output is correct
28 Correct 203 ms 16260 KB Output is correct
29 Correct 101 ms 16416 KB Output is correct
30 Correct 113 ms 16428 KB Output is correct
31 Correct 85 ms 14976 KB Output is correct
32 Correct 75 ms 15092 KB Output is correct
33 Correct 90 ms 15524 KB Output is correct
34 Correct 195 ms 18148 KB Output is correct
35 Correct 140 ms 17948 KB Output is correct
36 Correct 105 ms 18028 KB Output is correct
37 Correct 194 ms 17812 KB Output is correct
38 Correct 162 ms 17996 KB Output is correct
39 Correct 135 ms 17896 KB Output is correct
40 Correct 268 ms 17996 KB Output is correct
41 Correct 122 ms 17944 KB Output is correct
42 Correct 167 ms 17888 KB Output is correct
43 Correct 101 ms 17956 KB Output is correct
44 Correct 245 ms 17868 KB Output is correct
45 Correct 202 ms 17848 KB Output is correct
46 Correct 105 ms 17936 KB Output is correct
47 Correct 73 ms 18496 KB Output is correct
48 Correct 76 ms 18716 KB Output is correct
49 Correct 101 ms 18552 KB Output is correct
50 Correct 5 ms 8660 KB Output is correct
51 Correct 5 ms 8688 KB Output is correct
52 Correct 4 ms 8660 KB Output is correct
53 Correct 786 ms 52976 KB Output is correct
54 Correct 334 ms 60792 KB Output is correct
55 Correct 885 ms 59628 KB Output is correct
56 Correct 285 ms 32672 KB Output is correct
57 Correct 761 ms 60540 KB Output is correct
58 Correct 417 ms 30488 KB Output is correct
59 Correct 374 ms 30488 KB Output is correct
60 Correct 351 ms 45292 KB Output is correct
61 Correct 508 ms 61216 KB Output is correct
62 Correct 644 ms 61116 KB Output is correct
63 Correct 553 ms 61208 KB Output is correct
64 Correct 1109 ms 52868 KB Output is correct
65 Correct 1072 ms 52992 KB Output is correct
66 Correct 448 ms 45044 KB Output is correct
67 Correct 696 ms 44960 KB Output is correct
68 Correct 710 ms 50356 KB Output is correct
69 Correct 716 ms 59620 KB Output is correct
70 Correct 695 ms 59468 KB Output is correct
71 Correct 317 ms 57964 KB Output is correct
72 Correct 358 ms 57960 KB Output is correct
73 Correct 486 ms 57964 KB Output is correct
74 Correct 749 ms 59776 KB Output is correct
75 Correct 1047 ms 59716 KB Output is correct
76 Correct 685 ms 59400 KB Output is correct
77 Correct 688 ms 59628 KB Output is correct
78 Correct 569 ms 59592 KB Output is correct
79 Correct 506 ms 59184 KB Output is correct
80 Correct 663 ms 59828 KB Output is correct
81 Correct 1771 ms 59596 KB Output is correct
82 Correct 992 ms 59408 KB Output is correct
83 Correct 833 ms 59428 KB Output is correct
84 Correct 878 ms 59812 KB Output is correct
85 Correct 1283 ms 59800 KB Output is correct
86 Correct 495 ms 59520 KB Output is correct
87 Correct 769 ms 59728 KB Output is correct
88 Correct 1484 ms 59508 KB Output is correct
89 Correct 588 ms 59852 KB Output is correct
90 Correct 541 ms 59464 KB Output is correct
91 Correct 5 ms 8660 KB Output is correct
92 Correct 6 ms 8660 KB Output is correct
93 Correct 4 ms 8676 KB Output is correct
94 Correct 683 ms 57632 KB Output is correct
95 Correct 354 ms 60564 KB Output is correct
96 Correct 600 ms 60116 KB Output is correct
97 Correct 516 ms 32224 KB Output is correct
98 Correct 314 ms 32676 KB Output is correct
99 Correct 627 ms 60716 KB Output is correct
100 Correct 361 ms 30716 KB Output is correct
101 Correct 342 ms 30444 KB Output is correct
102 Correct 606 ms 45636 KB Output is correct
103 Correct 660 ms 61204 KB Output is correct
104 Correct 813 ms 61220 KB Output is correct
105 Correct 523 ms 61260 KB Output is correct
106 Correct 813 ms 57932 KB Output is correct
107 Correct 943 ms 53168 KB Output is correct
108 Correct 452 ms 44852 KB Output is correct
109 Correct 820 ms 52184 KB Output is correct
110 Correct 728 ms 49624 KB Output is correct
111 Correct 554 ms 61260 KB Output is correct
112 Correct 934 ms 60584 KB Output is correct
113 Correct 419 ms 58224 KB Output is correct
114 Correct 382 ms 58152 KB Output is correct
115 Correct 476 ms 58044 KB Output is correct
116 Correct 1334 ms 60240 KB Output is correct
117 Correct 636 ms 60460 KB Output is correct
118 Correct 818 ms 60236 KB Output is correct
119 Correct 511 ms 60084 KB Output is correct
120 Correct 1315 ms 60460 KB Output is correct
121 Correct 863 ms 52988 KB Output is correct
122 Correct 331 ms 60516 KB Output is correct
123 Correct 1118 ms 59620 KB Output is correct
124 Correct 4 ms 8660 KB Output is correct
125 Correct 5 ms 8660 KB Output is correct
126 Correct 4 ms 8660 KB Output is correct