結果
問題 | No.2825 Sum of Scores of Sets of Specified Sections |
ユーザー |
![]() |
提出日時 | 2024-07-26 22:05:53 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 97 ms / 3,000 ms |
コード長 | 4,379 bytes |
コンパイル時間 | 2,004 ms |
コンパイル使用メモリ | 205,352 KB |
最終ジャッジ日時 | 2025-02-23 18:40:23 |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 32 |
ソースコード
//#define NDEBUG #pragma warning(disable : 4146) #include <bits/stdc++.h> namespace n91 { using i32 = std::int32_t; using i64 = std::int64_t; using u32 = std::uint32_t; using u64 = std::uint64_t; using isize = std::ptrdiff_t; using usize = std::size_t; using f64 = double; struct rep { struct itr { usize i; constexpr itr(const usize i) noexcept : i(i) {} void operator++() noexcept { ++i; } constexpr usize operator*() const noexcept { return i; } constexpr bool operator!=(const itr x) const noexcept { return i != x.i; } }; const itr f, l; constexpr rep(const usize f, const usize l) noexcept : f(std::min(f, l)), l(l) {} constexpr auto begin() const noexcept { return f; } constexpr auto end() const noexcept { return l; } }; struct revrep { struct itr { usize i; constexpr itr(const usize i) noexcept : i(i) {} void operator++() noexcept { --i; } constexpr usize operator*() const noexcept { return i; } constexpr bool operator!=(const itr x) const noexcept { return i != x.i; } }; const itr f, l; constexpr revrep(const usize f, const usize l) noexcept : f(l - 1), l(std::min(f, l) - 1) {} constexpr auto begin() const noexcept { return f; } constexpr auto end() const noexcept { return l; } }; template <class T> auto md_vec(const usize n, const T& value) { return std::vector<T>(n, value); } template <class... Args> auto md_vec(const usize n, Args... args) { return std::vector<decltype(md_vec(args...))>(n, md_vec(args...)); } template <class T> constexpr T difference(const T& a, const T& b) noexcept { return a < b ? b - a : a - b; } template <class T> void chmin(T& a, const T& b) noexcept { if (b < a) a = b; } template <class T> void chmax(T& a, const T& b) noexcept { if (a < b) a = b; } template <class F> class rec_lambda { F f; public: rec_lambda(F&& f_) : f(std::forward<F>(f_)) {} template <class... Args> auto operator()(Args &&... args) const { return f(*this, std::forward<Args>(args)...); } }; template <class T> T scan() { T ret; std::cin >> ret; return ret; } constexpr char eoln = '\n'; i64 floor_div(const i64 n, const i64 d) { assert(d != 0); return n / d - static_cast<i64>((n ^ d) < 0 && n % d != 0); } i64 ceil_div(const i64 n, const i64 d) { assert(d != 0); return n / d + static_cast<i64>((n ^ d) >= 0 && n % d != 0); } #ifdef N91_LOCAL #define OJ_LOCAL(a, b) b #else #define OJ_LOCAL(a, b) a #endif } // namespace n91 #include <atcoder/modint> using mint = atcoder::modint998244353; #include <cassert> #include <utility> #include <vector> namespace noshi91 { template <class T> T determinant(std::vector<std::vector<T>> a) { if (a.empty()) return 1; const int n = a.size(); for (const auto& r : a) assert(r.size() == n); T ret = 1; for (int i = 0; i < n; i++) { for (int j = i; j < n; j++) { if (a[j][i] != 0) { if (i != j) { std::swap(a[i], a[j]); ret = -ret; } break; } } if (a[i][i] == 0) return 0; ret *= a[i][i]; const T inv = 1 / a[i][i]; for (int j = i; j < n; j++) a[i][j] *= inv; for (int j = i + 1; j < n; j++) { const T c = a[j][i]; for (int k = i; k < n; k++) a[j][k] -= c * a[i][k]; } } return ret; } } // namespace noshi91 namespace n91 { void main_() { const usize h = scan<usize>(); const usize w = scan<usize>(); auto A = md_vec(h, w, mint(0)), B = A; for (auto& r : A) for (auto& e : r) e = scan<u32>(); for (auto& r : B) for (auto& e : r) e = scan<u32>(); auto C = md_vec(h, h, mint(0)); for (const usize i : rep(0, h)) { for (const usize j : rep(0, h)) { for (const usize k : rep(0, w)) { C[i][j] += A[i][k] * B[j][k]; } } } for (const usize i : rep(0, h)) C[i][i] += 1; std::cout << (noshi91::determinant(C) - 1).val() << eoln; } } // namespace n91 int main() { //* std::ios::sync_with_stdio(false); std::cin.tie(nullptr); //*/ std::cout << std::fixed << std::setprecision(20); n91::main_(); return 0; }