cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A272818 Numbers such that (sum + product) of all their prime factors equals (sum + product) of all exponents in their prime factorization.

Original entry on oeis.org

1, 4, 27, 72, 96, 108, 486, 800, 1280, 3125, 6272, 10976, 12500, 14336, 21600, 30375, 36000, 48600, 51840, 54675, 69120, 84375, 121500, 134456, 169344, 174960, 192000, 225000, 240000, 247808, 337500, 340736, 395136, 435456, 451584, 703125, 750141, 781250, 787320, 823543, 857304, 885735
Offset: 1

Views

Author

Giuseppe Coppoletta, May 08 2016

Keywords

Comments

For p prime, p^p satisfy the condition, hence A051674 (and also A048102) is a subsequence. Moreover, if p and q are primes and i and j are positive integers, if p^i * q^j verify the condition, then the same is true for p^j * q^i. So 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. In addition, any number having no more than two distinct prime factors (apart their multiplicity) is a term iff it belongs also to A272858.

Examples

			885735 = 3^11 * 5 is included because (3+5) + 3*5 = (11+1) + 11*1.
2^10 * 3^6 * 19^2 is included because (2+3+19)+ 2*3*19 = (10+6+2)+ 10*6*2.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[10^6], Total@ First@ # + Times @@ First@ # == Total@ Last@ # + Times @@ Last@ # &@ Transpose@ FactorInteger@ # &] (* Michael De Vlieger, May 08 2016 *)
  • PARI
    spp(v) = vecsum(v) + prod(k=1, #v, v[k]);
    isok(n) = my(f = factor(n)); spp(f[,1]) == spp(f[,2]); \\ Michel Marcus, May 08 2016
  • Sage
    def d(n):
        v = factor(n)
        d1 = sum(w[0] for w in v) + prod(w[0] for w in v)
        d2 = sum(w[1] for w in v) + prod(w[1] for w in v)
        return d1 == d2
    [k for k in (1..10000) if d(k)]
    

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

Original entry on oeis.org

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

Views

Author

Giuseppe Coppoletta, May 08 2016

Keywords

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.

Examples

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

Crossrefs

Programs

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

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.
Showing 1-2 of 2 results.