A299757 -id:A299757 - cp's OEIS frontend

A302242 Total weight of the n-th multiset multisystem. Totally additive with a(prime(n)) = Omega(n).

0, 0, 1, 0, 1, 1, 2, 0, 2, 1, 1, 1, 2, 2, 2, 0, 1, 2, 3, 1, 3, 1, 2, 1, 2, 2, 3, 2, 2, 2, 1, 0, 2, 1, 3, 2, 3, 3, 3, 1, 1, 3, 2, 1, 3, 2, 2, 1, 4, 2, 2, 2, 4, 3, 2, 2, 4, 2, 1, 2, 3, 1, 4, 0, 3, 2, 1, 1, 3, 3, 3, 2, 2, 3, 3, 3, 3, 3, 2, 1, 4, 1, 1, 3, 2, 2, 3, 1, 4

Offset: 1

Gus Wiseman, Apr 03 2018

nonn

A multiset multisystem is a finite multiset of finite multisets of positive integers. The n-th multiset multisystem is constructed by factoring n into prime numbers and then factoring each prime index into prime numbers and taking their prime indices. This produces a unique multiset multisystem for each n, and every possible multiset multisystem is so constructed as n ranges over all positive integers.

			Sequence of finite multisets of finite multisets of positive integers begins: (), (()), ((1)), (()()), ((2)), (()(1)), ((11)), (()()()), ((1)(1)), (()(2)), ((3)), (()()(1)), ((12)), (()(11)), ((1)(2)), (()()()()), ((4)), (()(1)(1)), ((111)), (()()(2)).

Alois P. Heinz, Table of n, a(n) for n = 1..65536

Cf. A001222, A003963, A007716, A034691, A056239, A061775, A063834, A096443, A249620, A255397, A255906, A275024, A279789, A281113, A299757, A302243.

with(numtheory):
a:= n-> add(add(j[2], j=ifactors(pi(i[1]))[2])*i[2], i=ifactors(n)[2]):
seq(a(n), n=1..100);  # Alois P. Heinz, Sep 07 2018

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Total[PrimeOmega/@primeMS[n]],{n,100}]

a(n,f=factor(n))=sum(i=1,#f~, bigomega(primepi(f[i,1]))*f[i,2]) \\ Charles R Greathouse IV, Nov 10 2021

A300061 Heinz numbers of integer partitions of even numbers.

Original entry on oeis.org

1, 3, 4, 7, 9, 10, 12, 13, 16, 19, 21, 22, 25, 27, 28, 29, 30, 34, 36, 37, 39, 40, 43, 46, 48, 49, 52, 53, 55, 57, 61, 62, 63, 64, 66, 70, 71, 75, 76, 79, 81, 82, 84, 85, 87, 88, 89, 90, 91, 94, 100, 101, 102, 107, 108, 111, 112, 113, 115, 116, 117, 118, 120

Offset: 1

Gus Wiseman, Feb 23 2018

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			75 is the Heinz number of (3,3,2), which has even weight, so 75 belongs to the sequence.
Sequence of even-weight partitions begins: () (2) (1,1) (4) (2,2) (3,1) (2,1,1) (6) (1,1,1,1) (8) (4,2) (5,1) (3,3) (2,2,2) (4,1,1).

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Complement of A300063.

Cf. A000041, A000720, A001222, A056239, A063834, A100118, A112798, A122111, A215366, A296150, A299202, A299757, A300056, A300060, A304662.

a:= proc(n) option remember; local k; for k from 1+
     `if`(n=1, 0, a(n-1)) while add(numtheory[pi]
      (i[1])*i[2], i=ifactors(k)[2])::odd do od; k
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, May 22 2018

Select[Range[200],EvenQ[Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]]&]

A052330 Let S_k denote the first 2^k terms of this sequence and let b_k be the smallest positive integer that is not in S_k; then the numbers b_k*S_k are the next 2^k terms.

Original entry on oeis.org

1, 2, 3, 6, 4, 8, 12, 24, 5, 10, 15, 30, 20, 40, 60, 120, 7, 14, 21, 42, 28, 56, 84, 168, 35, 70, 105, 210, 140, 280, 420, 840, 9, 18, 27, 54, 36, 72, 108, 216, 45, 90, 135, 270, 180, 360, 540, 1080, 63, 126, 189, 378, 252, 504, 756, 1512, 315, 630, 945, 1890

Offset: 0

Christian G. Bower, Dec 15 1999

Inverse of sequence A064358 considered as a permutation of the positive integers. - Howard A. Landman, Sep 25 2001
This sequence is not exactly a permutation because it has offset 0 but doesn't contain 0. A052331 is its exact inverse, which has offset 1 and contains 0. See also A064358.
Are there any other values of n besides 4 and 36 with a(n) = n? - Thomas Ordowski, Apr 01 2005
4 = 100 = 4^1 * 3^0 * 2^0, 36 = 100100 = 9^1 * 7^0 * 5^0 * 4^1 * 3^0 * 2^0. - Thomas Ordowski, May 26 2005
Ordering of positive integers by increasing "Fermi-Dirac representation", which is a representation of the "Fermi-Dirac factorization", term implying that each prime power with a power of two as exponent may appear at most once in the "Fermi-Dirac factorization" of n. (Cf. comment in A050376; see also the OEIS Wiki page.) - Daniel Forgues, Feb 11 2011
The subsequence consisting of the squarefree terms is A019565. - Peter Munn, Mar 28 2018
Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. The FDH-number of a strict integer partition (y_1,...,y_k) is f(y_1)*...*f(y_k). A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793. Then a(n) is the number whose binary indices are the parts of the strict integer partition with FDH-number n. - Gus Wiseman, Aug 19 2019
The set of indices of odd-valued terms has asymptotic density 0. In this sense (using the order they appear in this permutation) 100% of numbers are even. - Peter Munn, Aug 26 2019

			Terms following 5 are 10, 15, 30, 20, 40, 60, 120; this is followed by 7 as 6 has already occurred. - _Philippe Deléham_, Jun 03 2015
From _Antti Karttunen_, Apr 13 2018, after also _Philippe Deléham_'s Jun 03 2015 example: (Start)
This sequence can be regarded also as an irregular triangle with rows of lengths 1, 1, 2, 4, 8, 16, ..., that is, it can be represented as a binary tree, where each left hand child contains A300841(k), and each right hand child contains 2*A300841(k), when their parent contains k:
                                     1
                                     |
                  ...................2...................
                 3                                       6
       4......../ \........8                  12......../ \........24
      / \                 / \                 / \                 / \
     /   \               /   \               /   \               /   \
    /     \             /     \             /     \             /     \
   5       10         15       30         20       40         60      120
  7 14   21  42     28  56   84  168    35  70  105  210   140 280  420 840
  etc.
Compare also to trees like A005940 and A283477, and sequences A207901 and A302783.
(End)

Antti Karttunen, Table of n, a(n) for n = 0..8191 (Terms 0..1023 from T. D. Noe)
Michael De Vlieger, Fan style binary tree showing a(n), n = 0..16383, with a color function representing primes in red, perfect powers of primes in gold, squarefree composites in greens, and numbers neither squarefree nor prime powers in blue or purple. Purple represents powerful numbers that are not prime powers. This is a 15 level version of Karttunen's diagram in the example.
OEIS Wiki, Ordering of positive integers by increasing "Fermi-Dirac representation"
Index entries for sequences that are permutations of the natural numbers

Subsequences: A019565 (squarefree terms), A050376 (the left edge from 2 onward), A336882 (odd terms).
Cf. A050030, A052331 (inverse), A096111, A096113, A096114, A096115, A096116, A096118, A096119, A182979, A207901, A300841, A302023, A302783.
Cf. A000120, A029931, A048793, A064547, A070939, A213925, A299755, A299757, A327041.

a = {1}; Do[a = Join[a, a*Min[Complement[Range[Max[a] + 1], a]]], {n, 1, 6}]; a (* Ivan Neretin, May 09 2015 *)

up_to_e = 13; \\ Good for computing up to n = (2^13)-1
v050376 = vector(up_to_e);
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to_e,break));
A050376(n) = v050376[n];
A052330(n) = { my(p=1,i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); }; \\ Antti Karttunen, Apr 12 2018

a(0)=1; a(n+2^k)=a(n)*b(k) for n < 2^k, k = 0, 1, ... where b is A050376. - Thomas Ordowski, Mar 04 2005
The binary representation of n, n = Sum_{i=0..1+floor(log_2(n))} n_i * 2^i, n_i in {0,1}, is taken as the "Fermi-Dirac representation" (A182979) of a(n), a(n) = Product_{i=0..1+floor(log_2(n))} (b_i)^(n_i) where b_i is A050376(i), i.e., the i-th "Fermi-Dirac prime" (prime power with exponent being a power of 2). - Daniel Forgues, Feb 12 2011
From Antti Karttunen, Apr 12 & 17 2018: (Start)
a(0) = 1; a(2n) = A300841(a(n)), a(2n+1) = 2*A300841(a(n)).
a(n) = A207901(A006068(n)) = A302783(A003188(n)) = A302781(A302845(n)).
(End)

Entry revised Mar 17 2005 by N. J. A. Sloane, based on comments from several people, especially David Wasserman and Thomas Ordowski

A300063 Heinz numbers of integer partitions of odd numbers.

Original entry on oeis.org

2, 5, 6, 8, 11, 14, 15, 17, 18, 20, 23, 24, 26, 31, 32, 33, 35, 38, 41, 42, 44, 45, 47, 50, 51, 54, 56, 58, 59, 60, 65, 67, 68, 69, 72, 73, 74, 77, 78, 80, 83, 86, 92, 93, 95, 96, 97, 98, 99, 103, 104, 105, 106, 109, 110, 114, 119, 122, 123, 124, 125, 126, 127

Offset: 1

Gus Wiseman, Feb 23 2018

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			15 is the Heinz number of (3,2), which has odd weight, so 15 belongs to the sequence.
Sequence of odd-weight partitions begins: (1) (3) (2,1) (1,1,1) (5) (4,1) (3,2) (7) (2,2,1) (3,1,1) (9) (2,1,1,1) (6,1).

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Complement of A300061.

Cf. A000041, A000720, A001222, A056239, A063834, A112798, A215366, A296150, A299757, A300056, A300060, A304662.

a:= proc(n) option remember; local k; for k from 1+
     `if`(n=1, 0, a(n-1)) while add(numtheory[pi]
      (i[1])*i[2], i=ifactors(k)[2])::even do od; k
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, May 22 2018

Select[Range[200],OddQ[Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]]&]

A299755 Triangle read by rows in which row n is the strict integer partition with FDH number n.

Original entry on oeis.org

1, 2, 3, 4, 2, 1, 5, 3, 1, 6, 4, 1, 7, 3, 2, 8, 5, 1, 4, 2, 9, 10, 6, 1, 11, 4, 3, 5, 2, 7, 1, 12, 3, 2, 1, 13, 8, 1, 6, 2, 5, 3, 14, 4, 2, 1, 15, 9, 1, 7, 2, 10, 1, 5, 4, 6, 3, 16, 11, 1, 8, 2, 4, 3, 1, 17, 5, 2, 1, 18, 7, 3, 6, 4, 12, 1, 19, 9, 2, 20, 13, 1

Offset: 1

Gus Wiseman, Feb 18 2018

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. Every positive integer n has a unique factorization of the form n = f(s_1)*...*f(s_k) where the s_i are strictly increasing positive integers. This determines a unique strict integer partition (s_k...s_1) whose FDH number is then defined to be n.

			Sequence of strict integer partitions begins: () (1) (2) (3) (4) (2,1) (5) (3,1) (6) (4,1) (7) (3,2) (8) (5,1) (4,2) (9) (10) (6,1) (11) (4,3) (5,2) (7,1) (12) (3,2,1) (13) (8,1) (6,2) (5,3) (14) (4,2,1) (15).

Row lengths are A064547.

Cf. A005117, A050376, A112798, A213925, A246867, A296150, A299756, A299757, A299759.

FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
nn=200;FDprimeList=Array[FDfactor,nn,1,Union];
FDrules=MapIndexed[(#1->#2[[1]])&,FDprimeList];
Join@@Table[Reverse[FDfactor[n]/.FDrules],{n,nn}]

A300272 Sorted list of Heinz numbers of odd partitions.

Original entry on oeis.org

2, 5, 8, 11, 17, 20, 23, 31, 32, 41, 44, 47, 50, 59, 67, 68, 73, 80, 83, 92, 97, 103, 109, 110, 124, 125, 127, 128, 137, 149, 157, 164, 167, 170, 176, 179, 188, 191, 197, 200, 211, 227, 230, 233, 236, 241, 242, 257, 268, 269, 272, 275, 277, 283, 292, 307, 310

Offset: 1

Gus Wiseman, Mar 01 2018

nonn

An odd partition is an integer partition of an odd number into an odd number of parts, all of which are odd.

Any product of three members of this sequence is also in the sequence.

			Sequence of odd partitions begins: (1), (3), (111), (5), (7), (311), (9), (11), (11111), (13), (511), (15), (331), (17), (19), (711), (21), (31111), (23), (911), (25), (27), (29), (531), (1111), (333), (31), (1111111).

Cf. A000009, A031368, A066208, A078408, A094457, A215366, A299757, A300063.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],OddQ[Total[primeMS[#]]]&&And@@OddQ/@primeMS[#]&]

A302243 Total weight of the n-th twice-odd-factored multiset partition.

Original entry on oeis.org

0, 1, 1, 2, 2, 1, 2, 2, 1, 3, 3, 2, 2, 3, 2, 1, 2, 3, 3, 3, 1, 2, 3, 2, 4, 2, 4, 2, 4, 1, 3, 4, 3, 1, 3, 3, 2, 3, 3, 2, 4, 1, 2, 3, 4, 4, 2, 4, 2, 3, 2, 3, 4, 3, 1, 4, 3, 3, 4, 3, 2, 2, 3, 1, 3, 5, 5, 4, 2, 2, 3, 3, 3, 5, 2, 4, 3, 2, 1, 5, 4, 2, 3, 2, 4, 5, 4, 4

Offset: 0

Gus Wiseman, Apr 03 2018

nonn

A multiset partition is a finite multiset of finite nonempty multisets of positive integers. The n-th twice-odd-factored multiset partition is constructed by factoring 2n + 1 into prime numbers and then factoring each prime index into prime numbers and taking their prime indices.

			Sequence of multiset partitions begins: (), ((1)), ((2)), ((11)), ((1)(1)), ((3)), ((12)), ((1)(2)), ((4)), ((111)), ((1)(11)), ((22)), ((2)(2)), ((1)(1)(1)), ((13)), ((5)), ((1)(3)), ((2)(11)), ((112)), ((1)(12)), ((6)).

Cf. A003963, A007716, A034691, A048673, A056239, A061775, A063834 A064216, A096443, A249620, A255397, A255906, A275024, A279789, A281113, A299757, A302242.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Sum[PrimeOmega[k],{k,primeMS[2n-1]}],{n,100}]

a(n) = A302242(2n + 1).

A246867 Triangle T(n,k) in which n-th row lists in increasing order all partitions lambda of n into distinct parts encoded as Product_{i:lambda} prime(i); n>=0, 1<=k<=A000009(n).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 14, 15, 13, 21, 22, 30, 17, 26, 33, 35, 42, 19, 34, 39, 55, 66, 70, 23, 38, 51, 65, 77, 78, 105, 110, 29, 46, 57, 85, 91, 102, 130, 154, 165, 210, 31, 58, 69, 95, 114, 119, 143, 170, 182, 195, 231, 330, 37, 62, 87, 115, 133, 138, 187

Offset: 0

Alois P. Heinz, Sep 05 2014

The concatenation of all rows (with offset 1) gives a permutation of the squarefree numbers A005117. The missing positive numbers are in A013929.

			The partitions of n=5 into distinct parts are {[5], [4,1], [3,2]}, encodings give {prime(5), prime(4)*prime(1), prime(3)*prime(2)} = {11, 7*2, 5*3} => row 5 = [11, 14, 15].
For n=0 the empty partition [] gives the empty product 1.
Triangle T(n,k) begins:
   1;
   2;
   3;
   5,  6;
   7, 10;
  11, 14, 15;
  13, 21, 22, 30;
  17, 26, 33, 35, 42;
  19, 34, 39, 55, 66,  70;
  23, 38, 51, 65, 77,  78, 105, 110;
  29, 46, 57, 85, 91, 102, 130, 154, 165, 210;
  ...
Corresponding triangle of strict integer partitions begins:
                  0
                 (1)
                 (2)
               (3) (21)
               (4) (31)
             (5) (41) (32)
          (6) (42) (51) (321)
        (7) (61) (52) (43) (421)
     (8) (71) (62) (53) (521) (431)
(9) (81) (72) (63) (54) (621) (432) (531). - _Gus Wiseman_, Feb 23 2018

Alois P. Heinz, Rows n = 0..42, flattened

Column k=1 gives: A008578(n+1).
Last elements of rows give: A246868.
Row sums give A147655.
Row lengths are: A000009.
Cf. A005117, A118462, A215366 (the same for all partitions), A258323, A299755, A299757, A299759.

b:= proc(n, i) option remember; `if`(n=0, [1], `if`(i<1, [], [seq(
      map(p->p*ithprime(i)^j, b(n-i*j, i-1))[], j=0..min(1, n/i))]))
    end:
T:= n-> sort(b(n$2))[]:
seq(T(n), n=0..14);

b[n_, i_] := b[n, i] = If[n==0, {1}, If[i<1, {}, Flatten[Table[Map[ #*Prime[i]^j&, b[n-i*j, i-1]], {j, 0, Min[1, n/i]}]]]]; T[n_] := Sort[b[n, n]]; Table[T[n], {n, 0, 14}] // Flatten (* Jean-François Alcover, Dec 18 2016, after Alois P. Heinz *)

A305829 Factor n into distinct Fermi-Dirac primes (A050376), normalize by replacing every instance of the k-th Fermi-Dirac prime with k, then multiply everything together.

Original entry on oeis.org

1, 1, 2, 3, 4, 2, 5, 3, 6, 4, 7, 6, 8, 5, 8, 9, 10, 6, 11, 12, 10, 7, 12, 6, 13, 8, 12, 15, 14, 8, 15, 9, 14, 10, 20, 18, 16, 11, 16, 12, 17, 10, 18, 21, 24, 12, 19, 18, 20, 13, 20, 24, 21, 12, 28, 15, 22, 14, 22, 24, 23, 15, 30, 27, 32, 14, 24, 30, 24, 20, 25

Offset: 1

Gus Wiseman, Jun 10 2018

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. Every positive integer n has a unique factorization of the form n = f(s_1)*...*f(s_k) where the s_i are strictly increasing positive integers. Then a(n) = s_1 * ... * s_k.

Multiplicative because for coprime m and n the Fermi-Dirac factorizations of m and n are disjoint and their union is the Fermi-Dirac factorization of m * n. - Andrew Howroyd, Aug 02 2018

Andrew Howroyd, Table of n, a(n) for n = 1..1000
Gus Wiseman, Tree of x -> a(x) for n = 1...75

Cf. A003963, A050376, A064547, A213925, A279065, A279614, A299755, A299756, A299757, A305830, A305831, A305832.

nn=100;
FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
FDprimeList=Array[FDfactor,nn,1,Union];FDrules=MapIndexed[(#1->#2[[1]])&,FDprimeList];
Table[Times@@(FDfactor[n]/.FDrules),{n,nn}]

\\ here isfd is membership test for A050376.
isfd(n)={my(e=isprimepower(n)); e && e == 1<Andrew Howroyd, Aug 02 2018

A318995 Totally additive with a(prime(n)) = n - 1.

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 3, 0, 2, 2, 4, 1, 5, 3, 3, 0, 6, 2, 7, 2, 4, 4, 8, 1, 4, 5, 3, 3, 9, 3, 10, 0, 5, 6, 5, 2, 11, 7, 6, 2, 12, 4, 13, 4, 4, 8, 14, 1, 6, 4, 7, 5, 15, 3, 6, 3, 8, 9, 16, 3, 17, 10, 5, 0, 7, 5, 18, 6, 9, 5, 19, 2, 20, 11, 5, 7, 7, 6, 21, 2, 4, 12

Offset: 1

Gus Wiseman, Sep 07 2018

nonn

Andrew Howroyd, Table of n, a(n) for n = 1..10000

Cf. A000040, A001222, A003963, A056239, A275024, A299757, A302242, A302243, A318994.

a:= n-> add((-1+numtheory[pi](i[1]))*i[2], i=ifactors(n)[2]):
seq(a(n), n=1..100);  # Alois P. Heinz, Sep 07 2018

Table[Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>(PrimePi[p]-1)*k]//Total,{n,200}]

a(n)={my(f=factor(n)); sum(i=1, #f~, my([p, e]=f[i, ]); (primepi(p)-1)*e)} \\ Andrew Howroyd, Sep 07 2018

cp's OEIS Frontend

A302242 Total weight of the n-th multiset multisystem. Totally additive with a(prime(n)) = Omega(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A300061 Heinz numbers of integer partitions of even numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A052330 Let S_k denote the first 2^k terms of this sequence and let b_k be the smallest positive integer that is not in S_k; then the numbers b_k*S_k are the next 2^k terms.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A300063 Heinz numbers of integer partitions of odd numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A299755 Triangle read by rows in which row n is the strict integer partition with FDH number n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A300272 Sorted list of Heinz numbers of odd partitions.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A302243 Total weight of the n-th twice-odd-factored multiset partition.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A246867 Triangle T(n,k) in which n-th row lists in increasing order all partitions lambda of n into distinct parts encoded as Product_{i:lambda} prime(i); n>=0, 1<=k<=A000009(n).

Views

Author

Keywords