結果
問題 |
No.620 ぐるぐるぐるりん
|
ユーザー |
![]() |
提出日時 | 2020-05-03 12:35:36 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
WA
|
実行時間 | - |
コード長 | 994 bytes |
コンパイル時間 | 83 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 11,008 KB |
最終ジャッジ日時 | 2024-06-11 19:24:14 |
合計ジャッジ時間 | 4,755 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | WA * 2 RE * 1 |
other | WA * 20 RE * 8 |
ソースコード
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) power = [r**i for i in range(T + 1)] g = complex(gx, gy) - power[T] s = sum(power[i] for i in range(T)) cx = g.real / s.real cy = g.imag / s.imag for i in range(T): u = cx * power[T - 1 - i].real v = cy * power[T - 1 - i].imag yield (u, v) """def check(p, w, v, gx, gy, uvs): z = 1 r = complex(1 + v, w) for dx, dy in uvs: z *= r z += complex(dx, dy) p -= dx**2 + dy**2 return abs(z - complex(gx, gy)) < 1e-6 and p > 0 uvs = ((1, 0), (1, 0), (1, 0), (0, 1), (0, 1)) check(5.5, 0.5, -0.5, -1.125, 2.125, uvs)""" N = int(readline()) probs = readlines() for prob in probs: uvs = solve(*prob.split()) print('\n'.join('{} {}'.format(x, y) for x, y in uvs))