login
The values of some algorithm.
2

%I #19 Feb 23 2025 01:59:53

%S 1,1,1,8,30,97,374,2185,7399,60475,303535,2332720,16630294,41419087,

%T 599216396

%N The values of some algorithm.

%C The algorithm: Take n consecutive primes starting with the a(n)-th prime: {p(a(n)), ..., p(a(n)+n-1)}. Calculate the absolute differences abs(p(i)-p(j)) for all relevant pairs i,j. The number of distinct entries of this n X n difference matrix equals binomial(n,2) = n(n-1)/2, the maximum possible. Also a(n) is the earliest index such that this diversity of differences reaches the maximum possible, binomial(n,2).

%C The diagonal [{p(i)-p(i)}] of difference matrices consists only of zeros and can be disregarded.

%C The complete diversity of k-1 consecutive prime differences [as in A079007] is a necessary but not sufficient condition for providing binomial(k,2) distinct entries in the corresponding k X k difference matrix of k consecutive primes. Consecutive prime differences are para-diagonal entries in the difference matrix. So the conditions here are stronger than in A079007.

%C Subscripts at which sequences like A098726, A098206-A098212 reach first their possible maximum, i.e., binomial(k,2) with the corresponding k.

%F a(n) = PrimePi(A099640(n)). - _Amiram Eldar_, Feb 23 2025

%e At n = 1, 2, 3, the maxima are binomial(n,2) = 0, 1, 3 reached at a(n) = 1, 1, 1 respectively.

%e For n = 7: a(7) = 374, primes = {p(374), ..., p(80)} = {2551, 2557, 2579, 2591, 2593, 2609, 2617}. Building the 7 X 7 matrix of values of abs(p(i)-p(j)), the number of its distinct positive entries equals binomial(7,2)=21, namely: {2, 6, 8, 12, 14, 16, 18, 22, 24, 26, 28, 30, 34, 36, 38, 40, 42, 52, 58, 60, 66}.

%e For n = 12: a(12) = 2332720, list of 12 primes = {p(n), ..., p(n+11)} = {38238461, ..., 38238737}. 12 X 12 matrix = {abs((p(i)-p(j))}, number of distinct entries = binomial(12,2) = 66, that of {2, 6, 8, ..., 266, 274, 276}.

%t a[n_] := Module[{k = 1, v = Prime[Range[n]]}, While[CountDistinct@ Flatten@ Abs@ Outer[Plus, v, -v] - 1 != Binomial[n, 2], k++; v = Join[v[[2 ;; -1]], {NextPrime[v[[-1]]]}]]; k]; Array[a, 10] (* _Amiram Eldar_, Feb 23 2025 *)

%o (PARI) s(v) = {my(d = List()); for(i = 1, #v, for(j = 1, i-1, listput(d, abs(v[i] - v[j])))); #Set(d);}

%o a(n) = {my(k = 1, v = primes(n), t = n*(n-1)/2, j = 1); while(s(v) != t, k++; v[j] = nextprime(1 + if(j==1, v[n], v[j-1])); j++; if(j > n, j -= n)); k;} \\ _Amiram Eldar_, Feb 23 2025

%Y Cf. A079007, A080370, A098726, A098206-A098216, A099640.

%K nonn,more,changed

%O 1,4

%A _Labos Elemer_, Oct 05 2004

%E Edited by _Jon E. Schoenfield_, Oct 27 2019

%E a(14)-a(15) from _Amiram Eldar_, Feb 23 2025