結果
問題 | No.2362 Inversion Number of Mod of Linear |
ユーザー | 遭難者 |
提出日時 | 2023-06-08 13:45:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 470 ms / 2,000 ms |
コード長 | 1,429 bytes |
コンパイル時間 | 281 ms |
コンパイル使用メモリ | 87,200 KB |
実行使用メモリ | 80,620 KB |
最終ジャッジ日時 | 2023-09-05 03:28:58 |
合計ジャッジ時間 | 3,658 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 74 ms
71,372 KB |
testcase_01 | AC | 73 ms
71,200 KB |
testcase_02 | AC | 70 ms
71,220 KB |
testcase_03 | AC | 327 ms
80,340 KB |
testcase_04 | AC | 140 ms
78,184 KB |
testcase_05 | AC | 133 ms
78,348 KB |
testcase_06 | AC | 470 ms
80,620 KB |
testcase_07 | AC | 272 ms
79,292 KB |
testcase_08 | AC | 201 ms
78,836 KB |
testcase_09 | AC | 190 ms
78,288 KB |
ソースコード
import sys input = sys.stdin.buffer.readline from math import gcd N = [0 for i in range(50)] M = [0 for i in range(50)] A = [0 for i in range(50)] B = [0 for i in range(50)] f = [0 for i in range(50)] g = [0 for i in range(50)] h = [0 for i in range(50)] def q(n, m, a, b, c=0) -> None: if m == 0: f[c] = g[c] = h[c] = 0 return a1, a2 = a // m, a % m b1, b2 = b // m, b % m N[c], M[c], A[c], B[c] = n, m, a2, b2 y = (a2 * n + b2) // m q(y, a2, m, m + a2 - b2 - 1, c + 1) nn = n * (n - 1) // 2 f[c] = n * y - f[c + 1] g[c] = n * y * y - f[c + 1] - 2 * h[c + 1] h[c] = nn * y + (f[c + 1] - g[c + 1]) // 2 g[c] += n * (2 * n - 1) * (n - 1) * a1 * a1 // 6 g[c] += 2 * nn * a1 * b1 g[c] += b1 * (b1) * n g[c] += 2 * a1 * h[c] g[c] += 2 * b1 * f[c] f[c] += a1 * nn + b1 * n h[c] += nn * (a1 * (2 * n - 1) + 3 * b1) // 3 return def main() -> None: n, m, a, b = map(int, input().split()) gccd = gcd(m, a) m, a, b = m // gccd, a // gccd, b // gccd ans = 0 q(n, m, a, b) ans -= (n - 1) * f[0] ans += 2 * h[0] t = a * (n - 1) // m y = m - a * (n - 1) + t * m q(n, m, a, y) ans += f[0] ans += h[0] ans -= n * (n + 1) // 2 * t k1, k2 = n // m, n % m ans -= (k1 + 1) * (k1 + 2) // 2 * k2 ans -= k1 * (k1 + 1) // 2 * (m - k2) print(ans) return for i in range(int(input())): main()