A048769 -id:A048769 - cp's OEIS frontend

A048767 If n = Product (p_j^k_j) then a(n) = Product ( prime(k_j)^pi(p_j) ) where pi is A000720.

1, 2, 4, 3, 8, 8, 16, 5, 9, 16, 32, 12, 64, 32, 32, 7, 128, 18, 256, 24, 64, 64, 512, 20, 27, 128, 25, 48, 1024, 64, 2048, 11, 128, 256, 128, 27, 4096, 512, 256, 40, 8192, 128, 16384, 96, 72, 1024, 32768, 28, 81, 54, 512, 192, 65536, 50, 256, 80, 1024, 2048

Offset: 1

Naohiro Nomoto

If the prime power factors p^e of n are replaced by prime(e)^pi(p), then the prime terms q in the sequence pertain to 2^m with m > 1, since pi(2) = 1. - Michael De Vlieger, Apr 25 2017

Also the Heinz number of the integer partition obtained by applying the map described in A217605 (which interchanges the parts with their multiplicities) to the integer partition with Heinz number n, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The image of this map (which is the union of this sequence) is A130091. - Gus Wiseman, May 04 2019

			For n=6, 6 = (2^1)*(3^1), a(6) = ([first prime]^pi(2))*([first prime]^pi(3)) = (2^1)*(2^2) = 8.
From _Gus Wiseman_, May 04 2019: (Start)
For n = 1..20, the prime indices of n together with the prime indices of a(n) are the following:
   1: {} {}
   2: {1} {1}
   3: {2} {1,1}
   4: {1,1} {2}
   5: {3} {1,1,1}
   6: {1,2} {1,1,1}
   7: {4} {1,1,1,1}
   8: {1,1,1} {3}
   9: {2,2} {2,2}
  10: {1,3} {1,1,1,1}
  11: {5} {1,1,1,1,1}
  12: {1,1,2} {1,1,2}
  13: {6} {1,1,1,1,1,1}
  14: {1,4} {1,1,1,1,1}
  15: {2,3} {1,1,1,1,1}
  16: {1,1,1,1} {4}
  17: {7} {1,1,1,1,1,1,1}
  18: {1,2,2} {1,2,2}
  19: {8} {1,1,1,1,1,1,1,1}
  20: {1,1,3} {1,1,1,2}
(End)

Alois P. Heinz, Table of n, a(n) for n = 1..10000
Stephan Wagner, The Number of Fixed Points of Wilf's Partition Involution, The Electronic Journal of Combinatorics, 20(4) (2013), #P13.

Cf. A008477, A048768, A048769, A064988.

Cf. A056239, A098859, A109298, A112798, A118914, A130091, A217605, A320348, A325326, A325337.

A048767 := proc(n)
    local a,p,e,f;
    a := 1 ;
    for f in ifactors(n)[2] do
        p := op(1,f) ;
        e := op(2,f) ;
        a := a*ithprime(e)^numtheory[pi](p) ;
    end do:
    a ;
end proc: # R. J. Mathar, Nov 08 2012

Table[{p, k} = Transpose@ FactorInteger[n]; Times @@ (Prime[k]^PrimePi[p]), {n, 58}] (* Ivan Neretin, Jun 02 2016 *)
Array[Apply[Times, FactorInteger[#] /. {p_, e_} /; e >= 0 :> Prime[e]^PrimePi[p]] &, 65] (* Michael De Vlieger, Apr 25 2017 *)

a(1)=1 prepended by Alois P. Heinz, Jul 26 2015

A048768 Numbers n such that A048767(n) = n.

Original entry on oeis.org

1, 2, 9, 12, 18, 40, 112, 125, 250, 352, 360, 675, 832, 1008, 1125, 1350, 1500, 2176, 2250, 2401, 3168, 3969, 4802, 4864, 7488, 7938, 11776, 14000, 19584, 21609, 28812, 29403, 29696, 43218, 43776, 44000, 58806, 63488, 75600, 96040, 104000, 105984, 123201, 126000

Offset: 1

Naohiro Nomoto

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions that are fixed points under the map described in A217605 (which interchanges the parts with their multiplicities). The enumeration of these partitions by sum is given by A217605. - Gus Wiseman, May 04 2019

			12 = (2^2)*(3^1) = (2nd prime)^pi(2) * (first prime)^pi(3).
From _Gus Wiseman_, May 04 2019: (Start)
The sequence of terms together with their prime indices begins:
     1: {}
     2: {1}
     9: {2,2}
    12: {1,1,2}
    18: {1,2,2}
    40: {1,1,1,3}
   112: {1,1,1,1,4}
   125: {3,3,3}
   250: {1,3,3,3}
   352: {1,1,1,1,1,5}
   360: {1,1,1,2,2,3}
   675: {2,2,2,3,3}
   832: {1,1,1,1,1,1,6}
  1008: {1,1,1,1,2,2,4}
  1125: {2,2,3,3,3}
  1350: {1,2,2,2,3,3}
  1500: {1,1,2,3,3,3}
  2176: {1,1,1,1,1,1,1,7}
  2250: {1,2,2,3,3,3}
  2401: {4,4,4,4}
(End)

Amiram Eldar, Table of n, a(n) for n = 1..128

A subsequence of A130091.

Cf. A008478, A048767, A048769, A056239, A098859, A109297, A109298, A112798, A118914, A217605, A325326, A325368.

wt[n_]:=Times@@Cases[FactorInteger[n],{p_,k_}:>Prime[k]^PrimePi[p]];
Select[Range[1000],wt[#]==#&] (* Gus Wiseman, May 04 2019 *)

is(n) = {my(f = factor(n), p = f[, 1], e = f[, 2]); #Set(e) == #e && prod(i = 1, #e, prime(e[i])^primepi(p[i])) == n;} \\ Amiram Eldar, Oct 20 2023

a(1) inserted and more terms added by Amiram Eldar, Oct 20 2023

cp's OEIS Frontend

A048767 If n = Product (p_j^k_j) then a(n) = Product ( prime(k_j)^pi(p_j) ) where pi is A000720.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions