A006171 -id:A006171 - cp's OEIS frontend

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

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

A107742 G.f.: Product_{j>=1} Product_{i>=1} (1 + x^(i*j)).

Original entry on oeis.org

1, 1, 2, 4, 6, 10, 17, 25, 38, 59, 86, 125, 184, 260, 369, 524, 726, 1005, 1391, 1894, 2576, 3493, 4687, 6272, 8373, 11090, 14647, 19294, 25265, 32991, 42974, 55705, 72025, 92895, 119349, 152965, 195592, 249280, 316991, 402215, 508932, 642598, 809739, 1017850, 1276959, 1599015, 1997943, 2491874, 3102477, 3855165, 4782408, 5922954

Offset: 0

Vladeta Jovovic, Jun 11 2005

From Gus Wiseman, Sep 13 2022: (Start)
Also the number of multiset partitions of integer partitions of n into intervals, where an interval is a set of positive integers with all differences of adjacent elements equal to 1. For example, the a(1) = 1 through a(4) = 6 multiset partitions are:
  {{1}}  {{2}}      {{3}}          {{4}}
         {{1},{1}}  {{1,2}}        {{1},{3}}
                    {{1},{2}}      {{2},{2}}
                    {{1},{1},{1}}  {{1},{1,2}}
                                   {{1},{1},{2}}
                                   {{1},{1},{1},{1}}
Intervals are counted by A001227, ranked by A073485.
The initial version is A007294.
The strict version is A327731.
The version for gapless multisets instead of intervals is A356941.
The case of strict partitions is A356957.
Also the number of multiset partitions of integer partitions of n into distinct constant blocks. For example, the a(1) = 1 through a(4) = 6 multiset partitions are:
  {{1}}  {{2}}    {{3}}        {{4}}
         {{1,1}}  {{1,1,1}}    {{2,2}}
                  {{1},{2}}    {{1},{3}}
                  {{1},{1,1}}  {{1,1,1,1}}
                               {{2},{1,1}}
                               {{1},{1,1,1}}
Constant multisets are counted by A000005, ranked by A000961.
The non-strict version is A006171.
The unlabeled version is A089259.
The non-constant block version is A261049.
The version for twice-partitions is A279786, factorizations A296131.
Also the number of multiset partitions of integer partitions of n into constant blocks of odd length. For example, a(1) = 1 through a(4) = 6 multiset partitions are:
  {{1}}  {{2}}      {{3}}          {{4}}
         {{1},{1}}  {{1,1,1}}      {{1},{3}}
                    {{1},{2}}      {{2},{2}}
                    {{1},{1},{1}}  {{1},{1,1,1}}
                                   {{1},{1},{2}}
                                   {{1},{1},{1},{1}}
The strict version is A327731 (also).
(End)

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
Lida Ahmadi, Ricardo Gómez Aíza, and Mark Daniel Ward, A unified treatment of families of partition functions, La Matematica (2024). Preprint available as arXiv:2303.02240 [math.CO], 2023.
N. J. A. Sloane, Transforms

Cf. A006171, A109386, A219554, A280473, A280486, A288007.
Product_{k>=1} (1 + x^k)^sigma_m(k): this sequence (m=0), A192065 (m=1), A288414 (m=2), A288415 (m=3), A301548 (m=4), A301549 (m=5), A301550 (m=6), A301551 (m=7), A301552 (m=8).
A000041 counts integer partitions, strict A000009.
A000110 counts set partitions.
A072233 counts partitions by sum and length.
Cf. A001055, A001970, A011782, A055887, A061260, A270995, A279784, A294617, A304969.

nmax = 50; CoefficientList[Series[Product[(1+x^(i*j)), {i, 1, nmax}, {j, 1, nmax/i}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 04 2017 *)
nmax = 50; CoefficientList[Series[Product[(1+x^k)^DivisorSigma[0, k], {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Mar 23 2018 *)
nmax = 50; s = 1 + x; Do[s *= Sum[Binomial[DivisorSigma[0, k], j]*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]];, {k, 2, nmax}]; Take[CoefficientList[s, x], nmax + 1] (* Vaclav Kotesovec, Aug 28 2018 *)
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
chQ[y_]:=Length[y]<=1||Union[Differences[y]]=={1};
Table[Length[Select[Join@@mps/@IntegerPartitions[n],And@@chQ/@#&]],{n,0,5}] (* Gus Wiseman, Sep 13 2022 *)

a(n)=polcoeff(prod(k=1,n,prod(j=1,n\k,1+x^(j*k)+x*O(x^n))),n) /* Paul D. Hanna */

N=66;  x='x+O('x^N); gf=1/prod(j=0,N, eta(x^(2*j+1))); gf=prod(j=1,N,(1+x^j)^numdiv(j)); Vec(gf) /* Joerg Arndt, May 03 2008 */

{a(n)=if(n==0,1,polcoeff(exp(sum(m=1,n,sigma(m)*x^m/(1-x^(2*m)+x*O(x^n))/m)),n))} /* Paul D. Hanna, Mar 28 2009 */

Euler transform of A001227.
Weigh transform of A000005.
G.f. satisfies: log(A(x)) = Sum_{n>=1} A109386(n)/n*x^n, where A109386(n) = Sum_{d|n} d*Sum_{m|d} (m mod 2). - Paul D. Hanna, Jun 26 2005
G.f.: A(x) = exp( Sum_{n>=1} sigma(n)*x^n/(1-x^(2n)) /n ). - Paul D. Hanna, Mar 28 2009
G.f.: Product_{n>=1} Q(x^n) where Q(x) is the g.f. of A000009. - Joerg Arndt, Feb 27 2014
a(0) = 1, a(n) = (1/n)*Sum_{k=1..n} A109386(k)*a(n-k) for n > 0. - Seiichi Manyama, Jun 04 2017
Conjecture: log(a(n)) ~ Pi*sqrt(n*log(n)/6). - Vaclav Kotesovec, Aug 29 2018

More terms from Paul D. Hanna, Jun 26 2005

A061256 Euler transform of sigma(n), cf. A000203.

Original entry on oeis.org

1, 1, 4, 8, 21, 39, 92, 170, 360, 667, 1316, 2393, 4541, 8100, 14824, 26071, 46422, 80314, 139978, 238641, 408201, 686799, 1156062, 1920992, 3189144, 5238848, 8589850, 13963467, 22641585, 36447544, 58507590, 93334008, 148449417, 234829969, 370345918

Offset: 0

Vladeta Jovovic, Apr 21 2001

This is also the number of ordered triples of permutations f, g, h in Symm(n) which all commute, divided by n!. This was conjectured by Franklin T. Adams-Watters, Jan 16 2006, and proved by J. R. Britnell in 2012.
According to a message on a blog page by "Allan" (see Secret Blogging Seminar link) it appears that a(n) = number of conjugacy classes of commutative ordered pairs in Symm(n).
John McKay (email to N. J. A. Sloane, Apr 23 2013) observes that A061256 and A006908 coincide for a surprising number of terms, and asks for an explanation. - N. J. A. Sloane, May 19 2013

			1 + x + 4*x^2 + 8*x^3 + 21*x^4 + 39*x^5 + 92*x^6 + 170*x^7 + 360*x^8 + ...

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe)
Lida Ahmadi, Ricardo Gómez Aíza, and Mark Daniel Ward, A unified treatment of families of partition functions, La Matematica (2024). Preprint available as arXiv:2303.02240 [math.CO], 2023.
J. R. Britnell, A formal identity involving commuting triples of permutations, arXiv:1203.5079 [math.CO], 2012.
J. R. Britnell, A formal identity involving commuting triples of permutations, Preprint 2012. - _N. J. A. Sloane_, Jun 13 2012
J. R. Britnell, A formal identity involving commuting triples of permutations, Journal of Combinatorial Theory, Series A, Volume 120, Issue 4, May 2013.
E. Marberg, How to compute the Frobenius-Schur indicator of a unipotent character of a finite Coxeter system, arXiv preprint arXiv:1202.1311 [math.RT], 2012. - _N. J. A. Sloane_, Jun 10 2012
Secret Blogging Seminar, A peculiar numerical coincidence.
N. J. A. Sloane, Transforms
Tad White, Counting Free Abelian Actions, arXiv:1304.2830 [math.CO], 2013.

Product_{k>=1} 1/(1 - x^k)^sigma_m(k): A006171 (m=0), this sequence (m=1), A275585 (m=2), A288391 (m=3), A301542 (m=4), A301543 (m=5), A301544 (m=6), A301545 (m=7), A301546 (m=8), A301547 (m=9).

Cf. A000203, A001001, A001970, A053529, A061255, A061257, A006908, A192065.

with(numtheory):
a:= proc(n) option remember; `if`(n=0, 1, add(add(
      d*sigma(d), d=divisors(j)) *a(n-j), j=1..n)/n)
    end:
seq(a(n), n=0..40);  # Alois P. Heinz, Jun 08 2017

nn = 30; b = Table[DivisorSigma[1, n], {n, nn}]; CoefficientList[Series[Product[1/(1 - x^m)^b[[m]], {m, nn}], {x, 0, nn}], x] (* T. D. Noe, Jun 18 2012 *)
nmax = 40; CoefficientList[Series[Product[1/QPochhammer[x^k]^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Nov 29 2015 *)

N=66; x='x+O('x^N); gf=1/prod(j=1,N, eta(x^j)^j); Vec(gf) /* Joerg Arndt, May 03 2008 */

{a(n)=if(n==0,1,polcoeff(exp(sum(m=1,n,sigma(m)*x^m/(1-x^m+x*O(x^n))^2/m)),n))} /* Paul D. Hanna, Mar 28 2009 */

a(n) = A072169(n) / n!.
G.f.: Product_{k=1..infinity} (1 - x^k)^(-sigma(k)). a(n)=1/n*Sum_{k=1..n} a(n-k)*b(k), n>1, a(0)=1, b(k)=Sum_{d|k} d*sigma(d), cf. A001001.
G.f.: exp( Sum_{n>=1} sigma(n)*x^n/(1-x^n)^2 /n ). [Paul D. Hanna, Mar 28 2009]
G.f.: exp( Sum_{n>=1} sigma_2(n)*x^n/(1-x^n)/n ). [Vladeta Jovovic, Mar 28 2009]
G.f.: prod(n>=1, E(x^n)^n ) where E(x) = prod(k>=1, 1-x^k). [Joerg Arndt, Apr 12 2013]
a(n) ~ exp((3*Pi)^(2/3) * Zeta(3)^(1/3) * n^(2/3)/2 - Pi^(4/3) * n^(1/3) / (4 * 3^(2/3) * Zeta(3)^(1/3)) - 1/24 - Pi^2/(288*Zeta(3))) * A^(1/2) * Zeta(3)^(11/72) / (2^(11/24) * 3^(47/72) * Pi^(11/72) * n^(47/72)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Mar 23 2018

Entry revised by N. J. A. Sloane, Jun 13 2012

A381636 Numbers whose prime indices cannot be partitioned into constant blocks with distinct sums.

Original entry on oeis.org

12, 60, 63, 84, 120, 126, 132, 156, 204, 228, 252, 276, 300, 315, 325, 348, 372, 420, 444, 492, 504, 516, 560, 564, 588, 630, 636, 650, 660, 693, 708, 720, 732, 780, 804, 819, 840, 852, 876, 924, 931, 948, 975, 996, 1008, 1020, 1068, 1071, 1092, 1140, 1164

Offset: 1

Gus Wiseman, Mar 10 2025

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Also numbers that cannot be written as a product of prime powers > 1 with distinct sums of prime indices (A056239).
Contains no squarefree numbers.
Conjecture: These are the zeros of A382876.

			The prime indices of 300 are {1,1,2,3,3}, with partitions into constant blocks:
  {{2},{1,1},{3,3}}
  {{1},{1},{2},{3,3}}
  {{2},{3},{3},{1,1}}
  {{1},{1},{2},{3},{3}}
but none of these has distinct block-sums, so 300 is in the sequence.
The terms together with their prime indices begin:
   12: {1,1,2}
   60: {1,1,2,3}
   63: {2,2,4}
   84: {1,1,2,4}
  120: {1,1,1,2,3}
  126: {1,2,2,4}
  132: {1,1,2,5}
  156: {1,1,2,6}
  204: {1,1,2,7}
  228: {1,1,2,8}
  252: {1,1,2,2,4}
  276: {1,1,2,9}
  300: {1,1,2,3,3}

More on multiset partitions into constant blocks: A006171, A279784, A295935.
These are the positions of 0 in A381635, after taking block-sums A381716.
Partitions of this type are counted by A381717.
For strict instead of constant blocks we have A381806, zeros of A381633.
For equal instead of distinct block-sums we have A381871.
A000688 counts multiset partitions into constant, see A381455 (upper), A381453 (lower).
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower).
A050361 counts multiset partitions into distinct constant blocks, after sums A381715.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
Cf. A000720, A001222, A005117, A050320, A059404, A213242, A293243, A299202, A300385, A381078, A381454, A381634.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
pfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[pfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],PrimePowerQ]}]];
Select[Range[100],Select[pfacs[#],UnsameQ@@hwt/@#&]=={}&]

A060640 If n = Product p_i^e_i then a(n) = Product (1 + 2p_i + 3p_i^2 + ... + (e_i+1)*p_i^e_i).

Original entry on oeis.org

1, 5, 7, 17, 11, 35, 15, 49, 34, 55, 23, 119, 27, 75, 77, 129, 35, 170, 39, 187, 105, 115, 47, 343, 86, 135, 142, 255, 59, 385, 63, 321, 161, 175, 165, 578, 75, 195, 189, 539, 83, 525, 87, 391, 374, 235, 95, 903, 162, 430, 245, 459, 107, 710, 253, 735, 273, 295, 119

Offset: 1

N. J. A. Sloane, Apr 17 2001

Equals row sums of triangle A143313. - Gary W. Adamson, Aug 06 2008
Equals row sums of triangle A127099. - Gary W. Adamson, Jul 27 2008
Sum of the divisors d2 from the ordered pairs of divisors of n, (d1,d2) with d1<=d2, such that d1|d2. - Wesley Ivan Hurt, Mar 22 2022

			a(4) = a(2^2) = 1 + (2)*(2) + (3)*(2^2) = 17;
a(6) = a(2)*a(3) = (1 + (2)*(2))*(1+(2)*(3)) = (5)*(7) = 35.
a(6) = tau(1) + 2*tau(2) + 3*tau(3) + 6*tau(6) = 1 + 2*2 + 3*2 + 6*4 = 35.

D. M. Burton, Elementary Number Theory, Allyn and Bacon Inc., Boston, MA, 1976, p. 120.

Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)

Cf. A000005, A000203, A001001, A006171, A038040 (Mobius transform), A049060, A057660, A057723, A327960 (Dirichlet inverse).

Cf. also triangles A027750, A127099, A143313.

a060640 n = sum [d * a000005 d | d <- a027750_row n]
-- Reinhard Zumkeller, Feb 29 2012

A060640 := proc(n) local ans, i, j; ans := 1: for i from 1 to nops(ifactors(n)[2]) do ans := ans*(1+sum((j+1)*ifactors(n)[2][i][1]^j,j=1..ifactors(n)[2][i][2])): od: RETURN(ans) end:

a[n_] := Total[#*DivisorSigma[1, n/#] & /@ Divisors[n]];
a /@ Range[59] (* Jean-François Alcover, May 19 2011, after Vladeta Jovovic *)
f[p_, e_] := ((e + 1)*p^(e + 2) - (e + 2)*p^(e + 1) + 1)/(p - 1)^2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Apr 10 2022 *)

j=[]; for(n=1,200,j=concat(j,sumdiv(n,d,n/d*sigma(d)))); j

a(n)=if(n<1,0,direuler(p=2,n,1/(1-X)/(1-p*X)^2)[n]) /* Ralf Stephan */

N=66; default(seriesprecision,N); x=z+O(z^(N+1))
c=sum(j=1,N,j*x^j); t=1/prod(j=1,N, eta(x^(j)));
t=log(t);t=serconvol(t,c);
Vec(t) /* Joerg Arndt, May 03 2008 */

{ for (n=1, 1000, write("b060640.txt", n, " ", direuler(p=2, n, 1/(1 - X)/(1 - p*X)^2)[n]); ) } /* Harry J. Smith, Jul 08 2009 */

def A060640(n) :
    sigma = sloane.A000203
    return add(sigma(k)*(n/k) for k in divisors(n))
[A060640(i) for i in (1..59)] # Peter Luschny, Sep 15 2012

a(n) = Sum_{d|n} d*tau(d), where tau(d) is the number of divisors of d, cf. A000005. a(n) = Sum_{d|n} d*sigma(n/d), where sigma(n)=sum of divisors of n, cf. A000203. - Vladeta Jovovic, Apr 23 2001
Multiplicative with a(p^e) = ((e+1)*p^{e+2} - (e+2)*p^{e+1} + 1) / (p-1)^2. Dirichlet g.f.: zeta(s)*zeta(s-1)^2. - Franklin T. Adams-Watters, Aug 03 2006
L.g.f.: Sum(A060640(n)*x^n/n) = -log( Product_{j>=1} P(x^j) ) where P(x) = Product_{k>=1} (1-x^k). - Joerg Arndt, May 03 2008
G.f.: Sum_{k>=1} k*tau(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Sep 06 2018
Sum_{k=1..n} a(k) ~ n^2/24 * ((4*gamma - 1)*Pi^2 + 2*Pi^2 * log(n) + 12*Zeta'(2)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Feb 01 2019

More terms from James Sellers, Vladeta Jovovic and Matthew Conroy, Apr 17 2001

A381635 Number of ways to partition the prime indices of n into constant blocks with distinct sums.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Mar 09 2025

nonn

First differs from A381716 at a(1728) = 5, A381716(1728) = 4.

Also the number of factorizations on n into prime powers > 1 with distinct sums of prime indices (A056239).

			The a(432) = 3 multiset partitions:
  {{2,2,2},{1,1,1,1}}
  {{1},{1,1,1},{2,2,2}}
  {{1},{2},{2,2},{1,1,1}}
Note {{2},{2,2},{1,1,1,1}} is not included, as it does not have distinct block-sums.

Without distinct block-sums we have A000688, after sums A381455 (upper), A381453 (lower).
For distinct blocks instead of sums we have A050361, after sums A381715.
For strict instead of constant we have A381633 (zeros A381806), after sums A381634.
Positions of 0 are A381636.
Taking block-sums (and sorting) gives A381716.
Other multiset partitions of prime indices:
More on multiset partitions into constant blocks: A006171, A279784, A295935.
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower).
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A001222, A002846, A005117, A213242, A213385, A293511, A299202, A300385, A317142, A381870.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
pfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[pfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],PrimePowerQ]}]];
Table[Length[Select[pfacs[n],UnsameQ@@hwt/@#&]],{n,100}]

A061255 Euler transform of Euler totient function phi(n), cf. A000010.

Original entry on oeis.org

1, 1, 2, 4, 7, 13, 21, 37, 60, 98, 157, 251, 392, 612, 943, 1439, 2187, 3293, 4930, 7330, 10839, 15935, 23315, 33933, 49170, 70914, 101861, 145713, 207638, 294796, 417061, 588019, 826351, 1157651, 1616849, 2251623, 3126775, 4330271, 5981190

Offset: 0

Vladeta Jovovic, Apr 21 2001

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe)
Alessandro De Luca, Gabriele Fici, and Andrea Frosini, Digital Convexity and Combinatorics on Words, arXiv:2502.19926 [math.CO], 2025. See p. 11.
N. J. A. Sloane, Transforms

Cf. A000010, A057660, A006171, A001970, A061256, A061257, A299069, A318975.

nn = 20; b = Table[EulerPhi[n], {n, nn}]; CoefficientList[Series[Product[1/(1 - x^m)^b[[m]], {m, nn}], {x, 0, nn}], x] (* T. D. Noe, Jun 19 2012 *)

G.f.: Product_{k>=1} (1 - x^k)^(-phi(k)).
a(n) = 1/n*Sum_{k=1..n} a(n-k)*b(k), n>1, a(0)=1, b(k) = Sum_{d|k} d*phi(d), cf. A057660.
Logarithmic derivative yields A057660 (equivalent to above formula). - Paul D. Hanna, Sep 05 2012
a(n) ~ exp(3^(4/3) * Zeta(3)^(1/3) * n^(2/3) / (2^(1/3) * Pi^(2/3)) - 1/6) * A^2 * Zeta(3)^(1/9) / (2^(4/9) * 3^(7/18) * Pi^(8/9) * n^(11/18)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Mar 23 2018
G.f.: exp(Sum_{k>=1} (sigma_2(k^2)/sigma_1(k^2)) * x^k/k). - Ilya Gutkovskiy, Apr 22 2019

A381716 Number of multisets that can be obtained by taking the sum of each block of a multiset partition of the prime indices of n into constant blocks with distinct sums.

Original entry on oeis.org

Offset: 1

Gus Wiseman, Mar 10 2025

nonn

First differs from A381635 at a(1728) = 4, A381635(1728) = 5.

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The prime indices of 1728 are {1,1,1,1,1,1,2,2,2}, with multiset partitions into constant multisets with distinct sums:
  {{1,1,1,1,1,1},{2,2},{2}}
  {{1,1,1,1,1},{1},{2,2,2}}
  {{1,1,1,1,1},{1},{2,2},{2}}
  {{1,1,1,1},{1,1},{2,2,2}}
  {{1,1,1},{1,1},{1},{2,2,2}}
with block-sums: {1,5,6}, {2,4,6}, {1,2,3,6}, {1,2,4,5}, so a(1728) = 4.

Without distinct sums we have A000688, after sums A381455 (upper), A381453 (lower).
More on multiset partitions into constant blocks: A006171, A279784, A295935.
For strict instead of constant we have A381633, before sums A381634.
Before taking sums we had A381635.
Positions of 0 are A381636.
For distinct blocks instead of sums we have A381715.
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower).
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
Cf. A000720, A001222, A002846, A005117, A050361, A213242, A265947, A293511, A299202, A300385, A362421.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
mce[y_]:=Table[ConstantArray[y[[1]],#]&/@ptn,{ptn,IntegerPartitions[Length[y]]}];
Table[Length[Union[Sort[Total/@#]&/@Select[Join@@@Tuples[mce/@Split[prix[n]]],UnsameQ@@Total/@#&]]],{n,100}]

A381453 Number of multisets that can be obtained by choosing a constant integer partition of each prime index of n and taking the multiset union.

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 3, 1, 3, 2, 2, 2, 4, 3, 4, 1, 2, 3, 4, 2, 6, 2, 3, 2, 3, 4, 4, 3, 4, 4, 2, 1, 4, 2, 6, 3, 6, 4, 8, 2, 2, 6, 4, 2, 6, 3, 4, 2, 6, 3, 4, 4, 5, 4, 4, 3, 8, 4, 2, 4, 6, 2, 8, 1, 8, 4, 2, 2, 6, 6, 6, 3, 4, 6, 6, 4, 6, 8, 4, 2, 5, 2, 2, 6, 4, 4, 8

Offset: 1

Gus Wiseman, Mar 08 2025

nonn

First differs from A355733 and A355735 at a(21) = 6, A355733(21) = A355735(21) = 5.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Multisets of constant multisets are generally not transitive. For example, we have arrows: {{1,1},{2}}: {1,1,2} -> {2,2} and {{2,2}}: {2,2} -> {4}, but there is no multiset of constant multisets {1,1,2} -> {4}.

			The a(21) = 6 multisets are: {2,4}, {1,1,4}, {2,2,2}, {1,1,2,2}, {2,1,1,1,1}, {1,1,1,1,1,1}.
The a(n) partitions for n = 1, 3, 7, 13, 53, 21 (G = 16):
  ()  (2)   (4)     (6)       (G)                 (42)
      (11)  (22)    (33)      (88)                (411)
            (1111)  (222)     (4444)              (222)
                    (111111)  (22222222)          (2211)
                              (1111111111111111)  (21111)
                                                  (111111)

Robert Price, Table of n, a(n) for n = 1..300

Positions of 1 are A000079.
The strict case is A008966.
Before sorting we had A355731.
Choosing divisors instead of constant multisets gives A355733.
The upper version is A381455, before taking sums A000688.
Multiset partitions of prime indices:
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For set multipartitions (A050320) see A381078 (upper), A381454 (lower).
- For set systems (A050326, zeros A293243) see A381441 (upper).
- For sets of constant multisets (A050361) see A381715.
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For set systems with distinct sums (A381633, zeros A381806) see A381634.
- For sets of constant multisets with distinct sums (A381635, zeros A381636) see A381716.
More on multiset partitions into constant blocks: A006171, A279784, A295935.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A000961, A001222, A002577, A018818, A213242, A213385, A213427, A275870, A299200, A300273, A300385.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Union[Sort/@Join@@@Tuples[Select[IntegerPartitions[#],SameQ@@#&]&/@prix[n]]]],{n,nn}]

a(A002110(n)) = A381807(n).

cp's OEIS Frontend

A288098 Convolution inverse of A006171.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

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

A107742 G.f.: Product_{j>=1} Product_{i>=1} (1 + x^(i*j)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

Extensions

A061256 Euler transform of sigma(n), cf. A000203.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A381636 Numbers whose prime indices cannot be partitioned into constant blocks with distinct sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A060640 If n = Product p_i^e_i then a(n) = Product (1 + 2*p_i + 3*p_i^2 + ... + (e_i+1)*p_i^e_i).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

A060640 If n = Product p_i^e_i then a(n) = Product (1 + 2p_i + 3p_i^2 + ... + (e_i+1)*p_i^e_i).