結果
問題 | No.2651 [Cherry 6th Tune B] $\mathbb{C}$omplex комбинат |
ユーザー | tassei903 |
提出日時 | 2024-02-23 22:10:11 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,338 ms / 2,500 ms |
コード長 | 1,346 bytes |
コンパイル時間 | 742 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 127,772 KB |
最終ジャッジ日時 | 2024-04-29 23:30:14 |
合計ジャッジ時間 | 42,281 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 45 ms
52,352 KB |
testcase_01 | AC | 48 ms
52,352 KB |
testcase_02 | AC | 1,145 ms
78,628 KB |
testcase_03 | AC | 1,226 ms
89,284 KB |
testcase_04 | AC | 1,303 ms
93,176 KB |
testcase_05 | AC | 1,295 ms
119,476 KB |
testcase_06 | AC | 1,338 ms
121,768 KB |
testcase_07 | AC | 1,287 ms
120,556 KB |
testcase_08 | AC | 1,273 ms
120,832 KB |
testcase_09 | AC | 1,276 ms
123,656 KB |
testcase_10 | AC | 1,294 ms
127,772 KB |
testcase_11 | AC | 1,272 ms
124,288 KB |
testcase_12 | AC | 1,286 ms
123,692 KB |
testcase_13 | AC | 1,276 ms
126,204 KB |
testcase_14 | AC | 1,278 ms
120,992 KB |
testcase_15 | AC | 1,276 ms
124,212 KB |
testcase_16 | AC | 1,278 ms
125,112 KB |
testcase_17 | AC | 1,264 ms
124,432 KB |
testcase_18 | AC | 1,264 ms
124,408 KB |
testcase_19 | AC | 1,295 ms
126,884 KB |
testcase_20 | AC | 1,272 ms
79,520 KB |
testcase_21 | AC | 427 ms
87,808 KB |
testcase_22 | AC | 687 ms
99,840 KB |
testcase_23 | AC | 907 ms
110,500 KB |
testcase_24 | AC | 1,074 ms
117,924 KB |
testcase_25 | AC | 164 ms
79,104 KB |
testcase_26 | AC | 626 ms
96,032 KB |
testcase_27 | AC | 1,174 ms
120,448 KB |
testcase_28 | AC | 888 ms
108,160 KB |
testcase_29 | AC | 437 ms
92,800 KB |
testcase_30 | AC | 176 ms
78,976 KB |
testcase_31 | AC | 998 ms
118,072 KB |
testcase_32 | AC | 241 ms
82,048 KB |
testcase_33 | AC | 913 ms
114,496 KB |
testcase_34 | AC | 591 ms
97,920 KB |
testcase_35 | AC | 1,095 ms
121,928 KB |
testcase_36 | AC | 1,010 ms
78,764 KB |
testcase_37 | AC | 1,181 ms
127,104 KB |
testcase_38 | AC | 1,255 ms
127,232 KB |
ソースコード
import sys input = lambda :sys.stdin.readline()[:-1] ni = lambda :int(input()) na = lambda :list(map(int,input().split())) yes = lambda :print("yes");Yes = lambda :print("Yes") no = lambda :print("no");No = lambda :print("No") ####################################################################### mod = 998244353 def mul(a, b): return ((a[0] * b[0] - a[1] * b[1]) % mod, (a[0] * b[1] + a[1] * b[0]) % mod) def add(a, b): return ((a[0] + b[0]) % mod, (a[1] + b[1]) % mod) def norm(a): return (a[0] * a[0] + a[1] * a[1]) % mod def div(a, b): return mul(mul(a, (b[0], -b[1])) , (pow(norm(b), mod - 2, mod), 0)) for i in range(ni()): n = ni() p = [na() for _ in range(n)] normP = [norm(x) for x in p] B = 0 for i in range(n): B = (B + normP[i]) % mod C = 0 for i in range(n): C = (C + pow(normP[i], mod - 2, mod)) % mod D = (2 * B * C) % mod D -= n * 2 D *= pow(2, mod - 2, mod) for i in range(n): p[i] = mul(p[i], p[i]) p[i] = div(p[i], (normP[i], 0)) sum = (0, 0) inv_p = [div((1, 0), x) for x in p] sumP = (0, 0) for i in range(n): sumP = add(sumP, p[i]) sumP_inv = (0, 0) for i in range(n): sumP_inv = add(sumP_inv, inv_p[i]) ans = mul(sumP, sumP_inv) ans = add(ans, (-n, 0)) print((D-ans[0])%mod)