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

A005843 The nonnegative even numbers: a(n) = 2n.

Original entry on oeis.org

0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120
Offset: 0

Views

Author

Keywords

Comments

-2, -4, -6, -8, -10, -12, -14, ... are the trivial zeros of the Riemann zeta function. - Vivek Suri (vsuri(AT)jhu.edu), Jan 24 2008
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-2) is the number of 2-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
A134452(a(n)) = 0; A134451(a(n)) = 2 for n > 0. - Reinhard Zumkeller, Oct 27 2007
Omitting the initial zero gives the number of prime divisors with multiplicity of product of terms of n-th row of A077553. - Ray Chandler, Aug 21 2003
A059841(a(n))=1, A000035(a(n))=0. - Reinhard Zumkeller, Sep 29 2008
(APSO) Alternating partial sums of (a-b+c-d+e-f+g...) = (a+b+c+d+e+f+g...) - 2*(b+d+f...), it appears that APSO(A005843) = A052928 = A002378 - 2*(A116471), with A116471=2*A008794. - Eric Desbiaux, Oct 28 2008
A056753(a(n)) = 1. - Reinhard Zumkeller, Aug 23 2009
Twice the nonnegative numbers. - Juri-Stepan Gerasimov, Dec 12 2009
The number of hydrogen atoms in straight-chain (C(n)H(2n+2)), branched (C(n)H(2n+2), n > 3), and cyclic, n-carbon alkanes (C(n)H(2n), n > 2). - Paul Muljadi, Feb 18 2010
For n >= 1; a(n) = the smallest numbers m with the number of steps n of iterations of {r - (smallest prime divisor of r)} needed to reach 0 starting at r = m. See A175126 and A175127. A175126(a(n)) = A175126(A175127(n)) = n. Example (a(4)=8): 8-2=6, 6-2=4, 4-2=2, 2-2=0; iterations has 4 steps and number 8 is the smallest number with such result. - Jaroslav Krizek, Feb 15 2010
For n >= 1, a(n) = numbers k such that arithmetic mean of the first k positive integers is not integer. A040001(a(n)) > 1. See A145051 and A040001. - Jaroslav Krizek, May 28 2010
Union of A179082 and A179083. - Reinhard Zumkeller, Jun 28 2010
a(k) is the (Moore lower bound on and the) order of the (k,4)-cage: the smallest k-regular graph having girth four: the complete bipartite graph with k vertices in each part. - Jason Kimberley, Oct 30 2011
For n > 0: A048272(a(n)) <= 0. - Reinhard Zumkeller, Jan 21 2012
Let n be the number of pancakes that have to be divided equally between n+1 children. a(n) is the minimal number of radial cuts needed to accomplish the task. - Ivan N. Ianakiev, Sep 18 2013
For n > 0, a(n) is the largest number k such that (k!-n)/(k-n) is an integer. - Derek Orr, Jul 02 2014
a(n) when n > 2 is also the number of permutations simultaneously avoiding 213, 231 and 321 in the classical sense which can be realized as labels on an increasing strict binary tree with 2n-1 nodes. See A245904 for more information on increasing strict binary trees. - Manda Riehl Aug 07 2014
It appears that for n > 2, a(n) = A020482(n) + A002373(n), where all sequences are infinite. This is consistent with Goldbach's conjecture, which states that every even number > 2 can be expressed as the sum of two prime numbers. - Bob Selcoe, Mar 08 2015
Number of partitions of 4n into exactly 2 parts. - Colin Barker, Mar 23 2015
Number of neighbors in von Neumann neighborhood. - Dmitry Zaitsev, Nov 30 2015
Unique solution b( ) of the complementary equation a(n) = a(n-1)^2 - a(n-2)*b(n-1), where a(0) = 1, a(1) = 3, and a( ) and b( ) are increasing complementary sequences. - Clark Kimberling, Nov 21 2017
Also the maximum number of non-attacking bishops on an (n+1) X (n+1) board (n>0). (Cf. A000027 for rooks and queens (n>3), A008794 for kings or A030978 for knights.) - Martin Renner, Jan 26 2020
Integer k is even positive iff phi(2k) > phi(k), where phi is Euler's totient (A000010) [see reference De Koninck & Mercier]. - Bernard Schott, Dec 10 2020
Number of 3-permutations of n elements avoiding the patterns 132, 213, 312 and also number of 3-permutations avoiding the patterns 213, 231, 321. See Bonichon and Sun. - Michel Marcus, Aug 20 2022
a(n) gives the y-value of the integral solution (x,y) of the Pellian equation x^2 - (n^2 + 1)*y^2 = 1. The x-value is given by 2*n^2 + 1 (see Tattersall). - Stefano Spezia, Jul 24 2025

Examples

			G.f. = 2*x + 4*x^2 + 6*x^3 + 8*x^4 + 10*x^5 + 12*x^6 + 14*x^7 + 16*x^8 + ...
		

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 2.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 28.
  • J.-M. De Koninck and A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 529a pp. 71 and 257, Ellipses, 2004, Paris.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 256.

Crossrefs

a(n)=2*A001477(n). - Juri-Stepan Gerasimov, Dec 12 2009
Moore lower bound on the order of a (k,g) cage: A198300 (square); rows: A000027 (k=2), A027383 (k=3), A062318 (k=4), A061547 (k=5), A198306 (k=6), A198307 (k=7), A198308 (k=8), A198309 (k=9), A198310 (k=10), A094626 (k=11); columns: A020725 (g=3), this sequence (g=4), A002522 (g=5), A051890 (g=6), A188377 (g=7). - Jason Kimberley, Oct 30 2011
Cf. A231200 (boustrophedon transform).

Programs

Formula

G.f.: 2*x/(1-x)^2.
E.g.f.: 2*x*exp(x). - Geoffrey Critzer, Aug 25 2012
G.f. with interpolated zeros: 2x^2/((1-x)^2 * (1+x)^2); e.g.f. with interpolated zeros: x*sinh(x). - Geoffrey Critzer, Aug 25 2012
Inverse binomial transform of A036289, n*2^n. - Joshua Zucker, Jan 13 2006
a(0) = 0, a(1) = 2, a(n) = 2a(n-1) - a(n-2). - Jaume Oliver Lafont, May 07 2008
a(n) = Sum_{k=1..n} floor(6n/4^k + 1/2). - Vladimir Shevelev, Jun 04 2009
a(n) = A034856(n+1) - A000124(n) = A000217(n) + A005408(n) - A000124(n) = A005408(n) - 1. - Jaroslav Krizek, Sep 05 2009
a(n) = Sum_{k>=0} A030308(n,k)*A000079(k+1). - Philippe Deléham, Oct 17 2011
Digit sequence 22 read in base n-1. - Jason Kimberley, Oct 30 2011
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Vincenzo Librandi, Dec 23 2011
a(n) = 2*n = Product_{k=1..2*n-1} 2*sin(Pi*k/(2*n)), n >= 0 (undefined product := 1). See an Oct 09 2013 formula contribution in A000027 with a reference. - Wolfdieter Lang, Oct 10 2013
From Ilya Gutkovskiy, Aug 19 2016: (Start)
Convolution of A007395 and A057427.
Sum_{n>=1} (-1)^(n+1)/a(n) = log(2)/2 = (1/2)*A002162 = (1/10)*A016655. (End)
From Bernard Schott, Dec 10 2020: (Start)
Sum_{n>=1} 1/a(n)^2 = Pi^2/24 = A222171.
Sum_{n>=1} (-1)^(n+1)/a(n)^2 = Pi^2/48 = A245058. (End)

A016742 Even squares: a(n) = (2*n)^2.

Original entry on oeis.org

0, 4, 16, 36, 64, 100, 144, 196, 256, 324, 400, 484, 576, 676, 784, 900, 1024, 1156, 1296, 1444, 1600, 1764, 1936, 2116, 2304, 2500, 2704, 2916, 3136, 3364, 3600, 3844, 4096, 4356, 4624, 4900, 5184, 5476, 5776, 6084, 6400, 6724, 7056, 7396, 7744, 8100, 8464
Offset: 0

Views

Author

Keywords

Comments

4 times the squares.
Number of edges in the complete bipartite graph of order 5n, K_{n,4n} - Roberto E. Martinez II, Jan 07 2002
It is conjectured (I think) that a regular Hadamard matrix of order n exists iff n is an even square (cf. Seberry and Yamada, Th. 10.11). A Hadamard matrix is regular if the sum of the entries in each row is the same. - N. J. A. Sloane, Nov 13 2008
Sequence arises from reading the line from 0, in the direction 0, 16, ... and the line from 4, in the direction 4, 36, ... in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008
The entries from a(1) on can be interpreted as pair sums of (2, 2), (8, 8), (18, 18), (32, 32) etc. that arise from a re-arrangement of the subshell orbitals in the periodic table of elements. 8 becomes the maximum number of electrons in the (2s,2p) or (3s,3p) orbitals, 18 the maximum number of electrons in (4s,3d,4p) or (5s,3d,5p) shells, for example. - Julio Antonio Gutiérrez Samanez, Jul 20 2008
The first two terms of the sequence (n=1, 2) give the numbers of chemical elements using only n types of atomic orbitals, i.e., there are a(1)=4 elements (H,He,Li,Be) where electrons reside only on s-orbitals, there are a(2)=16 elements (B,C,N,O,F,Ne,Na,Mg,Al,Si,P,S,Cl,Ar,K,Ca) where electrons reside only on s- and p-orbitals. However, after that, there is 37 (which is one more than a(3)=36) elements (from Sc, Scandium, atomic number 21 to La, Lanthanum, atomic number 57) where electrons reside only on s-, p- and d-orbitals. This is because Lanthanum (with the electron configuration [Xe]5d^1 6s^2) is an exception to the Aufbau principle, which would predict that its electron configuration is [Xe]4f^1 6s^2. - Antti Karttunen, Aug 14 2008.
Number of cycles of length 3 in the king's graph associated with an (n+1) X (n+1) chessboard. - Anton Voropaev (anton.n.voropaev(AT)gmail.com), Feb 01 2009
a(n+1) is the molecular topological index of the n-star graph S_n. - Eric W. Weisstein, Jul 11 2011
a(n) is the sum of two consecutives odd numbers 2*n^2-1 and 2*n^2+1 and the difference of two squares (n^2+1)^2 - (n^2-1)^2. - Pierre CAMI, Jan 02 2012
For n > 3, a(n) is the area of the irregular quadrilateral created by the points ((n-4)*(n-3)/2,(n-3)*(n-2)/2), ((n-2)*(n-1)/2,(n-1)*n/2), ((n+1)*(n+2)/2,n*(n+1)/2), and ((n+3)*(n+4)/2,(n+2)*(n+3)/2). - J. M. Bergot, May 27 2014
Number of terms less than 10^k: 1, 2, 5, 16, 50, 159, 500, 1582, 5000, 15812, 50000, 158114, 500000, ... - Muniru A Asiru, Jan 28 2018
Right-hand side of the binomial coefficient identity Sum_{k = 0..2*n} (-1)^(k+1)* binomial(2*n,k)*binomial(2*n + k,k)*(2*n - k) = a(n). - Peter Bala, Jan 12 2022

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd ed., 1994, p. 99.
  • Seberry, Jennifer and Yamada, Mieko; Hadamard matrices, sequences and block designs, in Dinitz and Stinson, eds., Contemporary design theory, pp. 431-560, Wiley-Intersci. Ser. Discrete Math. Optim., Wiley, New York, 1992.
  • W. D. Wallis, Anne Penfold Street and Jennifer Seberry Wallis, Combinatorics: Room squares, sum-free sets, Hadamard matrices, Lecture Notes in Mathematics, Vol. 292, Springer-Verlag, Berlin-New York, 1972. iv+508 pp.

Crossrefs

Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.
Cf. sequences listed in A254963.
Other n X n king graph cycle counts: A288918 (4-cycles), A288919 (5-cycles), A288920 (6-cycles).
Cf. A016813.

Programs

Formula

O.g.f.: 4*x*(1+x)/(1-x)^3. - R. J. Mathar, Jul 28 2008
a(n) = A000290(n)*4 = A001105(n)*2. - Omar E. Pol, May 21 2008
a(n) = A155955(n,2) for n > 1. - Reinhard Zumkeller, Jan 31 2009
Sum_{n>=1} 1/a(n) = (1/4)*Pi^2/6 = Pi^2/24. - Ant King, Nov 04 2009
a(n) = a(n-1) + 8*n - 4 (with a(0)=0). - Vincenzo Librandi, Nov 19 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) with a(0) = 0, a(1) = 4, a(2) = 16. - Philippe Deléham, Mar 26 2013
a(n) = A118729(8n+3). - Philippe Deléham, Mar 26 2013
Pi = 2*Product_{n>=1} (1 + 1/(a(n)-1)). - Adriano Caroli, Aug 04 2013
Pi = Sum_{n>=0} 8/(a(2n+1)-1). - Adriano Caroli, Aug 06 2013
E.g.f.: exp(x)*(4x^2 + 4x). - Geoffrey Critzer, Oct 07 2013
a(n) = A000384(n) + A014105(n). - Bruce J. Nicholson, Nov 11 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/48 (A245058). - Amiram Eldar, Oct 10 2020
From Amiram Eldar, Jan 25 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = sinh(Pi/2)/(Pi/2) (A308716).
Product_{n>=1} (1 - 1/a(n)) = sin(Pi/2)/(Pi/2) = 2/Pi (A060294). (End)
a(n) = A016754(n) - A016813(n). - Leo Tavares, Feb 24 2022

Extensions

More terms from Sabir Abdus-Samee (sabdulsamee(AT)prepaidlegal.com), Mar 13 2006

A139098 a(n) = 8*n^2.

Original entry on oeis.org

0, 8, 32, 72, 128, 200, 288, 392, 512, 648, 800, 968, 1152, 1352, 1568, 1800, 2048, 2312, 2592, 2888, 3200, 3528, 3872, 4232, 4608, 5000, 5408, 5832, 6272, 6728, 7200, 7688, 8192, 8712, 9248, 9800, 10368, 10952, 11552, 12168, 12800, 13448, 14112, 14792, 15488, 16200
Offset: 0

Views

Author

Omar E. Pol, Apr 25 2008

Keywords

Comments

Opposite numbers to the centered 16-gonal numbers (A069129) in the square spiral whose vertices are the triangular numbers (A000217).
8 times the squares. - Omar E. Pol, Dec 09 2008
a(n-1) is the molecular topological index of the n-wheel graph W_n. - Eric W. Weisstein, Jul 11 2011
An n X n pandiagonal magic square has a(n) orientations. - Kausthub Gudipati, Sep 15 2011
Area of a square with diagonal 4n. - Wesley Ivan Hurt, Jun 19 2014
Sum of all the parts in the partitions of 4n into exactly two parts. - Wesley Ivan Hurt, Jul 23 2014
Equivalently: integers k such that k$ / (k/2-1)! and k$ / (k/2)! are both squares when A000178 (k) = k$ = 1!*2!*...*k! is the superfactorial of k (see A348692 for further information). - Bernard Schott, Dec 02 2021

Crossrefs

Programs

Formula

a(n) = 8*A000290(n) = 4*A001105(n) = 2*A016742(n). - Omar E. Pol, Dec 13 2008
G.f.: -8*x*(1+x)/(x-1)^3. - R. J. Mathar, Nov 27 2015
From Amiram Eldar, Feb 03 2021: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/48 (A245058).
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/96.
Product_{n>=1} (1 + 1/a(n)) = sqrt(8)*sinh(Pi/sqrt(8))/Pi.
Product_{n>=1} (1 - 1/a(n)) = sqrt(8)*sin(Pi/sqrt(8))/Pi. (End)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Wesley Ivan Hurt, Dec 03 2021
From Elmo R. Oliveira, Dec 01 2024: (Start)
E.g.f.: 8*x*(1 + x)*exp(x).
a(n) = n*A008590(n) = A001105(2*n). (End)

A050449 a(n) = Sum_{d|n, d == 1 (mod 4)} d.

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 1, 1, 10, 6, 1, 1, 14, 1, 6, 1, 18, 10, 1, 6, 22, 1, 1, 1, 31, 14, 10, 1, 30, 6, 1, 1, 34, 18, 6, 10, 38, 1, 14, 6, 42, 22, 1, 1, 60, 1, 1, 1, 50, 31, 18, 14, 54, 10, 6, 1, 58, 30, 1, 6, 62, 1, 31, 1, 84, 34, 1, 18, 70, 6, 1, 10, 74, 38, 31, 1
Offset: 1

Views

Author

N. J. A. Sloane, Dec 23 1999

Keywords

Comments

Not multiplicative: a(3)*a(7) != a(21), for example. - R. J. Mathar, Dec 20 2011

Crossrefs

Cf. Sum_{d|n, d==1 (mod k)} d: A000593 (k=2), A078181 (k=3), this sequence (k=4), A284097 (k=5), A284098 (k=6), A284099 (k=7), A284100 (k=8).

Programs

  • Maple
    A050449 := proc(n)
            a := 0 ;
            for d in numtheory[divisors](n) do
                    if d mod 4 = 1 then
                            a := a+d ;
                    end if;
            end do:
            a;
    end proc:
    seq(A050449(n),n=1..40) ; # R. J. Mathar, Dec 20 2011
  • Mathematica
    a[n_] := DivisorSum[n, Boole[Mod[#, 4] == 1]*#&]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Jan 30 2018 *)
  • PARI
    a(n) = sumdiv(n, d, d*((d % 4) == 1)); \\ Michel Marcus, Jan 30 2018

Formula

G.f.: Sum_{n>=0} (4*n+1)*x^(4*n+1)/(1-x^(4*n+1)). - Vladeta Jovovic, Nov 14 2002
a(n) = A000593(n) - A050452(n). - Reinhard Zumkeller, Apr 18 2006
G.f.: Sum_{n >= 1} x^n*(1 + 3*x^(4*n))/(1 - x^(4*n))^2. - Peter Bala, Dec 19 2021
Sum_{k=1..n} a(k) = c * n^2 + O(n*log(n)), where c = Pi^2/48 = 0.205616... (A245058). - Amiram Eldar, Nov 26 2023

Extensions

More terms from Vladeta Jovovic, Nov 14 2002
More terms from Reinhard Zumkeller, Apr 18 2006

A050452 a(n) = Sum_{d|n, d == 3 (mod 4)} d.

Original entry on oeis.org

0, 0, 3, 0, 0, 3, 7, 0, 3, 0, 11, 3, 0, 7, 18, 0, 0, 3, 19, 0, 10, 11, 23, 3, 0, 0, 30, 7, 0, 18, 31, 0, 14, 0, 42, 3, 0, 19, 42, 0, 0, 10, 43, 11, 18, 23, 47, 3, 7, 0, 54, 0, 0, 30, 66, 7, 22, 0, 59, 18, 0, 31, 73, 0, 0, 14, 67, 0, 26, 42, 71, 3, 0, 0, 93, 19, 18
Offset: 1

Views

Author

N. J. A. Sloane, Dec 23 1999

Keywords

Crossrefs

Cf. Sum_{d|n, d=k-1 mod k} d: A000593 (k=2), A078182 (k=3), this sequence (k=4).

Programs

  • Maple
    A050452 := proc(n)
            a := 0 ;
            for d in numtheory[divisors](n) do
                    if d mod 4 = 3 then
                            a := a+d ;
                    end if;
            end do:
            a;
    end proc:
    seq(A050452(n),n=1..40) ; # R. J. Mathar, Dec 20 2011
  • Mathematica
    Table[Total[Select[Divisors[n],Mod[#,4]==3&]],{n,80}] (* Harvey P. Dale, Jul 07 2013 *)
  • PARI
    a(n) = sumdiv(n, d, d*((d % 4) == 3)); \\ Amiram Eldar, Nov 26 2023

Formula

a(n) = A000593(n) - A050449(n). - Reinhard Zumkeller, Apr 18 2006
G.f.: Sum_{k>=1} (4*k - 1)*x^(4*k-1)/(1 - x^(4*k-1)). - Ilya Gutkovskiy, Mar 21 2017
Sum_{k=1..n} a(k) = c * n^2 + O(n*log(n)), where c = Pi^2/48 = 0.205616... (A245058). - Amiram Eldar, Nov 26 2023

A096033 Difference between leg and hypotenuse in primitive Pythagorean triangles.

Original entry on oeis.org

1, 2, 8, 9, 18, 25, 32, 49, 50, 72, 81, 98, 121, 128, 162, 169, 200, 225, 242, 288, 289, 338, 361, 392, 441, 450, 512, 529, 578, 625, 648, 722, 729, 800, 841, 882, 961, 968, 1058, 1089, 1152, 1225, 1250, 1352, 1369, 1458, 1521, 1568, 1681, 1682, 1800, 1849
Offset: 1

Views

Author

Lekraj Beedassy, Jun 16 2004

Keywords

Comments

Consists of the odd squares and the halves of the even squares. - Andrew Weimholt, Sep 07 2010
Question: Do we have a(n) mod 2 = A004641(n)? - David A. Corneth, Jan 02 2019

References

  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 170.

Crossrefs

Programs

  • Mathematica
    nmax = 100;
    Union[2 Range[nmax]^2, (2 Range[0, Ceiling[nmax/Sqrt[2]]] + 1)^2] (* Jean-François Alcover, Jan 01 2019 *)
  • PARI
    upto(n) = vecsort(concat(vector((sqrtint(n)+1)\2, i, (2*i-1)^2), vector(sqrtint(n\2), i, 2*i^2))) \\ David A. Corneth, Jan 02 2019

Formula

Union of A001105 (integers of form 2*n^2) and A016754 (the odd squares).
Sum_{n>=1} 1/a(n) = 5*Pi^2/24 = 10 * A245058. - Amiram Eldar, Feb 14 2021

Extensions

Corrected and extended by Matthew Vandermast and Ray Chandler, Jun 17 2004
Erroneous comment deleted by Andrew Weimholt, Sep 07 2010

A284372 a(n) = Sum_{d|n, d = 0, 1, or 11 mod 12} d.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 12, 13, 14, 1, 1, 1, 1, 1, 1, 1, 1, 12, 24, 37, 26, 14, 1, 1, 1, 1, 1, 1, 12, 1, 36, 49, 38, 1, 14, 1, 1, 1, 1, 12, 1, 24, 48, 85, 50, 26, 1, 14, 1, 1, 12, 1, 1, 1, 60, 73, 62, 1, 1, 1, 14, 12, 1, 1, 24, 36, 72, 145, 74, 38, 26, 1
Offset: 1

Views

Author

Seiichi Manyama, Mar 25 2017

Keywords

Examples

			From _Peter Bala_, Dec 11 2020: (Start)
n = 24: n is not of the form m*(6*m +- 5), so e(n) = 0 and a(24) = a(23) + a(13) - a(10)  = 24 + 14 - 1  = 37;
n = 39: n = m*(6*m - 5) for m = 3, so e(n) = 39 and a(39) = 39 + a(38) + a(28) - a(25) - a(5) = 39 + 1 + 1 - 26 - 1 = 14;
n = 76: n = m*(6*m - 5) for m = 4, so e(n) = -76 and a(4) = -76 + a(75) + a(65) - a(62) - a(42) + a(37) + a(7) = -76 + 26 + 14  - 1 - 1 + 38 + 1 = 1. (End)
		

Crossrefs

Cf. A210964 (1/f(-x, -x^11)), A245058.
Cf. Sum_{d|n, d = 0, 1, or k-1 mod k} d: A000203 (k=3), A113184(k=4), A284361 (k=5), A284362 (k=6), A284363 (k=7), this sequence (k=12).

Programs

  • Mathematica
    Table[Sum[If[Mod[d, 12]<2 || Mod[d, 12]==11, d, 0], {d, Divisors[n]}], {n, 80}] (* Indranil Ghosh, Mar 25 2017 *)
    sd12[n_]:=Total[Select[Divisors[n],MemberQ[{0,1,11},Mod[#,12]]&]]; Array[sd12,80] (* Harvey P. Dale, Aug 29 2024 *)
  • PARI
    a(n) = sumdiv(n, d, ((d + 1) % 12 < 3) * d); \\ Amiram Eldar, Apr 12 2024

Formula

From Peter Bala, Dec 11 2020: (Start)
O.g.f.: Sum_{k >= 1, k == 0, 1 or 11 (mod 12)} k*x^k/(1 - x^k).
Define a(n) = 0 for n < 1. Then a(n) = e(n) + a(n-1) + a(n-11) - a(n-14) - a(n-34) + + - -, where [1, 11, 14, 34, ...] is the sequence of generalized 14-gonal numbers A195818, and e(n) = (-1)^(m+1)*n if n is a generalized 14-gonal number of the form m*(6*m+-5); otherwise e(n) = 0. Examples of this recurrence are given below. (End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/48 = -A245058 = 0.205616... . - Amiram Eldar, Apr 12 2024

A378995 Numerator of sigma(2*n)/(2*n).

Original entry on oeis.org

3, 7, 2, 15, 9, 7, 12, 31, 13, 21, 18, 5, 21, 2, 12, 63, 27, 91, 30, 9, 16, 21, 36, 31, 93, 49, 20, 15, 45, 14, 48, 127, 24, 63, 72, 65, 57, 35, 28, 93, 63, 8, 66, 45, 13, 42, 72, 21, 171, 217, 36, 105, 81, 70, 108, 31, 40, 105, 90, 3, 93, 56, 52, 255, 126, 28, 102, 135, 48, 12, 108, 403, 111, 133, 62, 75, 144, 98, 120
Offset: 1

Views

Author

Antti Karttunen, Dec 13 2024

Keywords

Crossrefs

Even bisection of A017665.
Topmost row of array A341605.
Cf. A062731, A245058, A378994, A378996 (denominators).
Cf. also A341525.

Programs

  • Mathematica
    a[n_] := Numerator[DivisorSigma[-1, 2*n]]; Array[a, 100] (* Amiram Eldar, Dec 14 2024 *)
  • PARI
    A378995(n) = numerator(sigma(2*n)/(2*n));

Formula

a(n) = A017665(2*n).
a(n) = A062731(n) / A378994(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A378996(k) = 5*Pi^2/24 (= 10 * A245058). - Amiram Eldar, Dec 14 2024

A143348 a(n) = -(-1)^n times sum of divisors of n.

Original entry on oeis.org

1, -3, 4, -7, 6, -12, 8, -15, 13, -18, 12, -28, 14, -24, 24, -31, 18, -39, 20, -42, 32, -36, 24, -60, 31, -42, 40, -56, 30, -72, 32, -63, 48, -54, 48, -91, 38, -60, 56, -90, 42, -96, 44, -84, 78, -72, 48, -124, 57, -93, 72, -98, 54, -120, 72, -120, 80, -90, 60, -168, 62, -96, 104, -127, 84, -144, 68, -126, 96
Offset: 1

Views

Author

Michael Somos, Aug 09 2008

Keywords

Examples

			q - 3*q^2 + 4*q^3 - 7*q^4 + 6*q^5 - 12*q^6 + 8*q^7 - 15*q^8 + 13*q^9 + ...
		

Crossrefs

-(-1)^n * A000203(n) = a(n). A143337(n) = 24 * a(n) unless n=0.
Cf. A245058 (-Pi^2/48).

Programs

  • Mathematica
    Table[-(-1)^n*DivisorSigma[1, n], {n, 69}] (* Michael De Vlieger, Aug 19 2017 *)
  • PARI
    {a(n) = if( n<1, 0, -(-1)^n * sigma(n))}

Formula

a(n) is multiplicative with a(2^e) = 1 - 2^(e+1) if e > 0, a(p^e) = (p^(e+1) - 1) / (p - 1) if p > 2.
G.f.: Sum_{k>0} k * -(-x)^k / (1 - (-x)^k) = Sum_{k>0} -(-x)^k / (1 - (-x)^k)^2.
Dirichlet g.f.: zeta(s) * zeta(s-1) * (1 - 6 / 2^s + 4 / 4^s).
Sum_{k=1..n} a(k) ~ -(Pi^2/48) * n^2. - Amiram Eldar, Oct 28 2023

A362984 Decimal expansion of the asymptotic mean of the abundancy index of the powerful numbers (A001694).

Original entry on oeis.org

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

Views

Author

Amiram Eldar, May 12 2023

Keywords

Comments

The abundancy index of a positive integer k is A000203(k)/k = A017665(k)/A017666(k).
The asymptotic mean of the abundancy index over all the positive integers is lim_{m->oo} (1/m) * Sum_{k=1..m} A000203(k)/k = Pi^2/6 = zeta(2) = 1.644934... (A013661).

Examples

			2.14968690301526765128219042105109416145987653275100999873...
		

Crossrefs

Similar constants (the asymptotic mean of the abundancy index of other sequences): A013661 (all positive integers), A082020 (cubefree), A111003 (odd), A157292 (5-free), A157294 (7-free), A157296 (9-free), A240976 (squares), A245058 (even), A306633 (squarefree), A362985 (cubefull).

Programs

  • Mathematica
    $MaxExtraPrecision = 1000; m = 1000; c = LinearRecurrence[{2, -3, 4, -6, 7, -7, 7, -6, 5, -3, 2, -1}, {0, 0, 0, 4, 5, 6, 0, -12, -9, -5, 0, 22}, m]; RealDigits[(2^4 + 2^2 + 2^(3/2) - 1)/(2^4 - 2)*(3^4 + 3^2 + 3^(3/2) - 1)/(3^4 - 3) * Exp[NSum[Indexed[c, n]*(PrimeZetaP[n/2] - 1/2^(n/2) - 1/3^(n/2))/n, {n, 4, m}, NSumTerms -> m, WorkingPrecision -> m]], 10, 120][[1]]
  • PARI
    prodeulerrat((p^8 + p^4 + p^3 - 1)/(p^8 - p^2), 1/2)

Formula

Equals lim_{m->oo} (1/m) * Sum_{k=1..m} A180114(k)/A001694(k).
Equals Product_{p prime} (p^4 + p^2 + p^(3/2) - 1)/(p^4 - p) = Product_{p prime} (1 + (p^2 + p^(3/2) + p - 1)/(p^4 - p)) (Jakimczuk and Lalín, 2022).
Showing 1-10 of 12 results. Next