login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(1)=2; thereafter, a(n) is the smallest number not occurring earlier such that Kronecker(a(k), a(n)) = 1 for the next n indices k = n+1, n+2, ..., 2n.
1

%I #16 Oct 02 2017 02:12:10

%S 2,7,3,4,11,25,9,19,16,49,17,59,67,121,81,169,43,115,64,191,293,361,

%T 289,587,269,841,863,961,1031,1369,441,1681,1867,2209,1849,2809,65,

%U 529,256,643,3481,3517,1639,1877,3721,4489,5041,5329,5591,6241,3557,6889,7921

%N a(1)=2; thereafter, a(n) is the smallest number not occurring earlier such that Kronecker(a(k), a(n)) = 1 for the next n indices k = n+1, n+2, ..., 2n.

%C Kronecker(i,j) is an extension of the Jacobi symbol to all integers. The sequence with the condition Kronecker(a(i), a(n)) = -1/+1 is given by A247665.

%e a(1) = 2 because the next term is 7 and k(2,7) = 1;

%e a(2) = 7 because the next two terms are (3,4) => k(7,3) = 1 and k(7,4) = 1;

%e a(3) = 3 because the next three terms are (4,11,25) => k(3,4)= 1, k(3,11) = 1 and k(3,25) = 1.

%o (PARI) m=55; v=vector(m); u=vectorsmall(1000*m); for(n=1, m, for(i=2, 10^9, if(!u[i], for(j=(n+1)\2, n-1, if(kronecker(v[j], i)==-1 || kronecker(v[j], i)==0, next(2))); v[n]= i; u[i]=1; break))); v

%Y Cf. A247665, A249692.

%K nonn

%O 1,1

%A _Michel Lagneau_, Nov 05 2014