login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A272858 Numbers m such that Product(1 + p_i) = Product(1 + e_i), where m = Product((p_i)^e_i). 3
1, 4, 27, 72, 96, 108, 486, 800, 1280, 3125, 6272, 10976, 12500, 14336, 21600, 28800, 30375, 34560, 36000, 38880, 48600, 54675, 84375, 92160, 96000, 121500, 134456, 153600, 169344, 217728, 218700, 225000, 247808, 262440, 296352, 300000, 337500, 340736, 387072, 395136, 489888, 666792, 703125, 750141, 781250, 823543, 857304, 885735 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
A048102 is clearly a subsequence, as for any prime p, p^p satisfy the herein condition. Similarly, A122406 is also a subsequence. More generally, if a number is a term, then any permutation of the exponents in its prime factorization (i.e., any permutation of its prime signature) gives also a term.
The condition defining this sequence coincides with the condition in A272859 at least for the terms of A114129.
LINKS
Giuseppe Coppoletta and Giovanni Resta, Table of n, a(n) for n = 1..6058 (terms < 10^18, first 100 terms from G. Coppoletta)
FORMULA
If N is a positive integer and N = Product_{i=1..k} (p_i)^e_i is its prime factorization, then N is in A272858 iff Product_{i=1..k} (1 + p_i) = Product_{i=1..k} (1 + e_i).
For a number with three different prime factors N = p1^e1 * p2^e2 * p3^e3, the defining condition can be expressed as: p1 + p2 + p3 + p1*p2 + p1*p3 + p2*p3 + p1*p2*p3 = e1 + e2 + e3 + e1*e2 + e1*e3 + e2*e3 + e1*e2*e3.
EXAMPLE
92160 is included because 92160 = 2^11 * 3^2 * 5 and (2+1)*(3+1)*(5+1) = (11+1)*(2+1)*(1+1).
MATHEMATICA
ok[n_] := Block[{p, e}, {p, e} = Transpose@ FactorInteger@ n; Times @@ (1+p) == Times @@ (1+e)]; Select[Range[10^6], ok] (* Giovanni Resta, May 08 2016 *)
PROG
(Sage)
def d(n):
v = factor(n)
d1 = prod(1 + w[0] for w in v)
d2 = prod(1 + w[1] for w in v)
return d1 == d2
[k for k in (1..10000) if d(k)]
(PARI) is(n)=my(f=factor(n)); prod(i=1, #f~, f[i, 1]+1)==prod(i=1, #f~, f[i, 2]) \\ Charles R Greathouse IV, Sep 08 2016
CROSSREFS
Sequence in context: A356814 A180576 A158186 * A272818 A356433 A272859
KEYWORD
nonn
AUTHOR
Giuseppe Coppoletta, May 08 2016
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 17:25 EDT 2024. Contains 371254 sequences. (Running on oeis4.)