%I #44 Jun 28 2024 18:10:28
%S 1,2,4,5,6,9,10,85,105,136,186,262,820,1161,2626,2926,4924,10396,
%T 11656,19689,27637,33736,36046,42886,42901,53866,55189,82741,95266,
%U 103762,106822,127401,135460,251506,366796,375220,413326,466966,531445,553456,568876
%N Numbers m such that m - 3 divides m^m + 3.
%C Numbers m such that (m^m + 3)/(m - 3) is an integer.
%C Most but not all terms are congruent to 4 modulo 6. - _Robert G. Wilson v_, Dec 19 2014
%C Note that m^m == 3^m (mod m-3). - _Robert Israel_, Dec 19 2014
%H Robert G. Wilson v, <a href="/A252041/b252041.txt">Table of n, a(n) for n = 1..331</a>
%e 2 is in this sequence because (2^2 + 3)/(2 - 3) = -7 is an integer.
%e 4 is in this sequence because (4^4 + 3)/(4 - 3) = 259 is an integer.
%e 7 is not in the sequence because (7^7 + 3)/4 = 411773/2, which is not an integer.
%p select(t -> 3 &^t + 3 mod (t-3) = 0, [1,2,$4..10^6]); # _Robert Israel_, Dec 19 2014
%t fQ[n_] := Mod[PowerMod[n, n, n - 3] + 3, n - 3] == 0; Select[Range@ 1000000, fQ] (* _Michael De Vlieger_, Dec 13 2014; modified by _Robert G. Wilson v_, Dec 19 2014 *)
%o (Magma) [n: n in [4..50000] | Denominator((n^n+3)/(n-3)) eq 1];
%o (PARI) isok(n) = (n != 3) && (Mod(n, n-3)^n == -3); \\ _Michel Marcus_, Dec 13 2014
%Y Cf. ...............Numbers n such that x divides y, where:
%Y ...x......y....k = 0.....k = 1.....k = 2......k = 3.......
%Y ..n-k..n^n-k..A000027...A087156...A242787....A242788......
%Y ..n-k..n^n+k..A000027..see below..A249751..this sequence..
%Y ..n+k..n^n-k..A000027...A004275...A251603....A251862......
%Y ..n+k..n^n+k..A000027...A004273...A213382....A242800......
%Y (For x=n-1 and y=n^n+1, the only terms are 0, 2 and 3. - _David L. Harden_, Dec 28 2014}
%K nonn
%O 1,2
%A _Juri-Stepan Gerasimov_, Dec 12 2014
%E More terms from _Michel Marcus_, Dec 13 2014