Submission #296841

# Submission time Handle Problem Language Result Execution time Memory
296841 2020-09-11T01:17:11 Z eggag32 Library (JOI18_library) C++17
19 / 100
2000 ms 4472 KB
#pragma GCC optimize ("O3")
#pragma GCC target ("sse4")
#include <bits/stdc++.h>
#include "library.h"
using namespace std;
typedef long long ll;
typedef long double ld;
typedef vector<int> vi;
typedef pair<int, int> pi;
#define debug(x) cerr << #x << ": " << x << endl
#define debug2(x, y) debug(x), debug(y)
#define repn(i, a, b) for(int i = (int)(a); i < (int)(b); i++)
#define rep(i, a) for(int i = 0; i < (int)(a); i++)
#define all(v) v.begin(), v.end() 
#define mp make_pair
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define fi first
#define se second
#define sq(x) ((x) * (x))
const int mxN = 1005;

template<class T> T gcd(T a, T b){ return ((b == 0) ? a : gcd(b, a % b)); }

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

int pos[mxN][mxN];
vi g[mxN];

void Solve(int n){
	if(n == 1){
		vi ans = {1};
		Answer(ans);
		return;
	}
	int sz = n / 10; //the size of our guess bucket
	int num = (n * (n - 1)) / 2;
	int rem = 20000;
	memset(pos, -1, sizeof(pos));
	while(num > rem){
		vi vis(n, 0), cr;
		rep(k, sz){
			while(true){
				int ind = rng() % n;
				if(!vis[ind]){
					vis[ind] = 1;
					cr.pb(ind);
					break;
				}
			}
		}
		int ans = Query(vis);
		if(ans == sz){
			rep(i, cr.size()) rep(j, cr.size()) if(i != j && pos[cr[i]][cr[j]] == -1){
				num--;
				pos[cr[i]][cr[j]] = 0;
				pos[cr[j]][cr[i]] = 0;
			}		
		}
		rem--;
	}
	rep(i, n) rep(j, n) if(i != j && pos[i][j] == -1){
		vi nw(n, 0);
		nw[i] = nw[j] = 1;
		int ans = Query(nw);
		if(ans == 1) pos[i][j] = pos[j][i] = 1;
		else pos[i][j] = pos[j][i] = 0;
	}
	//now we know all the comparisons so we can just reconstruct the answer
	int st = -1;
	rep(i, n){
		rep(j, n){
			if(pos[i][j] == 1) g[i].pb(j);
		}
		if(g[i].size() == 1) st = i;
	}
	vi ans;
	ans.pb(st + 1);
	int nxt = g[st][0];
	repn(i, 1, n){
		ans.pb(nxt + 1);
		rep(j, g[nxt].size()){
			if((ans[(int)ans.size() - 2] - 1) != g[nxt][j]){
				nxt = g[nxt][j];
				break;
			}
		}
	}
	Answer(ans);
	return;
}
# Verdict Execution time Memory Grader output
1 Correct 233 ms 4352 KB # of queries: 18336
2 Correct 248 ms 4472 KB # of queries: 18145
3 Correct 315 ms 4352 KB # of queries: 19900
4 Correct 245 ms 4472 KB # of queries: 19900
5 Correct 292 ms 4352 KB # of queries: 19900
6 Correct 305 ms 4352 KB # of queries: 19900
7 Correct 309 ms 4448 KB # of queries: 19900
8 Correct 271 ms 4344 KB # of queries: 18528
9 Correct 268 ms 4472 KB # of queries: 19701
10 Correct 115 ms 4352 KB # of queries: 8256
11 Correct 1 ms 384 KB # of queries: 0
12 Correct 3 ms 4352 KB # of queries: 1
13 Correct 3 ms 4224 KB # of queries: 3
14 Correct 3 ms 4352 KB # of queries: 6
15 Correct 4 ms 4352 KB # of queries: 105
16 Correct 8 ms 4352 KB # of queries: 351
# Verdict Execution time Memory Grader output
1 Correct 233 ms 4352 KB # of queries: 18336
2 Correct 248 ms 4472 KB # of queries: 18145
3 Correct 315 ms 4352 KB # of queries: 19900
4 Correct 245 ms 4472 KB # of queries: 19900
5 Correct 292 ms 4352 KB # of queries: 19900
6 Correct 305 ms 4352 KB # of queries: 19900
7 Correct 309 ms 4448 KB # of queries: 19900
8 Correct 271 ms 4344 KB # of queries: 18528
9 Correct 268 ms 4472 KB # of queries: 19701
10 Correct 115 ms 4352 KB # of queries: 8256
11 Correct 1 ms 384 KB # of queries: 0
12 Correct 3 ms 4352 KB # of queries: 1
13 Correct 3 ms 4224 KB # of queries: 3
14 Correct 3 ms 4352 KB # of queries: 6
15 Correct 4 ms 4352 KB # of queries: 105
16 Correct 8 ms 4352 KB # of queries: 351
17 Execution timed out 3040 ms 4352 KB Time limit exceeded
18 Halted 0 ms 0 KB -