結果
問題 | No.2443 特殊線形群の標準表現 |
ユーザー |
![]() |
提出日時 | 2023-08-25 22:24:52 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 707 ms / 3,000 ms |
コード長 | 1,810 bytes |
コンパイル時間 | 4,718 ms |
コンパイル使用メモリ | 244,008 KB |
実行使用メモリ | 47,872 KB |
最終ジャッジ日時 | 2024-12-24 09:39:58 |
合計ジャッジ時間 | 10,812 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 21 |
ソースコード
#include <bits/stdc++.h> using namespace std; #include <atcoder/all> using namespace atcoder; #define rep(i, n) for(int i=0;i<(n);++i) #define rep1(i, n) for(int i=1;i<=(n);i++) #define ll long long using mint = modint; using P = pair<ll,ll>; using lb = long double; using T = tuple<ll, ll, ll>; #ifdef LOCAL # include <debug_print.hpp> # define dbg(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else # define dbg(...) (static_cast<void>(0)) #endif int main() { int n, b, q; cin >> n >> b >> q; mint::set_mod(b); vector<vector<vector<int>>> a(n, vector<vector<int>>(2, vector<int>(2))); rep(i,n){ rep(j,2)rep(k,2){ cin >> a[i][j][k]; } } vector<vector<vector<mint>>> sum(n+1, vector<vector<mint>>(2, vector<mint>(2))); vector<vector<vector<mint>>> inv(n+1, vector<vector<mint>>(2, vector<mint>(2))); sum[0] = {{1,0},{0,1}}; // dbg(sum[0]); inv[0] = {{1,0},{0,1}}; rep(i,n){ rep(k,2)rep(j,2)rep(l,2){ sum[i+1][k][j] += a[i][k][l] * sum[i][l][j]; } vector<vector<int>> inva = a[i]; swap(inva[0][0], inva[1][1]); inva[0][1] *= -1; inva[1][0] *= -1; rep(k,2)rep(j,2)rep(l,2){ inv[i+1][k][j] += inv[i][k][l] * inva[l][j]; } } rep(qi,q){ int l, r, x, y; cin >> l >> r >> x >> y; // --l; vector<vector<mint>> tmp = sum[r]; vector<vector<mint>> tinv = inv[l]; vector<vector<mint>> now(2, vector<mint>(2)); rep(k,2)rep(j,2)rep(l,2){ now[k][j] += tmp[k][l] * tinv[l][j]; } mint ax, ay = 0; ax = now[0][0]*x + now[0][1]*y; ay = now[1][0]*x + now[1][1]*y; cout << ax.val() << " " << ay.val() << endl; } return 0; }