A272818 -id:A272818 - cp's OEIS frontend

A272858 Numbers m such that Product(1 + p_i) = Product(1 + e_i), where m = Product((p_i)^e_i).

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

Offset: 1

Giuseppe Coppoletta, May 08 2016

nonn

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.

			92160 is included because 92160 = 2^11 * 3^2 * 5 and (2+1)*(3+1)*(5+1) = (11+1)*(2+1)*(1+1).

Giuseppe Coppoletta and Giovanni Resta, Table of n, a(n) for n = 1..6058 (terms < 10^18, first 100 terms from G. Coppoletta)

Cf. A048102, A054411, A054412, A071174, A114129, A122406, A272818, A272859.

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 *)

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

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)]

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.

A272859 Numbers m such that sigma(Product(p_j)) = sigma(Product(e_j)), where m = Product((p_i)^e_i) and sigma = A000203.

Original entry on oeis.org

1, 4, 27, 72, 108, 192, 800, 1458, 3125, 5120, 6144, 6272, 10976, 12500, 21600, 30375, 36000, 48600, 54675, 77760, 84375, 114688, 116640, 121500, 134456, 138240, 169344, 173056, 225000, 229376, 247808, 337500, 354294, 384000, 395136, 600000, 653184, 655360, 703125, 750141, 823543, 857304, 913952, 979776

Offset: 1

Giuseppe Coppoletta, May 08 2016

nonn

A048102 is clearly a subsequence, as for any prime p, p^p satisfy the herein condition. Moreover, due to the multiplicativity of the arithmetic function sigma, 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 A272858 at least for the terms of A114129.

			173056 is included because 173056 = 2^10 * 13^2 and sigma(2*13) = sigma(10*2).
653184 is included because 653184 = 2^7 * 3^6 * 7 and sigma(2*3*7) = sigma(7*6*1).

Giuseppe Coppoletta and Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 100 terms from G. Coppoletta)

Cf. A048102, A054411, A054412, A071174, A114129, A122406, A272818, A272858.

Select[Range[10^6], First@ # == Last@ # &@ Map[DivisorSigma[1, Times @@ #] &, Transpose@ FactorInteger@ #] &] (* Michael De Vlieger, May 12 2016 *)

A272859 = []
for n in (1..10000):
    v = factor(n)
    if prod(1 + w[0] for w in v) == sigma(prod(w[1] for w in v)): A272859.append(n)
print(A272859)

cp's OEIS Frontend

A272858 Numbers m such that Product(1 + p_i) = Product(1 + e_i), where m = Product((p_i)^e_i).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula