OFFSET
1,1
COMMENTS
Contains all positive multiples of 9 (A008591), because 4^n-1 == 0 (mod 27) for these and (4^n-1)/3 is a multiple of 3^2 then. Contains also all positive multiples of 10 (A008592), because 4^n-1 == 0 (mod 125) for these and (4^n-1)/3 is a multiple of 5^2 then. Contains all positive multiples of 21 (A008603), because 4^n-1 == 0 (mod 147) for these and (4^n-1)/3 is a multiple of 7^2 then. - R. J. Mathar, Aug 02 2015
Complement of A259178. - Omar E. Pol, Aug 03 2015
REFERENCES
James R. Buddenhagen, Posting to Math Fun Mailing List, Jul 22 2015.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..169
EXAMPLE
(4^9-1)/3 = 3^2*7*19*73 is not squarefree, so 9 is in the sequence. - R. J. Mathar, Aug 02 2015
MATHEMATICA
Select[Range[120], !SquareFreeQ[(4^#-1)/3]&] (* Ivan N. Ianakiev, Jul 23 2015 *)
PROG
(Magma) [n: n in [1..120]| not IsSquarefree((4^n-1) div 3)]; // Vincenzo Librandi, Jul 27 2015
(PARI) isok(k) = !issquarefree((4^k-1)/3); \\ Michel Marcus, Feb 25 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Jul 22 2015
EXTENSIONS
a(24)-a(31) from Ivan N. Ianakiev, Jul 23 2015
a(32)-a(45) from Chai Wah Wu, Jul 26 2015
a(46)-a(57) from Lars Blomberg, Aug 06 2017
STATUS
approved