結果

問題 No.2443 特殊線形群の標準表現
コンテスト
ユーザー hiromi_ayase
提出日時 2023-08-25 23:06:52
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,399 ms / 3,000 ms
コード長 1,359 bytes
コンパイル時間 5,823 ms
コンパイル使用メモリ 325,440 KB
実行使用メモリ 57,088 KB
最終ジャッジ日時 2024-12-24 10:49:49
合計ジャッジ時間 15,420 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#include <atcoder/all>
using namespace std;
using i32 = int;
using u32 = unsigned int;
using i64 = long long;
using u64 = unsigned long long;
#define FAST_IO                \
  ios::sync_with_stdio(false); \
  cin.tie(0);
const i64 INF = 1001001001001001001;
using Modint = atcoder::static_modint<998244353>;

using mint = atcoder::modint;

int main() {
  FAST_IO

  int n,b,q;
  cin >> n >> b >> q;

  mint::set_mod(b);

  using M = vector<vector<mint>>;
  auto op = [](M x, M y) {
    M z(2, vector<mint>(2, 0));
    for (auto i : views::iota(0, 2)) {
      for (auto j : views::iota(0, 2)) {
        for (auto k : views::iota(0, 2)) {
          z[i][j] += y[i][k] * x[k][j];
        }
      }
    }
    return z;
  };
  auto e = []() {
    M I(2, vector<mint>(2, 0));
    I[0][0] = I[1][1] = 1;
    return I;
  };

  atcoder::segtree<M, op, e> st(n);

  for (auto i : views::iota(0, n)) {
    int a1,a2,a3,a4;
    cin >> a1 >> a2 >> a3 >> a4;
    M a(2, vector<mint>(2, 0));
    
    a[0][0] = a1;
    a[0][1] = a2;
    a[1][0] = a3;
    a[1][1] = a4;
    st.set(i, a);
  }

  for (auto i : views::iota(0,q)) {
    int l,r,x,y;
    cin >> l >> r >> x >> y;
    M v(2, vector<mint>(2, 0));
    v[0][0] = x;
    v[1][0] = y;
    auto ret = op(v, st.prod(l, r));

    cout << ret[0][0].val() << " " << ret[1][0].val() << endl;
  }
}
0