結果

問題 No.1932 動く点 P / Moving Point P
ユーザー lilictakalilictaka
提出日時 2022-05-21 13:07:16
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,991 bytes
コンパイル時間 443 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 283,276 KB
最終ジャッジ日時 2024-09-20 11:36:03
合計ジャッジ時間 11,441 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 144 ms
82,176 KB
testcase_01 TLE -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math
from decimal import Decimal
class SegTree:
    """
    init(init_val, ide_ele): 配列init_valで初期化 O(N)
    update(k, x): k番目の値をxに更新 O(logN)
    query(l, r): 区間[l, r)をsegfuncしたものを返す O(logN)
    """
    def __init__(self, init_val, segfunc, ide_ele):
        """
        init_val: 配列の初期値
        segfunc: 区間にしたい操作
        ide_ele: 単位元
        n: 要素数
        num: n以上の最小の2のべき乗
        tree: セグメント木(1-index)
        """
        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] = self.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 value(self,index):
        return self.query(index,index+1)


def conv(p,q,r):
    rad = math.radians(r)
    COS = Decimal(math.cos(rad))
    SIN = Decimal(math.sin(rad))
    rev = [[COS,-SIN, -p * COS + q * SIN + p],
    [SIN,COS,- p * SIN -q * COS + q],
    [0,0,1]
    ]
    return rev
def multply(A,B):
    C = [[0] * len(B[0]) for _ in range(len(A))]
    for i in range(len(A)):
        for k in range(len(B[0])):
            for m in range(len(B)):
                C[i][k] += A[i][m] * B[m][k]
    return C
def remultply(A,B):
    C = [[0] * len(A[0]) for _ in range(len(B))]
    for i in range(len(B)):
        for k in range(len(A[0])):
            for m in range(len(A)):
                C[i][k] += B[i][m] * A[m][k]
    return C
board = []
N = int(input())
for _ in range(N):
    p,q,r = map(Decimal,input().split())
    mat = conv(p,q,r)
    board.append(mat)
e = [[1,0,0],[0,1,0],[0,0,1]]
seg = SegTree(board,remultply,e)
Q = int(input())
ANS = []
for _ in range(Q):
    s,t,x,y = map(Decimal,input().split())
    s = int(s)
    t = int(t)
    s-=1
    t-=1
    ans = multply(seg.query(s,t+1),[[x],[y],[1]])
    ANS.append((ans[0][0],ans[1][0]))
for x,y in ANS:
    print(float(x),float(y))
0