|
%I #16 Sep 08 2022 08:45:51
%S 17,19,167,162287,159463,354751,2486784401,12589253,71549927,
%T 83695120256591,556123,10902870290089,25793,362561,26757560905578361,
%U 6015993258685545623,55044966589,1182622403647894934149,9973660056412561,38783069052902219,4576344458074395243073
%N a(n) is the largest prime factor of (n-1)^n - n^(n-1).
%C A007925(n) = (n-1)^n - n^(n-1).
%C A007925(n) = -1 for n = 1..3, which has no prime factors, so the offset is 4.
%H Vincenzo Librandi, <a href="/A174379/b174379.txt">Table of n, a(n) for n = 4..50</a>
%e 19 is a term because 4^5 - 5^4 = 399 = 19 * 21. - _Bernard Schott_, Jan 14 2019
%t Table[FactorInteger[(n - 1)^n - n^(n - 1)][[-1, 1]], {n, 4, 20}] (* _Vincenzo Librandi_, Jan 14 2019 *)
%o (PARI) a(n) = vecmax(factor((n-1)^n - n^(n-1))[,1]);
%o for (n=4, 30, print1(a(n), ", ")) \\ _Michel Marcus_, Jan 14 2019
%o (Magma) [Maximum(PrimeDivisors((n-1)^n - n^(n-1))): n in [4..25]]; // _Vincenzo Librandi_, Jan 14 2019
%Y Cf. A007925.
%K nonn
%O 4,1
%A Torbjorn Alm (talm(AT)tele2.se), Mar 17 2010
%E a(4) corrected, a(22)-a(24) added, and other edits by _Jon E. Schoenfield_, Jan 13 2019
|
