結果
問題 |
No.3230 Mutual Corresponding System
|
ユーザー |
![]() |
提出日時 | 2025-08-08 23:40:07 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 121 ms / 2,000 ms |
コード長 | 2,896 bytes |
コンパイル時間 | 2,986 ms |
コンパイル使用メモリ | 279,076 KB |
実行使用メモリ | 7,716 KB |
最終ジャッジ日時 | 2025-08-08 23:40:12 |
合計ジャッジ時間 | 5,415 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include <bits/stdc++.h> //#include <atcoder/all> //using namespace std; //using namespace atcoder; //using mint = modint1000000007; //const int mod = 1000000007; //using mint = modint998244353; //const int mod = 998244353; const int INF = 1e9; //const long long LINF = 1e18; #define rep(i, n) for (int i = 0; i < (n); ++i) #define rep2(i,l,r)for(int i=(l);i<(r);++i) #define rrep(i, n) for (int i = (n) - 1; i >= 0; --i) #define rrep2(i,l,r)for(int i=(r) - 1;i>=(l);--i) #define all(x) (x).begin(),(x).end() #define allR(x) (x).rbegin(),(x).rend() #define P pair<int,int> template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; } template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; } template<class T, T(*add)(T, T), T(*zero)(), T(*mul)(T, T), T(*one)(), int N> struct Matrix : std::array<std::array<T, N>, N> { using M = Matrix<T, add, zero, mul, one, N>; void make_identity() { for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { (*this)[i][j] = zero(); } } for (int i = 0; i < N; i++) { (*this)[i][i] = one(); } } M& operator+=(const M& rhs) { for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { (*this)[i][j] = add((*this)[i][j], rhs[i][j]); } } return *this; } M& operator*=(const M& rhs) { static M temp; for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { temp[i][j] = zero(); } } for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { for (int k = 0; k < N; k++) { temp[i][k] = add(temp[i][k], mul((*this)[i][j], rhs[j][k])); } } } *this = temp; return *this; } template <class I> void inplace_pow(I k) { assert(k >= 0); static M temp; temp = *this; make_identity(); while (k) { if (k & 1) *this *= temp; k >>= 1; if (k) temp *= temp; } } friend std::ostream& operator<<(std::ostream& os, const M& A) { for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { os << A[i][j] << " \n"[j + 1 == N]; } } return os; } }; // tropical semiring long long add(long long x, long long y) { return std::max(x, y); } long long zero() { return -1002003004005006007; } long long mul(long long x, long long y) { return x + y; } long long one() { return 0; } using M = Matrix<long long, add, zero, mul, one, 100>; int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); int n; std::cin >> n; long long m; std::cin >> m; auto mat = Matrix<long long, add, zero, mul, one, 100>(); rep(i, 100) rep(j, 100) mat[i][j] = zero(); std::vector<int>t(n); rep(i, n) std::cin >> t[i]; rep(i, n) rep(j, n) std::cin >> mat[i][j]; mat.inplace_pow(m); rep(j, n) { long long ans = zero(); rep(i, n) ans = add(ans, mul(t[i], mat[i][j])); std::cout << ans << " \n"[j + 1 == n]; } return 0; }