OFFSET
1,2
COMMENTS
This is the union of Sophie Germain primes and Sophie Germain nonprimes, so it might be called "Sophie Germain numbers".
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
J. S. Gerasimov, Sophie Germain nonprimes [title corrected], SeqFan mailing list, Jan 15 2013.
MAPLE
A210495 := proc(n)
option remember;
local a;
if n = 1 then
1 ;
else
for a from procname(n-1)+1 do
if isprime(numtheory[tau](a)*a+1) then
return a;
end if;
end do:
end if;
end proc: # R. J. Mathar, Jan 27 2013
MATHEMATICA
Select[Range[200], PrimeQ[# DivisorSigma[0, #]+1]&] (* Harvey P. Dale, Aug 26 2013 *)
PROG
(PARI) is(n)=isprime(numdiv(n)*n+1) \\ Charles R Greathouse IV, Jan 24 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jan 24 2013
STATUS
approved