OFFSET
1,2
COMMENTS
From Robert Israel, Aug 13 2020: (Start)
Except for 1, all terms are divisible by 3, but not 5 or 9.
All terms > 3 are divisible by at least one of 7, 19 and 73.
Are all terms squarefree? (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..500 (first 49 terms from Robert Price)
MAPLE
filter:= n -> 4&^(n^2)-1 mod (n^3) = 0:
select(filter, [1, seq(i, i=3..10^6, 6)]); # Robert Israel, Aug 13 2020
MATHEMATICA
k=2; Do[ p=Prime[k]; If[ IntegerQ[ (PowerMod[ p+1, n^2, n^3 ] - 1 )/n^3 ], Print[ {k, p, n} ]], {n, 1, 200000} ]
k=2; Do[ p=Prime[k]; If[ IntegerQ[ (PowerMod[ p+1, n^2, n^3 ] - 1 )/n^3 ], Print[ {k, p, n} ]], {n, 1000000} ] (* Robert G. Wilson v, Apr 06 2007 *)
Join[{1}, Select[Range[3000000], PowerMod[4, #^2, #^3] == 1 &]] (* Robert Price, Mar 31 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Apr 03 2007
EXTENSIONS
More terms from Robert G. Wilson v, Apr 06 2007
STATUS
approved