A088902 Numbers n such that n = product (p_k)^(c_k) and set of its (c_k k's)'s is a self-conjugate partition, where p_k is k-th prime and c_k > 0.
1, 2, 6, 9, 20, 30, 56, 75, 84, 125, 176, 210, 264, 350, 416, 441, 624, 660, 735, 1088, 1100, 1386, 1560, 1632, 1715, 2310, 2401, 2432, 2600, 3267, 3276, 3648, 4080, 5390, 5445, 5460, 5888, 6800, 7546, 7722, 8568, 8832, 9120, 12705, 12740, 12870, 13689
Offset: 1
Examples
20 is in the sequence because 20 = 2^2 * 5^1 = (p_1)^2 *(p_3)^1, (two 1's, one 3's) = (1,1,3) is a self-conjugate partition of 5. From _Gus Wiseman_, Jun 28 2022: (Start) The terms together with their prime indices begin: 1: () 2: (1) 6: (2,1) 9: (2,2) 20: (3,1,1) 30: (3,2,1) 56: (4,1,1,1) 75: (3,3,2) 84: (4,2,1,1) 125: (3,3,3) 176: (5,1,1,1,1) 210: (4,3,2,1) 264: (5,2,1,1,1) (End)
Links
- Alois P. Heinz, Table of n, a(n) for n = 1..400
Crossrefs
Fixed points of A122111.
A002110 (primorial numbers) is a subsequence.
After a(1) and a(2), a subsequence of A241913.
These partitions are counted by A000700.
The same count comes from A258116.
These are the positions of zeros in A352491.
Heinz number (rank) and partition:
- A056239 = sum of partition.
Programs
-
Maple
with(numtheory): c := proc (n) local B, C: B := proc (n) local pf: pf := op(2, ifactors(n)): [seq(seq(pi(op(1, op(i, pf))), j = 1 .. op(2, op(i, pf))), i = 1 .. nops(pf))] end proc: C := proc (P) local a: a := proc (j) local c, i: c := 0: for i to nops(P) do if j <= P[i] then c := c+1 else end if end do: c end proc: [seq(a(k), k = 1 .. max(P))] end proc: mul(ithprime(C(B(n))[q]), q = 1 .. nops(C(B(n)))) end proc: SC := {}: for i to 14000 do if c(i) = i then SC := `union`(SC, {i}) else end if end do: SC; # Emeric Deutsch, May 09 2015
-
Mathematica
Select[Range[14000], Function[n, n == If[n == 1, 1, Module[{l = #, m = 0}, Times @@ Power @@@ Table[l -= m; l = DeleteCases[l, 0]; {Prime@ Length@ l, m = Min@ l}, Length@ Union@ l]] &@ Catenate[ConstantArray[PrimePi@ #1, #2] & @@@ FactorInteger@ n]]]] (* Michael De Vlieger, Aug 27 2016, after JungHwan Min at A122111 *)
Extensions
More terms from David Wasserman, Aug 26 2005
Comments