Submission #532665

#	Time	Username	Problem	Language	Result	Execution time	Memory
532665		sidon	Genetics (BOI18_genetics)	C++17	46 / 100	2075 ms	48832 KiB

genetics

#pragma GCC optimize("O3", "unroll-loops")
#pragma GCC target("avx2", "popcnt", "sse4.2")
#include <bits/stdc++.h>
using namespace std;
 
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
 
const int Z = 4100;
using num = unsigned long long;
struct Bitset {
	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 b[N], 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;
		b[i] = i;
	}

	shuffle(o, o + N, rng);
 
 
	for(const int &i : o) {
		bool ok = 1;
		a[i].flip();


		shuffle(b, b + N, rng);

 
		for(const int &j : b) 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: 0 / 28
• Subtask 4: 0 / 26