結果

問題 No.2443 特殊線形群の標準表現
コンテスト
ユーザー ゼット
提出日時 2023-08-25 22:06:06
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 923 ms / 3,000 ms
コード長 1,096 bytes
コンパイル時間 336 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 192,988 KB
最終ジャッジ日時 2024-12-24 08:53:37
合計ジャッジ時間 8,414 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

N,B,Q=map(int,input().split())
A1=[]
A2=[]
T=[]
T2=[]
for i in range(N):
  v=[]
  v2=[]
  for j in range(2):
    a,b=map(int,input().split())
    a%=B
    b%=B
    v.append([a,b])
    v2.append([0,0])
  v2[1][1]=v[0][0]
  v2[0][0]=v[1][1]
  v2[0][1]=-v[0][1]
  v2[1][0]=-v[1][0]
  A1.append(v)
  A2.append(v2)
  if i==0:
    T.append(v.copy())
    T2.append(v2.copy())
for w in range(1,N):
  r=[[0]*2 for i in range(2)]
  r2=[[0]*2 for i in range(2)]
  for i in range(2):
    for j in range(2):
      for k in range(2):
        r[i][j]+=A1[w][i][k]*T[w-1][k][j]
        r[i][j]%=B
        r2[i][j]+=T2[w-1][i][k]*A2[w][k][j]
        r2[i][j]%=B
  T.append(r.copy())
  T2.append(r2.copy())
for _ in range(Q):
  l,r,x,y=map(int,input().split())
  x%=B
  y%=B
  if r==0:
    print(x,y)
    continue
  p=[[0]*2 for i in range(2)]
  for i in range(2):
    for j in range(2):
      if l==0:
        p[i][j]=T[r-1][i][j]
      for k in range(2):
        if l>0:
          p[i][j]+=T[r-1][i][k]*T2[l-1][k][j]
          p[i][j]%=B
  z=p[0][0]*x+p[0][1]*y
  w=p[1][0]*x+p[1][1]*y
  z%=B
  w%=B
  print(z,w)
0