OFFSET
1,1
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000
MAPLE
with(numtheory): for n from 1 to 110 do div:=divisors(n): A:=[seq(div[j]+1, j=1..tau(n))]: B:={}: for i from 1 to tau(n) do if isprime(A[i])=true then B:=B union {A[i]} else B:=B: fi: od: C:=convert(B, list): a[n]:=C[nops(C)] od: seq(a[n], n=1..110); # Emeric Deutsch, Aug 05 2005
PROG
(PARI) A108077(n) = fordiv(n, d, if(isprime(1+(n/d)), return(1+(n/d)))); \\ Antti Karttunen, Dec 19 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Philippe Deléham, Jul 20 2005
EXTENSIONS
More terms from Emeric Deutsch, Aug 05 2005
STATUS
approved