結果
問題 | No.2443 特殊線形群の標準表現 |
ユーザー |
|
提出日時 | 2023-08-25 22:14:48 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,502 ms / 3,000 ms |
コード長 | 954 bytes |
コンパイル時間 | 5,299 ms |
コンパイル使用メモリ | 263,296 KB |
最終ジャッジ日時 | 2025-02-16 14:05:23 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 21 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using ll = long long; #define rep(i,n) for(int i=0;i<(int)(n);i++) using mint = atcoder::modint998244353; int MOD; using mat=vector<vector<ll>>; mat op(mat a,mat b){ mat rt(2,vector<ll>(2,0)); rep(i,2){ rep(k,2){ rep(j,2){ rt.at(i).at(j)=(rt.at(i).at(j)+a.at(i).at(k)*b.at(k).at(j))%MOD; } } } return rt; } mat e(){ return {{1,0},{0,1}}; } using seg = atcoder::segtree<mat,op,e>; int main(){ int n,b,q; cin>>n>>b>>q; MOD=b; vector<mat> a(n,mat(2,vector<ll>(2))); rep(i,n){ for(auto &v:a.at(i)){ for(auto &x:v){ cin>>x; x=(x%MOD+MOD)%MOD; } } } reverse(a.begin(),a.end()); seg sg(a); rep(lp,q){ int l,r,x,y; cin>>l>>r>>x>>y; l=n-l,r=n-r; x=(x%MOD+MOD)%MOD,y=(y%MOD+MOD)%MOD; mat vv=sg.prod(r,l); cout<<((vv.at(0).at(0)*x)%MOD+(vv.at(0).at(1)*y)%MOD)%MOD<<" "<<((vv.at(1).at(0)*x)%MOD+(vv.at(1).at(1)*y)%MOD)%MOD<<endl; } }