結果
問題 | No.2177 Recurring ab |
ユーザー | aaaaaaaaaa2230 |
提出日時 | 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)