結果

問題 No.1932 動く点 P / Moving Point P
ユーザー Navier_BoltzmannNavier_Boltzmann
提出日時 2022-05-12 13:01:15
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 3,059 bytes
コンパイル時間 501 ms
コンパイル使用メモリ 82,132 KB
実行使用メモリ 187,012 KB
最終ジャッジ日時 2024-07-20 06:18:49
合計ジャッジ時間 32,048 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
53,120 KB
testcase_01 WA -
testcase_02 WA -
testcase_03 AC 3,897 ms
97,112 KB
testcase_04 AC 4,893 ms
112,548 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 TLE -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import *
def segfunc(x, y):
    return matrix_mul(y,x)

ide_ele = [[1,0,0],[0,1,0],[0,0,1]]


class SegTree:
    """
    Segment Tree
    """

    def __init__(self, init_val, segfunc, ide_ele):
        """
        初期化

        init_val: 配列の初期値
        """
        n = len(init_val)
        self.segfunc = segfunc
        self.ide_ele = ide_ele
        self.num = 1 << (n - 1).bit_length()
        self.tree = [ide_ele] * 2 * self.num
        # 配列の値を葉にセット
        for i in range(n):
            self.tree[self.num + i] = init_val[i]
        # 構築していく
        for i in range(self.num - 1, 0, -1):
            self.tree[i] = segfunc(self.tree[2 * i], self.tree[2 * i + 1])

    def update(self, k, x):
        """
        k番目の値をxに更新

        k: index(0-index)
        x: update value
        """
        k += self.num
        self.tree[k] = x
        while k > 1:
            self.tree[k >> 1] = self.segfunc(self.tree[k], self.tree[k ^ 1])
            k >>= 1

    def query(self, l, r):
        """
        [l, r)のsegfuncしたものを得る

        l: index(0-index)
        r: index(0-index)
        """
        res = self.ide_ele

        l += self.num
        r += self.num
        while l < r:
            if l & 1:
                res = self.segfunc(res, self.tree[l])
                l += 1
            if r & 1:
                res = self.segfunc(res, self.tree[r - 1])
            l >>= 1
            r >>= 1
        return res


def matrix_mul(A,B,mod = None):
    nA = len(A)
    mA = len(A[0])
    mB = len(B[0])
    
    tmp = [[0]*mB for _ in range(nA)]
    
    if mod is None:
        for i in range(nA):
            for j in range(mB):
                tmp[i][j] = sum(A[i][k]*B[k][j] for k in range(mA))
        return tmp
    
    for i in range(nA):
        for j in range(mB):
            tmp[i][j] = sum(A[i][k]*B[k][j]%mod for k in range(mA))%mod
    return tmp
    

def matrix_pow(A,n,mod = None):
    nbit = list(str(bin(n))[2:])
    nbit = [int(i) for i in nbit]
    N = len(A)
    C = [[0]*N for _ in range(N)]
    B = A
    for i in range(N):
        C[i][i] = 1
        
    if mod is None:
        
        for i in range(len(nbit)):
            if nbit[-1-i] == 1:
                C = matrix_mul(C,B)
            
            B = matrix_mul(B,B)
        return C
        
    for i in range(len(nbit)):
        if nbit[-1-i] == 1:
            C = matrix_mul(C,B,mod)
        
        B = matrix_mul(B,B,mod)
    
    return C
    
def R(x,y,r):
    
    c = cos(pi*r/180)
    s = sin(pi*r/180)
    
    M = [[c,-s,-x*c+y*s + x],
         [s,c,-x*s-y*c + y],
         [0,0,1]]
    return M

N = int(input())
X = []
for _ in range(N):
    p,q,r = map(float,input().split())
    X.append(R(p,q,r))
    
T = SegTree(X,segfunc,ide_ele)

Q = int(input())
for _ in range(Q):
    s,t,x,y = map(float,input().split())
    s = int(s)
    t = int(t)
    s -= 1
    Y = [[x],[y],[1]]
    tmp = matrix_mul(T.query(s,t),Y)
    print(tmp[0][0],tmp[1][0])
0