A275024 -id:A275024 - cp's OEIS frontend

A112798 Table where n-th row is factorization of n, with each prime p_i replaced by i.

1, 2, 1, 1, 3, 1, 2, 4, 1, 1, 1, 2, 2, 1, 3, 5, 1, 1, 2, 6, 1, 4, 2, 3, 1, 1, 1, 1, 7, 1, 2, 2, 8, 1, 1, 3, 2, 4, 1, 5, 9, 1, 1, 1, 2, 3, 3, 1, 6, 2, 2, 2, 1, 1, 4, 10, 1, 2, 3, 11, 1, 1, 1, 1, 1, 2, 5, 1, 7, 3, 4, 1, 1, 2, 2, 12, 1, 8, 2, 6, 1, 1, 1, 3, 13, 1, 2, 4, 14, 1, 1, 5, 2, 2, 3, 1, 9, 15, 1, 1, 1, 1

Offset: 2

Franklin T. Adams-Watters, Jan 22 2006

This is an enumeration of all partitions.
Technically this is an enumeration of all multisets (finite weakly increasing sequences of positive integers) rather than integer partitions. - Gus Wiseman, Dec 12 2016
A000040(a(n)) is a prime factor of A082288(n). - Reinhard Zumkeller, Feb 03 2008
Row n is the partition with Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1..r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. For a given n, the 2nd Maple program yields row n; for example, we obtain at once B(2436) = [1,1,2,4,10]. - Emeric Deutsch, Jun 04 2015
From Emeric Deutsch, May 05 2015: (Start)
Number of entries in row n is bigomega(n) (i.e., the number of prime factors of n, multiplicities included).
Product of entries in row n = A003963(n).
Row n contains the Matula numbers of the rooted trees obtained from the rooted tree with Matula number n by deleting the edges emanating from the root. Example: row 8 is 1,1,1; indeed the rooted tree with Matula number 8 is \|/ and deleting the edges emanating from the root we obtain three one-vertex trees, having Matula numbers 1, 1, 1. Example: row 7 is 4; indeed, the rooted tree with Matula number 7 is Y and deleting the edges emanating from the root we obtain the rooted tree V, having Matula number 4.
The Matula (or Matula-Goebel) number of a rooted tree can be defined in the following recursive manner: to the one-vertex tree there corresponds the number 1; to a tree T with root degree 1 there corresponds the t-th prime number, where t is the Matula-Goebel number of the tree obtained from T by deleting the edge emanating from the root; to a tree T with root degree m >= 2 there corresponds the product of the Matula-Goebel numbers of the m branches of T. (End)

			Row 20 is 1,1,3 because the prime factors of 20, namely 2,2,5 are the 1st, 1st, 3rd primes.
Table begins:
  1;
  2;
  1, 1;
  3;
  1, 2;
  4;
  1, 1, 1;
  ...
The sequence of all finite multisets of positive integers begins: (), (1), (2), (11), (3), (12), (4), (111), (22), (13), (5), (112), (6), (14), (23), (1111), (7), (122), (8), (113), (24), (15), (9), (1112), (33), (16), (222), (114). - _Gus Wiseman_, Dec 12 2016

Alois P. Heinz, Rows n = 2..3275, flattened
Emeric Deutsch, Rooted tree statistics from Matula numbers, arXiv:1111.4288 [math.CO], 2011.
Emeric Deutsch, Rooted tree statistics from Matula numbers, Discrete Appl. Math., 160, 2012, 2314-2322.
F. Goebel, On a 1-1-correspondence between rooted trees and natural numbers, J. Combin. Theory, B 29 (1980), 141-143.
I. Gutman and A. Ivic, On Matula numbers, Discrete Math., 150, 1996, 131-142.
I. Gutman and Yeong-Nan Yeh, Deducing properties of trees from their Matula numbers, Publ. Inst. Math., 53 (67), 1993, 17-22.
D. W. Matula, A natural rooted tree enumeration by prime factorization, SIAM Rev. 10 (1968) 273.
Index entries for sequences related to Matula-Goebel numbers
Index entries for sequences computed from indices in prime factorization

Row lengths are A001222. Cf. A000040, A027746, A000720, A036036.
Cf. A056239 (row sums).
Cf. A003963 (row products).
Cf. A049084, A082288, A215366, A241918, A265146, A275024.

a112798 n k = a112798_tabf !! (n-2) !! (n-1)
a112798_row n = a112798_tabf !! (n-2)
a112798_tabf = map (map a049084) $ tail a027746_tabf
-- Reinhard Zumkeller, Aug 04 2014

T:= n-> sort([seq(numtheory[pi](i[1])$i[2], i=ifactors(n)[2])])[]:
seq(T(n), n=2..50);  # Alois P. Heinz, Aug 09 2012
with(numtheory): B := proc (n) local nn, j, m: nn := op(2, ifactors(n)); for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: # Emeric Deutsch, Jun 04 2015. (This is equivalent to the first Maple program.)

PrimePi /@ Flatten[Table[#1, {#2}] & @@@ FactorInteger@ #] & /@ Range@ 60 // Flatten // Rest (* Michael De Vlieger, May 09 2015 *)

row(n)=my(v=List(),f=factor(n)); for(i=1,#f~,for(j=1,f[i,2], listput(v,primepi(f[i,1])))); Vec(v) \\ Charles R Greathouse IV, Nov 09 2021

T(n,k) = A000720(A027746(n,k)); A027746(n,k) = A000040(T(n,k)).

Also T(n,k) = A049084(A027746(n,k)). - Reinhard Zumkeller, Aug 04 2014

A001970 Functional determinants; partitions of partitions; Euler transform applied twice to all 1's sequence.

Original entry on oeis.org

1, 1, 3, 6, 14, 27, 58, 111, 223, 424, 817, 1527, 2870, 5279, 9710, 17622, 31877, 57100, 101887, 180406, 318106, 557453, 972796, 1688797, 2920123, 5026410, 8619551, 14722230, 25057499, 42494975, 71832114, 121024876, 203286806, 340435588, 568496753, 946695386

Offset: 0

N. J. A. Sloane

a(n) = number of partitions of n, when for each k there are p(k) different copies of part k. E.g., let the parts be 1, 2a, 2b, 3a, 3b, 3c, 4a, 4b, 4c, 4d, 4e, ... Then the a(4) = 14 partitions of 4 are: 4 = 4a = 4b = ... = 4e = 3a+1 = 3b+1 = 3c+1 = 2a+2a = 2a+2b = 2b+2b = 2a+1 = 2b+1 = 1+1+1+1.
Equivalently (Cayley), a(n) = number of 2-dimensional partitions of n. E.g., for n = 4 we have:
  4   31   3   22   2   211   21   2    2   1111   111   11    11   1
           1        2         1    11   1          1     11    1    1
                                        1                      1    1
                                                                    1
Also total number of different species of singularity for conjugate functions with n letters (Sylvester).
According to [Belmans], this sequence gives "[t]he number of Segre symbols for the intersection of two quadrics in a fixed dimension". - Eric M. Schmidt, Sep 02 2017
From Gus Wiseman, Jul 30 2022: (Start)
Also the number of non-isomorphic multiset partitions of weight n with all constant blocks. The strict case is A089259. For example, non-isomorphic representatives of the a(1) = 1 through a(3) = 6 multiset partitions are:
  {{1}}  {{1,1}}    {{1,1,1}}
         {{1},{1}}  {{1},{1,1}}
         {{1},{2}}  {{1},{2,2}}
                    {{1},{1},{1}}
                    {{1},{2},{2}}
                    {{1},{2},{3}}
A000688 counts factorizations into prime powers.
A007716 counts non-isomorphic multiset partitions by weight.
A279784 counts twice-partitions of type PPR, factorizations A295935.
Constant partitions are ranked by prime-powers: A000961, A023894, A054685, A246655, A355743.
(End)

			G.f. = 1 + x + 3*x^2 + 6*x^3 + 15*x^4 + 28*x^5 + 66*x^6 + 122*x^7 + ...
a(3) = 6 because we have (111) = (111) = (11)(1) = (1)(1)(1), (12) = (12) = (1)(2), (3) = (3).
The a(4)=14 multiset partitions whose total sum of parts is 4 are:
((4)),
((13)), ((1)(3)),
((22)), ((2)(2)),
((112)), ((1)(12)), ((2)(11)), ((1)(1)(2)),
((1111)), ((1)(111)), ((11)(11)), ((1)(1)(11)), ((1)(1)(1)(1)). - _Gus Wiseman_, Dec 19 2016

A. Cayley, Recherches sur les matrices dont les termes sont des fonctions linéaires d'une seule indéterminée, J. Reine angew. Math., 50 (1855), 313-317; Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 2, p. 219.
V. A. Liskovets, Counting rooted initially connected directed graphs. Vesci Akad. Nauk. BSSR, ser. fiz.-mat., No 5, 23-32 (1969), MR44 #3927.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
J. J. Sylvester, An Enumeration of the Contacts of Lines and Surfaces of the Second Order, Phil. Mag. 1 (1851), 119-140. Reprinted in Collected Papers, Vol. 1. See p. 239, where one finds a(n)-2, but with errors.
J. J. Sylvester, Note on the 'Enumeration of the Contacts of Lines and Surfaces of the Second Order', Phil. Mag., Vol. VII (1854), pp. 331-334. Reprinted in Collected Papers, Vol. 2, pp. 30-33.

Reinhard Zumkeller, Table of n, a(n) for n = 0..5000 (first 500 terms from T. D. Noe)
Pieter Belmans, Segre symbols, 2016.
Philip Boalch, Counting the fission trees and nonabelian Hodge graphs, arXiv:2410.23358 [math.AG], 2024. See pp. 10, 16.
P. J. Cameron, Some sequences of integers, Discrete Math., 75 (1989), 89-102; also in "Graph Theory and Combinatorics 1988", ed. B. Bollobas, Annals of Discrete Math., 43 (1989), 89-102.
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 148
R. Kaneiwa, An asymptotic formula for Cayley's double partition function p(2; n), Tokyo J. Math. 2, 137-158 (1979).
L. Kaylor and D. Offner, Counting matrices over a finite field with all eigenvalues in the field, Involve, a Journal of Mathematics, Vol. 7 (2014), No. 5, 627-645. [DOI]
M. Kozek, F. Luca, P. Pollack, and C. Pomerance, Harmonious pairs, 2014.
M. Kozek, F. Luca, P. Pollack, and C. Pomerance, Harmonious numbers, IJNT, to appear.
XiKun Li, JunLi Li, Bin Liu and CongFeng Qiao, The parametric symmetry and numbers of the entangled class of 2 × M × N system, Science China Physics, Mechanics & Astronomy, Volume 54, Number 8, 1471-1475, DOI: 10.1007/s11433-011-4395-9.
Jessie Pitsillides, Segre Characteristic Equivalence, arXiv:2506.12065 [math.GM], 2025.
Paul Pollack and Carl Pomerance, Some problems of Erdős on the sum-of-divisors function, For Richard Guy on his 99th birthday: May his sequence be unbounded, Trans. Amer. Math. Soc. Ser. B, Vol. 3 (2016), pp. 1-26; Errata.
N. J. A. Sloane, Transforms.
N. J. A. Sloane and Thomas Wieder, The Number of Hierarchical Orderings, arXiv:math/0307064 [math.CO], 2003; Order 21 (2004), 83-89.
J. J. Sylvester, The collected mathematical papers of James Joseph Sylvester, vol. 2, vol. 3, vol. 4.
Index entries for sequences related to rooted trees

Related to A001383 via generating function.
The multiplicative version (factorizations) is A050336.
The ordered version (sequences of partitions) is A055887.
Row-sums of A061260.
Main diagonal of A055885.
We have A271619(n) <= a(n) <= A063834(n).
Column k=3 of A290353.
The strict case is A316980.
Cf. A000041, A061259, A006171, A061255, A061256, A061257, A089292, A000219.
Cf. A089300.
Cf. A001055, A072233, A112798, A275024.

Following Vladeta Jovovic:
a001970 n = a001970_list !! (n-1)
a001970_list = 1 : f 1 [1] where
   f x ys = y : f (x + 1) (y : ys) where
            y = sum (zipWith (*) ys a061259_list) `div` x
-- Reinhard Zumkeller, Oct 31 2015

with(combstruct); SetSetSetU := [T, {T=Set(S), S=Set(U,card >= 1), U=Set(Z,card >=1)},unlabeled];
# second Maple program:
with(numtheory): with(combinat):
a:= proc(n) option remember; `if`(n=0, 1, add(add(d*
      numbpart(d), d=divisors(j))*a(n-j), j=1..n)/n)
    end:
seq(a(n), n=0..35);  # Alois P. Heinz, Dec 19 2016

m = 32; f[x_] = Product[1/(1-x^k)^PartitionsP[k], {k, 1, m}]; CoefficientList[ Series[f[x], {x, 0, m-1}], x] (* Jean-François Alcover, Jul 19 2011, after g.f. *)

{a(n) = if( n<0, 0, polcoeff( 1 / prod(k=1, n, 1 - numbpart(k) * x^k + x * O(x^n)), n))}; /* Michael Somos, Dec 20 2016 */

from sympy.core.cache import cacheit
from sympy import npartitions, divisors
@cacheit
def a(n): return 1 if n == 0 else sum([sum([d*npartitions(d) for d in divisors(j)])*a(n - j) for j in range(1, n + 1)]) / n
[a(n) for n in range(51)]  # Indranil Ghosh, Aug 19 2017, after Maple code
# (Sage) # uses[EulerTransform from A166861]
b = BinaryRecurrenceSequence(0, 1, 1)
a = EulerTransform(EulerTransform(b))
print([a(n) for n in range(36)]) # Peter Luschny, Nov 17 2022

G.f.: Product_{k >= 1} 1/(1-x^k)^p(k), where p(k) = number of partitions of k = A000041. [Cayley]
a(n) = (1/n)*Sum_{k = 1..n} a(n-k)*b(k), n > 1, a(0) = 1, b(k) = Sum_{d|k} d*numbpart(d), where numbpart(d) = number of partitions of d, cf. A061259. - Vladeta Jovovic, Apr 21 2001
Logarithmic derivative yields A061259 (equivalent to above formula from Vladeta Jovovic). - Paul D. Hanna, Sep 05 2012
a(n) = Sum_{k=1..A000041(n)} A001055(A215366(n,k)) = number of factorizations of Heinz numbers of integer partitions of n. - Gus Wiseman, Dec 19 2016
a(n) = |{m>=1 : n = Sum_{k=1..A001222(m)} A056239(A112798(m,k)+1)}| = number of normalized twice-prime-factored multiset partitions (see A275024) whose total sum of parts is n. - Gus Wiseman, Dec 19 2016

Additional comments from Valery A. Liskovets

Sylvester references from Barry Cipra, Oct 07 2003

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

Original entry on oeis.org

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

A302696 Numbers whose prime indices (with repetition) are pairwise coprime. Nonprime Heinz numbers of integer partitions with pairwise coprime parts.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 14, 15, 16, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 38, 40, 44, 46, 48, 51, 52, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 74, 76, 77, 80, 82, 85, 86, 88, 92, 93, 94, 95, 96, 102, 104, 106, 110, 112, 116, 118, 119, 120, 122, 123, 124, 128, 132

Offset: 1

Gus Wiseman, Apr 11 2018

nonn

A prime index of n is a number m such that prime(m) divides n. Two or more numbers are coprime if no pair has a common divisor other than 1. A single number is not considered coprime unless it is equal to 1.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Number 36 = prime(1)*prime(1)*prime(2)*prime(2) is not included in the sequence, because the pair of prime indices {2,2} is not coprime. - Gus Wiseman, Dec 06 2021

			Sequence of integer partitions with pairwise coprime parts begins: (), (1), (11), (21), (111), (31), (211), (41), (32), (1111), (311), (51), (2111), (61), (411), (321).
Missing from this list are: (2), (3), (4), (22), (5), (6), (7), (221), (8), (42), (9), (33), (222).

Cf. A000837, A000961, A001222, A005117, A007359, A051424, A275024, A289508, A289509, A298748, A302568, A302569, A302697, A302698, A327512, A327513.

filter:= proc(n) local F;
   F:= ifactors(n)[2];
   if nops(F)=1 then if F[1][1] = 2 then return true else return false fi fi;
   if ormap(t -> t[2]>1 and t[1] <> 2, F) then return false fi;
   F:= map(t -> numtheory:-pi(t[1]), F);
   ilcm(op(F))=convert(F,`*`)
end proc:
select(filter, [$1..200]); # Robert Israel, Sep 10 2020

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[200],Or[#===1,CoprimeQ@@primeMS[#]]&]

isA302696(n) = if(isprimepower(n),!(n%2), if(!issquarefree(n>>valuation(n,2)), 0, my(pis=apply(primepi,factor(n)[,1])); (lcm(pis)==factorback(pis)))); \\ Antti Karttunen, Dec 06 2021

Clarification (with repetition) added to the definition by Antti Karttunen, Dec 06 2021

A302569 Numbers that are either prime or whose prime indices are pairwise coprime. Heinz numbers of integer partitions with pairwise coprime parts.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 22, 23, 24, 26, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 40, 41, 43, 44, 46, 47, 48, 51, 52, 53, 55, 56, 58, 59, 60, 61, 62, 64, 66, 67, 68, 69, 70, 71, 73, 74, 76, 77, 79, 80, 82, 83, 85, 86, 88, 89

Offset: 1

Gus Wiseman, Apr 10 2018

nonn

A prime index of n is a number m such that prime(m) divides n.

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

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of multiset systems.
02: {{}}
03: {{1}}
04: {{},{}}
05: {{2}}
06: {{},{1}}
07: {{1,1}}
08: {{},{},{}}
10: {{},{2}}
11: {{3}}
12: {{},{},{1}}
13: {{1,2}}
14: {{},{1,1}}
15: {{1},{2}}
16: {{},{},{},{}}
17: {{4}}
19: {{1,1,1}}
20: {{},{},{2}}
22: {{},{3}}
23: {{2,2}}
24: {{},{},{},{1}}
26: {{},{1,2}}
28: {{},{},{1,1}}
29: {{1,3}}
30: {{},{1},{2}}
31: {{5}}
32: {{},{},{},{},{}}

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

Subsequence of A122132.

Cf. A000961, A001222, A005117, A007359, A007716, A051424, A056239, A076610, A101268, A275024, A302505, A302568.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[200],Or[PrimeQ[#],CoprimeQ@@primeMS[#]]&]

is(n)=if(n<9, return(n>1)); n>>=valuation(n,2); if(n<9, return(1)); my(f=factor(n)); if(vecmax(f[,2])>1, return(0)); if(#f~==1, return(1)); my(v=apply(primepi, f[,1]),P=vecprod(v)); for(i=1,#v, if(gcd(v[i],P/v[i])>1, return(0))); 1 \\ Charles R Greathouse IV, Nov 11 2021

A062457 a(n) = prime(n)^n.

Original entry on oeis.org

2, 9, 125, 2401, 161051, 4826809, 410338673, 16983563041, 1801152661463, 420707233300201, 25408476896404831, 6582952005840035281, 925103102315013629321, 73885357344138503765449, 12063348350820368238715343, 3876269050118516845397872321

Offset: 1

Labos Elemer, Jul 09 2001

Heinz numbers of square integer partitions, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 14 2018
Main diagonal of A182944. - Omar E. Pol, Sep 12 2018
Second diagonal of A319075. - Omar E. Pol, Sep 13 2018

T. D. Noe, Table of n, a(n) for n = 1..100

Cf. A000040, A000312, A000720, A000961, A001222, A051674, A062006, A076146, A093358, A201614, A275024, A302242, A302493, A302498, A319075.

[NthPrime(n)^n : n in [1..20]]; // Wesley Ivan Hurt, Jan 18 2016

A062457:=n->ithprime(n)^n: seq(A062457(n), n=1..20); # Wesley Ivan Hurt, Jan 18 2016

Table[Prime[n]^n, {n,20}] (* Harvey P. Dale, Jun 22 2011 *)

a(n) = prime(n)^n; \\ Harry J. Smith, Aug 07 2009

a(n) = A062006(n) - 1. - Wesley Ivan Hurt, Jan 18 2016
From Amiram Eldar, Nov 16 2020: (Start)
Sum_{n>=1} 1/a(n) = A093358.
Sum_{n>=1} (-1)^(n+1)/a(n) = A201614. (End)

A302478 Products of prime numbers of squarefree index.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 20, 22, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 36, 39, 40, 41, 43, 44, 45, 47, 48, 50, 51, 52, 54, 55, 58, 59, 60, 62, 64, 65, 66, 67, 68, 72, 73, 75, 78, 79, 80, 81, 82, 83, 85, 86, 87, 88, 90, 93, 94

Offset: 1

Gus Wiseman, Apr 08 2018

nonn

A prime index of n is a number m such that prime(m) divides n.

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of set multisystems.
01:  {}
02:  {{}}
03:  {{1}}
04:  {{},{}}
05:  {{2}}
06:  {{},{1}}
08:  {{},{},{}}
09:  {{1},{1}}
10:  {{},{2}}
11:  {{3}}
12:  {{},{},{1}}
13:  {{1,2}}
15:  {{1},{2}}
16:  {{},{},{},{}}
17:  {{4}}
18:  {{},{1},{1}}
20:  {{},{},{2}}
22:  {{},{3}}
24:  {{},{},{},{1}}
25:  {{2},{2}}
26:  {{},{1,2}}
27:  {{1},{1},{1}}
29:  {{1,3}}
30:  {{},{1},{2}}
31:  {{5}}
32:  {{},{},{},{},{}}

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

Cf. A000961, A001222, A003963, A005117, A007716, A050320, A056239, A063834, A076610, A270995, A275024, A281113, A296119, A301753, A302242, A302243, A302491.

Select[Range[100],Or[#===1,And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]

ok(n)={!#select(p->!issquarefree(primepi(p)), factor(n)[,1])} \\ Andrew Howroyd, Aug 26 2018

A303386 Number of aperiodic factorizations of n > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 7, 1, 2, 2, 4, 1, 5, 1, 6, 2, 2, 2, 7, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 1, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 3, 2, 1, 11, 2, 2, 2, 7, 1, 11, 2, 4, 2, 2, 2, 19, 1, 4, 4, 7, 1, 5, 1, 7, 5

Offset: 2

Gus Wiseman, Apr 23 2018

nonn

An aperiodic factorization of n is a finite multiset of positive integers greater than 1 whose product is n and whose multiplicities are relatively prime.

			The a(36) = 7 aperiodic factorizations are (2*2*9), (2*3*6), (2*18), (3*3*4), (3*12), (4*9), and (36). Missing from this list are (2*2*3*3) and (6*6).

Antti Karttunen, Table of n, a(n) for n = 2..65537

Cf. A000740, a(2^n) = A000837(n), A001055, A007716, A007916, A045778, A052409, A052410, A100953, A162247, A275024, A275870, A281113, A281116, A301700, A302242, A303431.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[n],GCD@@Length/@Split[#]===1&]],{n,2,100}]

A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
A303386(n) = if(1==n,n,my(r); sumdiv(A052409(n),d, ispower(n,d,&r); moebius(d)*A001055(r))); \\ Antti Karttunen, Sep 25 2018

a(n) = Sum_{d|A052409(n)} mu(d) * A001055(n^(1/d)), where mu = A008683.

More terms from Antti Karttunen, Sep 25 2018

A303431 Aperiodic tree numbers. Matula-Goebel numbers of aperiodic rooted trees.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 12, 13, 15, 18, 20, 22, 24, 26, 29, 30, 31, 33, 37, 39, 40, 41, 44, 45, 47, 48, 50, 52, 54, 55, 58, 60, 61, 62, 65, 66, 71, 72, 74, 75, 78, 79, 80, 82, 87, 88, 89, 90, 93, 94, 96, 99, 101, 104, 108, 109, 110, 111, 113, 116, 117, 120, 122

Offset: 1

Gus Wiseman, Apr 23 2018

nonn

A positive integer is an aperiodic tree number iff either it is equal to 1 or it belongs to A007916 (numbers that are not perfect powers, or numbers whose prime multiplicities are relatively prime) and all of its prime indices are also aperiodic tree numbers, where a prime index of n is a number m such that prime(m) divides n.

			Sequence of aperiodic rooted trees begins:
01 o
02 (o)
03 ((o))
05 (((o)))
06 (o(o))
10 (o((o)))
11 ((((o))))
12 (oo(o))
13 ((o(o)))
15 ((o)((o)))
18 (o(o)(o))
20 (oo((o)))
22 (o(((o))))
24 (ooo(o))
26 (o(o(o)))
29 ((o((o))))
30 (o(o)((o)))
31 (((((o)))))
33 ((o)(((o))))

Cf. A000081, A000740, A000837, A007097, A007916, A052409, A052410, A061775, A111299, A214577, A275024, A276625, A290760, A290822, A291442, A298533, A298536, A301700, A302242, A303386.

zapQ[1]:=True;zapQ[n_]:=And[GCD@@FactorInteger[n][[All,2]]===1,And@@zapQ/@PrimePi/@FactorInteger[n][[All,1]]];
Select[Range[100],zapQ]

A302494 Products of distinct primes of squarefree index.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 15, 17, 22, 26, 29, 30, 31, 33, 34, 39, 41, 43, 47, 51, 55, 58, 59, 62, 65, 66, 67, 73, 78, 79, 82, 83, 85, 86, 87, 93, 94, 101, 102, 109, 110, 113, 118, 123, 127, 129, 130, 134, 137, 139, 141, 143, 145, 146, 149, 155, 157, 158, 163

Offset: 1

Gus Wiseman, Apr 08 2018

nonn

A prime index of n is a number m such that prime(m) divides n.

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of set systems.
01: {}
02: {{}}
03: {{1}}
05: {{2}}
06: {{},{1}}
10: {{},{2}}
11: {{3}}
13: {{1,2}}
15: {{1},{2}}
17: {{4}}
22: {{},{3}}
26: {{},{1,2}}
29: {{1,3}}
30: {{},{1},{2}}
31: {{5}}
33: {{1},{3}}
34: {{},{4}}
39: {{1},{1,2}}

Cf. A000961, A001222, A003963, A005117, A007716, A056239, A270995, A275024, A276625, A277098, A279785, A281113, A296120, A301754, A302242, A302243.

Select[Range[100],Or[#===1,SquareFreeQ[#]&&And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]

is(n) = if(bigomega(n)!=omega(n), return(0), my(f=factor(n)[, 1]~); for(k=1, #f, if(!issquarefree(primepi(f[k])) && primepi(f[k])!=1, return(0)))); 1 \\ Felix Fröhlich, Apr 10 2018

cp's OEIS Frontend

A112798 Table where n-th row is factorization of n, with each prime p_i replaced by i.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A001970 Functional determinants; partitions of partitions; Euler transform applied twice to all 1's sequence.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

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

A302696 Numbers whose prime indices (with repetition) are pairwise coprime. Nonprime Heinz numbers of integer partitions with pairwise coprime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A302569 Numbers that are either prime or whose prime indices are pairwise coprime. Heinz numbers of integer partitions with pairwise coprime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A062457 a(n) = prime(n)^n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI