結果

問題 No.3175 転移迷宮 (Easy)
ユーザー tnakao0123
提出日時 2025-06-07 18:16:37
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 85 ms / 3,000 ms
コード長 2,608 bytes
コンパイル時間 583 ms
コンパイル使用メモリ 44,416 KB
実行使用メモリ 6,272 KB
最終ジャッジ日時 2025-06-07 18:16:41
合計ジャッジ時間 4,175 ms
ジャッジサーバーID
(参考情報)
judge1 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 44
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:69:8: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   69 |   scanf("%d", &tn);
      |   ~~~~~^~~~~~~~~~~
main.cpp:73:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   73 |     scanf("%d%d", &n, &k);
      |     ~~~~~^~~~~~~~~~~~~~~~

ソースコード

diff #

/* -*- coding: utf-8 -*-
 *
 * 3175.cc:  No.3175 霆「遘サ霑キ螳ョ (Easy) - yukicoder
 */

#include<cstdio>
#include<algorithm>

using namespace std;

/* constant */

const int MAX_N = 500;
const int MOD = 998244353;

/* typedef */

template<const int MOD>
struct MI {
  int v;
  MI(): v() {}
  MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; }
  MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; }

  explicit operator int() const { return v; }
  
  MI operator+(const MI m) const { return MI(v + m.v); }
  MI operator-(const MI m) const { return MI(v + MOD - m.v); }
  MI operator-() const { return MI(MOD - v); }
  MI operator*(const MI m) const { return MI((long long)v * m.v); }

  MI &operator+=(const MI m) { return (*this = *this + m); }
  MI &operator-=(const MI m) { return (*this = *this - m); }
  MI &operator*=(const MI m) { return (*this = *this * m); }

  bool operator==(const MI m) const { return v == m.v; }
  bool operator!=(const MI m) const { return v != m.v; }

  MI pow(int n) const {  // a^n % MOD
    MI pm = 1, a = *this;
    while (n > 0) {
      if (n & 1) pm *= a;
      a *= a;
      n >>= 1;
    }
    return pm;
  }

  MI inv() const { return pow(MOD - 2); }
  MI operator/(const MI m) const { return *this * m.inv(); }
  MI &operator/=(const MI m) { return (*this = *this / m); }
};

using mi = MI<MOD>;
using vec = mi[MAX_N];
using mat = vec[MAX_N];

/* global variables */

mat ma;
vec va;

/* subroutines */

/* main */

int main() {
  int tn;
  scanf("%d", &tn);

  while (tn--) {
    int n, k;
    scanf("%d%d", &n, &k);

    // xi = 1 + sum_i(xi)/(2k+1)

    for (int i = 0; i < n; i++) {
      fill(ma[i], ma[i] + n, 0);
      for (int j = i - k; j <= i + k; j++)
	if (j >= 0 && j < n) ma[i][j] = -1;
      ma[i][i] = k * 2;
      va[i] = k * 2 + 1;
    }

    for (int i = 0; i < n; i++) {
      if (ma[i][i] == 0) {
	for (int i0 = i + 1; i0 < n; i0++)
	  if (ma[i0][i] != 0) {
	    for (int j = i; j < n; j++) swap(ma[i][j], ma[i0][j]);
	    swap(va[i], va[i0]);
	    break;
	  }
      }

      mi invmii = ma[i][i].inv();
      for (int j = i; j < n; j++) ma[i][j] *= invmii;
      va[i] *= invmii;
      
      for (int i0 = 0; i0 < n; i0++)
	if (i0 != i && ma[i0][i] != 0) {
	  mi d = ma[i0][i];
	  for (int j = i; j < n; j++) ma[i0][j] -= d * ma[i][j];
	  va[i0] -= d * va[i];
	}
    }

    //for (int i = 0; i < n; i++) {
    //  for (int j = 0; j < n; j++) printf(" %d", ma[i][j]);
    //  printf(" = %d\n", va[i]);
    //} putchar('\n');

    for (int i = 0; i < n; i++)
      printf("%d%c", (int)va[i], (i + 1 < n) ? ' ' : '\n');
  }

  return 0;
}
0