OFFSET
1,1
COMMENTS
m such that m+3 divides (-3)^m - 3. - Robert Israel, Dec 14 2014
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..96
EXAMPLE
3 is in this sequence because 3 + 3 = 6 divides 3^3 - 3 = 24.
MAPLE
select(t ->((-3) &^ (t) - 3) mod (t+3) = 0, [$1..10^6]); # Robert Israel, Dec 14 2014
MATHEMATICA
a251862[n_] := Select[Range[n], Mod[PowerMod[#, #, # + 3] - 3, # + 3] == 0 &]; a251862[10^6] (* Michael De Vlieger, Dec 14 2014, after Robert G. Wilson v at A252041 *)
PROG
(Magma) [n: n in [2..10000] | Denominator((n^n-3)/(n+3)) eq 1];
(PARI) isok(n) = Mod(n, n+3)^n == 3; \\ Michel Marcus, Dec 10 2014
(Sage)
[n for n in range(10^4) if (n + 3).divides((-3)^n - 3)] # Peter Luschny, Jan 17 2015
(Python)
A251862_list = [n for n in range(10**6) if pow(-3, n, n+3) == 3] # Chai Wah Wu, Jan 19 2015
CROSSREFS
Cf. ...............Numbers n such that x divides y, where:
...x.....y......k=0.......k=1.......k=2........k=3........
(For x=n-1 and y=n^n+1, the only terms are 0, 2 and 3. - David L. Harden, Jan 14 2015)
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Dec 10 2014
EXTENSIONS
More terms from Michel Marcus, Dec 10 2014
STATUS
approved