結果
| 問題 | 
                            No.1918 Simple Math ?
                             | 
                    
| コンテスト | |
| ユーザー | 
                             とりゐ
                         | 
                    
| 提出日時 | 2022-03-26 10:45:34 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 279 ms / 2,000 ms | 
| コード長 | 482 bytes | 
| コンパイル時間 | 210 ms | 
| コンパイル使用メモリ | 82,176 KB | 
| 実行使用メモリ | 136,832 KB | 
| 最終ジャッジ日時 | 2024-06-29 02:08:01 | 
| 合計ジャッジ時間 | 5,429 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge1 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 1 | 
| other | AC * 33 | 
ソースコード
def floor_sqrt(n):
  ng,ok=n+1,-1
  while ng-ok>1:
    mid=(ok+ng)//2
    if mid*mid<=n:
      ok=mid
    else:
      ng=mid
  return ok
def solve():
  a,n=map(int,input().split())
  m=floor_sqrt(a*n)
  c=[0]*(a+1)
  for i in range(1,a+1):
    c[i]=c[i-1]+(i*i-1)%a
    c[i]%=mod
  
  num=m*(2*m*m+3*m-5)%mod*inv6%mod
  p,q=m//a,m%a
  num-=c[a]*p+c[q]
  num=num*pow(a,mod-2,mod)%mod
  print((n*m-num)%mod)
mod=10**9+7
inv6=pow(6,mod-2,mod)
for _ in range(int(input())):
  solve()
            
            
            
        
            
とりゐ