結果
問題 |
No.1750 ラムドスウイルスの感染拡大-hard
|
ユーザー |
![]() |
提出日時 | 2021-11-19 22:04:02 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 236 ms / 2,000 ms |
コード長 | 2,451 bytes |
コンパイル時間 | 2,601 ms |
コンパイル使用メモリ | 209,764 KB |
最終ジャッジ日時 | 2025-01-25 20:22:28 |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> using namespace std; struct iofast_t { iofast_t() { ios::sync_with_stdio(false); cin.tie(nullptr); } } iofast; struct uns_t {} uns; template <typename Element, typename Head, typename ...Args> auto vec(Element init, Head arg, Args ...args) { if constexpr (sizeof...(Args) == 0) return std::vector(arg, init); else return std::vector(arg, vec(init, args...)); } template <typename Element, typename Head, typename ...Args> auto vec(uns_t, Head arg, Args ...args) { return vec(Element(), arg, args...); } template <typename T, typename Compare = less<T>> T &chmin(T &l, T r, Compare &&f = less<T>()) { return l = min(l, r, f); } template <typename T, typename Compare = less<T>> T &chmax(T &l, T r, Compare &&f = less<T>()) { return l = max(l, r, f); } #include <atcoder/modint> using mint = atcoder::modint998244353; using matrix = vector<vector<mint>>; auto melem(const matrix &mat) { matrix m = mat; for (int i = 0; i < size(m); ++i) { for (int j = 0; j < size(m[0]); ++j) { m[i][j] = (i == j); } } return m; } auto mmul(matrix l, matrix r) { int h = size(l); int w = size(r[0]); int s = size(l[0]); matrix m; m.resize(h); for (auto &r : m) { r.resize(w); } for (int i = 0; i < h; ++i) { for (int j = 0; j < w; ++j) { m[i][j] = 0; for (int k = 0; k < s; ++k) { m[i][j] += l[i][k] * r[k][j]; } } } return m; } auto mpow(matrix l, int64_t x) { matrix r = melem(l); while (0 < x) { if (x & 1) { r = mmul(r, l); } l = mmul(l, l); x >>= 1; } return r; } int main() { int64_t n, m, t; cin >> n >> m >> t; auto g = vec<int>(uns, n, 0); for (int i = 0; i < m; ++i) { int s, t; cin >> s >> t; g[s].push_back(t); g[t].push_back(s); } matrix dp; dp.resize(n); for (auto &r : dp) { r.resize(n); } for (auto &u : dp) { for (auto &v : u) { v = 0; } } for (int i = 0; i < n; ++i) { for (auto v : g[i]) { dp[i][v] = 1; } } matrix mat; mat.resize(n); for (auto &v : mat) { v.resize(1); v[0] = 0; } mat[0][0] = 1; auto ans = mmul(mpow(dp, t), mat); cout << ans[0][0].val() << endl; }