結果
| 問題 | No.2177 Recurring ab | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2023-01-06 22:47:29 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 213 ms / 2,000 ms | 
| コード長 | 790 bytes | 
| コンパイル時間 | 191 ms | 
| コンパイル使用メモリ | 82,176 KB | 
| 実行使用メモリ | 63,420 KB | 
| 最終ジャッジ日時 | 2024-11-30 19:16:38 | 
| 合計ジャッジ時間 | 3,182 ms | 
| ジャッジサーバーID (参考情報) | judge5 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 18 | 
ソースコード
n = int(input())
# x = (p*a+b)/(p**2-1) > 1/n
# (p*a+b)*n > p**2-1
# 0 > p**2 - p*a*n -1 - b*n
# 0 > (p-a*n/2)**2 -1 - b*n - (a*n/2)**2
M = 10**9+1
def f(p,a,b):
    return p**2-p*a*n-1-b*n < 0
ans = 0
for a in range(10):
    for b in range(10):
        if a == b:
            continue
        if n >= 10**4:
            for p in range(max(a,b)+1,100):
                if f(p,a,b):
                    ans += 1
            ll = 100
            l = ll
            r = M
            while l+1 < r:
                m = (r+l)//2
                if f(m,a,b):
                    l = m
                else:
                    r = m
            ans += l-ll+1
        else:
            for p in range(max(a,b)+1,2*10**5):
                if f(p,a,b):
                    ans += 1
print(ans)
            
            
            
        