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”).

Starting with a(0) = 1, smallest squarefree number k such that, for all a(m), m < n, k + a(m) is not squarefree.
5

%I #18 Aug 22 2023 08:00:28

%S 1,3,15,17,233,291,577,723,1455,3615,8117,8835,9505,30833,128773,

%T 130827,239595,273435,426891,654135,676297,926117,1455533,1662533,

%U 2389517,2762427,2820927,7994449,8098527,14319073,16766835,20506733,27606617,31627817,43558023,55566015

%N Starting with a(0) = 1, smallest squarefree number k such that, for all a(m), m < n, k + a(m) is not squarefree.

%e 17 belongs to this sequence as 17 + 1, 17 + 3, 17 + 15 all are divisible by some square.

%t a[0] = 1; a[n_] := a[n] = Module[{t = Array[a, n, 0], k = a[n - 1] + 1}, While[! SquareFreeQ[k] || AnyTrue[t, SquareFreeQ[k + #] &], k++]; k]; Array[a, 20, 0] (* _Amiram Eldar_, Aug 21 2023 *)

%o (PARI) v=vector(60); v[1]=1; print1("1,");for(n=2,60, for(k=1,10^15,if(issquarefree(k),s=0;for(l=1,n-1,if(issquarefree(k+v[l]),break);s=s+1)); if(s==n-1,print1(k",");v[n]=k;break)))

%Y Cf. A013929, A077224, A080793, A085902.

%K nonn

%O 0,2

%A _Amarnath Murthy_, Nov 03 2002

%E Edited by _Ralf Stephan_, Mar 25 2003

%E More terms from _Sam Alexander_, Dec 12 2003

%E a(22) corrected and a(33)-a(35) added by _Amiram Eldar_, Aug 21 2023