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-6 of 6 results.

A007494 Numbers that are congruent to 0 or 2 mod 3.

Original entry on oeis.org

0, 2, 3, 5, 6, 8, 9, 11, 12, 14, 15, 17, 18, 20, 21, 23, 24, 26, 27, 29, 30, 32, 33, 35, 36, 38, 39, 41, 42, 44, 45, 47, 48, 50, 51, 53, 54, 56, 57, 59, 60, 62, 63, 65, 66, 68, 69, 71, 72, 74, 75, 77, 78, 80, 81, 83, 84, 86, 87, 89, 90, 92, 93, 95, 96, 98, 99, 101, 102, 104, 105, 107
Offset: 0

Views

Author

Christopher Lam Cham Kee (Topher(AT)CyberDude.Com)

Keywords

Comments

The map n -> a(n) (where a(n) = 3n/2 if n even or (3n+1)/2 if n odd) was studied by Mahler, in connection with "Z-numbers" and later by Flatto. One question was whether, iterating from an initial integer, one eventually encountered an iterate = 1 (mod 4). - Jeff Lagarias, Sep 23 2002
Partial sums of 0,2,1,2,1,2,1,2,1,... . - Paul Barry, Aug 18 2007
a(n) = numbers k such that antiharmonic mean of the first k positive integers is not integer. A169609(a(n-1)) = 3. See A146535 and A169609. Complement of A016777. - Jaroslav Krizek, May 28 2010
Range of A173732. - Reinhard Zumkeller, Apr 29 2012
Number of partitions of 6n into two odd parts. - Wesley Ivan Hurt, Nov 15 2014
Numbers m such that 3 divides A000217(m). - Bruno Berselli, Aug 04 2017
Maximal length of a snake like polyomino that fits in a 2 X n rectangle. - Alain Goupil, Feb 12 2020

References

  • L. Flatto, Z-numbers and beta-transformations, in Symbolic dynamics and its applications (New Haven, CT, 1991), 181-201, Contemp. Math., 135, Amer. Math. Soc., Providence, RI, 1992.

Crossrefs

Complement of A016777.
Range of A002517.
Cf. A274406. [Bruno Berselli, Jun 26 2016]

Programs

Formula

a(n) = 3*n/2 if n even, otherwise (3*n+1)/2.
If u(1)=0, u(n) = n + floor(u(n-1)/3), then a(n-1) = u(n). - Benoit Cloitre, Nov 26 2002
G.f.: x*(x+2)/((1-x)^2*(1+x)). - Ralf Stephan, Apr 13 2002
a(n) = 3*floor(n/2) + 2*(n mod 2) = A032766(n) + A000035(n). - Reinhard Zumkeller, Apr 04 2005
a(n) = (6*n+1)/4 - (-1)^n/4; a(n) = Sum_{k=0..n-1} (1 + (-1)^(k/2)*cos(k*Pi/2)). - Paul Barry, Aug 18 2007
A145389(a(n)) <> 1. - Reinhard Zumkeller, Oct 10 2008
a(n) = A002943(n) - A173511(n). - Reinhard Zumkeller, Feb 20 2010
a(n) = 3*n - a(n-1) - 1 (with a(0)=0). - Vincenzo Librandi, Nov 18 2010
a(n) = Sum_{k>=0} A030308(n,k)*A042950(k). - Philippe Deléham, Oct 17 2011
a(n) = n + ceiling(n/2). - Arkadiusz Wesolowski, Sep 18 2012
a(n) = 2n - floor(n/2) = floor((3n+1)/2) = n + (n + (n mod 2))/2. - Wesley Ivan Hurt, Oct 19 2013
a(n) = A000217(n+1) - A099392(n+1). - Bui Quang Tuan, Mar 27 2015
a(n) = n + floor(n/2) + (n mod 2). - Bruno Berselli, Apr 04 2016
a(n) = Sum_{i=1..n} numerator(2/i). - Wesley Ivan Hurt, Feb 26 2017
a(n) = Sum_{k=0..n-1} Sum_{i=0..k} C(i,k)+(-1)^(k-i). - Wesley Ivan Hurt, Sep 20 2017
E.g.f.: (3*exp(x)*x + sinh(x))/2. - Stefano Spezia, Feb 11 2020
Sum_{n>=1} (-1)^(n+1)/a(n) = log(3)/2 - Pi/(6*sqrt(3)). - Amiram Eldar, Dec 04 2021

A003105 Schur's 1926 partition theorem: number of partitions of n into parts 6n+1 or 6n-1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 3, 3, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20, 23, 26, 30, 34, 38, 42, 47, 53, 60, 67, 74, 82, 91, 102, 114, 126, 139, 153, 169, 187, 207, 228, 250, 274, 301, 331, 364, 399, 436, 476, 520, 569, 622, 679, 739, 804, 875, 953, 1038, 1128, 1224, 1327
Offset: 0

Views

Author

Keywords

Comments

There are many (at least 8) equivalent definitions of this sequence (besides the comments below, see also Schur, Alladi, Andrews). - N. J. A. Sloane, Jun 17 2011
Coefficients of replicable function number 72e. - N. J. A. Sloane, Jun 10 2015
Also number of partitions of n into odd parts in which no part appears more than twice, cf. A070048 and A096938. - Vladeta Jovovic, Jan 18 2005
Also number of partitions of n into distinct parts congruent to 1 or 2 modulo 3. (Follows from second g.f.) - N. Sato, Jul 20 2005
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Convolution of A262928 and A261612. - Vaclav Kotesovec, Jan 13 2017
Convolution of A109702 and A109701. - Vaclav Kotesovec, Jan 21 2017

Examples

			G.f: A(x) = 1 + x + x^2 + x^3 + x^4 + 2*x^5 + 2*x^6 + 3*x^7 + 3*x^8 + 3*x^9 + 4*x^10 + ...
T72e = 1/q + q^11 + q^23 + q^35 + q^47 + 2*q^59 + 2*q^71 + 3*q^83 + ...
The logarithm of the g.f. begins:
log(A(x)) = x + x^2/2 + x^3/3 + x^4/4 + 6*x^5/5 + x^6/6 + 8*x^7/7 + x^8/8 + x^9/9 + 6*x^10/10 + 12*x^11/11 + x^12/12 + ... + A186099(n)*x^n/n + ... . - _Paul D. Hanna_, Feb 17 2013
		

References

  • K. Alladi, Refinements of Rogers-Ramanujan type identities. In Special Functions, q-Series and Related Topics (Toronto, ON, 1995), 1-35, Fields Inst. Commun., 14, Amer. Math. Soc., Providence, RI, 1997.
  • G. E. Andrews, Schur's theorem, partitions with odd parts and the Al-Salam-Carlitz polynomials. In q-Series From a Contemporary Perspective (South Hadley, MA, 1998), 45-56, Contemp. Math., 254, Amer. Math. Soc., Providence, RI, 2000.
  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • I. Schur, Zur Additiven Zahlentheorie, Ges. Abh., Vol. 2, Springer, pp. 43-50.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a003105 n = p 1 n where
       p k m | m == 0 = 1 | m < k = 0 | otherwise = q k (m-k) + p (k+2) m
       q k m | m == 0 = 1 | m < k = 0 | otherwise = p (k+2) (m-k) + p (k+2) m
    -- Reinhard Zumkeller, Nov 12 2011
  • Maple
    with(combinat);
    A:=proc(n) local i, j, t3, t2, t1;
        t2:=0;
        t1:=firstpart(n);
        for j from 1 to numbpart(n)+2 do
            t3:=1;
            for i from 1 to nops(t1) do
                if (t1[i] mod 6) <> 1 and (t1[i] mod 6) <> 5 then t3:=0; fi;
            od;
            if t3=1 then t2:=t2+1; fi;
            if nops(t1) = 1 then RETURN(t2); fi;
            t1:=nextpart(t1);
        od;
    end;
    # brute-force Maple program from N. J. A. Sloane, Jun 17 2011
  • Mathematica
    max = 63; f[x_] := 1/Product[1 - x^k + x^(2k), {k, 0, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 01 2011, after Vladeta Jovovic *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x] / QPochhammer[ -x^3, x^3], {x, 0, n}]; (* Michael Somos, Jul 05 2014 *)
    nmax = 100; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[If[Mod[k, 3] != 0, Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}];], {k, 2, nmax}]; poly (* Vaclav Kotesovec, Jan 13 2017 *)
    nmax = 63; kmax = nmax/6;
    s = Flatten[{Range[0, kmax]*6 + 1}~Join~{Range[kmax]*6 - 1}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Jul 31 2020 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^3 + A) / (eta(x + A) * eta(x^6 + A)), n))}; /* Michael Somos, Jan 09 2005 */
    
  • PARI
    {S(n,x)=sumdiv(n,d,d*(1-x^d)^(n/d))}
    {a(n)=polcoeff(exp(sum(k=1,n,S(k,x)*x^k/k)+x*O(x^n)),n)}
    for(n=0,60,print1(a(n),", "))
    /* Paul D. Hanna, Feb 17 2013 */
    

Formula

G.f.: 1/Product_{k>=0} (1-x^(6*k+1))*(1-x^(6*k+5)) = Product_{k>=0} (1+x^(3*k+1))*(1+x^(3*k+2)) = 1/Product_{k>=0} (1-x^k+x^(2*k)). - Vladeta Jovovic, Jun 08 2003
Expansion of chi(-x^3) / chi(-x) in powers of x where chi() is a Ramanujan theta function. - Michael Somos, Mar 04 2012
Expansion of f(x, x^2) / f(-x^3) = f(-x^6) / f(-x, -x^5) in powers of x where f() is Ramanujan theta function. - Michael Somos, Jul 05 2014
Expansion of q^(1/12) * eta(q^2) * eta(q^3) / (eta(q) * eta(q^6)) in powers of q. - Michael Somos, Jan 09 2005
Euler transform of period 6 sequence [1, 0, 0, 0, 1, 0, ...]. - Michael Somos, Jan 09 2005
Given g.f. A(x), then B(q) = (A(q^12) / q)^4 satisfies 0 = f(B(q), B(q^2)) where f(u, v) = u*v^4 + (1 - u^3) * v^3 + 6*u^2*v^2 + (u^4 - u)*v + u^3. - Michael Somos, Jan 09 2005
The logarithmic derivative equals A186099. - Paul D. Hanna, Feb 17 2013
G.f.: exp( Sum_{n>=1} A186099(n) * x^n/n ) where A186099(n) = sum of divisors of n congruent to 1 or 5 mod 6. - Paul D. Hanna, Feb 17 2013
G.f.: exp( Sum_{n>=1} S(n,x) * x^n/n ) where S(n,x) = Sum_{d|n} d*(1-x^d)^(n/d). - Paul D. Hanna, Feb 17 2013
a(n) ~ Pi*sqrt(2) / sqrt(3*(12*n-1)) * BesselI(1, Pi*sqrt(12*n-1) / (3*sqrt(6))) ~ exp(Pi*sqrt(2*n)/3) / (2^(5/4) * sqrt(3) * n^(3/4)) * (1 - (9/(8*Pi) + Pi/36)/sqrt(2*n) + (5 - 135/(4*Pi^2) + Pi^2/81)/(64*n)). - Vaclav Kotesovec, Aug 23 2015, extended Jan 09 2017
a(n) = (1/n)*Sum_{k=1..n} A186099(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 21 2017

Extensions

More terms from Vladeta Jovovic, Jun 08 2003

A132463 Number of partitions of n into distinct parts congruent to 0 or 1 modulo 3.

Original entry on oeis.org

1, 1, 0, 1, 2, 1, 1, 3, 2, 2, 5, 4, 3, 7, 7, 5, 10, 11, 8, 14, 17, 13, 20, 25, 19, 27, 36, 29, 37, 50, 43, 51, 69, 61, 69, 94, 86, 93, 126, 120, 125, 167, 164, 167, 220, 222, 222, 287, 297, 294, 373, 393, 386, 481, 516, 505, 617, 672, 657, 788, 868, 850, 1002, 1114, 1094
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 22 2007

Keywords

Examples

			a(7)=3 because we have 7, 61 and 43.
		

Crossrefs

Programs

  • Maple
    g:=product((1+x^(3*k))*(1+x^(3*k-2)),k=1..30): gser:=series(g,x=0,100): seq(coeff(gser,x,n),n=0..65); # Emeric Deutsch, Aug 26 2007
  • Mathematica
    nmax = 100; CoefficientList[Series[Product[((1+x^(3*k))*(1+x^(3*k-2))), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 24 2015 *)

Formula

G.f.: Product(k>=1, (1+x^(3*k))*(1+x^(3*k-2)) ). - Emeric Deutsch, Aug 26 2007
a(n) ~ exp(Pi*sqrt(2*n)/3) / (2^(19/12) * sqrt(3) * n^(3/4)). - Vaclav Kotesovec, Aug 24 2015

Extensions

Prepended a(0) = 1, Joerg Arndt, Feb 22 2015

A035361 Number of partitions of n into parts 3k or 3k+2.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 3, 2, 5, 5, 6, 9, 11, 11, 18, 19, 23, 31, 36, 40, 56, 60, 73, 92, 105, 121, 155, 170, 204, 247, 280, 325, 397, 440, 521, 615, 695, 805, 954, 1063, 1244, 1442, 1626, 1873, 2176, 2431, 2813, 3221, 3623, 4146, 4751, 5309, 6086, 6905, 7746, 8807
Offset: 0

Views

Author

Keywords

Comments

Euler transform of period 3 sequence [ 0, 1, 1, ...]. - Kevin T. Acres, Apr 28 2018

Examples

			1 + x^2 + x^3 + x^4 + 2*x^5 + 3*x^6 + 2*x^7 + 5*x^8 + 5*x^9 + 6*x^10 + ...
		

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[1/((1 - x^(3*k))*(1 - x^(3*k-1))), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 23 2015 *)
    nmax = 55; kmax = nmax/3;
    s = Flatten[{Range[0, kmax]*3}~Join~{Range[0, kmax]*3 + 2}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Aug 02 2020 *)

Formula

a(n) ~ Gamma(2/3) * exp(2*Pi*sqrt(n)/3) / (4 * sqrt(3) * n^(13/12) * Pi^(1/3)). - Vaclav Kotesovec, Aug 23 2015

A374065 Expansion of Product_{k>=1} 1 / (1 + x^(3*k-2)).

Original entry on oeis.org

1, -1, 1, -1, 0, 0, 0, -1, 2, -2, 1, 0, -1, 0, 2, -3, 3, -1, -1, 1, 1, -4, 5, -3, 0, 2, 0, -4, 7, -6, 1, 3, -2, -3, 9, -10, 4, 3, -5, -1, 11, -15, 10, 1, -8, 3, 10, -20, 17, -3, -10, 9, 7, -24, 26, -10, -10, 15, 2, -27, 37, -21, -8, 22, -6, -28, 49, -36, -2, 30, -19, -24, 61, -56, 10, 35
Offset: 0

Views

Author

Ilya Gutkovskiy, Jun 27 2024

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 75; CoefficientList[Series[Product[1/(1 + x^(3 k - 2)), {k, 1, nmax}], {x, 0, nmax}], x]
    a[0] = 1; a[n_] := a[n] = (1/n) Sum[DivisorSum[k, (-1)^(k/#) # &, Mod[#, 3] == 1 &] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 75}]

Formula

a(0) = 1; a(n) = -Sum_{k=1..n} A261612(k) * a(n-k).
a(n) = Sum_{k=0..n} A081362(k) * A132462(n-k).
a(n) = Sum_{k=0..n} A109389(k) * A262928(n-k).

A374060 Expansion of Product_{k>=1} (1 - x^(3*k-1)) * (1 - x^(3*k)).

Original entry on oeis.org

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

Views

Author

Ilya Gutkovskiy, Jun 27 2024

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 85; CoefficientList[Series[Product[(1 - x^(3 k - 1)) (1 - x^(3 k)), {k, 1, nmax}], {x, 0, nmax}], x]
    a[0] = 1; a[n_] := a[n] = -(1/n) Sum[Plus @@ Select[Divisors[k], Mod[#, 3] != 1 &] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 85}]

Formula

a(0) = 1; a(n) = -(1/n) * Sum_{k=1..n} A082050(k) * a(n-k).
a(0) = 1; a(n) = -Sum_{k=1..n} A035361(k) * a(n-k).
a(n) = Sum_{k=0..n} A010815(k) * A035382(n-k).
Showing 1-6 of 6 results.