OFFSET
1,1
COMMENTS
Primes of the form m*tau(m) - sigma(m), listed in the order in which the values of m appear in A096847.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
Do[s=n*DivisorSigma[0, n]-DivisorSigma[1, n]; If[PrimeQ[s], Print[{n, s}]; ta[[u]]=n; tb[[u]]=s; u=u+1], {n, 1, 1000000}]; ta
s[n_] := Module[{f = FactorInteger[n], p, e}, p = f[[;; , 1]]; e = f[[;; , 2]]; n * Times @@ (e + 1) - Times @@ ((p^(e + 1) - 1)/(p - 1))]; q[n_] := PrimeQ[s[n]]; seq[lim_] := Module[{m1 = Floor[Sqrt[lim/2]], m2 = Floor[Sqrt[lim]/2]}, s /@ Join[{3}, Union[Select[2*Range[m1]^2, q], Select[4*Range[m2]^2, q]]]]; seq[200000] (* Amiram Eldar, Feb 14 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jul 15 2004
STATUS
approved