結果
問題 | 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;}}