|
|
A051903
|
|
Maximal exponent in prime factorization of n.
|
|
239
|
|
|
0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 3, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 5, 1, 2, 2, 2, 1, 1, 1, 3, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
Maximum number of invariant factors among abelian groups of order n. - Álvar Ibeas, Nov 01 2014
a(n) is the highest of the frequencies of the parts of the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1..r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(24) = 3; indeed, the partition having Heinz number 24 = 2*2*2*3 is [1,1,1,2], where the distinct parts 1 and 2 have frequencies 3 and 1, respectively. - Emeric Deutsch, Jun 04 2015
a(n) is the smallest k such that b^(phi(n)+k) == b^k (mod n) for all b.
The Euler phi function can be replaced by the Carmichael lambda function.
Problems:
(*) Are there composite numbers n > 4 such that n == a(n) (mod phi(n))? By Lehmer's totient conjecture, there are no such squarefree numbers.
(**) Are there odd numbers n such that a(n) > 1 and n == a(n) (mod lambda(n))? These are odd numbers n such that a(n) > 1 and b^n == b^a(n) (mod n) for all b.
(***) Are there odd numbers n such that a(n) > 1 and n == a(n) (mod ord_{n}(2))? These are odd numbers n such that a(n) > 1 and 2^n == 2^a(n) (mod n).
Note: if (***) do not exist, then (**) do not exist. (End)
Niven (1969) proved that the asymptotic mean of this sequence is 1 + Sum_{j>=2} 1 - (1/zeta(j)) (A033150). - Amiram Eldar, Jul 10 2020
|
|
LINKS
|
|
|
FORMULA
|
Conjecture: a(n) = a(A003557(n)) + 1. This relation together with a(1) = 0 defines the sequence. - Velin Yanev, Sep 02 2017
This conjecture seems very easily provable to me: if the factorization of n is p1^k1 * p2^k2 * ... * pm^km, then the factorization of the largest squarefree divisor of n is p1 * p2 * ... * pm. So the factorization of A003557(n) is p1^(k1-1) * p2^(k2-1) * ... * pm^(km-1) if exponents of zero are allowed, or with the product terms that have an exponent of zero removed if they're not (if that results in an empty product, consider it to be 1 as usual).
The formula then follows from the fact that provided all ki >= 1, Max(k1, k2, ..., km) = Max(k1-1, k2-1, ..., km-1) + 1, and Max(k1-1, k2-1, ..., km-1) is not altered by removing the ki-1 values that are 0, provided we treat the empty Max() as being 0. That proves the formula and the provisos about empty products and Max() correspond to a(1) = 0.
Also, for any n, applying the formula Max(k1, k2, ..., km) times to n = p1^k1 * p2^k2 * ... * pm^km reduces all the exponents to zero, i.e., to the case a(1) = 0, so that case and the formula generate the sequence. (End)
|
|
EXAMPLE
|
For n = 72 = 2^3*3^2, a(72) = max(exponents) = max(3,2) = 3.
|
|
MAPLE
|
a := 0 ;
for f in ifactors(n)[2] do
a := max(a, op(2, f)) ;
end do:
a ;
# second Maple program:
a:= n-> max(0, seq(i[2], i=ifactors(n)[2])):
|
|
MATHEMATICA
|
Table[If[n == 1, 0, Max @@ Last /@ FactorInteger[n]], {n, 100}] (* Ray Chandler, Jan 24 2006 *)
|
|
PROG
|
(Haskell)
a051903 1 = 0
(Python)
from sympy import factorint
return max(factorint(n).values()) if n > 1 else 0
(Scheme, with memoization-macro definec)
|
|
CROSSREFS
|
Average value is A033150 = 1.7052....
Cf. A002322, A005361, A008479, A028234, A051904, A052409, A067029, A091050, A129132, A327295, A328310, A329885.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|