結果
| 問題 | 
                            No.2501 Maximum Inversion Number
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2023-07-24 20:33:05 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 1,524 ms / 2,000 ms | 
| コード長 | 1,186 bytes | 
| コンパイル時間 | 248 ms | 
| コンパイル使用メモリ | 82,176 KB | 
| 実行使用メモリ | 137,208 KB | 
| 最終ジャッジ日時 | 2024-09-15 13:56:18 | 
| 合計ジャッジ時間 | 7,960 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 17 | 
ソースコード
import sys
from typing import List
def main(_: int, m: int, l: List[int], r: List[int]):
    if not (sum(l) <= m <= sum(r)):
        return -1
    
    if m == sum(l):
        return (m * m - sum(e * e for e in l)) // 2
    
    def clamp(v: int, l: int, r: int):
        return max(l, min(v, r))
    
    def get_sum(t: int):
        return sum(clamp(t, l[i], r[i]) for i in range(n))
    
    tl, tr = 0, 10 ** 9 + 10
    while abs(tr - tl) > 1:
        t = (tl + tr) // 2
        if get_sum(t) < m:
            tl = t
        else:
            tr = t
    t0, d = tl, m - get_sum(tl)
    x = 0
    for i in range(n):
        a, b = l[i], r[i]
        c = clamp(t0, a, b)
        if d > 0 and a <= t0 < b:
            c += 1
            d -= 1
        x += c * c
    
    return (m * m - x) // 2
if __name__ == '__main__':
    t = int(sys.stdin.readline().rstrip())
    answers = []
    for _ in range(t):
        n, m = map(int, sys.stdin.readline().rstrip().split())
        l = list(map(int, sys.stdin.readline().rstrip().split()))
        r = list(map(int, sys.stdin.readline().rstrip().split()))
        answers.append(main(n, m, l, r))
    print('\n'.join(map(str, answers)))