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 37 results. Next

A257990 The side-length of the Durfee square of the partition having Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 2, 2, 1, 2, 1, 1
Offset: 1

Views

Author

Emeric Deutsch, May 18 2015

Keywords

Comments

The Durfee square of a partition is the largest square that fits inside the Ferrers board of the partition.
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.
In the Maple program the subprogram B yields the partition with Heinz number n.
First appearance of k is a(prime(k)^k) = k. - Gus Wiseman, Apr 12 2019

Examples

			a(9)=2; indeed, 9 = 3*3 is the Heinz number of the partition [2,2] and, clearly its Durfee square has side-length =2.
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, Reading, Mass. 1976.
  • G. E. Andrews, K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
  • M. Bona, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Crossrefs

Positions of 1's are A093641. Positions of 2's are A325164. Positions of 3's are A307386.

Programs

  • Maple
    with(numtheory): a := proc (p) local B, S, i: 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: S := {}: for i to nops(B(p)) do if i <= B(p)[nops(B(p))+1-i] then S := `union`(S, {i}) else  end if end do: max(S) end proc: seq(a(n), n = 2 .. 146);
    # second Maple program:
    a:= proc(n) local l, t;
          l:= sort(map(i-> numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`);
          for t from nops(l) to 1 by -1 do if l[t]>=t then break fi od; t
        end:
    seq(a(n), n=1..120);  # Alois P. Heinz, May 10 2016
  • Mathematica
    a[n_] := a[n] = Module[{l, t}, l = Reverse[Sort[Flatten[Table[PrimePi[ f[[1]] ], {f, FactorInteger[n]}, {f[[2]]}]]]]; For[t = Length[l], t >= 1, t--, If[l[[t]] >= t, Break[]]]; t]; Table[a[n], {n, 1, 120}] (* Jean-François Alcover, Feb 17 2017, after Alois P. Heinz *)

Formula

For a partition (p_1 >= p_2 >= ... > = p_r) the side-length of its Durfee square is the largest i such that p_i >=i.

Extensions

a(1)=0 prepended by Alois P. Heinz, May 10 2016

A115720 Triangle T(n,k) is the number of partitions of n with Durfee square k.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 3, 0, 4, 1, 0, 5, 2, 0, 6, 5, 0, 7, 8, 0, 8, 14, 0, 9, 20, 1, 0, 10, 30, 2, 0, 11, 40, 5, 0, 12, 55, 10, 0, 13, 70, 18, 0, 14, 91, 30, 0, 15, 112, 49, 0, 16, 140, 74, 1, 0, 17, 168, 110, 2, 0, 18, 204, 158, 5, 0, 19, 240, 221, 10, 0, 20, 285, 302, 20, 0, 21, 330, 407
Offset: 0

Views

Author

Keywords

Comments

T(n,k) is number of partitions of n-k^2 into parts of 2 kinds with at most k of each kind.

Examples

			Triangle starts:
  1;
  0,  1;
  0,  2;
  0,  3;
  0,  4,  1;
  0,  5,  2;
  0,  6,  5;
  0,  7,  8;
  0,  8, 14;
  0,  9, 20,  1;
  0, 10, 30,  2;
From _Gus Wiseman_, Apr 12 2019: (Start)
Row n = 9 counts the following partitions:
  (9)          (54)       (333)
  (81)         (63)
  (711)        (72)
  (6111)       (432)
  (51111)      (441)
  (411111)     (522)
  (3111111)    (531)
  (21111111)   (621)
  (111111111)  (3222)
               (3321)
               (4221)
               (4311)
               (5211)
               (22221)
               (32211)
               (33111)
               (42111)
               (222111)
               (321111)
               (2211111)
(End)
		

Crossrefs

For a version without zeros see A115994. Row lengths are A003059. Row sums are A000041. Column k = 2 is A006918. Column k = 3 is A117485.
Related triangles are A096771, A325188, A325189, A325192, with Heinz-encoded versions A263297, A325169, A065770, A325178.

Programs

  • Maple
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    T:= (n, k)-> add(b(m, k)*b(n-k^2-m, k), m=0..n-k^2):
    seq(seq(T(n, k), k=0..floor(sqrt(n))), n=0..30); # Alois P. Heinz, Apr 09 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; T[n_, k_] := Sum[b[m, k]*b[n-k^2-m, k], {m, 0, n-k^2}]; Table[ T[n, k], {n, 0, 30}, {k, 0, Sqrt[n]}] // Flatten (* Jean-François Alcover, Dec 03 2015, after Alois P. Heinz *)
    durf[ptn_]:=Length[Select[Range[Length[ptn]],ptn[[#]]>=#&]];
    Table[Length[Select[IntegerPartitions[n],durf[#]==k&]],{n,0,10},{k,0,Sqrt[n]}] (* Gus Wiseman, Apr 12 2019 *)

Formula

T(n,k) = Sum_{i=0..n-k^2} P*(i,k)*P*(n-k^2-i), where P*(n,k) = P(n+k,k) is the number of partitions of n objects into at most k parts.

A065770 Number of prime cascades to reach 1, where a prime cascade (A065769) is multiplicative with a(p(m)^k) = p(m-1) * p(m)^(k-1).

Original entry on oeis.org

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

Views

Author

Henry Bottomley, Nov 19 2001

Keywords

Comments

It seems that a(n) <= A297113(n) for all n. Of the first 10000 positive natural numbers, 6454 are such that a(n) = A297113(n). - Antti Karttunen, Dec 31 2017
Also one plus the maximum number of unit steps East or South in the Young diagram of the integer partition with Heinz number n > 1, starting from the upper-left square and ending in a boundary square in the lower-right quadrant. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 06 2019

Examples

			a(50) = 4 since the cascade goes from 50 = 2^1 * 5^2 to 15 = 3^1 * 5^1 to 6 = 2^1 * 3^1 to 2 = 2^1 to 1.
From _Gus Wiseman_, Apr 06 2019: (Start)
The partition with Heinz number 7865 is (6,5,5,3), with diagram
  o o o o o o
  o o o o o
  o o o o o
  o o o
which has longest path from (1,1) to (5,3) of length 6, so a(7865) = 7.
(End)
		

Crossrefs

Cf. A065769.
Differs from A297113 for the first time at n=20, where a(20) = 3, while A297113(20) = 4.

Programs

  • Mathematica
    Table[If[n==1,0,Max@@Total/@Position[PadRight[ConstantArray[1,#]&/@Sort[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]],Greater]],1]-1],{n,100}] (* Gus Wiseman, Apr 06 2019 *)
  • Scheme
    (definec (A065770 n) (if (= 1 n) 0 (+ 1 (A065770 (A065769 n))))) ;; Antti Karttunen, Dec 31 2017

Formula

Inverse of primes, powers of 2 and primorials in sense that a(A000040(n))=n; a(A000079(n))=n; a(A002110(n))=n. If n>0: a(3^n)=n+1; a(2^n*3^k)=n+k; a(p(k)^n)=n+k-1; a(n!)=A022559(n).
a(1) = 0; and for n > 1, a(n) = 1 + A065769(n). - Antti Karttunen, Dec 31 2017

A352822 Number of fixed points y(i) = i, where y is the weakly increasing sequence of prime indices of n.

Original entry on oeis.org

0, 1, 0, 1, 0, 2, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 2, 0, 1, 0, 1, 0, 1, 1, 1, 0, 3, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 2, 1, 0, 1, 0, 2, 0, 1, 0, 2, 0, 2, 0, 1, 0, 1, 0, 1, 1, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 2, 0, 1, 1, 1, 0, 2, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 1, 1, 2, 0, 2, 0, 1, 0
Offset: 1

Views

Author

Gus Wiseman, Apr 05 2022

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 prime indices of 6500 are {1,1,3,3,3,6} with fixed points at positions {1,3,6}, so a(6500) = 3.
		

Crossrefs

* = unproved
Positions of first appearances are A002110.
The triangle version is A238352.
Positions of 0's are A352830, counted by A238394.
Positions of 1's are A352831, counted by A352832.
A version for compositions is A352512, complement A352513, triangle A238349.
The complement is A352823.
The reverse version is A352824, complement A352825.
A000700 counts self-conjugate partitions, ranked by A088902.
A001222 counts prime indices, distinct A001221.
*A001522 counts partitions with a fixed point, ranked by A352827.
A056239 adds up prime indices, row sums of A112798 and A296150.
*A064428 counts partitions without a fixed point, ranked by A352826.
A122111 represents partition conjugation using Heinz numbers.
A124010 gives prime signature, sorted A118914, conjugate rank A238745.
A115720 and A115994 count partitions by their Durfee square.
A238395 counts reversed partitions with a fixed point, ranked by A352872.

Programs

  • Maple
    f:= proc(n) local F,J,t;
      F:= sort(ifactors(n)[2],(s,t) -> s[1] numtheory:-pi(t[1])$t[2], F);
      nops(select(t -> J[t]=t, [$1..nops(J)]));
    end proc:
    map(f, [$1..200]); # Robert Israel, Apr 11 2023
  • Mathematica
    pq[y_]:=Length[Select[Range[Length[y]],#==y[[#]]&]];
    Table[pq[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]],{n,100}]
  • PARI
    A352822(n) = { my(f=factor(n),i=0,c=0); for(k=1,#f~,while(f[k,2], f[k,2]--; i++; c += (i==primepi(f[k,1])))); (c); }; \\ Antti Karttunen, Apr 11 2022

Formula

a(n) = A001222(n) - A352823(n). - Antti Karttunen, Apr 11 2022

Extensions

Data section extended up to 105 terms by Antti Karttunen, Apr 11 2022

A297113 a(1) = 0, a(2) = 1, after which, a(n) = a(n/2) if n is of the form 4k+2, and otherwise a(n) = 1+a(A252463(n)).

Original entry on oeis.org

0, 1, 2, 2, 3, 2, 4, 3, 3, 3, 5, 3, 6, 4, 3, 4, 7, 3, 8, 4, 4, 5, 9, 4, 4, 6, 4, 5, 10, 3, 11, 5, 5, 7, 4, 4, 12, 8, 6, 5, 13, 4, 14, 6, 4, 9, 15, 5, 5, 4, 7, 7, 16, 4, 5, 6, 8, 10, 17, 4, 18, 11, 5, 6, 6, 5, 19, 8, 9, 4, 20, 5, 21, 12, 4, 9, 5, 6, 22, 6, 5, 13, 23, 5, 7, 14, 10, 7, 24, 4, 6, 10, 11, 15, 8, 6, 25
Offset: 1

Views

Author

Antti Karttunen, Dec 26 2017

Keywords

Comments

From Gus Wiseman, Apr 06 2019: (Start)
Also the number of squares in the Young diagram of the integer partition with Heinz number n that are graph-distance 1 from the lower-right boundary, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). For example, the partition (6,5,5,3) with Heinz number 7865 has diagram
o o o o o o
o o o o o
o o o o o
o o o
with inner rim
o
o
o o
o o o
of size 7, so a(7867) = 7.
(End)

Crossrefs

Programs

  • Mathematica
    Table[If[n==1,0,PrimePi[FactorInteger[n][[-1,1]]]+PrimeOmega[n]-PrimeNu[n]],{n,100}] (* Gus Wiseman, Apr 06 2019 *)
  • PARI
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A297113(n) = if(n<=2,n-1,if(n%2,1+A297113(A064989(n)), !(n%4)+A297113(n/2)));
    
  • PARI
    \\ More complex way, after Moebius transform:
    A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n))));
    A297112(n) = sumdiv(n,d,moebius(n/d)*A156552(d));
    A297113(n) = if(1==n,0,1+valuation(A297112(n),2));
    
  • Scheme
    ;; With memoization-macro definec.
    (definec (A297113 n) (cond ((<= n 2) (- n 1)) ((= 2 (modulo n 4)) (A297113 (/ n 2))) (else (+ 1 (A297113 (A252463 n))))))

Formula

a(1) = 0, a(2) = 1, after which, a(n) = a(n/2) if n is of the form 4k+2, and otherwise a(n) = 1+a(A252463(n)) .
For n > 1, a(n) = A001511(A297112(n)), where A297112(n) = Sum_{d|n} moebius(n/d)*A156552(d).
a(n) = A252464(n) - A297155(n).
For n > 1, a(n) = 1+A033265(A156552(n)) = 1+A297167(n) = A046660(n) + A061395(n). - Last two sums added by Antti Karttunen, Sep 02 2018
Other identities. For all n >= 1:
a(A000040(n)) = n. [Each n occurs for the first time at the n-th prime.]

A071724 a(n) = 3*binomial(2n, n-1)/(n+2), n > 0, with a(0)=1.

Original entry on oeis.org

1, 1, 3, 9, 28, 90, 297, 1001, 3432, 11934, 41990, 149226, 534888, 1931540, 7020405, 25662825, 94287120, 347993910, 1289624490, 4796857230, 17902146600, 67016296620, 251577050010, 946844533674, 3572042254128, 13505406670700
Offset: 0

Views

Author

N. J. A. Sloane, Jun 06 2002

Keywords

Comments

Number of standard tableaux of shape (n+1,n-1) (n>=1). - Emeric Deutsch, May 30 2004
From Gus Wiseman, Apr 12 2019: (Start)
Also the number of integer partitions (of any positive integer) such that n is the maximum number of unit steps East or South in the Young diagram starting from the upper-left square and ending in a boundary square in the lower-right quadrant. Also the number of integer partitions fitting in a triangular partition of length n but not of length n - 1. For example, the a(0) = 1 through a(4) = 9 partitions are:
() (1) (2) (3)
(11) (22)
(21) (31)
(32)
(111)
(211)
(221)
(311)
(321)
(End)
The sequence (-1)^(n+1)*a(n), for n >= 1 and +1 for n = 0, is the so-called Z-sequence of the Riordan triangle A158909. For the notion of Z- and A-sequences for Riordan arrays see the W. Lang link under A006232 with details and references. - Wolfdieter Lang, Oct 22 2019

Crossrefs

Number of times n appears in A065770.
Column sums of A325189.
Row sums of A030237.

Programs

  • Magma
    [1] cat [3*Binomial(2*n,n-1)/(n+2): n in [1..29]]; // Vincenzo Librandi, Jul 12 2017
    
  • Maple
    A071724:= n-> 3*binomial(2*n, n-1)/(n+2); 1,seq(A071724(n), n=1..30); # G. C. Greubel, Mar 17 2021
  • Mathematica
    Join[{1}, Table[3Binomial[2n, n-1]/(n+2), {n,1,30}]] (* Vincenzo Librandi, Jul 12 2017 *)
    nn=7;
    otbmax[ptn_]:=Max@@MapIndexed[#1+#2[[1]]-1&,Append[ptn,0]];
    allip=Join@@Table[IntegerPartitions[n],{n,0,nn*(nn+1)/2}];
    Table[Length[Select[allip,otbmax[#]==n&]],{n,0,nn}] (* Gus Wiseman, Apr 12 2019 *)
  • PARI
    a(n)=if(n<1,n==0,3*(2*n)!/(n+2)!/(n-1)!)
    
  • Sage
    [1]+[3*n*catalan_number(n)/(n+2) for n in (1..30)] # G. C. Greubel, Mar 17 2021

Formula

a(n) = A000245(n), n>0.
G.f.: (C(x)-1)*(1-x)/x = (1 + x^2 * C(x)^3)*C(x), where C(x) is g.f. for Catalan numbers, A000108.
G.f.: ((1-sqrt(1-4*x))/(2*x)-1)*(1-x)/x = A(x) satisfies x^2*A(x)^2 + (x-1)*(2*x-1)*A(x) + (x-1)^2 = 0.
G.f.: 1 + x*C(x)^3, where C(x) is g.f. for the Catalan numbers (A000108). Sequence without the first term is the 3-fold convolution of the Catalan sequence. - Emeric Deutsch, May 30 2004
a(n) is the n-th moment of the function defined on the segment (0, 4) of x axis: a(n) = Integral_{x=0..4} x^n*(-x^(1/2)*cos(3*arcsin((1/2)*x^(1/2)))/Pi) dx, n=0, 1... . - Karol A. Penson, Sep 29 2004
D-finite with recurrence -(n+2)*(n-1)*a(n) + 2*n*(2*n-1)*a(n-1) = 0. - R. J. Mathar, Jul 10 2017
a(n) ~ c*2^(2*n)*n^(-3/2), where c = 3/sqrt(Pi). - Stefano Spezia, Sep 23 2022
From Amiram Eldar, Sep 29 2022: (Start)
Sum_{n>=0} 1/a(n) = 14*(Pi/(3*sqrt(3)) + 1)/9.
Sum_{n>=0} (-1)^n/a(n) = 18/25 - 164*log(phi)/(75*sqrt(5)), where phi is the golden ratio (A001622). (End)

A096771 Triangle read by rows: T(n,m) is the number of partitions of n that (just) fit inside an m X m box, but not in an (m-1) X (m-1) box. Partitions of n with Max(max part, length) = m.

Original entry on oeis.org

1, 0, 2, 0, 1, 2, 0, 1, 2, 2, 0, 0, 3, 2, 2, 0, 0, 3, 4, 2, 2, 0, 0, 2, 5, 4, 2, 2, 0, 0, 1, 7, 6, 4, 2, 2, 0, 0, 1, 6, 9, 6, 4, 2, 2, 0, 0, 0, 7, 11, 10, 6, 4, 2, 2, 0, 0, 0, 5, 14, 13, 10, 6, 4, 2, 2, 0, 0, 0, 5, 15, 19, 14, 10, 6, 4, 2, 2, 0, 0, 0, 3, 17, 22, 21, 14, 10, 6, 4, 2, 2, 0, 0, 0, 2, 17, 29
Offset: 1

Views

Author

Wouter Meeussen, Aug 21 2004

Keywords

Comments

Row sums are A000041. Columns are finite and sum to A051924. The final floor(n/2) terms of each row are the reverse of the initial terms of 2*A000041.

Examples

			T(5,3)=3, counting 32, 311 and 221.
From _Gus Wiseman_, Apr 12 2019: (Start)
Triangle begins:
  1
  0  2
  0  1  2
  0  1  2  2
  0  0  3  2  2
  0  0  3  4  2  2
  0  0  2  5  4  2  2
  0  0  1  7  6  4  2  2
  0  0  1  6  9  6  4  2  2
  0  0  0  7 11 10  6  4  2  2
  0  0  0  5 14 13 10  6  4  2  2
  0  0  0  5 15 19 14 10  6  4  2  2
  0  0  0  3 17 22 21 14 10  6  4  2  2
  0  0  0  2 17 29 27 22 14 10  6  4  2  2
  0  0  0  1 17 33 36 29 22 14 10  6  4  2  2
  0  0  0  1 15 39 45 41 30 22 14 10  6  4  2  2
  0  0  0  0 14 41 57 52 43 30 22 14 10  6  4  2  2
  0  0  0  0 11 47 67 69 57 44 30 22 14 10  6  4  2  2
  0  0  0  0  9 46 81 85 76 59 44 30 22 14 10  6  4  2  2
(End)
		

Crossrefs

A version with reflected rows is A338621.
Related triangles are A115720, A325188, A325189, A325192, A325200, with Heinz-encoded versions A257990, A325169, A065770, A325178, A325195.

Programs

  • Mathematica
    Table[Count[Partitions[n], q_List /; Max[Length[q], Max[q]]===k], {n, 16}, {k, n}]
  • PARI
    row(n)={my(r=vector(n)); forpart(p=n, r[max(#p,p[#p])]++); r} \\ Andrew Howroyd, Jan 12 2024

Formula

Sum_{k>=1} k*T(n,k) = A368985(n). - Andrew Howroyd, Jan 12 2024

A325188 Regular triangle read by rows where T(n,k) is the number of integer partitions of n with origin-to-boundary graph-distance equal to k.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 0, 2, 1, 0, 0, 2, 3, 0, 0, 0, 2, 5, 0, 0, 0, 0, 2, 8, 1, 0, 0, 0, 0, 2, 9, 4, 0, 0, 0, 0, 0, 2, 12, 8, 0, 0, 0, 0, 0, 0, 2, 13, 15, 0, 0, 0, 0, 0, 0, 0, 2, 16, 23, 1, 0, 0, 0, 0, 0, 0, 0, 2, 17, 32, 5, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Apr 11 2019

Keywords

Comments

The origin-to-boundary graph-distance of a Young diagram is the minimum number of unit steps right or down from the upper-left square to a nonsquare in the lower-right quadrant. It is also the side-length of the maximum triangular partition contained inside the diagram.

Examples

			Triangle begins:
  1
  0  1
  0  2  0
  0  2  1  0
  0  2  3  0  0
  0  2  5  0  0  0
  0  2  8  1  0  0  0
  0  2  9  4  0  0  0  0
  0  2 12  8  0  0  0  0  0
  0  2 13 15  0  0  0  0  0  0
  0  2 16 23  1  0  0  0  0  0  0
  0  2 17 32  5  0  0  0  0  0  0  0
  0  2 20 43 12  0  0  0  0  0  0  0  0
  0  2 21 54 24  0  0  0  0  0  0  0  0  0
  0  2 24 67 42  0  0  0  0  0  0  0  0  0  0
  0  2 25 82 66  1  0  0  0  0  0  0  0  0  0  0
		

Crossrefs

Programs

  • Mathematica
    otb[ptn_]:=Min@@MapIndexed[#1+#2[[1]]-1&,Append[ptn,0]];
    Table[Length[Select[IntegerPartitions[n],otb[#]==k&]],{n,0,15},{k,0,n}]
  • PARI
    row(n)={my(r=vector(n+1)); forpart(p=n, my(w=#p); for(i=1, #p, w=min(w,#p-i+p[i])); r[w+1]++); r} \\ Andrew Howroyd, Jan 12 2024

Formula

Sum_{k=1..n} k*T(n,k) = A368986(n).

A325189 Regular triangle read by rows where T(n,k) is the number of integer partitions of n with maximum origin-to-boundary graph-distance equal to k.

Original entry on oeis.org

1, 0, 1, 0, 0, 2, 0, 0, 1, 2, 0, 0, 0, 3, 2, 0, 0, 0, 3, 2, 2, 0, 0, 0, 1, 6, 2, 2, 0, 0, 0, 0, 7, 4, 2, 2, 0, 0, 0, 0, 6, 8, 4, 2, 2, 0, 0, 0, 0, 4, 12, 6, 4, 2, 2, 0, 0, 0, 0, 1, 15, 12, 6, 4, 2, 2, 0, 0, 0, 0, 0, 17, 15, 10, 6, 4, 2, 2
Offset: 0

Views

Author

Gus Wiseman, Apr 11 2019

Keywords

Comments

The maximum origin-to-boundary graph-distance of an integer partition is one plus the maximum number of unit steps East or South in the Young diagram that can be followed, starting from the upper-left square, to reach a boundary square in the lower-right quadrant. It is also the side-length of the minimum triangular partition containing the diagram.

Examples

			Triangle begins:
  1
  0  1
  0  0  2
  0  0  1  2
  0  0  0  3  2
  0  0  0  3  2  2
  0  0  0  1  6  2  2
  0  0  0  0  7  4  2  2
  0  0  0  0  6  8  4  2  2
  0  0  0  0  4 12  6  4  2  2
  0  0  0  0  1 15 12  6  4  2  2
  0  0  0  0  0 17 15 10  6  4  2  2
  0  0  0  0  0 14 23 16 10  6  4  2  2
  0  0  0  0  0 10 30 23 14 10  6  4  2  2
  0  0  0  0  0  5 39 29 24 14 10  6  4  2  2
  0  0  0  0  0  1 42 42 31 22 14 10  6  4  2  2
Row 9 counts the following partitions:
  (432)   (54)     (63)      (72)       (81)        (9)
  (3321)  (333)    (621)     (711)      (21111111)  (111111111)
  (4221)  (441)    (6111)    (2211111)
  (4311)  (522)    (222111)  (3111111)
          (531)    (321111)
          (3222)   (411111)
          (5211)
          (22221)
          (32211)
          (33111)
          (42111)
          (51111)
		

Crossrefs

Programs

  • Mathematica
    otbmax[ptn_]:=Max@@MapIndexed[#1+#2[[1]]-1&,Append[ptn,0]];
    Table[Length[Select[IntegerPartitions[n],otbmax[#]==k&]],{n,0,15},{k,0,n}]
  • PARI
    row(n)={my(r=vector(n+1)); forpart(p=n, my(w=0); for(i=1, #p, w=max(w,#p-i+p[i])); r[w+1]++); r} \\ Andrew Howroyd, Jan 12 2024

Formula

Sum_{k=1..n} k*T(n,k) = A366157(n). - Andrew Howroyd, Jan 12 2024

A325225 Lesser of the number of prime factors of n counted with multiplicity and the maximum prime index of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 3, 2, 2, 1, 2, 2, 2, 2, 3, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 3, 1, 3, 1, 3, 3, 2, 1, 2, 2, 3, 2, 3, 1, 2, 2, 4, 2, 2, 1, 3, 1, 2, 3, 1, 2, 3, 1, 3, 2, 3, 1, 2, 1, 2, 3, 3, 2, 3, 1, 3, 2, 2, 1, 4, 2, 2, 2, 4, 1, 3, 2, 3, 2, 2, 2, 2, 1, 3, 3, 3, 1, 3, 1, 4, 3
Offset: 1

Views

Author

Gus Wiseman, Apr 12 2019

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

			88 has 4 prime indices {1,1,1,5}, the maximum of which is 5, so a(88) = min(4,5) = 4.
		

Crossrefs

Positions of 1's are A174090. Positions of 2's are A325229.

Programs

  • Mathematica
    Table[Min[PrimeOmega[n],PrimePi[FactorInteger[n][[-1,1]]]],{n,100}]
  • PARI
    A061395(n) = if(1==n, 0, primepi(vecmax(factor(n)[, 1])));
    A325225(n) = min(bigomega(n), A061395(n)); \\ Antti Karttunen, Apr 14 2019

Formula

a(n) = min(A001222(n), A061395(n)).

Extensions

More terms from Antti Karttunen, Apr 14 2019
Showing 1-10 of 37 results. Next