Submission #532650

#	Time	Username	Problem	Language	Result	Execution time	Memory
532650		sidon	Genetics (BOI18_genetics)	C++17	74 / 100	2045 ms	57856 KiB

genetics

#pragma GCC optimize("O3", "unroll-loops")
#pragma GCC target("avx", "avx2", "popcnt")
#include <bits/stdc++.h>
using namespace std;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

const int Z = 4100;

struct Bitset {
	using num = unsigned long long;
	static const int B = 1<<6;

#define setBit(X, I) (X ^= (num(-1) ^ X) & (num(1)<<(I)))
#define blk(X) ((X)>>6)

	num v[blk(Z<<2)+1] {};

	void assign(const int &i) {
		setBit(v[blk(i)], (i & 63));
	}
	inline int countAnd(const Bitset &x) const {
		int res = 0;
		for(int i = blk(Z<<2); i >= 0; --i)
			res += __builtin_popcountll(v[i] & x.v[i]);

		return res;
	}
	void flip() {
		for(num &i : v) i = ~i;
	}
};
 
int N, M, K;
Bitset a[Z];
 
int g[100];
 
bool chk[Z][Z], cnt[Z][Z];

int main() {
	ios::sync_with_stdio(0), cin.tie(0);
 
	g['C'] = 1;
	g['G'] = 2;
	g['T'] = 3;
 
	cin >> N >> M >> K;

	int o[N];

	for(int i = 0; i < N; ++i) {
		string inp; cin >> inp;
		for(int j = 0; j < M; ++j)
			a[i].assign(g[int(inp[j])] * Z + j);
		o[i] = i;
	}

	shuffle(o, o + N, rng);

	for(const int &i : o) {
		bool ok = 1;
		a[i].flip();
 
		for(int j = 0; j < N; ++j) if(i != j) {

			bool c = (chk[i][j] ? cnt[i][j] : a[i].countAnd(a[j]) == K);
 
			chk[j][i] = 1;
			cnt[j][i] = c;
 
			if(!c) {
				ok = 0;
				break;
			}
		}
 
		if(ok) {
			cout << i + 1;
			return 0;
		}
		a[i].flip();
	}
}

• Subtask 1: 27 / 27
• Subtask 2: 19 / 19
• Subtask 3: 28 / 28
• Subtask 4: 0 / 26