結果
問題 | No.752 mod数列 |
ユーザー |
|
提出日時 | 2024-07-05 05:42:59 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 380 ms / 2,000 ms |
コード長 | 496 bytes |
コンパイル時間 | 141 ms |
コンパイル使用メモリ | 82,288 KB |
実行使用メモリ | 79,416 KB |
最終ジャッジ日時 | 2024-07-05 05:43:09 |
合計ジャッジ時間 | 9,269 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
from bisect import bisect_right def f(r): i = bisect_right(R,r)-1 res = r*p - acc[i] r = min(r,p) res -= (r+R[i]+1)*(r-R[i])//2*V[i] return res p,Q = map(int,input().split()) M = int(p**0.5) V = [p//i for i in range(1,M)] + list(range(1,p//M+1)[::-1]) R = [0]+[p//i for i in V] V += [0] acc = [0]*(len(V)+1) for i in range(1,len(V)): acc[i] = acc[i-1] + (R[i]+R[i-1]+1)*(R[i]-R[i-1])//2*V[i-1] for _ in range(Q): l,r = map(int,input().split()) print(f(r)-f(l-1))