|
|
A067062
|
|
Numbers n such that sigma(n+2) - sigma(n) = prime(n+1) - prime(n).
|
|
1
|
|
|
3, 5, 17, 40, 41, 203, 1949, 2309, 2711, 2789, 2801, 3299, 3329, 3359, 3917, 4157, 4217, 4259, 4637, 5009, 5021, 5231, 6449, 6497, 7757, 8087, 8219, 8627, 9419, 9929, 10007, 10937, 11777, 12071, 14321, 15647, 15971, 16061, 16901, 18131, 18251
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
If n and n+2, prime(n) and prime(n+1) are twin primes, then n is in the sequence. But some values of n are composite.
|
|
LINKS
|
|
|
MAPLE
|
with(numtheory); A067062:=n->`if`(sigma(n+2)-sigma(n)=ithprime(n+1)-ithprime(n), n, NULL); seq(A067062(n), n=1..20000); # Wesley Ivan Hurt, Nov 27 2013
|
|
MATHEMATICA
|
okQ[n_]:=DivisorSigma[1, n+2]-DivisorSigma[1, n]==Prime[n+1]-Prime[n]; Select[Range[20000], okQ] (* Harvey P. Dale, May 24 2011 *)
|
|
PROG
|
(PARI) { default(primelimit, 4294965247); n=0; for (m=1, 10^10, if (sigma(m+2)-sigma(m) == prime(m+1)-prime(m), write("b067062.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, May 02 2010
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|