login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A085020 a(n) = Sum_{d|n, (d+1) prime} (d + 1). 4
2, 5, 2, 10, 2, 12, 2, 10, 2, 16, 2, 30, 2, 5, 2, 27, 2, 31, 2, 21, 2, 28, 2, 30, 2, 5, 2, 39, 2, 54, 2, 27, 2, 5, 2, 86, 2, 5, 2, 62, 2, 55, 2, 33, 2, 52, 2, 47, 2, 16, 2, 63, 2, 31, 2, 39, 2, 64, 2, 133, 2, 5, 2, 27, 2, 102, 2, 10, 2, 87, 2, 159, 2, 5, 2, 10, 2, 91, 2, 79, 2, 88, 2, 102, 2, 5 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

Michel Marcus, Table of n, a(n) for n = 1..10000

EXAMPLE

a(18) = 31 because the divisors of 18 are [1, 2, 3, 6, 9, 18] and 2 + 3 + 7 + 19 = 31.

MAPLE

T := proc(n, k) local i; numtheory[divisors](n); select(isprime, map(i->i+k, %)); add(i, i=%) end: seq(T(n+1, 1), n=0..20); # Peter Luschny, May 04 2009

MATHEMATICA

a[n_] := Sum[If[PrimeQ[d+1], d+1, 0], {d, Divisors[n]}]; Array[a, 100] (* Jean-Fran├žois Alcover, Jun 04 2019 *)

PROG

(PARI) a(n) = sumdiv(n, d, if (isprime(q=d+1), q)); \\ Michel Marcus, Aug 14 2017

CROSSREFS

Cf. A067513.

Cf. A008472. [Peter Luschny, May 04 2009]

Sequence in context: A119764 A178310 A210697 * A276609 A292591 A240760

Adjacent sequences:  A085017 A085018 A085019 * A085021 A085022 A085023

KEYWORD

easy,nonn

AUTHOR

Jason Earls, Jun 18 2003

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 22 10:13 EST 2022. Contains 350481 sequences. (Running on oeis4.)