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

A005044 Alcuin's sequence: expansion of x^3/((1-x^2)*(1-x^3)*(1-x^4)).

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 1, 2, 1, 3, 2, 4, 3, 5, 4, 7, 5, 8, 7, 10, 8, 12, 10, 14, 12, 16, 14, 19, 16, 21, 19, 24, 21, 27, 24, 30, 27, 33, 30, 37, 33, 40, 37, 44, 40, 48, 44, 52, 48, 56, 52, 61, 56, 65, 61, 70, 65, 75, 70, 80, 75, 85, 80, 91, 85, 96, 91, 102, 96, 108, 102, 114, 108, 120
Offset: 0

Views

Author

Keywords

Comments

a(n) is the number of triangles with integer sides and perimeter n.
Also a(n) is the number of triangles with distinct integer sides and perimeter n+6, i.e., number of triples (a, b, c) such that 1 < a < b < c < a+b, a+b+c = n+6. - Roger Cuculière
With a different offset (i.e., without the three leading zeros, as in A266755), the number of ways in which n empty casks, n casks half-full of wine and n full casks can be distributed to 3 persons in such a way that each one gets the same number of casks and the same amount of wine [Alcuin]. E.g., for n=2 one can give 2 people one full and one empty and the 3rd gets two half-full. (Comment corrected by Franklin T. Adams-Watters, Oct 23 2006)
For m >= 2, the sequence {a(n) mod m} is periodic with period 12*m. - Martin J. Erickson (erickson(AT)truman.edu), Jun 06 2008
Number of partitions of n into parts 2, 3, and 4, with at least one part 3. - Joerg Arndt, Feb 03 2013
For several values of p and q the sequence (A005044(n+p) - A005044(n-q)) leads to known sequences, see the crossrefs. - Johannes W. Meijer, Oct 12 2013
For n>=3, number of partitions of n-3 into parts 2, 3, and 4. - David Neil McGrath, Aug 30 2014
Also, a(n) is the number of partitions mu of n of length 3 such that mu_1-mu_2 is even and mu_2-mu_3 is even (see below example). - John M. Campbell, Jan 29 2016
For n > 1, number of triangles with odd side lengths and perimeter 2*n-3. - Wesley Ivan Hurt, May 13 2019
Number of partitions of n+1 into 4 parts whose largest two parts are equal. - Wesley Ivan Hurt, Jan 06 2021
For n>=3, number of weak partitions of n-3 (that is, allowing parts of size 0) into three parts with no part exceeding (n-3)/2. Also, number of weak partitions of n-3 into three parts, all of the same parity as n-3. - Kevin Long, Feb 20 2021
Also, a(n) is the number of incongruent acute triangles formed from the vertices of a regular n-gon. - Frank M Jackson, Nov 04 2022

Examples

			There are 4 triangles of perimeter 11, with sides 1,5,5; 2,4,5; 3,3,5; 3,4,4. So a(11) = 4.
G.f. = x^3 + x^5 + x^6 + 2*x^7 + x^8 + 3*x^9 + 2*x^10 + 4*x^11 + 3*x^12 + ...
From _John M. Campbell_, Jan 29 2016: (Start)
Letting n = 15, there are a(n)=7 partitions mu |- 15 of length 3 such that mu_1-mu_2 is even and mu_2-mu_3 is even:
(13,1,1) |- 15
(11,3,1) |- 15
(9,5,1) |- 15
(9,3,3) |- 15
(7,7,1) |- 15
(7,5,3) |- 15
(5,5,5) |- 15
(End)
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 74, Problem 7.
  • I. Niven and H. S. Zuckerman, An Introduction to the Theory of Numbers. Wiley, NY, Chap.10, Section 10.2, Problems 5 and 6, pp. 451-2.
  • D. Olivastro: Ancient Puzzles. Classic Brainteasers and Other Timeless Mathematical Games of the Last 10 Centuries. New York: Bantam Books, 1993. See p. 158.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • A. M. Yaglom and I. M. Yaglom: Challenging Mathematical Problems with Elementary Solutions. Vol. I. Combinatorial Analysis and Probability Theory. New York: Dover Publications, Inc., 1987, p. 8, #30 (First published: San Francisco: Holden-Day, Inc., 1964)

Crossrefs

See A266755 for a version without the three leading zeros.
Both bisections give (essentially) A001399.
(See the comments.) Cf. A008615 (p=1, q=3, offset=0), A008624 (3, 3, 0), A008679 (3, -1, 0), A026922 (1, 5, 1), A028242 (5, 7, 0), A030451 (6, 6, 0), A051274 (3, 5, 0), A052938 (8, 4, 0), A059169 (0, 6, 1), A106466 (5, 4, 0), A130722 (2, 7, 0)
Cf. this sequence (k=3), A288165 (k=4), A288166 (k=5).
Number of k-gons that can be formed with perimeter n: this sequence (k=3), A062890 (k=4), A069906 (k=5), A069907 (k=6), A288253 (k=7), A288254 (k=8), A288255 (k=9), A288256 (k=10).

Programs

  • Haskell
    a005044 = p [2,3,4] . (subtract 3) where
    p _ 0 = 1
    p [] _ = 0
    p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Feb 28 2013
  • Maple
    A005044 := n-> floor((1/48)*(n^2+3*n+21+(-1)^(n-1)*3*n)): seq(A005044(n), n=0..73);
    A005044 := -1/(z**2+1)/(z**2+z+1)/(z+1)**2/(z-1)**3; # Simon Plouffe in his 1992 dissertation
  • Mathematica
    a[n_] := Round[If[EvenQ[n], n^2, (n + 3)^2]/48] (* Peter Bertok, Jan 09 2002 *)
    CoefficientList[Series[x^3/((1 - x^2)*(1 - x^3)*(1 - x^4)), {x, 0, 105}], x] (* Robert G. Wilson v, Jun 02 2004 *)
    me[n_] := Module[{i, j, sum = 0}, For[i = Ceiling[(n - 3)/3], i <= Floor[(n - 3)/2], i = i + 1, For[j = Ceiling[(n - i - 3)/2], j <= i, j = j + 1, sum = sum + 1] ]; Return[sum]; ] mine = Table[me[n], {n, 1, 11}]; (* Srikanth (sriperso(AT)gmail.com), Aug 02 2008 *)
    LinearRecurrence[{0,1,1,1,-1,-1,-1,0,1},{0,0,0,1,0,1,1,2,1},80] (* Harvey P. Dale, Sep 22 2014 *)
    Table[Length@Select[IntegerPartitions[n, {3}], Max[#]*180 < 90 n &], {n, 1, 100}] (* Frank M Jackson, Nov 04 2022 *)
  • PARI
    a(n) = round(n^2 / 12) - (n\2)^2 \ 4
    
  • PARI
    a(n) = (n^2 + 6*n * (n%2) + 24) \ 48
    
  • PARI
    a(n)=if(n%2,n+3,n)^2\/48 \\ Charles R Greathouse IV, May 02 2016
    
  • PARI
    concat(vector(3), Vec((x^3)/((1-x^2)*(1-x^3)*(1-x^4)) + O(x^70))) \\ Felix Fröhlich, Jun 07 2017
    

Formula

a(n) = a(n-6) + A059169(n) = A070093(n) + A070101(n) + A024155(n).
For odd indices we have a(2*n-3) = a(2*n). For even indices, a(2*n) = nearest integer to n^2/12 = A001399(n).
For all n, a(n) = round(n^2/12) - floor(n/4)*floor((n+2)/4) = a(-3-n) = A069905(n) - A002265(n)*A002265(n+2).
For n = 0..11 (mod 12), a(n) is respectively n^2/48, (n^2 + 6*n - 7)/48, (n^2 - 4)/48, (n^2 + 6*n + 21)/48, (n^2 - 16)/48, (n^2 + 6*n - 7)/48, (n^2 + 12)/48, (n^2 + 6*n + 5)/48, (n^2 - 16)/48, (n^2 + 6*n + 9)/48, (n^2 - 4)/48, (n^2 + 6*n + 5)/48.
Euler transform of length 4 sequence [ 0, 1, 1, 1]. - Michael Somos, Sep 04 2006
a(-3 - n) = a(n). - Michael Somos, Sep 04 2006
a(n) = sum(ceiling((n-3)/3) <= i <= floor((n-3)/2), sum(ceiling((n-i-3)/2) <= j <= i, 1 ) ) for n >= 1. - Srikanth K S, Aug 02 2008
a(n) = a(n-2) + a(n-3) + a(n-4) - a(n-5) - a(n-6) - a(n-7) + a(n-9) for n >= 9. - David Neil McGrath, Aug 30 2014
a(n+3) = a(n) if n is odd; a(n+3) = a(n) + floor(n/4) + 1 if n is even. Sketch of proof: There is an obvious injective map from perimeter-n triangles to perimeter-(n+3) triangles defined by f(a,b,c) = (a+1,b+1,c+1). It is easy to show f is surjective for odd n, while for n=2k the image of f is only missing the triangles (a,k+2-a,k+1) for 1 <= a <= floor(k/2)+1. - James East, May 01 2016
a(n) = round(n^2/48) if n is even; a(n) = round((n+3)^2/48) if n is odd. - James East, May 01 2016
a(n) = (6*n^2 + 18*n - 9*(-1)^n*(2*n + 3) - 36*sin(Pi*n/2) - 36*cos(Pi*n/2) + 64*cos(2*Pi*n/3) - 1)/288. - Ilya Gutkovskiy, May 01 2016
a(n) = A325691(n-3) + A000035(n) for n>=3. The bijection between partition(n,[2,3,4]) and not-over-half partition(n,3,n/2) + partition(n,2,n/2) can be built by a Ferrers(part)[0+3,1,2] map. And the last partition(n,2,n/2) is unique [n/2,n/2] if n is even, it is given by A000035. - Yuchun Ji, Sep 24 2020
a(4n+3) = a(4n) + n+1, a(4n+4) = a(4n+1) = A000212(n+1), a(4n+5) = a(4n+2) + n+1, a(4n+6) = a(4n+3) = A007980(n). - Yuchun Ji, Oct 10 2020
a(n)-a(n-4) = A008615(n-1). - R. J. Mathar, Jun 23 2021
a(n)-a(n-2) = A008679(n-3). - R. J. Mathar, Jun 23 2021

Extensions

Additional comments from Reinhard Zumkeller, May 11 2002
Yaglom reference and mod formulas from Antreas P. Hatzipolakis (xpolakis(AT)otenet.gr), May 27 2000
The reference to Alcuin of York (735-804) was provided by Hermann Kremer (hermann.kremer(AT)onlinehome.de), Jun 18 2004

A266755 Expansion of 1/((1-x^2)*(1-x^3)*(1-x^4)).

Original entry on oeis.org

1, 0, 1, 1, 2, 1, 3, 2, 4, 3, 5, 4, 7, 5, 8, 7, 10, 8, 12, 10, 14, 12, 16, 14, 19, 16, 21, 19, 24, 21, 27, 24, 30, 27, 33, 30, 37, 33, 40, 37, 44, 40, 48, 44, 52, 48, 56, 52, 61, 56, 65, 61, 70, 65, 75, 70, 80, 75, 85, 80, 91, 85, 96, 91, 102, 96, 108, 102, 114, 108, 120, 114, 127, 120, 133, 127, 140, 133, 147, 140, 154, 147, 161, 154, 169
Offset: 0

Views

Author

N. J. A. Sloane, Jan 10 2016

Keywords

Comments

This is the same as A005044 but without the three leading zeros. There are so many situations where one wants this sequence rather than A005044 that it seems appropriate for it to have its own entry.
But see A005044 (still the main entry) for numerous applications and references.
Also, Molien series for invariants of finite Coxeter group D_3.
The Molien series for the finite Coxeter group of type D_k (k >= 3) has g.f. = 1/Product_i (1-x^(1+m_i)) where the m_i are [1,3,5,...,2k-3,k-1]. If k is even only even powers of x appear, and we bisect the sequence.
Also, Molien series for invariants of finite Coxeter group A_3. The Molien series for the finite Coxeter group of type A_k (k >= 1) has g.f. = 1/Product_{i=2..k+1} (1-x^i). Note that this is the root system A_k not the alternating group Alt_k.
a(n) is the number of partitions of n into parts 2, 3, and 4. - Joerg Arndt, Apr 16 2017
From Gus Wiseman, May 23 2021: (Start)
Also the number of integer partitions of n into at most n/2 parts, none greater than 3. The case of any maximum is A110618. The case of any length is A001399. The Heinz numbers of these partitions are given by A344293.
For example, the a(2) = 1 through a(13) = 5 partitions are:
2 3 22 32 33 322 332 333 3322 3332 3333 33322
31 222 331 2222 3222 3331 32222 33222 33331
321 3221 3321 22222 33221 33321 322222
3311 32221 33311 222222 332221
33211 322221 333211
332211
333111
(End)

Examples

			G.f. = 1 + x^2 + x^3 + 2*x^4 + x^5 + 3*x^6 + 2*x^7 + 4*x^8 + ... - _Michael Somos_, Jan 29 2022
		

References

  • J. E. Humphreys, Reflection Groups and Coxeter Groups, Cambridge, 1990. See Table 3.1, page 59.

Crossrefs

Molien series for finite Coxeter groups A_1 through A_12 are A059841, A103221, A266755, A008667, A037145, A001996, and A266776-A266781.
Molien series for finite Coxeter groups D_3 through D_12 are A266755, A266769, A266768, A003402, and A266770-A266775.
A variant of A005044.
Cf. A001400 (partial sums).
Cf. A308065.
Number of partitions of n whose Heinz number is in A344293.
A001399 counts partitions with all parts <= 3, ranked by A051037.
A025065 counts partitions of n with >= n/2 parts, ranked by A344296.
A035363 counts partitions of n with n/2 parts, ranked by A340387.
A110618 counts partitions of n into at most n/2 parts, ranked by A344291.

Programs

  • Magma
    I:=[1,0,1,1,2,1,3,2,4]; [n le 9 select I[n] else Self(n-2)+ Self(n-3)+Self(n-4)-Self(n-5)-Self(n-6)-Self(n-7)+Self(n-9): n in [1..100]]; // Vincenzo Librandi, Jan 11 2016
    
  • Mathematica
    CoefficientList[Series[1/((1-x^2)(1-x^3)(1-x^4)), {x, 0, 100}], x] (* JungHwan Min, Jan 10 2016 *)
    LinearRecurrence[{0,1,1,1,-1,-1,-1,0,1}, {1,0,1,1,2,1,3,2,4}, 100] (* Vincenzo Librandi, Jan 11 2016 *)
    Table[Length[Select[IntegerPartitions[n],Length[#]<=n/2&&Max@@#<=3&]],{n,0,30}] (* Gus Wiseman, May 23 2021 *)
    a[ n_] := Round[(n + 3*(2 - Mod[n,2]))^2/48]; (* Michael Somos, Jan 29 2022 *)
  • PARI
    Vec(1/((1-x^2)*(1-x^3)*(1-x^4)) + O(x^100)) \\ Michel Marcus, Jan 11 2016
    
  • PARI
    {a(n) = round((n + 3*(2-n%2))^2/48)}; /* Michael Somos, Jan 29 2022 */
    
  • Sage
    (1/((1-x^2)*(1-x^3)*(1-x^4))).series(x, 100).coefficients(x, sparse=False) # G. C. Greubel, Jun 13 2019

Formula

a(n) = a(n-2) + a(n-3) + a(n-4) - a(n-5) - a(n-6) - a(n-7) + a(n-9) for n>8. - Vincenzo Librandi, Jan 11 2016
a(n) = a(-9-n) for all n in Z. a(n) = a(n+3) for all n in 2Z. - Michael Somos, Jan 29 2022
E.g.f.: exp(-x)*(81 - 18*x + exp(2*x)*(107 + 60*x + 6*x^2) + 64*exp(x/2)*cos(sqrt(3)*x/2) + 36*exp(x)*(cos(x) - sin(x)))/288. - Stefano Spezia, Mar 05 2023
For n >= 3, if n is even, a(n) = a(n-3) + floor(n/4) + 1, otherwise a(n) = a(n-3). - Robert FERREOL, Feb 05 2024
a(n) = floor((n^2+9*n+(3*n+9)*(-1)^n+39)/48). - Hoang Xuan Thanh, Jun 03 2025

A344296 Numbers with at least as many prime factors (counted with multiplicity) as half their sum of prime indices.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 27, 28, 30, 32, 36, 40, 48, 54, 56, 60, 64, 72, 80, 81, 84, 88, 90, 96, 100, 108, 112, 120, 128, 144, 160, 162, 168, 176, 180, 192, 200, 208, 216, 224, 240, 243, 252, 256, 264, 270, 280, 288, 300, 320, 324, 336, 352
Offset: 1

Views

Author

Gus Wiseman, May 16 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.
These are the Heinz numbers of certain partitions counted by A025065, but different from palindromic partitions, which have Heinz numbers A265640.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}            30: {1,2,3}
      2: {1}           32: {1,1,1,1,1}
      3: {2}           36: {1,1,2,2}
      4: {1,1}         40: {1,1,1,3}
      6: {1,2}         48: {1,1,1,1,2}
      8: {1,1,1}       54: {1,2,2,2}
      9: {2,2}         56: {1,1,1,4}
     10: {1,3}         60: {1,1,2,3}
     12: {1,1,2}       64: {1,1,1,1,1,1}
     16: {1,1,1,1}     72: {1,1,1,2,2}
     18: {1,2,2}       80: {1,1,1,1,3}
     20: {1,1,3}       81: {2,2,2,2}
     24: {1,1,1,2}     84: {1,1,2,4}
     27: {2,2,2}       88: {1,1,1,5}
     28: {1,1,4}       90: {1,2,2,3}
		

Crossrefs

The case with difference at least 1 is A322136.
The case of equality is A340387, counted by A000041 or A035363.
The opposite version is A344291, counted by A110618.
The conjugate version is A344414, with even-weight case A344416.
A025065 counts palindromic partitions, ranked by A265640.
A056239 adds up prime indices, row sums of A112798.
A300061 lists numbers whose sum of prime indices is even.

Programs

  • Mathematica
    Select[Range[100],PrimeOmega[#]>=Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]/2&]

Formula

A056239(a(n)) <= 2*A001222(a(n)).
a(n) = A322136(n)/4.

A008642 Quarter-squares repeated.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 6, 6, 9, 9, 12, 12, 16, 16, 20, 20, 25, 25, 30, 30, 36, 36, 42, 42, 49, 49, 56, 56, 64, 64, 72, 72, 81, 81, 90, 90, 100, 100, 110, 110, 121, 121, 132, 132, 144, 144, 156, 156, 169, 169, 182, 182, 196, 196, 210, 210, 225, 225
Offset: 0

Views

Author

Keywords

Comments

The area of the largest rectangle whose perimeter is not greater than n. - Dmitry Kamenetsky, Aug 30 2006
Also number of partitions of n into parts 1, 2 or 4. - Reinhard Zumkeller, Aug 12 2011
Let us consider a rectangle composed of unit squares. Then count how many squares are necessary to surround this rectangle by a layer whose width is 1 unit. And repeat this surrounding ad libitum. This sequence, prepended by 4 zeros and with offset 0, gives the number of rectangles that need 2*n unit squares in one of their surrounding layers. - Michel Marcus, Sep 19 2015
a(n) is the number of nonnegative integer solutions (x,y,z) for n-2 <= 2*x + 3*y + 4*z <= n. For example, the two solutions for 1 <= 2*x + 3*y + 4*z <= 3 are (1,0,0) and (0,1,0). - Ran Pan, Oct 07 2015
Conjecture: Consider the number of compositions of n>=4*k+8 into odd parts, where the order of the parts 1,3,..,2k+1 does not count. Then, as k approaches infinity, a(n-4*k-8) is equal to the number of these restricted compositions minus A000009(n), the number of strict partitions of n. - Gregory L. Simay, Aug 12 2016
From Gus Wiseman, May 17 2019: (Start)
Also the number of length-3 integer partitions of n + 4 whose largest part is greater than the sum of the other two. These are unordered triples that cannot be the sides of a triangle. For example, the a(1) = 1 through a(10) = 9 partitions are (A = 10, B = 11, C = 12):
(311) (411) (421) (521) (522) (622) (632) (732) (733) (833)
(511) (611) (531) (631) (641) (741) (742) (842)
(621) (721) (722) (822) (751) (851)
(711) (811) (731) (831) (832) (932)
(821) (921) (841) (941)
(911) (A11) (922) (A22)
(931) (A31)
(A21) (B21)
(B11) (C11)
(End)
This sequence, prepended by four 0's and with offset 0, is the number of partitions of n into four parts whose smallest two parts are equal. - Wesley Ivan Hurt, Jan 05 2021
This sequence, prepended by four 0's and with offset 0, is the number of incongruent obtuse triangles formed from the vertices of a regular n-gon. - Frank M Jackson, Nov 27 2022

References

  • D. J. Benson, Polynomial Invariants of Finite Groups, Cambridge, 1993, p. 105.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 112, D(n).

Crossrefs

Cf. A002620.
Cf. A001399, A005044 (triangles without self-intersections), A069905, A124278, A266223, A325686, A325689, A325690, A325691, A325695.

Programs

  • Magma
    [Floor(((n+1)*((-1)^n+n+6)+9)/16): n in [0..70]]; // Vincenzo Librandi, Apr 02 2014
    
  • Maple
    seq((7/8+(-1)^k/8 + k + k^2/4)$2, k=0..100); # Robert Israel, Oct 08 2015
  • Mathematica
    CoefficientList[Series[1/((1-x)(1-x^2)(1-x^4)), {x, 0, 70}], x] (* Vincenzo Librandi, Apr 02 2014 *)
    LinearRecurrence[{1,1,-1,1,-1,-1,1},{1,1,2,2,4,4,6}, 70] (* Harvey P. Dale, Jun 03 2015 *)
    Table[Floor[((n + 1) ((-1)^n + n + 6) + 9)/16], {n, 0, 70}] (* Michael De Vlieger, Aug 14 2016 *)
  • PARI
    Vec(1/((1-x)*(1-x^2)*(1-x^4)) + O(x^70)) \\ Michel Marcus, Mar 31 2014
    
  • PARI
    vector(70, n, n--; floor(((n+1)*((-1)^n+n+6)+9)/16)) \\ Altug Alkan, Oct 08 2015
    
  • Sage
    [floor(floor(n/2+2)^2/2)/2 for n in (0..70)] # Bruno Berselli, Mar 03 2016

Formula

G.f.: 1/((1-x)*(1-x^2)*(1-x^4)).
a(n) = (2*n^2 + 14*n + 21 + (2*n + 7)*(-1)^n)/32 + ((1 + (-1)^n)/2 - (1 - (-1)^n)*i/2)*i^n/8, with i = sqrt(-1).
a(n) = floor(((n+1)*((-1)^n+n+6)+9)/16). - Tani Akinari, Jun 16 2013
a(n) = Sum_{i=1..floor((n+6)/2)} floor((n+6-2*i-(n mod 2))/4). - Wesley Ivan Hurt, Mar 31 2014
a(0)=1, a(1)=1, a(2)=2, a(3)=2, a(4)=4, a(5)=4, a(6)=6; for n>6, a(n) = a(n-1) + a(n-2) - a(n-3) + a(n-4) - a(n-5) - a(n-6) + a(n-7). - Harvey P. Dale, Jun 03 2015
a(n) = floor(floor(n/2+2)^2/4) = floor(floor(n/2+2)^2/2)/2. - Bruno Berselli, Mar 03 2016
E.g.f.: ((14 + 7*x + x^2)*cosh(x) + 2*(cos(x) + sin(x)) + (7 + 9*x + x^2)*sinh(x))/16. - Stefano Spezia, Mar 05 2023
a(n) = floor((n + 4)/4)*floor((n + 6)/4). - Ridouane Oudra, Apr 01 2023

A344293 5-smooth numbers n whose sum of prime indices A056239(n) is at least twice the number of prime indices A001222(n).

Original entry on oeis.org

1, 3, 5, 9, 10, 15, 25, 27, 30, 45, 50, 75, 81, 90, 100, 125, 135, 150, 225, 243, 250, 270, 300, 375, 405, 450, 500, 625, 675, 729, 750, 810, 900, 1000, 1125, 1215, 1250, 1350, 1500, 1875, 2025, 2187, 2250, 2430, 2500, 2700, 3000, 3125, 3375, 3645, 3750, 4050
Offset: 1

Views

Author

Gus Wiseman, May 16 2021

Keywords

Comments

A number is 5-smooth if its prime divisors are all <= 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.

Examples

			The sequence of terms together with their prime indices begins:
       1: {}            125: {3,3,3}
       3: {2}           135: {2,2,2,3}
       5: {3}           150: {1,2,3,3}
       9: {2,2}         225: {2,2,3,3}
      10: {1,3}         243: {2,2,2,2,2}
      15: {2,3}         250: {1,3,3,3}
      25: {3,3}         270: {1,2,2,2,3}
      27: {2,2,2}       300: {1,1,2,3,3}
      30: {1,2,3}       375: {2,3,3,3}
      45: {2,2,3}       405: {2,2,2,2,3}
      50: {1,3,3}       450: {1,2,2,3,3}
      75: {2,3,3}       500: {1,1,3,3,3}
      81: {2,2,2,2}     625: {3,3,3,3}
      90: {1,2,2,3}     675: {2,2,2,3,3}
     100: {1,1,3,3}     729: {2,2,2,2,2,2}
		

Crossrefs

Allowing any number of parts and sum gives A051037, counted by A001399.
These are Heinz numbers of the partitions counted by A266755.
Allowing parts > 5 gives A344291, counted by A110618.
The non-3-smooth case is A344294, counted by A325691.
Requiring the sum of prime indices to be even gives A344295.
A000070 counts non-multigraphical partitions, ranked by A344292.
A025065 counts partitions of n with >= n/2 parts, ranked by A344296.
A035363 counts partitions of n with n/2 parts, ranked by A340387.
A056239 adds up prime indices, row sums of A112798.
A300061 ranks partitions of even numbers, with 5-smooth case A344297.

Programs

  • Mathematica
    Select[Range[1000],PrimeOmega[#]<=Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]/2&&Max@@First/@FactorInteger[#]<=5&]

Formula

Intersection of A051037 and A344291.

A344294 5-smooth but not 3-smooth numbers k such that A056239(k) >= 2*A001222(k).

Original entry on oeis.org

5, 10, 15, 25, 30, 45, 50, 75, 90, 100, 125, 135, 150, 225, 250, 270, 300, 375, 405, 450, 500, 625, 675, 750, 810, 900, 1000, 1125, 1215, 1250, 1350, 1500, 1875, 2025, 2250, 2430, 2500, 2700, 3000, 3125, 3375, 3645, 3750, 4050, 4500, 5000, 5625, 6075, 6250
Offset: 1

Views

Author

Gus Wiseman, May 16 2021

Keywords

Comments

A number is d-smooth iff its prime divisors are all <= d.
A prime index of k is a number m such that prime(m) divides k, and the multiset of prime indices of k is row k of A112798. This row has length A001222(k) and sum A056239(k).

Examples

			The sequence of terms together with their prime indices begins:
       5: {3}           270: {1,2,2,2,3}
      10: {1,3}         300: {1,1,2,3,3}
      15: {2,3}         375: {2,3,3,3}
      25: {3,3}         405: {2,2,2,2,3}
      30: {1,2,3}       450: {1,2,2,3,3}
      45: {2,2,3}       500: {1,1,3,3,3}
      50: {1,3,3}       625: {3,3,3,3}
      75: {2,3,3}       675: {2,2,2,3,3}
      90: {1,2,2,3}     750: {1,2,3,3,3}
     100: {1,1,3,3}     810: {1,2,2,2,2,3}
     125: {3,3,3}       900: {1,1,2,2,3,3}
     135: {2,2,2,3}    1000: {1,1,1,3,3,3}
     150: {1,2,3,3}    1125: {2,2,3,3,3}
     225: {2,2,3,3}    1215: {2,2,2,2,2,3}
     250: {1,3,3,3}    1250: {1,3,3,3,3}
		

Crossrefs

Allowing any number of parts and sum gives A080193, counted by A069905.
The partitions with these Heinz numbers are counted by A325691.
Relaxing the smoothness conditions gives A344291, counted by A110618.
Allowing 3-smoothness gives A344293, counted by A266755.
A025065 counts partitions of n with at least n/2 parts, ranked by A344296.
A035363 counts partitions of n whose length is n/2, ranked by A340387.
A051037 lists 5-smooth numbers (complement: A279622).
A056239 adds up prime indices, row sums of A112798.
A257993 gives the least gap of the partition with Heinz number n.
A300061 lists numbers with even sum of prime indices (5-smooth: A344297).
A342050/A342051 list Heinz numbers of partitions with even/odd least gap.

Programs

  • Mathematica
    Select[Range[1000],PrimeOmega[#]<=Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]/2&&Max@@First/@FactorInteger[#]==5&]

Formula

Intersection of A080193 and A344291.

A325688 Number of length-3 compositions of n such that every distinct consecutive subsequence has a different sum.

Original entry on oeis.org

0, 0, 0, 1, 0, 4, 5, 12, 12, 25, 24, 40, 41, 60, 60, 85, 84, 112, 113, 144, 144, 181, 180, 220, 221, 264, 264, 313, 312, 364, 365, 420, 420, 481, 480, 544, 545, 612, 612, 685, 684, 760, 761, 840, 840, 925, 924, 1012, 1013, 1104, 1104, 1201, 1200, 1300, 1301, 1404
Offset: 0

Views

Author

Gus Wiseman, May 15 2019

Keywords

Comments

A composition of n is a finite sequence of positive integers summing to n.
Confirmed recurrence relation from Colin Barker for n <= 5000. - Fausto A. C. Cariboni, Feb 13 2022

Examples

			The a(3) = 1 through a(8) = 12 compositions:
  (111)  (113)  (114)  (115)  (116)
         (122)  (132)  (124)  (125)
         (221)  (222)  (133)  (143)
         (311)  (231)  (142)  (152)
                (411)  (214)  (215)
                       (223)  (233)
                       (241)  (251)
                       (322)  (332)
                       (331)  (341)
                       (412)  (512)
                       (421)  (521)
                       (511)  (611)
		

Crossrefs

Column k = 3 of A325687.
Cf. A000217 (all length-3).

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n,{3}],UnsameQ@@Total/@Union[ReplaceList[#,{_,s__,_}:>{s}]]&]],{n,0,30}]

Formula

Conjectures from Colin Barker, May 16 2019: (Start)
G.f.: x^3*(1 + 2*x^2 + 4*x^3 + 5*x^4) / ((1 - x)^3*(1 + x)^2*(1 + x + x^2)).
a(n) = 2*a(n-2) + a(n-3) - a(n-4) - 2*a(n-5) + a(n-7) for n>7.
(End)

A325690 Number of length-3 integer partitions of n whose largest part is not the sum of the other two.

Original entry on oeis.org

0, 0, 0, 1, 0, 2, 2, 4, 3, 7, 6, 10, 9, 14, 13, 19, 17, 24, 23, 30, 28, 37, 35, 44, 42, 52, 50, 61, 58, 70, 68, 80, 77, 91, 88, 102, 99, 114, 111, 127, 123, 140, 137, 154, 150, 169, 165, 184, 180, 200, 196, 217, 212, 234, 230, 252, 247, 271, 266, 290, 285, 310
Offset: 0

Views

Author

Gus Wiseman, May 15 2019

Keywords

Comments

Confirmed recurrence relation from Colin Barker for n <= 10000. - Fausto A. C. Cariboni, Feb 19 2022

Examples

			The a(3) = 1 through a(13) = 14 partitions (A = 10, B = 11):
  (111)  (221)  (222)  (322)  (332)  (333)  (433)  (443)  (444)   (544)
         (311)  (411)  (331)  (521)  (432)  (442)  (533)  (543)   (553)
                       (421)  (611)  (441)  (622)  (542)  (552)   (643)
                       (511)         (522)  (631)  (551)  (732)   (652)
                                     (531)  (721)  (632)  (741)   (661)
                                     (621)  (811)  (641)  (822)   (733)
                                     (711)         (722)  (831)   (742)
                                                   (731)  (921)   (751)
                                                   (821)  (A11)   (832)
                                                   (911)          (841)
                                                                  (922)
                                                                  (931)
                                                                  (A21)
                                                                  (B11)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n,{3}],#[[1]]!=#[[2]]+#[[3]]&]],{n,0,30}]

Formula

Conjectures from Colin Barker, May 15 2019: (Start)
G.f.: x^3*(1 + x^2 + x^3 + x^4) / ((1 - x)^3*(1 + x)^2*(1 + x^2)*(1 + x + x^2)).
a(n) = a(n-2) + a(n-3) + a(n-4) - a(n-5) - a(n-6) - a(n-7) + a(n-9) for n>8.
(End)

A325686 Number of strict length-3 compositions x + y + z = n satisfying x + y != z and x != y + z.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 2, 6, 8, 18, 16, 30, 34, 48, 48, 72, 72, 96, 98, 126, 128, 162, 160, 198, 202, 240, 240, 288, 288, 336, 338, 390, 392, 450, 448, 510, 514, 576, 576, 648, 648, 720, 722, 798, 800, 882, 880, 966, 970, 1056, 1056, 1152, 1152, 1248, 1250, 1350, 1352
Offset: 0

Views

Author

Gus Wiseman, May 13 2019

Keywords

Comments

A strict composition of n is a finite sequence of distinct positive integers summing to n.
From Kevin O'Bryant, Jun 02 2025: (Start)
Also the number of Sidon sets in {0,1,...,n} with 4 elements that contain both 0 and n.
Also, the number of 3-tuples of positive integers with the 6 numbers x, y, z, x+y, y+z, x+y+z=n all distinct. (End)

Examples

			The a(6) = 2 through a(10) = 16 compositions:
  (132)  (124)  (125)  (126)  (127)
  (231)  (142)  (143)  (135)  (136)
         (214)  (152)  (153)  (154)
         (241)  (215)  (162)  (163)
         (412)  (251)  (216)  (172)
         (421)  (341)  (234)  (217)
                (512)  (243)  (253)
                (521)  (261)  (271)
                       (315)  (316)
                       (324)  (352)
                       (342)  (361)
                       (351)  (451)
                       (423)  (613)
                       (432)  (631)
                       (513)  (712)
                       (531)  (721)
                       (612)
                       (621)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Cases[Join@@Permutations/@IntegerPartitions[n,{3}],{x_,y_,z_}/;x!=y!=z&&x+y!=z &&x!=y+z]],{n,0,30}]

Formula

Conjectures from Colin Barker, May 14 2019: (Start)
G.f.: 2*x^6*(1 + 3*x + 3*x^2 + 5*x^3) / ((1 - x)^3*(1 + x)^2*(1 + x^2)*(1 + x + x^2)).
a(n) = a(n-2) + a(n-3) + a(n-4) - a(n-5) - a(n-6) - a(n-7) + a(n-9) for n>9. (End)
Above conjecture confirmed for n <= 5000. - Fausto A. C. Cariboni, Feb 17 2022

A325689 Number of length-3 compositions of n such that no part is the sum of the other two.

Original entry on oeis.org

0, 0, 0, 1, 0, 6, 4, 15, 12, 28, 24, 45, 40, 66, 60, 91, 84, 120, 112, 153, 144, 190, 180, 231, 220, 276, 264, 325, 312, 378, 364, 435, 420, 496, 480, 561, 544, 630, 612, 703, 684, 780, 760, 861, 840, 946, 924, 1035, 1012, 1128, 1104, 1225, 1200, 1326, 1300, 1431
Offset: 0

Views

Author

Gus Wiseman, May 15 2019

Keywords

Comments

A composition of n is a finite sequence of positive integers summing to n.
Confirmed recurrence relation from Colin Barker for n <= 5000. - Fausto A. C. Cariboni, Feb 15 2022

Examples

			The a(3) = 1 through a(8) = 12 compositions (empty columns not shown):
  (111)  (113)  (114)  (115)  (116)
         (122)  (141)  (124)  (125)
         (131)  (222)  (133)  (152)
         (212)  (411)  (142)  (161)
         (221)         (151)  (215)
         (311)         (214)  (233)
                       (223)  (251)
                       (232)  (323)
                       (241)  (332)
                       (313)  (512)
                       (322)  (521)
                       (331)  (611)
                       (412)
                       (421)
                       (511)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n,{3}],And@@Table[#[[i]]!=Total[Delete[#,i]],{i,3}]&]],{n,0,30}]

Formula

Conjectures from Colin Barker, May 16 2019: (Start)
G.f.: x^3*(1 - x + 4*x^2) / ((1 - x)^3*(1 + x)^2) for n>5.
a(n) = -(5 + 3*(-1)^n - 2*n) * (n-2) / 4 for n>0.
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5).
(End)
Showing 1-10 of 11 results. Next