Submission #990557

# Submission time Handle Problem Language Result Execution time Memory
990557 2024-05-30T16:07:00 Z jakobrs Ancient Machine 2 (JOI23_ancient2) C++17
97 / 100
52 ms 904 KB
#include "ancient2.h"

#include <bitset>
#include <memory>
#include <numeric>
#include <string>
#include <vector>

void gauss(std::bitset<2000> (&matrix)[1000]) {
  for (int col = 0; col < 1000; col++) {
    int pivot_row;
    for (int row = col; row < 1000; row++) {
      if (matrix[row][col]) {
        pivot_row = row;
        break;
      }
    }
    std::swap(matrix[col], matrix[pivot_row]);
    for (int row = col + 1; row < 1000; row++)
      if (matrix[row][col])
        matrix[row] ^= matrix[col];
  }
  for (int col = 999; col >= 0; col--) {
    for (int row = col - 1; row >= 0; row--)
      if (matrix[row][col])
        matrix[row] ^= matrix[col];
  }
}

std::string Solve(int N) {
  int m = 57 * 2;
  std::vector<int> a(m), b(m);
  std::vector<bool> results;
  int acc_phi = 0;

  auto matrix = std::make_unique<std::bitset<2000>[]>(1000);
  for (int i = 1; i < 58; i++) {
    for (int j = 0; j < i; j++) {
      a[j] = j + 1;
      a[j + i] = j + i + 1;
      b[j] = j + 1;
      b[j + i] = j + i + 1;
    }
    a[i - 1] = 0;
    a[i - 1 + i] = i;
    b[i - 1] = 0;
    b[i - 1 + i] = i;

    int phi = 0;
    for (int j = 1; j <= i; j++)
      if (std::gcd(i, j) == 1)
        phi += 1;

    for (int j = 0; j < phi; j++) {
      for (int k = 0; k < 1000; k++) {
        if (k % i == j) {
          matrix[results.size()][k] = 1;
        }
      }
      matrix[results.size()][1000 + results.size()] = 1;

      std::swap(b[j], b[j + i]);
      auto res = Query(m, a, b);
      results.push_back(res >= i);
      std::swap(b[j], b[j + i]);
    }
  }

  gauss(*(std::bitset<2000>(*)[1000])matrix.get());

  std::vector<bool> answer(1000);

  for (int i = 0; i < 1000; i++)
    if (results[i])
      for (int j = 0; j < 1000; j++)
        if (matrix[j][i + 1000])
          answer[j] = !answer[j];

  std::string s;

  for (int i = 0; i < 1000; i++)
    s.push_back(answer[i] ? '1' : '0');

  s.erase(s.begin() + N, s.end());

  return s;
}

Compilation message

ancient2.cpp: In function 'std::string Solve(int)':
ancient2.cpp:34:7: warning: unused variable 'acc_phi' [-Wunused-variable]
   34 |   int acc_phi = 0;
      |       ^~~~~~~
# Verdict Execution time Memory Grader output
1 Partially correct 50 ms 600 KB Output is partially correct
2 Partially correct 50 ms 664 KB Output is partially correct
3 Partially correct 30 ms 600 KB Output is partially correct
4 Partially correct 46 ms 660 KB Output is partially correct
5 Partially correct 32 ms 600 KB Output is partially correct
6 Partially correct 31 ms 692 KB Output is partially correct
7 Partially correct 33 ms 600 KB Output is partially correct
8 Partially correct 52 ms 600 KB Output is partially correct
9 Partially correct 31 ms 684 KB Output is partially correct
10 Partially correct 37 ms 664 KB Output is partially correct
11 Partially correct 45 ms 600 KB Output is partially correct
12 Partially correct 47 ms 692 KB Output is partially correct
13 Partially correct 33 ms 600 KB Output is partially correct
14 Partially correct 49 ms 600 KB Output is partially correct
15 Partially correct 31 ms 600 KB Output is partially correct
16 Partially correct 35 ms 856 KB Output is partially correct
17 Partially correct 50 ms 600 KB Output is partially correct
18 Partially correct 31 ms 692 KB Output is partially correct
19 Partially correct 40 ms 692 KB Output is partially correct
20 Partially correct 32 ms 856 KB Output is partially correct
21 Partially correct 36 ms 856 KB Output is partially correct
22 Partially correct 36 ms 848 KB Output is partially correct
23 Partially correct 30 ms 600 KB Output is partially correct
24 Partially correct 50 ms 688 KB Output is partially correct
25 Partially correct 32 ms 688 KB Output is partially correct
26 Partially correct 33 ms 600 KB Output is partially correct
27 Partially correct 42 ms 600 KB Output is partially correct
28 Partially correct 30 ms 692 KB Output is partially correct
29 Partially correct 42 ms 600 KB Output is partially correct
30 Partially correct 29 ms 600 KB Output is partially correct
31 Partially correct 32 ms 600 KB Output is partially correct
32 Partially correct 39 ms 856 KB Output is partially correct
33 Partially correct 36 ms 600 KB Output is partially correct
34 Partially correct 38 ms 668 KB Output is partially correct
35 Partially correct 41 ms 600 KB Output is partially correct
36 Partially correct 41 ms 600 KB Output is partially correct
37 Partially correct 29 ms 600 KB Output is partially correct
38 Partially correct 38 ms 600 KB Output is partially correct
39 Partially correct 50 ms 596 KB Output is partially correct
40 Partially correct 30 ms 668 KB Output is partially correct
41 Partially correct 35 ms 852 KB Output is partially correct
42 Partially correct 38 ms 600 KB Output is partially correct
43 Partially correct 31 ms 600 KB Output is partially correct
44 Partially correct 31 ms 668 KB Output is partially correct
45 Partially correct 32 ms 688 KB Output is partially correct
46 Partially correct 47 ms 688 KB Output is partially correct
47 Partially correct 30 ms 600 KB Output is partially correct
48 Partially correct 29 ms 600 KB Output is partially correct
49 Partially correct 42 ms 600 KB Output is partially correct
50 Partially correct 51 ms 856 KB Output is partially correct
51 Partially correct 43 ms 600 KB Output is partially correct
52 Partially correct 39 ms 600 KB Output is partially correct
53 Partially correct 39 ms 600 KB Output is partially correct
54 Partially correct 36 ms 600 KB Output is partially correct
55 Partially correct 37 ms 692 KB Output is partially correct
56 Partially correct 33 ms 600 KB Output is partially correct
57 Partially correct 34 ms 692 KB Output is partially correct
58 Partially correct 32 ms 600 KB Output is partially correct
59 Partially correct 35 ms 600 KB Output is partially correct
60 Partially correct 33 ms 600 KB Output is partially correct
61 Partially correct 36 ms 600 KB Output is partially correct
62 Partially correct 29 ms 600 KB Output is partially correct
63 Partially correct 44 ms 904 KB Output is partially correct