結果
| 問題 |
No.2651 [Cherry 6th Tune B] $\mathbb{C}$omplex комбинат
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-02-28 10:29:06 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 858 ms / 2,500 ms |
| コード長 | 1,643 bytes |
| コンパイル時間 | 5,515 ms |
| コンパイル使用メモリ | 315,768 KB |
| 実行使用メモリ | 6,820 KB |
| 最終ジャッジ日時 | 2024-09-29 12:19:04 |
| 合計ジャッジ時間 | 31,234 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 39 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
#define rep(i, n) for (int i = 0, i##_len = (n); i < i##_len; ++i)
#define all(v) begin(v), end(v)
using namespace std;
using namespace atcoder;
int main() {
ios_base::sync_with_stdio(false);
cin.tie(NULL);
using mint = modint998244353;
using P = pair<mint, mint>;
int Q;
cin >> Q;
while (Q--) {
int n;
cin >> n;
vector<P> z(n);
rep(i, n) {
int l, r;
cin >> l >> r;
z[i] = {l, r};
}
mint A = 0, B = 0, C = 0, D = 0, E = 0, F = 0, G = 0;
rep(i, n) {
auto [x, y] = z[i];
mint xxyy = x * x + y * y;
A += (2 * x * x) / xxyy;
B += (8 * x * y) / xxyy;
C += (2 * y * y) / xxyy;
D += x * x;
E += y * y;
F += x * x / xxyy / xxyy;
G += y * y / xxyy / xxyy;
}
mint ans = 0;
rep(i, n) {
auto [a, b] = z[i];
mint aabb = a * a + b * b;
ans += -A * (a * a) / aabb;
ans += A * (b * b) / aabb;
ans += -B * a * b / aabb;
ans += C * (a * a) / aabb;
ans += -C * (b * b) / aabb;
ans += D * (a * a) / aabb / aabb;
ans += D * (b * b) / aabb / aabb;
ans += E * (a * a) / aabb / aabb;
ans += E * (b * b) / aabb / aabb;
ans += F * (a * a);
ans += F * (b * b);
ans += G * (a * a);
ans += G * (b * b);
}
ans /= 2;
cout << ans.val() << endl;
}
return 0;
}