結果
問題 | No.2651 [Cherry 6th Tune B] omplex комбинат |
ユーザー |
|
提出日時 | 2024-02-23 23:24:33 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 653 ms / 2,500 ms |
コード長 | 1,375 bytes |
コンパイル時間 | 369 ms |
コンパイル使用メモリ | 82,560 KB |
実行使用メモリ | 114,504 KB |
最終ジャッジ日時 | 2024-09-29 08:57:40 |
合計ジャッジ時間 | 21,375 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 39 |
ソースコード
#int(input())#map(int, input().split())#list(map(int, input().split()))T = int(input())n = [0] * TX = [0] * Tfor i in range(T):n[i] = int(input())y = [0] * n[i]for j in range(n[i]):y[j] = list(map(int, input().split()))X[i] = ymod = 998244353b = []for k in range(T):# k = 1N = n[k]x = X[k]norm = [0] * Nfor i in range(N):norm[i] = x[i][0] ** 2 + x[i][1] ** 2norminv = [0] * Nfor i in range(N):norminv[i] = pow(norm[i], mod-2, mod)def f(x, y):re = x[0] * y[0] - x[1] * y[1]im = x[1] * y[0] + x[0] * y[1]return [re%mod, im%mod]s0 = [0, 0]for i in range(N):a = f(x[i], x[i])# print(x[i], a)for j in range(2):s0[j] = (s0[j] + a[j] * norminv[i]) % mods1 = [0, 0]for i in range(N):y = [x[i][0], -x[i][1]]a = f(y, y)# print(y, a)for j in range(2):s1[j] = (s1[j] + a[j] * norminv[i]) % modu0 = 0u1 = 0for i in range(N):u0 = (u0 + norm[i]) % modu1 = (u1 + norminv[i]) % modans = 0ans = (ans + u0*u1) % modans = (ans - (s0[0]*s1[0] - s0[1]*s1[1])) % mod# ans = (ans * pow(2, mod-2, mod)) % mod# print(u0, u1)# print(s0, s1)b.append(ans)for x in b:print(x)