OFFSET
1,1
COMMENTS
This sequence has positive density in N; the density is around 0.0025.
The sequence includes an infinite family of arithmetic progressions. Such AP's can be constructed to each term, with large differences [like e.g. square of primorials, A061742]. It is necessary to solve suitable systems of linear Diophantine equations. E.g.: subsequences of quadruples of terms = {44100k+29349, 44100k+29350, 44100k+29351, 44100k+29352} = {9(4900k+3261), 25(1764k+1174), 49(900k+599), 4(11025k+7338)}; starting terms in this sequence = {29349, 73449, 117549, ...}; difference = A002110(4)^2 = 210^2. - Labos Elemer, Nov 25 2002
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
f[n_] := Union[Transpose[FactorInteger[n]][[2]]][[ -1]]; a = 0; b = 1; c = 0; Do[d = f[n]; If[a > 1 && b > 1 && c > 1 && d > 1, Print[n - 3]]; a = b; b = c; c = d, {n, 4, 10^6}]
Flatten[Position[Partition[SquareFreeQ/@Range[60000], 4, 1], _?(Union[#] == {False}&), {1}, Heads->False]] (* Harvey P. Dale, May 24 2014 *)
PROG
(PARI) is(n)=for(i=n, n+3, if(!issquarefree(n), return(0))); 1 \\ Charles R Greathouse IV, Sep 14 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Sharon Sela (sharonsela(AT)hotmail.com), May 09 2002
EXTENSIONS
More terms from Robert G. Wilson v, May 09 2002
b-file from Charles R Greathouse IV, Jul 23 2010
STATUS
approved