Submission #871440

# Submission time Handle Problem Language Result Execution time Memory
871440 2023-11-10T19:58:03 Z rainboy Secret Permutation (RMI19_permutation) C++17
99.9181 / 100
825 ms 253000 KB
#include "permutation.h"
#include <iostream>
#include <cstring>
#include <vector>
 
using namespace std;
 
typedef vector<int> vi;
 
const int N = 495, M = 131072;
const long long INF = 0x3f3f3f3f3f3f3f3fLL;
 
int abs_(int a) { return a > 0 ? a : -a; }
int min(int a, int b) { return a < b ? a : b; }
 
unsigned int X = 12345;
 
int rand_() {
	return (X *= 3) >> 1;
}
 
char used[N]; int kk[N];
 
int check(short *aa, int n, int m) {
	int a = 0;
	for (int i = 0; i < m; i++)
		a = min(a, aa[i]);
	memset(used, 0, n * sizeof *used);
	for (int i = 0; i < m; i++) {
		if (aa[i] - a >= n || used[aa[i] - a])
			return 0;
		used[aa[i] - a] = 1;
	}
	return 1;
}
 
int rr[N], dd[N];
short aaa[M][N], aaa_[M][N];
 
void solve(int n) {
	for (int i = 0; i < n; i++)
		rr[i] = i;
	for (int i = 0; i < n; i++) {
		int j = rand_() % (i + 1), tmp;
		tmp = rr[i], rr[i] = rr[j], rr[j] = tmp;
	}
	vi ii(n);
	for (int i = 0; i < n; i++) {
		for (int h = 0; h < n; h++)
			ii[h] = rr[(i + h) % n] + 1;
		dd[i] = query(ii);
	}
	int sum = 0;
	for (int i = 0; i < n; i++)
		sum += dd[i];
	sum /= n - 1;
	for (int i = 0; i < n; i++)
		dd[i] = sum - dd[i];
	vi aa(n), aa_(n);
	int m, m_;
	m = 0;
	aaa[m][0] = 0, aaa[m][1] = dd[1], m++;
	int j = 2;
	while (1) {
		while (j < n) {
			m_ = 0;
			for (int h = 0; h < m; h++) {
				memcpy(aaa_[m_], aaa[h], j * sizeof *aaa[h]), aaa_[m_][j] = aaa[h][j - 1] + dd[j];
				if (check(aaa_[m_], n, j + 1))
					m_++;
				if (m_ >= M)
					goto out;
				memcpy(aaa_[m_], aaa[h], j * sizeof *aaa[h]), aaa_[m_][j] = aaa[h][j - 1] - dd[j];
				if (check(aaa_[m_], n, j + 1))
					m_++;
				if (m_ >= M)
					goto out;
			}
			m = m_;
			for (int h = 0; h < m_; h++)
				memcpy(aaa[h], aaa_[h], (j + 1) * sizeof *aaa_[h]);
			j++;
		}
out:
		if (j == n) {
			m_ = 0;
			for (int h = 0; h < m; h++)
				if (abs_(aaa[h][n - 1]) == dd[0])
					memcpy(aaa_[m_++], aaa[h], n * sizeof *aaa[h]);
			m = m_;
			for (int h = 0; h < m_; h++)
				memcpy(aaa[h], aaa_[h], n * sizeof *aaa_[h]);
			if (m == 1) {
				for (int i = 0; i < n; i++)
					aa[i] = aaa[0][i];
				break;
			}
		}
		long long z_ = INF; int i_ = -1;
		for (int i = 0; i < j; i++) {
			memset(kk, 0, n * sizeof *kk);
			for (int h = 0; h < m; h++)
				kk[abs_(aaa[h][i])]++;
			long long z = 0;
			for (int a = 0; a < n; a++)
				z += (long long) kk[a] * kk[a];
			if (z_ > z)
				z_ = z, i_ = i;
		}
		cerr << i_ << endl;
		for (int h = 0; h < n; h++)
			ii[h] = rr[h < i_ ? i_ - 1 - h : h] + 1;
		int a = query(ii) - (sum - dd[0]) + dd[i_];
		m_ = 0;
		for (int h = 0; h < m; h++)
			if (abs_(aaa[h][i_]) == a)
				memcpy(aaa_[m_++], aaa[h], n * sizeof *aaa[h]);
		m = m_;
		for (int h = 0; h < m_; h++)
			memcpy(aaa[h], aaa_[h], n * sizeof *aaa_[h]);
	}
	int a = 0;
	for (int i = 0; i < n; i++)
		a = min(a, aa[i]);
	for (int i = 0; i < n; i++)
		aa[i] -= a - 1;
	for (int i = 0; i < n; i++)
		aa_[rr[i]] = aa[i];
	for (int i = 0; i < n; i++)
		aa[i] = aa_[i];
	answer(aa);
}

Compilation message

stub.cpp: In function 'int query(int*)':
stub.cpp:15:9: warning: ignoring return value of 'int fscanf(FILE*, const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   15 |   fscanf(stdin, "%d", &x);
      |   ~~~~~~^~~~~~~~~~~~~~~~~
stub.cpp: In function 'int main(int, char**)':
stub.cpp:48:9: warning: ignoring return value of 'int fscanf(FILE*, const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   48 |   fscanf(stdin, "%d", &N);
      |   ~~~~~~^~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2500 KB Output is correct
2 Correct 0 ms 2496 KB Output is correct
3 Correct 1 ms 2496 KB Output is correct
4 Correct 0 ms 2496 KB Output is correct
5 Correct 1 ms 2504 KB Output is correct
6 Correct 1 ms 2500 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2500 KB Output is correct
2 Correct 0 ms 2496 KB Output is correct
3 Correct 1 ms 2496 KB Output is correct
4 Correct 0 ms 2496 KB Output is correct
5 Correct 1 ms 2504 KB Output is correct
6 Correct 1 ms 2500 KB Output is correct
7 Correct 1 ms 4548 KB Output is correct
8 Correct 1 ms 2500 KB Output is correct
9 Correct 1 ms 2500 KB Output is correct
10 Correct 2 ms 4544 KB Output is correct
11 Correct 1 ms 2500 KB Output is correct
12 Correct 1 ms 4552 KB Output is correct
13 Correct 1 ms 2496 KB Output is correct
14 Correct 1 ms 2500 KB Output is correct
15 Correct 1 ms 6592 KB Output is correct
16 Correct 1 ms 4544 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2500 KB Output is correct
2 Correct 0 ms 2496 KB Output is correct
3 Correct 1 ms 2496 KB Output is correct
4 Correct 0 ms 2496 KB Output is correct
5 Correct 1 ms 2504 KB Output is correct
6 Correct 1 ms 2500 KB Output is correct
7 Correct 1 ms 4548 KB Output is correct
8 Correct 1 ms 2500 KB Output is correct
9 Correct 1 ms 2500 KB Output is correct
10 Correct 2 ms 4544 KB Output is correct
11 Correct 1 ms 2500 KB Output is correct
12 Correct 1 ms 4552 KB Output is correct
13 Correct 1 ms 2496 KB Output is correct
14 Correct 1 ms 2500 KB Output is correct
15 Correct 1 ms 6592 KB Output is correct
16 Correct 1 ms 4544 KB Output is correct
17 Correct 685 ms 191432 KB Output is correct
18 Correct 10 ms 6584 KB Output is correct
19 Correct 410 ms 162476 KB Output is correct
20 Correct 176 ms 101052 KB Output is correct
21 Partially correct 512 ms 253000 KB Partially correct
22 Partially correct 269 ms 246704 KB Partially correct
23 Partially correct 245 ms 250812 KB Partially correct
24 Correct 825 ms 252844 KB Output is correct
25 Partially correct 408 ms 218112 KB Partially correct
26 Partially correct 191 ms 218112 KB Partially correct