結果

問題 No.2362 Inversion Number of Mod of Linear
ユーザー 遭難者遭難者
提出日時 2023-05-07 12:23:51
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 511 ms / 2,000 ms
コード長 1,397 bytes
コンパイル時間 265 ms
コンパイル使用メモリ 86,792 KB
実行使用メモリ 79,724 KB
最終ジャッジ日時 2023-09-05 03:29:15
合計ジャッジ時間 3,387 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 73 ms
71,264 KB
testcase_01 AC 72 ms
71,476 KB
testcase_02 AC 73 ms
71,260 KB
testcase_03 AC 373 ms
79,232 KB
testcase_04 AC 180 ms
79,384 KB
testcase_05 AC 187 ms
79,544 KB
testcase_06 AC 511 ms
79,724 KB
testcase_07 AC 325 ms
78,920 KB
testcase_08 AC 246 ms
79,104 KB
testcase_09 AC 242 ms
78,572 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
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()
0