結果

問題 No.620 ぐるぐるぐるりん
ユーザー maspymaspy
提出日時 2020-05-03 12:23:50
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 969 bytes
コンパイル時間 257 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-06-11 19:07:11
合計ジャッジ時間 3,769 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 1
other AC * 14 WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines

def solve(T, p, w, v, gx, gy):
    T = int(T)
    p, w, v, gx, gy = map(float, (p, w, v, gx, gy))
    r = complex(1 + v, w)
    g = complex(gx, gy) - r**T
    power = [0] * T
    n = 0 if abs(r) > 1 else T - 1
    r = r**(T - 1 - n)
    d = g / r
    for i in range(T):
        if i != n:
            yield (0, 0)
        else:
            yield (d.real, d.imag)

"""
def check(w, v, gx, gy, uvs):
    z = 1
    r = complex(1 + v, w)
    for dx, dy in uvs:
        z *= r
        z += complex(dx, dy)
    return abs(z - complex(gx, gy)) < 1e-6


uvs = ((1, 0), (1, 0), (1, 0), (0, 1), (0, 1))
check(0.5, -0.5, -1.125, 2.125, uvs)"""

N = int(readline())
probs = readlines()

for prob in probs:
    p, w, v, gx, gy = map(float, prob[2:].split())
    uvs = solve(*prob.split())
    print('\n'.join('{} {}'.format(x, y) for x, y in solve(*prob.split())))
0