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.

Previous Showing 11-20 of 115 results. Next

A240026 Number of partitions of n such that the successive differences of consecutive parts are nondecreasing.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 10, 12, 16, 21, 27, 32, 43, 50, 60, 75, 90, 103, 128, 146, 170, 203, 234, 264, 315, 355, 402, 467, 530, 589, 684, 764, 851, 969, 1083, 1195, 1360, 1504, 1659, 1863, 2063, 2258, 2531, 2779, 3039, 3379, 3709, 4032, 4474, 4880, 5304, 5846, 6373, 6891, 7578, 8227, 8894, 9727, 10550, 11357, 12405, 13404, 14419
Offset: 0

Views

Author

Joerg Arndt, Mar 31 2014

Keywords

Comments

Partitions (p(1), p(2), ..., p(m)) such that p(k-1) - p(k-2) <= p(k) - p(k-1) for all k >= 3.
The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (6,3,1) are (-3,-2). Then a(n) is the number of integer partitions of n whose differences are weakly increasing. The Heinz numbers of these partitions are given by A325360. Of course, the number of such integer partitions of n is also the number of reversed integer partitions of n whose differences are weakly increasing, which is the author's interpretation. - Gus Wiseman, May 03 2019

Examples

			There are a(10) = 27 such partitions of 10:
01:  [ 1 1 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 1 1 1 2 ]
03:  [ 1 1 1 1 1 1 1 3 ]
04:  [ 1 1 1 1 1 1 4 ]
05:  [ 1 1 1 1 1 2 3 ]
06:  [ 1 1 1 1 1 5 ]
07:  [ 1 1 1 1 2 4 ]
08:  [ 1 1 1 1 6 ]
09:  [ 1 1 1 2 5 ]
10:  [ 1 1 1 7 ]
11:  [ 1 1 2 6 ]
12:  [ 1 1 3 5 ]
13:  [ 1 1 8 ]
14:  [ 1 2 3 4 ]
15:  [ 1 2 7 ]
16:  [ 1 3 6 ]
17:  [ 1 9 ]
18:  [ 2 2 2 2 2 ]
19:  [ 2 2 2 4 ]
20:  [ 2 2 6 ]
21:  [ 2 3 5 ]
22:  [ 2 8 ]
23:  [ 3 3 4 ]
24:  [ 3 7 ]
25:  [ 4 6 ]
26:  [ 5 5 ]
27:  [ 10 ]
		

Crossrefs

Cf. A240027 (strictly increasing differences).
Cf. A179255 (distinct parts, nondecreasing), A179254 (distinct parts, strictly increasing).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],OrderedQ[Differences[#]]&]],{n,0,30}] (* Gus Wiseman, May 03 2019 *)
  • Ruby
    def partition(n, min, max)
      return [[]] if n == 0
      [max, n].min.downto(min).flat_map{|i| partition(n - i, min, i).map{|rest| [i, *rest]}}
    end
    def f(n)
      return 1 if n == 0
      cnt = 0
      partition(n, 1, n).each{|ary|
        ary0 = (1..ary.size - 1).map{|i| ary[i - 1] - ary[i]}
        cnt += 1 if ary0.sort == ary0.reverse
      }
      cnt
    end
    def A240026(n)
      (0..n).map{|i| f(i)}
    end
    p A240026(50) # Seiichi Manyama, Oct 13 2018

A051533 Numbers that are the sum of two positive triangular numbers.

Original entry on oeis.org

2, 4, 6, 7, 9, 11, 12, 13, 16, 18, 20, 21, 22, 24, 25, 27, 29, 30, 31, 34, 36, 37, 38, 39, 42, 43, 46, 48, 49, 51, 55, 56, 57, 58, 60, 61, 64, 65, 66, 67, 69, 70, 72, 73, 76, 79, 81, 83, 84, 87, 88, 90, 91, 92, 93, 94, 97, 99, 100, 101, 102, 106, 108
Offset: 1

Views

Author

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de)

Keywords

Comments

Numbers n such that 8n+2 is in A085989. - Robert Israel, Mar 06 2017

Examples

			666 is in the sequence because we can write 666 = 435 + 231 = binomial(22,2) + binomial(30,2).
		

Crossrefs

Cf. A000217, A020756 (sums of two triangular numbers), A001481 (sums of two squares), A007294, A051611 (complement).
Cf. A061336: minimal number of triangular numbers that sum up to n.
Cf. A085989.

Programs

  • Haskell
    a051533 n = a051533_list !! (n-1)
    a051533_list = filter ((> 0) . a053603) [1..]
    -- Reinhard Zumkeller, Jun 28 2013
    
  • Maple
    isA051533 := proc(n)
        local a,ta;
        for a from 1 do
            ta := A000217(a) ;
            if 2*ta > n then
                return false;
            end if;
            if isA000217(n-ta) then
                return true;
            end if;
        end do:
    end proc:
    for n from 1 to 200 do
        if isA051533(n) then
            printf("%d,",n) ;
        end if;
    end do: # R. J. Mathar, Dec 16 2015
  • Mathematica
    f[k_] := If[!
       Head[Reduce[m (m + 1) + n (n + 1) == 2 k && 0 < m && 0 < n, {m, n},
           Integers]] === Symbol, k, 0]; DeleteCases[Table[f[k], {k, 1, 108}], 0] (* Ant King, Nov 22 2010 *)
    nn=50; tri=Table[n(n+1)/2, {n,nn}]; Select[Union[Flatten[Table[tri[[i]]+tri[[j]], {i,nn}, {j,i,nn}]]], #<=tri[[-1]] &]
    With[{nn=70},Take[Union[Total/@Tuples[Accumulate[Range[nn]],2]],nn]] (* Harvey P. Dale, Jul 16 2015 *)
  • PARI
    is(n)=for(k=ceil((sqrt(4*n+1)-1)/2),(sqrt(8*n-7)-1)\2, if(ispolygonal(n-k*(k+1)/2, 3), return(1))); 0 \\ Charles R Greathouse IV, Jun 09 2015

Formula

A053603(a(n)) > 0. - Reinhard Zumkeller, Jun 28 2013
A061336(a(n)) = 2. - M. F. Hasler, Mar 06 2017

A068980 Number of partitions of n into nonzero tetrahedral numbers.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 11, 11, 12, 12, 15, 15, 16, 16, 19, 19, 22, 22, 25, 25, 28, 29, 32, 32, 35, 36, 42, 42, 45, 46, 52, 53, 56, 57, 63, 64, 70, 71, 77, 78, 84, 87, 94, 95, 101, 104, 115, 116, 122, 125, 136, 139, 146, 149, 160, 163, 175
Offset: 0

Views

Author

Keywords

Examples

			a(10) = 4 because we can write 10 = 10 = 4 + 4 + 1 + 1 = 4 + 1 + 1 + 1 + 1 + 1 = 1 + ... + 1.
		

Crossrefs

See also A007294 (partitions into triangular numbers), A000292 (tetrahedral numbers).

Programs

  • Mathematica
    nmax = 100; CoefficientList[Series[Product[1/(1-x^(k*(k+1)*(k+2)/6)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jun 09 2017 *)

Formula

G.f.: 1 / Product_{k>=3} (1 - z^binomial(k, 3)).
G.f.: Sum_{i>=0} x^(i*(i+1)*(i+2)/6) / Product_{j=1..i} (1 - x^(j*(j+1)*(j+2)/6)). - Ilya Gutkovskiy, Jun 08 2017

A325352 Heinz number of the differences plus one of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 5, 1, 6, 1, 7, 3, 8, 1, 6, 1, 10, 5, 11, 1, 12, 2, 13, 4, 14, 1, 9, 1, 16, 7, 17, 3, 12, 1, 19, 11, 20, 1, 15, 1, 22, 6, 23, 1, 24, 2, 10, 13, 26, 1, 12, 5, 28, 17, 29, 1, 18, 1, 31, 10, 32, 7, 21, 1, 34, 19, 15, 1, 24, 1, 37, 6, 38
Offset: 1

Views

Author

Gus Wiseman, Apr 23 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The only fixed point is 1 because otherwise the sequence decreases omega (A001222) by one.

Examples

			The partition (3,2,2,1) with Heinz number 90 has differences plus one (2,1,2) with Heinz number 18, so a(90) = 18.
		

Crossrefs

Positions of m's are A008578 (m = 1), A001248 (m = 2), A006094 (m = 3), A030078 (m = 4), A090076 (m = 5).

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    db[n_]:=Times@@Prime/@(1+Differences[primeMS[n]]);
    Table[db[n],{n,100}]

A002636 Number of ways of writing n as an unordered sum of at most 3 nonzero triangular numbers.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 2, 2, 1, 2, 2, 1, 3, 2, 1, 2, 3, 2, 2, 2, 1, 4, 3, 2, 2, 2, 2, 3, 3, 1, 4, 4, 2, 2, 3, 2, 3, 4, 2, 3, 3, 2, 4, 3, 2, 4, 4, 2, 4, 4, 1, 4, 5, 1, 2, 3, 4, 6, 4, 3, 2, 5, 2, 3, 3, 3, 6, 5, 2, 2, 5, 3, 5, 4, 2, 4, 5, 3, 4, 5, 2, 4, 6, 2, 6, 3, 3, 6, 3, 2, 3, 7, 3, 6, 6, 2, 4, 6, 3, 2
Offset: 0

Views

Author

N. J. A. Sloane, Sep 18 2001

Keywords

Comments

Fermat asserted that every number is the sum of three triangular numbers. This was proved by Gauss, who recorded in his Tagebuch entry for Jul 10 1796 that: EYPHKA! num = DELTA + DELTA + DELTA.
a(n) <= A167618(n). - Reinhard Zumkeller, Nov 07 2009
Equivalently, number of ways of writing n as an unordered sum of exactly 3 triangular numbers. - Jon E. Schoenfield, Mar 28 2021

Examples

			0 : empty sum
1 : 1
2 : 1+1
3 : 3 = 1+1+1
4 : 3+1
5 : 3+1+1
6 : 6 = 3+3
7 : 6+1 = 3+3+1
...
13 : 10 + 3 = 6 + 6 + 1, so a(13) = 2.
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 102, eq. (8).
  • D. H. Lehmer, Review of Loria article, Math. Comp. 2 (1947), 301-302.
  • G. Loria, Sulla scomposizione di un intero nella somma di numeri poligonali. (Italian) Atti Accad. Naz. Lincei. Rend. Cl. Sci. Fis. Mat. Nat. (8) 1, (1946). 7-15.
  • Mel Nathanson, Additive Number Theory: The Classical Bases, Graduate Texts in Mathematics, Volume 165, Springer-Verlag, 1996. See Chapter 1.
  • 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).

Crossrefs

Programs

  • Maple
    # reuses code in A000217
    A002636 := proc(n)
        local a,i,Ti, j,Tj, Tk ;
        a := 0 ;
        for i from 0 do
            Ti := A000217(i) ;
            if Ti > n then
                break ;
            end if;
            for j from i do
                Tj := A000217(j) ;
                if Ti+Tj > n then
                    break ;
                end if;
                Tk := n-Ti-Tj ;
                if Tk >= Tj and isA000217(Tk) then
                    a := a+1 ;
                end if;
                if Tk < Tj then
                    break ;
                end if;
            end do:
        end do:
        a ;
    end proc:
    seq(A002636(n),n=0..40) ; # R. J. Mathar, May 26 2025
  • Mathematica
    a = Table[ n(n + 1)/2, {n, 0, 15} ]; b = {0}; c = Table[ 0, {100} ]; Do[ b = Append[ b, a[ [ i ] ] + a[ [ j ] ] + a[ [ k ] ] ], {k, 1, 15}, {j, 1, k}, {i, 1, j} ]; b = Delete[ b, 1 ]; b = Sort[ b ]; l = Length[ b ]; Do[ If[ b[ [ n ] ] < 100, c[ [ b[ [ n ] ] + 1 ] ]++ ], {n, 1, l} ]; c
  • PARI
    first(n)=my(v=vector(n+1),A,B,C); for(a=0,n, A=a*(a+1)/2; if(A>n, break); for(b=0,a, B=A+b*(b+1)/2; if(B>n, break); for(c=0,b, C=B+c*(c+1)/2; if(C>n, break); v[C+1]++))); v \\ Charles R Greathouse IV, Jun 23 2017

Extensions

More terms from Robert G. Wilson v, Sep 20 2001
Entry revised by N. J. A. Sloane, Feb 25 2007

A053614 Numbers that are not the sum of distinct triangular numbers.

Original entry on oeis.org

2, 5, 8, 12, 23, 33
Offset: 1

Views

Author

Jud McCranie, Mar 19 2000

Keywords

Comments

The Mathematica program first computes A024940, the number of partitions of n into distinct triangular numbers. Then it finds those n having zero such partitions. It appears that A024940 grows exponentially, which would preclude additional terms in this sequence. - T. D. Noe, Jul 24 2006, Jan 05 2009

Examples

			a(2) = 5: the 7 partitions of 5 are 5, 4+1, 3+2, 3+1+1, 2+2+1, 2+1+1+1, 1+1+1+1+1. Among those the distinct ones are 5, 4+1, 3+2. None contains all distinct triangular numbers.
12 is a term as it is not a sum of 1, 3, 6 or 10 taken at most once.
		

References

  • Joe Roberts, Lure of the Integers, The Mathematical Association of America, 1992, page 184, entry 33.
  • David Wells in "The Penguin Dictionary of Curious and Interesting Numbers, Revised Edition, page 94, states that "33 is the largest number that is not the sum of distinct triangular numbers".

Crossrefs

Cf. A025524 (number of numbers not the sum of distinct n-th-order polygonal numbers)
Cf. A007419 (largest number not the sum of distinct n-th-order polygonal numbers)
Cf. A001422, A121405 (corresponding sequences for square and pentagonal numbers)

Programs

  • Mathematica
    nn=100; t=Rest[CoefficientList[Series[Product[(1+x^(k*(k+1)/2)), {k,nn}], {x,0,nn(nn+1)/2}], x]]; Flatten[Position[t,0]] (* T. D. Noe, Jul 24 2006 *)

Formula

Complement of A061208.

Extensions

Entry revised by N. J. A. Sloane, Jul 23 2006

A325362 Heinz numbers of integer partitions whose differences (with the last part taken to be 0) are weakly increasing.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 37, 38, 39, 41, 42, 43, 46, 47, 51, 53, 57, 58, 59, 61, 62, 65, 66, 67, 69, 71, 73, 74, 78, 79, 82, 83, 85, 86, 87, 89, 93, 94, 95, 97, 101, 102, 103, 106, 107, 109, 110, 111, 113
Offset: 1

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (x, y, z) are (y - x, z - y). We adhere to this standard for integer partitions also even though they are always weakly decreasing. For example, the differences of (6,3,1) (with the last part taken to be 0) are (-3,-2,-1).
The enumeration of these partitions by sum is given by A007294.
This sequence and A025487, considered as sets, are related by the partition conjugation function A122111(.), which maps the members of either set 1:1 onto the other set. - Peter Munn, Feb 10 2022

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
    2: {1}
    3: {2}
    5: {3}
    6: {1,2}
    7: {4}
   10: {1,3}
   11: {5}
   13: {6}
   14: {1,4}
   17: {7}
   19: {8}
   21: {2,4}
   22: {1,5}
   23: {9}
   26: {1,6}
   29: {10}
   30: {1,2,3}
   31: {11}
   33: {2,5}
		

Crossrefs

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],OrderedQ[Differences[Append[primeptn[#],0]]]&]

A061336 Smallest number of triangular numbers which sum to n.

Original entry on oeis.org

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

Views

Author

Henry Bottomley, Apr 25 2001

Keywords

Comments

a(n)=3 if n=5 or 8 mod 9, since triangular numbers are {0,1,3,6} mod 9.
From Bernard Schott, Jul 16 2022: (Start)
In September 1636, Fermat, in a letter to Mersenne, made the statement that every number is a sum of at most three triangular numbers. This was proved by Gauss, who noted this event in his diary on July 10 1796 with the notation:
EYPHKA! num = DELTA + DELTA + DELTA (where Y is in fact the Greek letter Upsilon and DELTA is the Greek letter of that name).
This proof was published in his book Disquisitiones Arithmeticae, Leipzig, 1801. (End)

Examples

			a(3)=1 since 3=3, a(4)=2 since 4=1+3, a(5)=3 since 5=1+1+3, with 1 and 3 being triangular.
		

References

  • Elena Deza and Michel Marie Deza, Fermat's polygonal number theorem, Figurate numbers, World Scientific Publishing (2012), Chapter 5, pp. 313-377.
  • C. F. Gauss, Disquisitiones Arithmeticae, Yale University Press, 1966, New Haven and London, p. 342, art. 293.

Crossrefs

Cf. A100878 (analog for A000326), A104246 (analog for A000292), A283365 (analog for A000332), A283370 (analog for A000389).

Programs

  • Mathematica
    t[n_]:=n*(n+1)/2; a[0]=0; a[n_]:=Block[ {k=1, tt= t/@ Range[Sqrt[2*n]]}, Off[IntegerPartitions::take]; While[{} == IntegerPartitions[n, {k}, tt, 1], k++]; k]; a/@ Range[0, 104] (* Giovanni Resta, Jun 09 2015 *)
  • PARI
    \\ see A283370 for generic code, working but not optimized for this case of triangular numbers. - M. F. Hasler, Mar 06 2017
    
  • PARI
    a(n)=my(m=n%9,f); if(m==5 || m==8, return(3)); f=factor(4*n+1); for(i=1,#f~, if(f[i,2]%2 && f[i,1]%4==3, return(3))); if(ispolygonal(n,3), n>0, 2) \\ Charles R Greathouse IV, Mar 17 2022

Formula

a(n) = 0 if n=0, otherwise 1 if n is in A000217, otherwise 2 if n is in A051533, otherwise 3 in which case n is in A020757.
a(n) <= 3 (proposed by Fermat and proved by Gauss). - Bernard Schott, Jul 16 2022

A325360 Heinz numbers of integer partitions whose differences are weakly increasing.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 71, 73
Offset: 1

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (x, y, z) are (y - x, z - y). We adhere to this standard for integer partitions also even though they are always weakly decreasing. For example, the differences of (6,3,1) are (-3,-2).
The enumeration of these partitions by sum is given by A240026.

Examples

			Most small numbers are in the sequence. However, the sequence of non-terms together with their prime indices begins:
   18: {1,2,2}
   36: {1,1,2,2}
   50: {1,3,3}
   54: {1,2,2,2}
   70: {1,3,4}
   72: {1,1,1,2,2}
   75: {2,3,3}
   90: {1,2,2,3}
   98: {1,4,4}
  100: {1,1,3,3}
		

Crossrefs

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],OrderedQ[Differences[primeptn[#]]]&]

A053604 Number of ways to write n as an ordered sum of 3 nonzero triangular numbers.

Original entry on oeis.org

0, 0, 0, 1, 0, 3, 0, 3, 3, 1, 6, 0, 6, 3, 6, 3, 3, 9, 1, 12, 0, 6, 9, 6, 6, 6, 9, 6, 12, 0, 10, 9, 12, 6, 9, 9, 3, 18, 3, 12, 12, 9, 9, 9, 12, 10, 12, 9, 9, 18, 6, 6, 27, 6, 12, 6, 9, 18, 15, 15, 6, 21, 9, 13, 12, 9, 18, 21, 9, 6, 21, 15, 15, 15, 12, 15, 18, 15, 9
Offset: 0

Views

Author

N. J. A. Sloane, Jan 20 2000

Keywords

Comments

Fermat asserted that every number is the sum of three triangular numbers. This was proved by Gauss, who recorded in his Tagebuch entry for Jul 10 1796 that: EYPHEKA! num = DELTA + DELTA + DELTA.

References

  • Mel Nathanson, Additive Number Theory: The Classical Bases, Graduate Texts in Mathematics, Volume 165, Springer-Verlag, 1996. See Chapter 1.

Crossrefs

Programs

  • Mathematica
    nmax = 100; m0 = 10; A053604 :=
    Table[a[n], {n, 0, nmax}]; Clear[counts];
    counts[m_] :=
    counts[m] = (Clear[a]; a[_] = 0;
       Do[s = i*(i + 1)/2 + j*(j + 1)/2 + k*(k + 1)/2;
        a[s] = a[s] + 1, {i, 1, m}, {j, 1, m}, {k, 1, m}];
       A053603); counts[m = m0]; counts[m = 2*m]; While[
    counts[m] != counts[m/2], m = 2*m]; A053604  (* G. C. Greubel, Dec 24 2016 *)

Formula

G.f.: ( Sum_{k>=1} x^(k*(k+1)/2) )^3. - Ilya Gutkovskiy, Dec 24 2016
Previous Showing 11-20 of 115 results. Next