OFFSET
2,1
COMMENTS
Theorem I. If n is prime and the prime factors of b^n - (b-1)^n are p_1, p_2, p_3, ..., p_i & i <=1, then p_i == 1 (mod n) for all i's. This is a result from Fermat's little theorem.
Theorem II. If n is not prime and n = d_1 * d_2 * d_3 * ... * d_j and the prime factors of b^n - (b-1)^n are p_1, p_2, p_3, ..., p_i & i < 1, then some prime factors, q_1, q_2, q_3, ..., q_k & k <= 1 are not primitive, i.e., they are prime factors of b^d_j - (b-1)^d_j. Excluding from the list of p_i those which are not primitive, then the rest are also == 1 (mod n). In fact, these two theorems may be generalized for a and b, (a,b) and (a^n +- b^n)/(a +- b).
LINKS
Chai Wah Wu, Table of n, a(n) for n = 2..598
EXAMPLE
3^9 - 2^9 = 19*1009.
3^17 - 2^17 = 129009091.
MATHEMATICA
PrimeFactors[n_] := Flatten[ Table[ #[[1]], {1}] & /@ FactorInteger[n]]; Do[ Print[ PrimeFactors[3^n - 2^n][[1]]], {n, 2, 60}] (* Robert G. Wilson v, Dec 05 2003 *)
PROG
(PARI) leastfactor(a, n) = { for(x=2, n, y = a^x-(a-1)^x; f = factor(y); v = component(component(f, 1), 1); print1(v", ") ) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, Dec 04 2003
STATUS
approved