結果

問題 No.144 エラトステネスのざる
ユーザー いた
提出日時 2024-12-04 13:26:08
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 319 bytes
コンパイル時間 304 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 65,280 KB
最終ジャッジ日時 2024-12-04 13:31:03
合計ジャッジ時間 2,457 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample RE * 3
other RE * 17
権限があれば一括ダウンロードができます

ソースコード

diff #

import sympy
N,p=map(float,input().split())
N=int(N)

def divisors(n):
  i=0
  ret=1
  while True:
    now=1
    while n%a[i]==0:
      n//=a[i]
      now+=1
    ret*=now
    i+=1
    if n<a[i]:
      return ret
a=list(sympy.primerange(2,10**6))
p=1-p
ans=0
for i in range(2,N+1):
  ans+=pow(p,divisors(i)-2)
print(ans)
0