A306812
Maximally idempotent integers with three or more factors.
Original entry on oeis.org
273, 455, 1729, 2109, 2255, 2387, 3367, 3515, 4433, 4697, 4921, 5673, 6643, 6935, 7667, 8103, 8723, 8729, 9139, 9455, 10235, 10787, 11543, 13237, 13505, 14497, 16211, 16385, 16523, 17507, 18031, 18907, 20033, 20801, 21437, 22649, 23579, 24583
Offset: 1
273 is the smallest maximally idempotent integer. Factorization is (3,7,13). Bipartite factorizations are (3,91), (7,39), (13,21). Lambda(273) = 12, (2*90),(6*38) and (12*20) are all divisible by 12, thus 273 is maximally idempotent. The same is true for 455 = 5*7*13. The next entry in the sequence, 1729=7*13*19, is a Carmichael number, but most Carmichael numbers are not maximally idempotent.
Subsequence of
A120944 (composite squarefree numbers). Subsequence of
A306330 (squarefree numbers that admit idempotent factorizations). Members of the sequence with >= 4 factors for a subsequence of
A306508 (squarefree integers with fully composite idempotent factorizations).
-
## This uses a custom library of number theory functions and the numbthy library.
## Hopefully the names of the functions make the process clear.
for n in range(2,max_n):
factor_list = numbthy.factor(n)
numFactors = len(factor_list)
if numFactors <= 2: # skip primes and semiprimes
continue
if not bsflib.is_composite_and_square_free_with_list(n,factor_list):
continue
ipList = bsflib.idempotentPartitions(n, factor_list)
if len(ipList) == 2**(numFactors-1)-1:
print(n)
A307537
a(n) is the smallest maximally idempotent integer with n factors, n >= 3.
Original entry on oeis.org
273, 63973, 72719023, 13006678091, 7817013532691
Offset: 3
273 is the smallest maximally idempotent integer. Factorization is (3,7,13). Bipartite factorizations are (3,91), (7,39), (13,21). Lambda(273) = 12, (2*90),(6*38) and (12*20) are all divisible by 12, thus 273 is maximally idempotent.
-
(* This program is not suitable to compute large terms. *)
okQ[n_] := Module[{partitions, p, q, lambda}, partitions = {p, q} /. {ToRules[Reduce[1= 3 && !IntegerQ[a[nu]], If[okQ[n], Print["a(", nu, ") = ", n]; a[nu] = n]]]]; (* Jean-François Alcover, Jun 20 2019 *)
-
# Partial Python code is shown below. It uses other routines:
# numbthy.factor(n) -- from the Python number theory library, returns a list of
# (p,e) pairs corresponding to the prime factors and their exponents in the factorizations of n
# partitions(n,factor_list) -- takes an integer n and the factor list from above,
# returns a list of all bipartite factorizations of n
# lambda_n -- calculates the carmichael lambda function
# returns True if all partitions of n are idempotent
def isMaximallyIdempotent(n):
factor_list = numbthy.factor(n)
partitions_of_n = partitions(n,factor_list)
lambda_n = carmichael_lambda_with_list(n,factor_list)
for (p,q) in partitions_of_n:
pseudo = (p-1)*(q-1)
if pseudo % lambda_n != 0:
return False
return True
M(7), now confirmed as being a(7), added by
Barry Fagin, Dec 04 2019
Showing 1-2 of 2 results.
Comments