A325359 -id:A325359 - cp's OEIS frontend

A307824 Heinz numbers of integer partitions whose augmented differences are all equal.

1, 2, 3, 4, 5, 7, 8, 11, 13, 15, 16, 17, 19, 23, 29, 31, 32, 37, 41, 43, 47, 53, 55, 59, 61, 64, 67, 71, 73, 79, 83, 89, 97, 101, 103, 105, 107, 109, 113, 119, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227

Offset: 1

Gus Wiseman, May 03 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).
The enumeration of these partitions by sum is given by A129654.

			The sequence of terms together with their prime indices begins:
    1: {}
    2: {1}
    3: {2}
    4: {1,1}
    5: {3}
    7: {4}
    8: {1,1,1}
   11: {5}
   13: {6}
   15: {2,3}
   16: {1,1,1,1}
   17: {7}
   19: {8}
   23: {9}
   29: {10}
   31: {11}
   32: {1,1,1,1,1}
   37: {12}
   41: {13}
   43: {14}

Gus Wiseman, Sequences counting and ranking integer partitions by the differences of their successive parts.

Cf. A049988, A056239, A093641, A112798, A129654, A325327, A325328, A325351, A325359, A325366, A325389, A325394, A325395, A325396.

primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
aug[y_]:=Table[If[i

A325407 Nonprime Heinz numbers of multiples of triangular partitions, or of finite arithmetic progressions with offset 0.

Original entry on oeis.org

1, 6, 21, 30, 65, 133, 210, 273, 319, 481, 731, 1007, 1403, 1495, 2059, 2310, 2449, 3293, 4141, 4601, 4921, 5187, 5311, 6943, 8201, 9211, 10921, 12283, 13213, 14993, 15247, 16517, 19847, 22213, 24139, 25853, 28141, 29341, 29539, 30030, 31753, 37211, 40741

Offset: 1

Gus Wiseman, May 03 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are numbers of the form Product_{k = 1...b} prime(k * c) for some b > 1 and c > 0.

			The sequence of terms together with their prime indices begins:
      1: {}
      6: {1,2}
     21: {2,4}
     30: {1,2,3}
     65: {3,6}
    133: {4,8}
    210: {1,2,3,4}
    273: {2,4,6}
    319: {5,10}
    481: {6,12}
    731: {7,14}
   1007: {8,16}
   1403: {9,18}
   1495: {3,6,9}
   2059: {10,20}
   2310: {1,2,3,4,5}
   2449: {11,22}
   3293: {12,24}
   4141: {13,26}
   4601: {14,28}

Cf. A007294, A007862, A049988, A056239, A112798, A325327, A325328, A325355, A325359, A325367, A325390.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[10000],!PrimeQ[#]&&SameQ@@Differences[Prepend[primeMS[#],0]]&]

A325355 One plus the number of steps applying A325351 (Heinz number of augmented differences of reversed prime indices) to reach a fixed point.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 2, 1, 1, 4, 1, 1, 1, 2, 1, 2, 1, 1, 3, 1, 1, 5, 1, 4, 2, 1, 1, 6, 1, 1, 4, 1, 1, 3, 1, 1, 1, 2, 2, 7, 1, 1, 2, 3, 1, 8, 1, 1, 3, 1, 1, 4, 1, 5, 5, 1, 1, 9, 4, 1, 2, 1, 1, 3, 1, 5, 6, 1, 1, 2, 1, 1, 4, 4, 1, 10, 1, 1, 3, 5, 1, 11, 1, 6, 1, 1, 2, 5, 2, 1, 7, 1, 1, 3

Offset: 1

Gus Wiseman, Apr 23 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).
The fixed points of A325351 are the Heinz numbers of hooks A093641.

			Repeatedly applying A325351 starting with 78 gives 78 -> 66 -> 42 -> 30 -> 18 -> 12, and 12 is a fixed point, so a(78) = 6.

Positions of 2's are A325359.

Cf. A056239, A093641, A112798, A130091, A289509, A325351, A325352, A325366, A325389, A325394, A325395, A325396.

primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
aug[y_]:=Table[If[i

augdiffs(n) = { my(diffs=List([]), f=factor(n), prevpi, pi=0, i=#f~); while(i, prevpi=pi; pi = primepi(f[i, 1]); if(prevpi, listput(diffs, 1+(prevpi-pi))); if(f[i, 2]>1, f[i, 2]--, i--)); if(pi, listput(diffs,pi)); Vec(diffs); };
A325351(n) = factorback(apply(prime,augdiffs(n)));
A325355(n) = { my(u=A325351(n)); if(u==n,1,1+A325355(u)); }; \\ Antti Karttunen, Nov 16 2019

More terms from Antti Karttunen, Nov 16 2019

A325358 Number of integer partitions of n whose augmented differences are strictly decreasing.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 3, 4, 4, 5, 6, 6, 7, 9, 10, 11, 13, 14, 15, 18, 20, 21, 24, 26, 28, 33, 36, 38, 43, 46, 49, 56, 60, 63, 71, 76, 80, 90, 96, 100, 112, 120, 125, 139, 149, 155, 171, 183, 190, 208, 223, 232, 252, 269, 280, 304, 325, 338, 364, 387, 403

Offset: 0

Gus Wiseman, May 01 2019

nonn

The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).

The Heinz numbers of these partitions are given by A325396.

			The a(1) = 1 through a(11) = 6 partitions:
  (1)  (2)  (3)   (4)   (5)   (6)   (7)    (8)    (9)    (10)   (11)
            (21)  (31)  (41)  (42)  (52)   (62)   (63)   (73)   (83)
                              (51)  (61)   (71)   (72)   (82)   (92)
                                    (421)  (521)  (81)   (91)   (101)
                                                  (621)  (631)  (731)
                                                         (721)  (821)

Fausto A. C. Cariboni, Table of n, a(n) for n = 0..1000
Gus Wiseman, Sequences counting and ranking integer partitions by the differences of their successive parts.

Cf. A049988, A240026, A320466, A325349, A325350, A325351, A325356, A325357, A325359, A325393.

Mathematica
```
aug[y_]:=Table[If[i
				
```

A069562 Numbers, m, whose odd part (largest odd divisor, A000265(m)) is a nontrivial square.

Original entry on oeis.org

9, 18, 25, 36, 49, 50, 72, 81, 98, 100, 121, 144, 162, 169, 196, 200, 225, 242, 288, 289, 324, 338, 361, 392, 400, 441, 450, 484, 529, 576, 578, 625, 648, 676, 722, 729, 784, 800, 841, 882, 900, 961, 968, 1058, 1089, 1152, 1156, 1225, 1250, 1296, 1352, 1369

Offset: 1

Benoit Cloitre, Apr 18 2002

Previous name: sum(d|n,6d/(2+mu(d))) is odd, where mu(.) is the Moebius function, A008683.
From Peter Munn, Jul 06 2020: (Start)
Numbers that have an odd number of odd nonsquarefree divisors.
[Proof of equivalence to the name, where m denotes a positive integer:
(1) These properties are equivalent: (a) m has an even number of odd squarefree divisors; (b) m has a nontrivial odd part.
(2) These properties are equivalent: (a) m has an odd number of odd divisors; (b) the odd part of m is square.
(3) m satisfies the condition at the start of this comment if and only if (1)(a) and (2)(a) are both true or both false.
(4) The trivial odd part, 1, is a square, so (1)(b) and (2)(b) cannot both be false, which (from (1), (2)) means (1)(a) and (2)(a) cannot both be false.
(5) From (3), (4), m satisfies the condition at the start of this comment if and only if (1)(a) and (2)(a) are true.
(6) m satisfies the condition in the name if and only if (1)(b) and (2)(b) are true, which (from (1), (2)) is equivalent to (1)(a) and (2)(a) being true, and hence from (5), to m satisfying the condition at the start of this comment.]
(End)
Numbers whose sum of non-unitary divisors (A048146) is odd. - Amiram Eldar, Sep 16 2024

			To determine the odd part of 18, remove all factors of 2, leaving 9. 9 is a nontrivial square, so 18 is in the sequence. - _Peter Munn_, Jul 06 2020

David A. Corneth, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Odd part.

A000265, A008683 are used in definitions of this sequence.
Lists of numbers whose odd part satisfies other conditions: A028982 (square), A028983 (nonsquare), A029747 (less than 6), A029750 (less than 8), A036349 (even number of prime factors), A038550 (prime), A070776 U {1} (power of a prime), A072502 (square of a prime), A091067 (has form 4k+3), A091072 (has form 4k+1), A093641 (noncomposite), A105441 (composite), A116451 (greater than 4), A116882 (less than or equal to even part), A116883 (greater than or equal to even part), A122132 (squarefree), A229829 (7-rough), A236206 (11-rough), A260488\{0} (has form 6k+1), A325359 (proper prime power), A335657 (odd number of prime factors), A336101 (prime power).
Cf. A048146, A091476.

Select[Range[1000], (odd = #/2^IntegerExponent[#, 2]) > 1 && IntegerQ @ Sqrt[odd] &] (* Amiram Eldar, Sep 29 2020 *)

upto(n) = { my(res = List()); forstep(i = 3, sqrtint(n), 2, for(j = 0, logint(n\i^2, 2), listput(res, i^2<David A. Corneth, Sep 28 2020

Sum_{n>=1} 1/a(n) = 2 * Sum_{k>=1} 1/(2*k+1)^2 = Pi^2/4 - 2 = A091476 - 2 = 0.467401... - Amiram Eldar, Feb 18 2021

New name from Peter Munn, Jul 06 2020

A325459 Sum of numbers of nontrivial divisors (greater than 1 and less than k) of k for k = 1..n.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 3, 3, 5, 6, 8, 8, 12, 12, 14, 16, 19, 19, 23, 23, 27, 29, 31, 31, 37, 38, 40, 42, 46, 46, 52, 52, 56, 58, 60, 62, 69, 69, 71, 73, 79, 79, 85, 85, 89, 93, 95, 95, 103, 104, 108, 110, 114, 114, 120, 122, 128, 130, 132, 132, 142

Offset: 0

Gus Wiseman, May 04 2019

nonn

Also the number of integer partitions of n that are not hooks but whose augmented differences are hooks (original name). The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and otherwise aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).

This sequence counts integer partitions with any number of ones and one part > 1 which appears at least twice. The Heinz numbers of these partitions are given by A325359.

			The a(4) = 1 through a(10) = 8 partitions:
  (22)  (221)  (33)    (331)    (44)      (333)      (55)
               (222)   (2221)   (2222)    (441)      (3331)
               (2211)  (22111)  (3311)    (22221)    (4411)
                                (22211)   (33111)    (22222)
                                (221111)  (222111)   (222211)
                                          (2211111)  (331111)
                                                     (2221111)
                                                     (22111111)

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

Cf. A049988, A093641, A325349, A325351, A325355, A325356, A325357, A325358.

Cf. A070824, A006218.

a:= proc(n) option remember; `if`(n<2, 0,
      numtheory[tau](n)-2+a(n-1))
    end:
seq(a(n), n=0..100);  # Alois P. Heinz, Oct 11 2019

Table[Length[Select[IntegerPartitions[n],MatchQ[#,{x_,y__,1...}/;x>1&&SameQ[x,y]]&]],{n,0,30}]
(* Second program: *)
a[n_] := a[n] = If[n<2, 0, DivisorSigma[0, n] - 2 + a[n-1]];
a /@ Range[0, 100] (* Jean-François Alcover, May 20 2021, after Alois P. Heinz *)

from math import isqrt
def A325459(n): return 0 if n == 0 else (lambda m: 2*(sum(n//k for k in range(1, m+1))-n)+(1-m)*(1+m))(isqrt(n)) # Chai Wah Wu, Oct 07 2021

From M. F. Hasler, Oct 11 2019: (Start)
a(n) = A006218(n) - 2*n + 1, in terms of partial sums of number of divisors.
a(n) = Sum_{k=1..n} A070824(k): partial sums of A070824 = number of nontrivial divisors. (End)

Name changed at the suggestion of Patrick James Smalley-Wall and Luc Rousseau by Gus Wiseman, Oct 11 2019

A325458 Triangle read by rows where T(n,k) is the number of integer partitions of n with largest hook of size k, i.e., with (largest part) + (number of parts) - 1 = k.

Original entry on oeis.org

1, 0, 1, 0, 0, 2, 0, 0, 0, 3, 0, 0, 0, 1, 4, 0, 0, 0, 0, 2, 5, 0, 0, 0, 0, 2, 3, 6, 0, 0, 0, 0, 0, 4, 4, 7, 0, 0, 0, 0, 0, 3, 6, 5, 8, 0, 0, 0, 0, 0, 1, 6, 8, 6, 9, 0, 0, 0, 0, 0, 0, 6, 9, 10, 7, 10, 0, 0, 0, 0, 0, 0, 2, 11, 12, 12, 8, 11

Offset: 0

Gus Wiseman, May 04 2019

Conjectured to be equal to A049597.

			Triangle begins:
  1
  0  1
  0  0  2
  0  0  0  3
  0  0  0  1  4
  0  0  0  0  2  5
  0  0  0  0  2  3  6
  0  0  0  0  0  4  4  7
  0  0  0  0  0  3  6  5  8
  0  0  0  0  0  1  6  8  6  9
  0  0  0  0  0  0  6  9 10  7 10
  0  0  0  0  0  0  2 11 12 12  8 11
  0  0  0  0  0  0  2  9 16 15 14  9 12
  0  0  0  0  0  0  0  7 16 21 18 16 10 13
  0  0  0  0  0  0  0  4 18 23 26 21 18 11 14
  0  0  0  0  0  0  0  3 12 29 30 31 24 20 12 15
  0  0  0  0  0  0  0  1 12 27 40 37 36 27 22 13 16
  0  0  0  0  0  0  0  0  8 26 42 51 44 41 30 24 14 17
  0  0  0  0  0  0  0  0  6 23 48 57 62 51 46 33 26 15 18
  0  0  0  0  0  0  0  0  2 21 44 70 72 73 58 51 36 28 16 19
Row n = 9 counts the following partitions:
  (333)  (54)     (63)      (72)       (9)
         (432)    (522)     (621)      (81)
         (441)    (531)     (5211)     (711)
         (3222)   (4221)    (42111)    (6111)
         (3321)   (4311)    (321111)   (51111)
         (22221)  (32211)   (2211111)  (411111)
                  (33111)              (3111111)
                  (222111)             (21111111)
                                       (111111111)

Row sums are A000041.
Column sums are 2^(k - 1) for k > 0.
Cf. A008284, A049597, A093641, A116608, A325242, A325349, A325351, A325355, A325359, A325459.

Table[Length[Select[IntegerPartitions[n],If[n==0,k==0,First[#]+Length[#]-1==k]&]],{n,0,19},{k,0,n}]

Franklin T. Adams-Watters has conjectured at A049597 that the k-th column gives the coefficients of the sum of Gaussian polynomials [k,m] for m = 0..k.

cp's OEIS Frontend

A307824 Heinz numbers of integer partitions whose augmented differences are all equal.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A325407 Nonprime Heinz numbers of multiples of triangular partitions, or of finite arithmetic progressions with offset 0.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A325355 One plus the number of steps applying A325351 (Heinz number of augmented differences of reversed prime indices) to reach a fixed point.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A325358 Number of integer partitions of n whose augmented differences are strictly decreasing.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A069562 Numbers, m, whose odd part (largest odd divisor, A000265(m)) is a nontrivial square.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A325459 Sum of numbers of nontrivial divisors (greater than 1 and less than k) of k for k = 1..n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Python

Formula

Extensions

A325458 Triangle read by rows where T(n,k) is the number of integer partitions of n with largest hook of size k, i.e., with (largest part) + (number of parts) - 1 = k.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula