結果
問題 |
No.3175 転移迷宮 (Easy)
|
ユーザー |
![]() |
提出日時 | 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); | ~~~~~^~~~~~~~~~~~~~~~
ソースコード
/* -*- 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; }