結果
問題 | No.2651 [Cherry 6th Tune B] omplex комбинат |
ユーザー |
![]() |
提出日時 | 2024-02-23 22:52:47 |
言語 | Crystal (1.14.0) |
結果 |
AC
|
実行時間 | 283 ms / 2,500 ms |
コード長 | 836 bytes |
コンパイル時間 | 13,631 ms |
コンパイル使用メモリ | 297,160 KB |
実行使用メモリ | 13,036 KB |
最終ジャッジ日時 | 2024-09-29 08:07:32 |
合計ジャッジ時間 | 21,838 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 39 |
ソースコード
MOD = 998244353i64read_line.to_i.times doputs solve()enddef solven = read_line.to_ixys = Array.new(n) { Tuple(Int64, Int64).from(read_line.split.map(&.to_i64)) }len = xys.map { |x, y| x ** 2 + y ** 2 }sum_len = len.sum % MODans = n.times.sum { |i| inv(len[i]) } % MODans *= sum_lenans %= MODsum_re = 0i64sum_im = 0i64n.times do |i|div = inv(len[i])sum_re += (xys[i][0] ** 2 - xys[i][1] ** 2) * divsum_re %= MODsum_im += xys[i][0] * xys[i][1] * 2 % MOD * divsum_im %= MODendans -= sum_re ** 2 + sum_im ** 2ans %= MODans += MOD if ans < 0return ansenddef inv(v : Int64)pow(v, MOD - 2)enddef pow(v : Int64, p)ret = 1i64while p > 0if (p & 1i64) != 0ret *= vret %= MODendv *= vv %= MODp >>= 1endretend