結果

問題 No.2177 Recurring ab
ユーザー aaaaaaaaaa2230aaaaaaaaaa2230
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0