結果
| 問題 | 
                            No.1383 Numbers of Product
                             | 
                    
| コンテスト | |
| ユーザー | 
                             convexineq
                         | 
                    
| 提出日時 | 2021-02-09 18:26:13 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                            (最新)
                                AC
                                 
                             
                            (最初)
                            
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 1,015 bytes | 
| コンパイル時間 | 284 ms | 
| コンパイル使用メモリ | 82,304 KB | 
| 実行使用メモリ | 243,852 KB | 
| 最終ジャッジ日時 | 2024-09-14 02:42:37 | 
| 合計ジャッジ時間 | 22,821 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge3 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 | 
| other | AC * 50 WA * 1 | 
ソースコード
from math import sqrt
n,k,m = map(int,input().split())
fk = float(k)
def val(a,s):
    v = float(a)
    for i in range(1,s):
        v *= a+i*fk
        if v > 2e18: return 2000000000000000000
    w = a
    for i in range(1,s):
        w *= a+i*k
    return w
from collections import Counter
d = Counter()
for s in range(3,100):
    if val(1,s) > n: break
    for a in range(1,10**6+1):
        x = val(a,s)
        if x <= n: d[x] += 1
        else: break
def check(x): # a(a+k) = x なる a ある?
    v = int((sqrt(float(k)**2 + 4.0*x) - k)/2)
    for a in range(max(1,v-3),v+5):
        if a >= 1 and a*(a+k)==x: return True
        if a*(a+k) > x: return False
    return False
ans = 0
for x,v in d.items():
    if v == m and not check(x):
        ans += 1
    elif v == m-1 and check(x):
        ans += 1
if m == 1:
    v = int((sqrt(float(k)**2 + 4.0*n) - k)/2)
    for a in range(max(1,v-3),v+5):
        if a*(a+k) > n: break
    ans += a-1
    for x in d:
        if check(x): ans -= 1
print(ans)
            
            
            
        
            
convexineq