fork(1) download
  1. a=->n,d=1{d<67?a[n,d+=2]-(2**d%d==2?d*(a[n/=d,d]+n*~n/2):0):0}
  2. p a[10**5],a[10**9]
Success #stdin #stdout 3s 7452KB
stdin
Standard input is empty
stdout
3712013581
370375739387672112