結果
問題 |
No.3175 転移迷宮 (Easy)
|
ユーザー |
|
提出日時 | 2025-06-12 14:02:22 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 215 ms / 3,000 ms |
コード長 | 1,896 bytes |
コンパイル時間 | 3,175 ms |
コンパイル使用メモリ | 281,156 KB |
実行使用メモリ | 6,272 KB |
最終ジャッジ日時 | 2025-06-12 14:03:19 |
合計ジャッジ時間 | 9,329 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 44 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; using pii = pair<int, int>; using pll = pair<ll, ll>; using vi = vector<int>; using vl = vector<ll>; #define rep3(i, a, b, c) for (ll i = (a); i < (b); i += (c)) #define rep2(i, a, b) rep3(i, a, b, 1) #define rep1(i, n) rep2(i, 0, n) #define rep0(n) rep1(aaaaa, n) #define ov4(a, b, c, d, name, ...) name #define rep(...) ov4(__VA_ARGS__, rep3, rep2, rep1, rep0)(__VA_ARGS__) #define per(i, a, b) for (ll i = (a) - 1; i >= (b); i--) #define fore(e, v) for (auto&& e : v) #define all(a) begin(a), end(a) #define sz(a) (int)(size(a)) #define lb(v, x) (lower_bound(all(v), x) - begin(v)) #define eb emplace_back template <typename T, typename S> bool chmin(T& a, const S& b) { return a > b ? a = b, 1 : 0; } template <typename T, typename S> bool chmax(T& a, const S& b) { return a < b ? a = b, 1 : 0; } const int INF = 1e9 + 100; const ll INFL = 3e18 + 100; #define i128 __int128_t struct _ { _() { cin.tie(0)->sync_with_stdio(0), cout.tie(0); } } __; void solve(); constexpr int MOD=998244353; int main(){ int T; cin>>T; while(T--)solve(); } ll modpow(ll a,ll n){ long long ret=1%MOD,t=a%MOD; while(n>0){ if(n&1)ret=ret*t%MOD; t=t*t%MOD; n/=2; } return ret; } void ch(ll &a,ll b){ a=(a+b)%MOD; } void solve(){ int N,K; cin>>N>>K; vector equ(N,vl(N+1)); ll m=modpow(K*2+1,MOD-2); rep(i,N){ equ[i][N]=1; ch(equ[i][i],-1); rep(j,-K,K+1){ if(i+j>=0&&i+j<N)ch(equ[i][i+j],m); } } rep(i,N){ ll c=modpow(equ[i][i],MOD-2); rep(k,N+1)equ[i][k]=equ[i][k]*c%MOD; rep(j,N){ if(i==j)continue; ll x=-equ[j][i]; rep(k,N+1){ ch(equ[j][k],equ[i][k]*x); } } } // rep(i,N){ // rep(j,N+1)cerr<<equ[i][j]<<' '; // cerr<<endl; // } rep(i,N){ if(i)cout<<' '; cout<<(MOD-equ[i][N])%MOD; } cout<<'\n'; }