結果
問題 |
No.2651 [Cherry 6th Tune B] $\mathbb{C}$omplex комбинат
|
ユーザー |
|
提出日時 | 2025-03-23 08:07:51 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 175 ms / 2,500 ms |
コード長 | 1,251 bytes |
コンパイル時間 | 3,735 ms |
コンパイル使用メモリ | 280,180 KB |
実行使用メモリ | 7,936 KB |
最終ジャッジ日時 | 2025-03-23 08:08:04 |
合計ジャッジ時間 | 12,718 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 39 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; #define rep(i, s, e) for (int i = (int)(s); i < (int)(e); ++i) #define all(a) (a).begin(),(a).end() #include <atcoder/modint> using namespace atcoder; using mint = modint998244353; int main() { cin.tie(nullptr); ios_base::sync_with_stdio(false); int T; cin >> T; rep(testcase, 0, T) { int N; cin >> N; vector<ll> X(N), Y(N); rep(i, 0, N) cin >> X[i] >> Y[i]; mint v = 0, inv = 0; rep(i, 0, N) { mint abssq = X[i]*X[i] + Y[i]*Y[i]; v += abssq; inv += 1 / abssq; } mint ans = v * inv; ans -= N; mint t1 = mint(X[0]*X[0] - Y[0]*Y[0]) / mint(X[0]*X[0] + Y[0]*Y[0]); mint t2 = mint(2 * X[0] * Y[0]) / mint(X[0]*X[0] + Y[0]*Y[0]); rep(i, 1, N) { mint abssq = X[i]*X[i] + Y[i]*Y[i]; mint invabssq = 1 / abssq; mint A = mint(X[i]*X[i] - Y[i]*Y[i]) * invabssq; mint B = mint(2 * X[i] * Y[i]) * invabssq; ans -= 2 * (A * t1 + B * t2); t1 += A; t2 += B; } cout << ans.val() << '\n'; } }