cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 53 results. Next

A316413 Heinz numbers of integer partitions whose length divides their sum.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 16, 17, 19, 21, 22, 23, 25, 27, 28, 29, 30, 31, 32, 34, 37, 39, 41, 43, 46, 47, 49, 53, 55, 57, 59, 61, 62, 64, 67, 68, 71, 73, 78, 79, 81, 82, 83, 84, 85, 87, 88, 89, 90, 91, 94, 97, 98, 99, 100, 101, 103, 105, 107, 109, 110
Offset: 1

Views

Author

Gus Wiseman, Jul 02 2018

Keywords

Comments

In other words, partitions whose average is an integer.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of partitions whose length divides their sum begins (1), (2), (11), (3), (4), (111), (22), (31), (5), (6), (1111), (7), (8), (42), (51), (9), (33), (222), (411).
		

Crossrefs

Programs

  • Maple
    isA326413 := proc(n)
        psigsu := A056239(n) ;
        psigle := numtheory[bigomega](n) ;
        if modp(psigsu,psigle) = 0 then
            true;
        else
            false;
        end if;
    end proc:
    n := 1:
    for i from 2 to 3000 do
        if isA326413(i) then
            printf("%d %d\n",n,i);
            n := n+1 ;
        end if;
    end do: # R. J. Mathar, Aug 09 2019
    # second Maple program:
    q:= n-> (l-> nops(l)>0 and irem(add(i, i=l), nops(l))=0)(map
            (i-> numtheory[pi](i[1])$i[2], ifactors(n)[2])):
    select(q, [$1..110])[];  # Alois P. Heinz, Nov 19 2021
  • Mathematica
    Select[Range[2,100],Divisible[Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]],PrimeOmega[#]]&]

A067538 Number of partitions of n in which the number of parts divides n.

Original entry on oeis.org

1, 2, 2, 4, 2, 8, 2, 11, 9, 14, 2, 46, 2, 24, 51, 66, 2, 126, 2, 202, 144, 69, 2, 632, 194, 116, 381, 756, 2, 1707, 2, 1417, 956, 316, 2043, 5295, 2, 511, 2293, 9151, 2, 10278, 2, 8409, 14671, 1280, 2, 36901, 8035, 21524, 11614, 25639, 2, 53138, 39810, 85004
Offset: 1

Views

Author

Naohiro Nomoto, Jan 27 2002

Keywords

Comments

Also sum of p(n,d) over the divisors d of n, where p(n,m) is the count of partitions of n in exactly m parts. - Wouter Meeussen, Jun 07 2009
From Gus Wiseman, Sep 24 2019: (Start)
Also the number of integer partitions of n whose maximum part divides n. The Heinz numbers of these partitions are given by A326836. For example, the a(1) = 1 through a(8) = 11 partitions are:
(1) (2) (3) (4) (5) (6) (7) (8)
(11) (111) (22) (11111) (33) (1111111) (44)
(211) (222) (422)
(1111) (321) (431)
(2211) (2222)
(3111) (4211)
(21111) (22211)
(111111) (41111)
(221111)
(2111111)
(11111111)
(End)

Examples

			a(3)=2 because 3 is a prime; a(4)=4 because the five partitions of 4 are {4}, {3, 1}, {2, 2}, {2, 1, 1}, {1, 1, 1, 1}, and the number of parts in each of them divides 4 except for {2, 1, 1}.
From _Gus Wiseman_, Sep 24 2019: (Start)
The a(1) = 1 through a(8) = 11 partitions whose length divides their sum are the following. The Heinz numbers of these partitions are given by A316413.
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (11111)  (33)      (1111111)  (44)
                    (31)             (42)                 (53)
                    (1111)           (51)                 (62)
                                     (222)                (71)
                                     (321)                (2222)
                                     (411)                (3221)
                                     (111111)             (3311)
                                                          (4211)
                                                          (5111)
                                                          (11111111)
(End)
		

Crossrefs

The strict case is A102627.
Partitions with integer geometric mean are A067539.

Programs

  • Mathematica
    Do[p = IntegerPartitions[n]; l = Length[p]; c = 0; k = 1; While[k < l + 1, If[ IntegerQ[ n/Length[ p[[k]] ]], c++ ]; k++ ]; Print[c], {n, 1, 57}, All]
    p[n_,k_]:=p[n,k]=p[n-1,k-1]+p[n-k,k];p[n_,k_]:=0/;k>n;p[n_,n_]:=1;p[n_,0]:=0
    Table[Plus @@ (p[n,# ]&/ @ Divisors[n]),{n,36}] (* Wouter Meeussen, Jun 07 2009 *)
    Table[Count[IntegerPartitions[n], q_ /; IntegerQ[Mean[q]]], {n, 50}]  (*Clark Kimberling, Apr 23 2019 *)
  • PARI
    a(n) = {my(nb = 0); forpart(p=n, if ((vecsum(Vec(p)) % #p) == 0, nb++);); nb;} \\ Michel Marcus, Jul 03 2018
    
  • Python
    # uses A008284_T
    from sympy import divisors
    def A067538(n): return sum(A008284_T(n,d) for d in divisors(n,generator=True)) # Chai Wah Wu, Sep 21 2023

Formula

a(p) = 2 for all primes p.

Extensions

Extended by Robert G. Wilson v, Oct 16 2002

A025065 Number of palindromic partitions of n.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 7, 7, 12, 12, 19, 19, 30, 30, 45, 45, 67, 67, 97, 97, 139, 139, 195, 195, 272, 272, 373, 373, 508, 508, 684, 684, 915, 915, 1212, 1212, 1597, 1597, 2087, 2087, 2714, 2714, 3506, 3506, 4508, 4508, 5763, 5763, 7338, 7338, 9296, 9296, 11732, 11732, 14742, 14742, 18460, 18460, 23025, 23025, 28629, 28629
Offset: 0

Views

Author

Keywords

Comments

That is, the number of partitions of n into parts which can be listed in palindromic order.
Alternatively, number of partitions of n into parts from the set {1,2,4,6,8,10,12,...}. - T. D. Noe, Aug 05 2005
Also, partial sums of A035363.
Also number of partitions of n with at most one part occurring an odd number of times. - Reinhard Zumkeller, Dec 18 2013
The first Mathematica program computes terms of A025065; the second computes the k palindromic partitions of user-chosen n. - Clark Kimberling, Jan 20 2014
a(n) is the number of partitions p of n+1 such that 2*max(p) > n+1. - Clark Kimberling, Apr 20 2014.
From Gus Wiseman, Nov 28 2018: (Start)
Also the number of integer partitions of n + 2 that are the vertex-degrees of some hypertree. For example, the a(6) = 7 partitions of 8 that are the vertex-degrees of some hypertree, together with a realizing hypertree are:
(41111): {{1,2},{1,3},{1,4},{1,5}}
(32111): {{1,2},{1,3},{1,4},{2,5}}
(22211): {{1,2},{1,3},{2,4},{3,5}}
(311111): {{1,2},{1,3},{1,4,5,6}}
(221111): {{1,2},{1,3},{2,4,5,6}}
(2111111): {{1,2},{1,3,4,5,6,7}}
(11111111): {{1,2,3,4,5,6,7,8}}
(End)
Conjecture: a(n) is the length of maximal initial segment of A308355(n-1) that is identical to row n of A128628, for n >= 2. - Clark Kimberling, May 24 2019
From Gus Wiseman, May 21 2021: (Start)
The Heinz numbers of palindromic partitions are given by A265640. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.
Also the number of integer partitions of n with a part greater than or equal to n/2. This is equivalent to Clark Kimberling's final comment above. The Heinz numbers of these partitions are given by A344414. For example, the a(1) = 1 through a(8) = 12 partitions are:
(1) (2) (3) (4) (5) (6) (7) (8)
(11) (21) (22) (32) (33) (43) (44)
(31) (41) (42) (52) (53)
(211) (311) (51) (61) (62)
(321) (421) (71)
(411) (511) (422)
(3111) (4111) (431)
(521)
(611)
(4211)
(5111)
(41111)
Also the number of integer partitions of n with at least n/2 parts. The Heinz numbers of these partitions are given by A344296. For example, the a(1) = 1 through a(8) = 12 partitions are:
(1) (2) (21) (22) (221) (222) (2221) (2222)
(11) (111) (31) (311) (321) (3211) (3221)
(211) (2111) (411) (4111) (3311)
(1111) (11111) (2211) (22111) (4211)
(3111) (31111) (5111)
(21111) (211111) (22211)
(111111) (1111111) (32111)
(41111)
(221111)
(311111)
(2111111)
(11111111)
(End)

Examples

			The partitions for the first few values of n are as follows:
n: partitions .......................... number
1: 1 ................................... 1
2: 2 11 ................................ 2
3: 3 111 ............................... 2
4: 4 22 121 1111 ....................... 4
5: 5 131 212 11111 ..................... 4
6: 6 141 33 222 1221 11211 111111 ...... 7
7: 7 151 313 11311 232 21112 1111111 ... 7
From _Reinhard Zumkeller_, Jan 23 2010: (Start)
Partitions into 1,2,4,6,... for the first values of n:
1: 1 ....................................... 1
2: 2 11 .................................... 2
3: 21 111 .................................. 2
4: 4 22 211 1111 ........................... 4
5: 41 221 2111 11111 ....................... 4
6: 6 42 4211 222 2211 21111 111111.......... 7
7: 61 421 42111 2221 22111 211111 1111111 .. 7. (End)
		

Crossrefs

Cf. A172033, A004277. - Reinhard Zumkeller, Jan 23 2010
The bisections are both A000070.
The ordered version (palindromic compositions) is A016116.
The complement is counted by A233771 and A210249.
The case of palindromic prime signature is A242414.
Palindromic partitions are ranked by A265640, with complement A229153.
The case of palindromic plane trees is A319436.
The multiplicative version (palindromic factorizations) is A344417.
A000569 counts graphical partitions.
A027187 counts partitions of even length, ranked by A028260.
A035363 counts partitions into even parts, ranked by A066207.
A058696 counts partitions of even numbers, ranked by A300061.
A110618 counts partitions with length <= half sum, ranked by A344291.

Programs

  • Haskell
    a025065 = p (1:[2,4..]) where
       p [] _ = 0
       p _  0 = 1
       p ks'@(k:ks) m | m < k     = 0
                      | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 12 2011
    
  • Haskell
    import Data.List (group)
    a025065 = length . filter (<= 1) .
                       map (sum . map ((`mod` 2) . length) . group) . ps 1
       where ps x 0 = [[]]
             ps x y = [t:ts | t <- [x..y], ts <- ps t (y - t)]
    -- Reinhard Zumkeller, Dec 18 2013
    
  • Mathematica
    Map[Length[Select[IntegerPartitions[#], Count[OddQ[Transpose[Tally[#]][[2]]], True] <= 1 &]] &, Range[40]] (* Peter J. C. Moses, Jan 20 2014 *)
    n = 8; Select[IntegerPartitions[n], Count[OddQ[Transpose[Tally[#]][[2]]], True] <= 1 &] (* Peter J. C. Moses, Jan 20 2014 *)
    CoefficientList[Series[1/((1 - x) Product[1 - x^(2 n), {n, 1, 50}]), {x, 0, 60}], x] (* Clark Kimberling, Mar 14 2014 *)
  • PARI
    N=66; q='q+O('q^N); Vec( 1/((1-q)*eta(q^2)) ) \\ Joerg Arndt, Mar 11 2014

Formula

a(n) = A000070(A004526(n)). - Reinhard Zumkeller, Jan 23 2010
G.f.: 1/((1-q)*prod(n>=1, 1-q^(2*n))). [Joerg Arndt, Mar 11 2014]
a(2*k+2) = a(2*k) + A000041(k + 1). - David A. Corneth, May 29 2021
a(n) ~ exp(Pi*sqrt(n/3)) / (2*Pi*sqrt(n)). - Vaclav Kotesovec, Nov 16 2021

Extensions

Edited by N. J. A. Sloane, Dec 29 2007
Prepended a(0)=1, added more terms, Joerg Arndt, Mar 11 2014

A298422 Number of rooted trees with n nodes in which all positive outdegrees are the same.

Original entry on oeis.org

1, 1, 2, 2, 3, 2, 5, 2, 6, 4, 9, 2, 20, 2, 26, 12, 53, 2, 120, 2, 223, 43, 454, 2, 1100, 11, 2182, 215, 4902, 2, 11446, 2, 24744, 1242, 56014, 58, 131258, 2, 293550, 7643, 676928, 2, 1582686, 2, 3627780, 49155, 8436382, 2, 19809464, 50, 46027323, 321202
Offset: 1

Views

Author

Gus Wiseman, Jan 19 2018

Keywords

Comments

Row sums of A298426.

Examples

			The a(9) = 6 trees: ((((((((o)))))))), (o(o(o(oo)))), (o((oo)(oo))), ((oo)(o(oo))), (ooo(oooo)), (oooooooo).
		

Crossrefs

Programs

  • Mathematica
    srut[n_]:=srut[n]=If[n===1,{{}},Select[Join@@Function[c,Union[Sort/@Tuples[srut/@c]]]/@Select[IntegerPartitions[n-1],Function[ptn,And@@(Divisible[#-1,Length[ptn]]&/@ptn)]],SameQ@@Length/@Cases[#,{},{0,Infinity}]&]];
    Table[srut[n]//Length,{n,20}]

Formula

a(n) = 2 <=> n in {A008864}. - Alois P. Heinz, Jan 20 2018

Extensions

a(44)-a(52) from Alois P. Heinz, Jan 20 2018

A349156 Number of integer partitions of n whose mean is not an integer.

Original entry on oeis.org

1, 0, 0, 1, 1, 5, 3, 13, 11, 21, 28, 54, 31, 99, 111, 125, 165, 295, 259, 488, 425, 648, 933, 1253, 943, 1764, 2320, 2629, 2962, 4563, 3897, 6840, 6932, 9187, 11994, 12840, 12682, 21635, 25504, 28892, 28187, 44581, 42896, 63259, 66766, 74463, 104278, 124752
Offset: 0

Views

Author

Gus Wiseman, Nov 14 2021

Keywords

Comments

Equivalently, partitions whose length does not divide their sum.
By conjugation, also the number of integer partitions of n with greatest part not dividing n.

Examples

			The a(3) = 1 through a(8) = 11 partitions:
  (21)  (211)  (32)    (2211)   (43)      (332)
               (41)    (3111)   (52)      (422)
               (221)   (21111)  (61)      (431)
               (311)            (322)     (521)
               (2111)           (331)     (611)
                                (421)     (22211)
                                (511)     (32111)
                                (2221)    (41111)
                                (3211)    (221111)
                                (4111)    (311111)
                                (22111)   (2111111)
                                (31111)
                                (211111)
		

Crossrefs

Below, "!" means either enumerative or set theoretical complement.
The version for nonempty subsets is !A051293.
The complement is counted by A067538, ranked by A316413.
The geometric version is !A067539, strict !A326625, ranked by !A326623.
The strict case is !A102627.
The version for prime factors is A175352, complement A078175.
The version for distinct prime factors is A176587, complement A078174.
The ordered version (compositions) is !A271654, ranked by !A096199.
The multiplicative version (factorizations) is !A326622, geometric !A326028.
The conjugate is ranked by !A326836.
The conjugate strict version is !A326850.
These partitions are ranked by A348551.
A000041 counts integer partitions.
A326567/A326568 give the mean of prime indices, conjugate A326839/A326840.
A236634 counts unbalanced partitions, complement of A047993.
A327472 counts partitions not containing their mean, complement of A237984.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!IntegerQ[Mean[#]]&]],{n,0,30}]

Formula

a(n > 0) = A000041(n) - A067538(n).

A168659 Number of partitions of n such that the number of parts is divisible by the greatest part. Also number of partitions of n such that the greatest part is divisible by the number of parts.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 6, 6, 8, 9, 14, 16, 22, 25, 33, 39, 51, 60, 79, 92, 116, 137, 174, 204, 254, 300, 368, 435, 530, 625, 760, 896, 1076, 1267, 1518, 1780, 2121, 2484, 2946, 3444, 4070, 4749, 5594, 6514, 7637, 8879, 10384, 12043, 14040, 16255
Offset: 1

Views

Author

Vladeta Jovovic, Dec 02 2009

Keywords

Examples

			a(5)=3 because in the partitions [1,1,1,1,1], [1,1,1,2], [1,1,3] the number of parts is divisible by the greatest part; not true for the partitions [1,2,2],[2,3], [1,4], and [5]. - _Emeric Deutsch_, Dec 04 2009
From _Gus Wiseman_, Feb 08 2021: (Start)
The a(1) = 1 through a(10) = 9 partitions of the first type:
  1  11  21   22    311    321     322      332       333        4222
         111  1111  2111   2211    331      2222      4221       4321
                    11111  111111  2221     4211      4311       4411
                                   4111     221111    51111      52111
                                   211111   311111    222111     222211
                                   1111111  11111111  321111     322111
                                                      21111111   331111
                                                      111111111  22111111
                                                                 1111111111
The a(1) = 1 through a(11) = 14 partitions of the second type (A=10, B=11):
  1   2   3    4    5     6     7      8      9       A       B
          21   22   41    42    43     44     63      64      65
                    311   321   61     62     81      82      83
                                322    332    333     622     A1
                                331    611    621     631     632
                                4111   4211   4221    4222    641
                                              4311    4321    911
                                              51111   4411    4322
                                                      52111   4331
                                                              4421
                                                              8111
                                                              52211
                                                              53111
                                                              611111
(End)
		

Crossrefs

Note: A-numbers of Heinz-number sequences are in parentheses below.
The case of equality is A047993 (A106529).
The Heinz numbers of these partitions are A340609/A340610.
If all parts (not just the greatest) are divisors we get A340693 (A340606).
The strict case in the second interpretation is A340828 (A340856).
A006141 = partitions whose length equals their minimum (A324522).
A067538 = partitions whose length/max divides their sum (A316413/A326836).
A200750 = partitions with length coprime to maximum (A340608).
Row sums of A350879.

Programs

  • Maple
    a := proc (n) local pn, ct, j: with(combinat): pn := partition(n): ct := 0: for j to numbpart(n) do if `mod`(nops(pn[j]), max(seq(pn[j][i], i = 1 .. nops(pn[j])))) = 0 then ct := ct+1 else end if end do: ct end proc: seq(a(n), n = 1 .. 50); # Emeric Deutsch, Dec 04 2009
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Divisible[Length[#],Max[#]]&]],{n,30}] (* Gus Wiseman, Feb 08 2021 *)
    nmax = 100; s = 0; Do[s += Normal[Series[Sum[x^((m+1)*k - 1) * Product[(1 - x^(m*k + j - 1))/(1 - x^j), {j, 1, k-1}], {k, 1, (1 + nmax)/(1 + m) + 1}], {x, 0, nmax}]], {m, 1, nmax}]; Rest[CoefficientList[s, x]] (* Vaclav Kotesovec, Oct 18 2024 *)

Formula

G.f.: Sum_{i>=1} Sum_{j>=1} x^((i+1)*j-1) * Product_{k=1..j-1} (1-x^(i*j+k-1))/(1-x^k). - Seiichi Manyama, Jan 24 2022
a(n) ~ c * exp(Pi*sqrt(2*n/3)) / n^(3/2), where c = 0.04628003... - Vaclav Kotesovec, Nov 16 2024

Extensions

Extended by Emeric Deutsch, Dec 04 2009

A064173 Number of partitions of n with positive rank.

Original entry on oeis.org

0, 1, 1, 2, 3, 5, 6, 10, 13, 19, 25, 35, 45, 62, 80, 106, 136, 178, 225, 291, 366, 466, 583, 735, 912, 1140, 1407, 1743, 2140, 2634, 3214, 3932, 4776, 5807, 7022, 8495, 10225, 12313, 14762, 17696, 21136, 25236, 30030, 35722, 42367, 50216, 59368, 70138, 82665
Offset: 1

Views

Author

Vladeta Jovovic, Sep 19 2001

Keywords

Comments

The rank of a partition is the largest summand minus the number of summands.
Also number of partitions of n with negative rank. - Omar E. Pol, Mar 05 2012
Column 1 of A208478. - Omar E. Pol, Mar 11 2012
Number of partitions p of n such that max(max(p), number of parts of p) is not a part of p. - Clark Kimberling, Feb 28 2014
The sequence enumerates the semigroup of partitions of positive rank for each number n. The semigroup is a subsemigroup of the monoid of partitions of nonnegative rank under the binary operation "*": Let A be the positive rank partition (a1,...,ak) where ak > k, and let B=(b1,...bj) with bj > j. Then let A*B be the partition (a1b1,...,a1bj,...,akb1,...,akbj), which has akbj > kj, thus having positive rank. For example, the partition (2,3,4) of 9 has rank 1, and its product with itself is (4,6,6,8,8,9,12,12,16) of 81, which has rank 7. A similar situation holds for partitions of negative rank--they are a subsemigroup of the monoid of nonpositive rank partitions. - Richard Locke Peterson, Jul 15 2018

Examples

			a(20) = p(18) - p(13) + p(5) = 385 - 101 + 7 = 291.
From _Gus Wiseman_, Feb 09 2021: (Start)
The a(2) = 1 through a(9) = 13 partitions of positive rank:
  (2)  (3)  (4)   (5)   (6)    (7)    (8)     (9)
            (31)  (32)  (33)   (43)   (44)    (54)
                  (41)  (42)   (52)   (53)    (63)
                        (51)   (61)   (62)    (72)
                        (411)  (421)  (71)    (81)
                               (511)  (422)   (432)
                                      (431)   (441)
                                      (521)   (522)
                                      (611)   (531)
                                      (5111)  (621)
                                              (711)
                                              (5211)
                                              (6111)
(End)
		

Crossrefs

Note: A-numbers of ranking sequences are in parentheses below.
The negative-rank version is also A064173 (A340788).
The case of odd positive rank is A101707 (A340604).
The case of even positive rank is A101708 (A340605).
These partitions are ranked by (A340787).
A063995/A105806 count partitions by rank.
A072233 counts partitions by sum and length.
A168659 counts partitions whose length is a multiple of the greatest part.
A200750 counts partitions whose length and greatest part are coprime.
- Rank -
A064174 counts partitions of nonnegative/nonpositive rank (A324562/A324521).
A101198 counts partitions of rank 1 (A325233).
A257541 gives the rank of the partition with Heinz number n.
A340601 counts partitions of even rank (A340602).
A340692 counts partitions of odd rank (A340603).
- Balance -
A047993 counts balanced partitions (A106529).
A340599 counts alt-balanced factorizations.
A340653 counts balanced factorizations.

Programs

  • Maple
    A064173 := proc(n)
        a := 0 ;
        for p in combinat[partition](n) do
            r := max(op(p))-nops(p) ;
            if r > 0 then
                a := a+1 ;
            end if;
        end do:
        a ;
    end proc:
    seq(A064173(n),n=0..40) ;# Emeric Deutsch, Dec 11 2004
  • Mathematica
    Table[Count[IntegerPartitions[n], q_ /; First[q] > Length[q]], {n, 24}] (* Clark Kimberling, Feb 12 2014 *)
    Table[Count[IntegerPartitions[n], p_ /; ! MemberQ[p, Max[Max[p], Length[p]]]], {n, 20}] (* Clark Kimberling, Feb 28 2014 *)
    P = PartitionsP;
    a[n_] := (P[n] - Sum[-(-1)^k (P[n - (3k^2 - k)/2] - P[n - (3k^2 + k)/2]), {k, 1, Floor[(1 + Sqrt[1 + 24n])/6]}])/2;
    a /@ Range[48] (* Jean-François Alcover, Jan 11 2020, after Wouter Meeussen in A047993 *)
  • PARI
    my(N=66, x='x+O('x^N)); concat(0, Vec(sum(k=1, N, x^k*prod(j=1, k, (1-x^(k+j-2))/(1-x^j))))) \\ Seiichi Manyama, Jan 25 2022

Formula

a(n) = (A000041(n) - A047993(n))/2.
a(n) = p(n-2) - p(n-7) + p(n-15) - ... - (-1)^k*p(n-(3*k^2+k)/2) + ..., where p() is A000041(). - Vladeta Jovovic, Aug 04 2004
G.f.: Product_{k>=1} (1/(1-q^k)) * Sum_{k>=1} ( (-1)^k * (-q^(3*k^2/2+k/2))) (conjectured). - Thomas Baruchel, May 12 2018
G.f.: Sum_{k>=1} x^k * Product_{j=1..k} (1-x^(k+j-2))/(1-x^j). - Seiichi Manyama, Jan 25 2022
a(n)+A064174(n) = A000041(n). - R. J. Mathar, Feb 22 2023

A340610 Numbers whose number of prime factors (A001222) divides their greatest prime index (A061395).

Original entry on oeis.org

2, 3, 5, 6, 7, 9, 11, 13, 14, 17, 19, 20, 21, 23, 26, 29, 30, 31, 35, 37, 38, 39, 41, 43, 45, 47, 49, 50, 52, 53, 56, 57, 58, 59, 61, 65, 67, 71, 73, 74, 75, 78, 79, 83, 84, 86, 87, 89, 91, 92, 95, 97, 101, 103, 106, 107, 109, 111, 113, 117, 122, 125, 126, 127
Offset: 1

Views

Author

Gus Wiseman, Jan 27 2021

Keywords

Comments

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.

Examples

			The sequence of terms together with their prime indices begins:
     2: {1}        29: {10}       56: {1,1,1,4}
     3: {2}        30: {1,2,3}    57: {2,8}
     5: {3}        31: {11}       58: {1,10}
     6: {1,2}      35: {3,4}      59: {17}
     7: {4}        37: {12}       61: {18}
     9: {2,2}      38: {1,8}      65: {3,6}
    11: {5}        39: {2,6}      67: {19}
    13: {6}        41: {13}       71: {20}
    14: {1,4}      43: {14}       73: {21}
    17: {7}        45: {2,2,3}    74: {1,12}
    19: {8}        47: {15}       75: {2,3,3}
    20: {1,1,3}    49: {4,4}      78: {1,2,6}
    21: {2,4}      50: {1,3,3}    79: {22}
    23: {9}        52: {1,1,6}    83: {23}
    26: {1,6}      53: {16}       84: {1,1,2,4}
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
The case of equality is A047993 (A106529).
The case where all parts are multiples, not just the maximum part, is A143773 (A316428), with strict case A340830, while the case of factorizations is A340853.
These are the Heinz numbers of certain partitions counted by A168659.
The reciprocal version is A340609.
The squarefree case is A340828 (A340856).
A001222 counts prime factors.
A006141 counts partitions whose length equals their minimum (A324522).
A056239 adds up prime indices.
A061395 selects the maximum prime index.
A067538 counts partitions whose length divides their sum (A316413).
A067538 counts partitions whose maximum divides their sum (A326836).
A112798 lists the prime indices of each positive integer.
A200750 counts partitions with length coprime to maximum (A340608).

Programs

  • Maple
    filter:= proc(n) local F,m,g,t;
      F:= ifactors(n)[2];
      m:= add(t[2],t=F);
      g:= numtheory:-pi(max(seq(t[1],t=F)));
      g mod m = 0;
    end proc:
    select(filter, [$2..1000]); # Robert Israel, Feb 08 2021
  • Mathematica
    Select[Range[2,100],Divisible[PrimePi[FactorInteger[#][[-1,1]]],PrimeOmega[#]]&]

Formula

A001222(a(n)) divides A061395(a(n)).

A062319 Number of divisors of n^n, or of A000312(n).

Original entry on oeis.org

1, 1, 3, 4, 9, 6, 49, 8, 25, 19, 121, 12, 325, 14, 225, 256, 65, 18, 703, 20, 861, 484, 529, 24, 1825, 51, 729, 82, 1653, 30, 29791, 32, 161, 1156, 1225, 1296, 5329, 38, 1521, 1600, 4961, 42, 79507, 44, 4005, 4186, 2209, 48, 9457, 99, 5151, 2704, 5565, 54
Offset: 0

Views

Author

Jason Earls, Jul 05 2001

Keywords

Comments

From Gus Wiseman, May 02 2021: (Start)
Conjecture: The number of divisors of n^n equals the number of pairwise coprime ordered n-tuples of divisors of n. Confirmed up to n = 30. For example, the a(1) = 1 through a(5) = 6 tuples are:
(1) (1,1) (1,1,1) (1,1,1,1) (1,1,1,1,1)
(1,2) (1,1,3) (1,1,1,2) (1,1,1,1,5)
(2,1) (1,3,1) (1,1,1,4) (1,1,1,5,1)
(3,1,1) (1,1,2,1) (1,1,5,1,1)
(1,1,4,1) (1,5,1,1,1)
(1,2,1,1) (5,1,1,1,1)
(1,4,1,1)
(2,1,1,1)
(4,1,1,1)
The unordered case (pairwise coprime n-multisets of divisors of n) is counted by A343654.
(End)

Examples

			From _Gus Wiseman_, May 02 2021: (Start)
The a(1) = 1 through a(5) = 6 divisors:
  1  1  1   1    1
     2  3   2    5
     4  9   4    25
        27  8    125
            16   625
            32   3125
            64
            128
            256
(End)
		

Crossrefs

Number of divisors of A000312(n).
Taking Omega instead of sigma gives A066959.
Positions of squares are A173339.
Diagonal n = k of the array A343656.
A000005 counts divisors.
A059481 counts k-multisets of elements of {1..n}.
A334997 counts length-k strict chains of divisors of n.
A343658 counts k-multisets of divisors.
Pairwise coprimality:
- A018892 counts coprime pairs of divisors.
- A084422 counts pairwise coprime subsets of {1..n}.
- A100565 counts pairwise coprime triples of divisors.
- A225520 counts pairwise coprime sets of divisors.
- A343652 counts maximal pairwise coprime sets of divisors.
- A343653 counts pairwise coprime non-singleton sets of divisors > 1.
- A343654 counts pairwise coprime sets of divisors > 1.

Programs

  • Magma
    [NumberOfDivisors(n^n): n in  [0..60]]; // Vincenzo Librandi, Nov 09 2014
    
  • Mathematica
    A062319[n_IntegerQ]:=DivisorSigma[0,n^n]; (* Enrique Pérez Herrero, Nov 09 2010 *)
    Join[{1},DivisorSigma[0,#^#]&/@Range[60]] (* Harvey P. Dale, Jun 06 2024 *)
  • PARI
    je=[]; for(n=0,200,je=concat(je,numdiv(n^n))); je
    
  • PARI
    { for (n=0, 1000, write("b062319.txt", n, " ", numdiv(n^n)); ) } \\ Harry J. Smith, Aug 04 2009
    
  • PARI
    a(n)=local(fm);fm=factor(n);prod(k=1,matsize(fm)[1],fm[k,2]*n+1) \\ Franklin T. Adams-Watters, May 03 2011
    
  • PARI
    a(n) = if(n==0, 1, sumdiv(n, d, n^omega(d))); \\ Seiichi Manyama, May 12 2021
    
  • Python
    from math import prod
    from sympy import factorint
    def A062319(n): return prod(n*d+1 for d in factorint(n).values()) # Chai Wah Wu, Jun 03 2021

Formula

a(n) = A000005(A000312(n)). - Enrique Pérez Herrero, Nov 09 2010
a(2^n) = A002064(n). - Gus Wiseman, May 02 2021
a(prime(n)) = prime(n) + 1. - Gus Wiseman, May 02 2021
a(n) = Product_{i=1..s} (1 + n * m_i) where (m_1,...,m_s) is the sequence of prime multiplicities (prime signature) of n. - Gus Wiseman, May 02 2021
a(n) = Sum_{d|n} n^omega(d) for n > 0. - Seiichi Manyama May 12 2021

A348551 Heinz numbers of integer partitions whose mean is not an integer.

Original entry on oeis.org

1, 6, 12, 14, 15, 18, 20, 24, 26, 33, 35, 36, 38, 40, 42, 44, 45, 48, 50, 51, 52, 54, 56, 58, 60, 63, 65, 66, 69, 70, 72, 74, 75, 76, 77, 80, 86, 92, 93, 95, 96, 102, 104, 106, 108, 112, 114, 117, 119, 120, 122, 123, 124, 126, 130, 132, 135, 136, 140, 141, 142
Offset: 1

Views

Author

Gus Wiseman, Nov 14 2021

Keywords

Comments

Equivalently, partitions whose length does not divide their sum.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The terms and their prime indices begin:
   1: {}
   6: {1,2}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  18: {1,2,2}
  20: {1,1,3}
  24: {1,1,1,2}
  26: {1,6}
  33: {2,5}
  35: {3,4}
  36: {1,1,2,2}
  38: {1,8}
  40: {1,1,1,3}
  42: {1,2,4}
  44: {1,1,5}
  45: {2,2,3}
  48: {1,1,1,1,2}
		

Crossrefs

A version counting nonempty subsets is A000079 - A051293.
A version counting factorizations is A001055 - A326622.
A version counting compositions is A011782 - A271654.
A version for prime factors is A175352, complement A078175.
A version for distinct prime factors A176587, complement A078174.
The complement is A316413, counted by A067538, strict A102627.
The geometric version is the complement of A326623.
The conjugate version is the complement of A326836.
These partitions are counted by A349156.
A000041 counts partitions.
A001222 counts prime factors with multiplicity.
A018818 counts partitions into divisors, ranked by A326841.
A143773 counts partitions into multiples of the length, ranked by A316428.
A236634 counts unbalanced partitions.
A047993 counts balanced partitions, ranked by A106529.
A056239 adds up prime indices, row sums of A112798.
A326567/A326568 gives the mean of prime indices, conjugate A326839/A326840.
A327472 counts partitions not containing their mean, complement A237984.

Programs

  • Maple
    q:= n-> (l-> nops(l)=0 or irem(add(i, i=l), nops(l))>0)(map
            (i-> numtheory[pi](i[1])$i[2], ifactors(n)[2])):
    select(q, [$1..142])[];  # Alois P. Heinz, Nov 19 2021
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],!IntegerQ[Mean[primeMS[#]]]&]
Showing 1-10 of 53 results. Next