OFFSET
1,2
COMMENTS
A001223(n) = A001223(A049084(a(n))); a(A001359(n)) = A001359(n+1); conjecture: a(n) > 0 for n > 1 (implies the twin prime conjecture). - Reinhard Zumkeller, Jan 26 2004
For n > 1, a(n) >= prime(n+1) and a(n) = prime(n+1) if prime(n+1) is a balanced prime (A006562). - Zak Seidov, Jun 03 2015
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Prime Difference Function
MATHEMATICA
a[1] = 0; a[n_] := For[p = Prime[n+1]; d = p - Prime[n], True, p = q, q = NextPrime[p]; If[d == q - p, Return[p]]]; (* Jean-François Alcover, Feb 24 2015 *)
PROG
(MATLAB) P = primes(5000); A = zeros(1, length(P)); D = P(2:end) - P(1:(length(P) - 1)); for i = 2:2:(max(D)); f = find(D == i); A(f(1:(length(f) - 1))) = P(f(2:end)); end; A(2:100) % David Wasserman, Jan 26 2004
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Jul 09 2003
EXTENSIONS
More terms from Reinhard Zumkeller and David Wasserman, Jan 26 2004
Edited by N. J. A. Sloane, Oct 21 2008 at the suggestion of R. J. Mathar
STATUS
approved
