#1168380 (C++23) No.3567 Modulo Grid

• 問題: No.3567 Modulo Grid
• ユーザー: risujiroh
• 提出日時: 2026-06-05 23:36:35
• 言語: C++23 (gcc 15.2.0 + boost 1.89.0)
• 結果: AC
• 実行時間: 27 ms / 2,000 ms
• コード長: 4,659 bytes
• コンパイル時間: 3,000 ms
• コンパイル使用メモリ: 362,356 KB
• 実行使用メモリ: 6,528 KB
• 最終ジャッジ日時: 2026-06-05 23:36:47
• 合計ジャッジ時間: 7,671 ms
• ジャッジサーバーID （参考情報）: judge3_0 / judge1_0

ファイルパターン	結果
sample	AC * 2
other	AC * 25

ソースコード

#if __INCLUDE_LEVEL__ == 0

#include __BASE_FILE__

Array2D<int> Solve(int h, int w) {
  if (h == 1 && w == 1)
    return Array2D<int>(h, w);

  int p = 2;
  while (h % p && w % p)
    ++p;

  int Bh = h;
  int Ch = 1;
  while (Bh % p == 0) {
    Bh /= p;
    Ch *= p;
  }

  int Bw = w;
  int Cw = 1;
  while (Bw % p == 0) {
    Bw /= p;
    Cw *= p;
  }

  Array2D<int> rem(Ch, Cw);
  for (int i : Rep(0, Ch))
    for (int j : Rep(0, Cw))
      rem[i][j] = i * Cw + j;
  {
    vector<pair<int, int>> order_h(Ch);
    for (int i : Rep(0, Ch)) {
      if (i == 0) {
        order_h[i].first = INF;
      } else {
        int t = i;
        while (t % p == 0) {
          t /= p;
          ++order_h[i].first;
        }
      }
      order_h[i].second = i;
    }
    ranges::sort(order_h, greater{});
    vector<pair<int, int>> order_w(Cw);
    for (int i : Rep(0, Cw)) {
      if (i == 0) {
        order_w[i].first = INF;
      } else {
        int t = i;
        while (t % p == 0) {
          t /= p;
          ++order_w[i].first;
        }
      }
      order_w[i].second = i;
    }
    ranges::sort(order_w, greater{});
    Array2D<int> new_rem(Ch, Cw);
    for (int i : Rep(0, Ch))
      for (int j : Rep(0, Cw))
        new_rem[i][j] = rem[order_h[i].second][order_w[j].second];
    rem = move(new_rem);
  }

  auto tmp = Solve(Bh, Bw);
  Array2D<int> ret(h, w);
  for (int i : Rep(0, h))
    for (int j : Rep(0, w)) {
      int v = tmp[i / Bh & 1 ? Bh - i % Bh - 1 : i % Bh][j / Bw & 1 ? Bw - j % Bw - 1 : j % Bw];
      ret[i][j] = int(atcoder::crt({rem[i / Bh][j / Bw], v}, {Ch * Cw, Bh * Bw}).first);
    }
  return ret;
}

void Solve() {
  int h, w;
  IN(h, w);
  auto ans = Solve(h, w);
  for (int i : Rep(0, h)) {
    ranges::replace(ans[i], 0, h * w);
    OUT(ans[i]);
  }
}

int main() {
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  Solve();
}

#elif __INCLUDE_LEVEL__ == 1

#include <bits/stdc++.h>

#include <atcoder/math.hpp>

template <class T>
class Array2D {
 public:
  Array2D() : Array2D(0, 0) {}
  Array2D(int h, int w, T v = {}) : h_(h), w_(w), data_(h * w, std::move(v)) {}

  auto operator[](int i) { return std::span<T>(data_.data() + i * w_, w_); }
  auto operator[](int i) const { return std::span<const T>(data_.data() + i * w_, w_); }

 private:
  int h_;
  int w_;
  std::vector<T> data_;
};

template <class T>
concept MyRange =
    std::ranges::range<T> &&
    !std::convertible_to<T, std::string_view> &&
    !std::convertible_to<T, std::filesystem::path>;

template <class T>
concept MyTuple = std::__is_tuple_like<T>::value && !MyRange<T>;

namespace std {

istream& operator>>(istream& is, MyRange auto&& r) {
  for (auto&& e : r) {
    is >> e;
  }
  return is;
}

istream& operator>>(istream& is, MyTuple auto&& t) {
  apply([&](auto&... xs) { (is >> ... >> xs); }, t);
  return is;
}

ostream& operator<<(ostream& os, MyRange auto&& r) {
  auto sep = "";
  for (auto&& e : r) {
    os << exchange(sep, " ") << forward<decltype(e)>(e);
  }
  return os;
}

ostream& operator<<(ostream& os, MyTuple auto&& t) {
  auto sep = "";
  apply([&](auto&... xs) { ((os << exchange(sep, " ") << xs), ...); }, t);
  return os;
}

}  // namespace std

template <class T>
class OneBased {
 public:
  explicit OneBased(T&& x) : ref_(std::forward<T>(x)) {}

  template <class... Ts>
  requires(sizeof...(Ts) > 1)
      OneBased(Ts&&... xs) : ref_(std::forward_as_tuple(std::forward<Ts>(xs)...)) {}

  friend std::istream& operator>>(std::istream& is, OneBased x) {
    if constexpr (MyRange<T>) {
      for (auto&& e : x.ref_) {
        is >> ::OneBased(e);
      }
    } else if constexpr (MyTuple<T>) {
      std::apply([&](auto&... xs) { (is >> ... >> ::OneBased(xs)); }, x.ref_);
    } else {
      is >> x.ref_;
      --x.ref_;
    }
    return is;
  }

  friend std::ostream& operator<<(std::ostream& os, OneBased x) {
    if constexpr (MyRange<T>) {
      auto f = [](auto&& e) { return ::OneBased(std::forward<decltype(e)>(e)); };
      os << (x.ref_ | std::views::transform(f));
    } else if constexpr (MyTuple<T>) {
      std::apply([&](auto&... xs) { os << std::tuple(::OneBased(xs)...); }, x.ref_);
    } else {
      os << ++x.ref_;
      --x.ref_;
    }
    return os;
  }

 private:
  T ref_;
};

template <class T>
OneBased(T&&) -> OneBased<T>;

template <class... Ts>
OneBased(Ts&&...) -> OneBased<std::tuple<Ts...>>;

using namespace std;

#define Rep(...) [](int l, int r) { return views::iota(min(l, r), r); }(__VA_ARGS__)
#define INF (INT_MAX / 2)
#define IN(...) (cin >> forward_as_tuple(__VA_ARGS__))
#define OUT(...) (cout << forward_as_tuple(__VA_ARGS__) << '\n')

#endif  // __INCLUDE_LEVEL__ == 1