login
a(n) is the common difference in the longest arithmetic progression of primes ending in prime(n). If there is more than one such arithmetic progression, the smallest difference is chosen.
4

%I #8 Aug 14 2024 08:33:53

%S 1,2,2,4,2,6,6,6,6,12,6,12,12,18,12,6,24,24,12,6,6,12,18,18,30,30,18,

%T 6,30,30,30,24,36,48,24,30,12,18,42,6,54,54,42,48,60,30,42,30,66,42,

%U 66,30,60,30,12,6,30,48,84,60,60,78,60,102,60,60,30,78,36,60,90,18,90,6,72,96,30,54

%N a(n) is the common difference in the longest arithmetic progression of primes ending in prime(n). If there is more than one such arithmetic progression, the smallest difference is chosen.

%C a(n) is the smallest common difference in an arithmetic progression of A373888(n) primes ending in prime(n).

%C a(n) is divisible by all primes < min(A373888(n) + 1, prime(n) - (A373888(n)-1)*a(n)).

%H Robert Israel, <a href="/A375386/b375386.txt">Table of n, a(n) for n = 2..10000</a>

%e a(4) = 2 because the 4th prime is 7 and the arithmetic progression of 3 primes ending in 7, namely 3, 5, 7, has common difference 2.

%p f:= proc(n) local s, i, m, dd, d, j;

%p m:= 1;

%p s:= ithprime(n);

%p for i from n-1 to 1 by -1 do

%p d:= s - ithprime(i);

%p if s - m*d < 2 then return dd fi;

%p for j from 2 while isprime(s-j*d) do od;

%p if j > m then m:= j; dd:= d fi;

%p od;

%p dd

%p end proc:

%p map(f, [$2..100]);

%Y Cf. A000040, A373888.

%K nonn,look

%O 2,2

%A _Robert Israel_, Aug 13 2024