OFFSET
1,2
COMMENTS
Sequence also is the sum of distinct prime divisors of n which are >= sqrt(n). At most one prime divisor of n is >= square root of n.
LINKS
Diana Mecum, Table of n, a(n) for n = 1..1000
MATHEMATICA
Do[l = Select[Select[Divisors[n], PrimeQ], # >= Sqrt[n]&]; If[Length[l] == 0, Print[0], Print[l[[1]]]], {n, 1, 50}] (* Ryan Propper, Jul 24 2005 *)
Array[Select[FactorInteger[#][[All, 1]], Function[p, p >= Sqrt@ #]] /. {{} -> {0}, {1} -> {0}} &, 83][[All, 1]] (* Michael De Vlieger, Dec 22 2017 *)
PROG
(PARI) a(n) = sumdiv(n, d, if (isprime(d) && (d^2 >= n), d)); \\ Michel Marcus, Dec 23 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Sep 07 2004
EXTENSIONS
More terms from Ryan Propper, Jul 24 2005
More terms from Stefan Steinerberger, Jan 21 2006
Further terms from Diana L. Mecum, Jun 15 2007
STATUS
approved