結果
問題 | No.2443 特殊線形群の標準表現 |
ユーザー | kmjp |
提出日時 | 2023-08-25 22:44:44 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 500 ms / 3,000 ms |
コード長 | 1,806 bytes |
コンパイル時間 | 2,165 ms |
コンパイル使用メモリ | 206,484 KB |
実行使用メモリ | 68,864 KB |
最終ジャッジ日時 | 2024-06-06 17:26:34 |
合計ジャッジ時間 | 6,964 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 20 ms
68,736 KB |
testcase_01 | AC | 19 ms
68,864 KB |
testcase_02 | AC | 19 ms
68,708 KB |
testcase_03 | AC | 18 ms
68,808 KB |
testcase_04 | AC | 19 ms
68,736 KB |
testcase_05 | AC | 19 ms
68,788 KB |
testcase_06 | AC | 19 ms
68,736 KB |
testcase_07 | AC | 19 ms
68,708 KB |
testcase_08 | AC | 20 ms
68,836 KB |
testcase_09 | AC | 19 ms
68,736 KB |
testcase_10 | AC | 18 ms
68,736 KB |
testcase_11 | AC | 18 ms
68,736 KB |
testcase_12 | AC | 19 ms
68,716 KB |
testcase_13 | AC | 23 ms
68,732 KB |
testcase_14 | AC | 63 ms
68,712 KB |
testcase_15 | AC | 480 ms
68,760 KB |
testcase_16 | AC | 500 ms
68,716 KB |
testcase_17 | AC | 489 ms
68,796 KB |
testcase_18 | AC | 485 ms
68,736 KB |
testcase_19 | AC | 494 ms
68,728 KB |
testcase_20 | AC | 384 ms
68,736 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef signed long long ll; #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define FORR2(x,y,arr) for(auto& [x,y]:arr) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) template<class T> bool chmax(T &a, const T &b) { if(a<b){a=b;return 1;}return 0;} template<class T> bool chmin(T &a, const T &b) { if(a>b){a=b;return 1;}return 0;} //------------------------------------------------------- ll mo; int N,Q; template<class V,int NV> class SegTree_1 { public: vector<V> val; V comp(V r,V l){ V m; m[0]=(l[0]*r[0]+l[1]*r[2])%mo; m[1]=(l[0]*r[1]+l[1]*r[3])%mo; m[2]=(l[2]*r[0]+l[3]*r[2])%mo; m[3]=(l[2]*r[1]+l[3]*r[3])%mo; return m; }; SegTree_1(){val=vector<V>(NV*2,{1,0,0,1});}; V getval(int x,int y,int l=0,int r=NV,int k=1) { // x<=i<y if(r<=x || y<=l) return {1,0,0,1}; if(x<=l && r<=y) return val[k]; return comp(getval(x,y,l,(l+r)/2,k*2),getval(x,y,(l+r)/2,r,k*2+1)); } void update(int entry, V v) { entry += NV; val[entry]=v; while(entry>1) entry>>=1, val[entry]=comp(val[entry*2],val[entry*2+1]); } }; SegTree_1<array<ll,4>,1<<20> st; void solve() { int i,j,k,l,r,x,y; string s; cin>>N>>mo>>Q; FOR(i,N) { array<ll,4> A; cin>>A[0]>>A[1]>>A[2]>>A[3]; st.update(i,A); } while(Q--) { ll L,R,X,Y; cin>>L>>R>>X>>Y; auto A=st.getval(L,R); ll X2=(A[0]*X+A[1]*Y)%mo; ll Y2=(A[2]*X+A[3]*Y)%mo; cout<<(X2+mo)%mo<<" "<<(Y2+mo)%mo<<endl; } } int main(int argc,char** argv){ string s;int i; if(argc==1) ios::sync_with_stdio(false), cin.tie(0); FOR(i,argc-1) s+=argv[i+1],s+='\n'; FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin); cout.tie(0); solve(); return 0; }