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.
%I A133871 #63 Jan 30 2025 11:24:43
%S A133871 2,4,6,10,12,20,24,34,44,64,78,116,148,208,286,410,556,808,1120,1620,
%T A133871 2308,3352,4784,6980,10064,14680,21296,31128,45276,66288,96712,141654,
%U A133871 207156,303716,444748,652612,956884,1404920,2062080,3029564,4450120
%N A133871 a(n) = the definite integral Integral_{0..1} Product_{j=1..n} 4*sin^2(Pi*j*x) dx.
%C A133871 This quantity arises in some examples associated to the dynamical Mertens's theorem for quasihyperbolic toral automorphisms.
%C A133871 The function being integrated to compute a_n vanishes on the set of points in the Farey sequence of level n. I am particularly interested in knowing how large the sequence is asymptotically.
%C A133871 a(n) = coefficient of x^(n*(n+1)/2) in the polynomial (-1)^n*Product_{k=1..n} (1-x^k)^2, and is the maximal such coefficient as well. - _Steven Finch_, Feb 03 2009
%H A133871 Vaclav Kotesovec, <a href="/A133871/b133871.txt">Table of n, a(n) for n = 1..1000</a> (terms 1..174 from Robert Israel)
%H A133871 Miklos Bóna, R. Gómez, M. D. Ward, <a href="http://www.birs.ca/workshops/2016/16w5048/report16w5048.pdf">Workshop in Analytic and Probabilistic Combinatorics BIRS-16w5048</a> 2016.
%H A133871 S. R. Finch, <a href="https://www.semanticscholar.org/paper/Signum-Equations-and-Extremal-Coefficients-Finch/66da606fec632a2e752a6f421641fb84148b810f#paperDetail">Signum equations and extremal coefficients</a>.
%H A133871 Steven R. Finch, <a href="/A000980/a000980.pdf">Signum equations and extremal coefficients</a>, February 7, 2009. [Cached copy, with permission of the author]
%H A133871 Jeffrey Gaither, Guy Louchard, Stephan Wagner, and Mark Daniel Ward, <a href="http://www.stat.purdue.edu/~mdw/papers/paper028.pdf">Resolution of T. Ward's Question and the Israel-Finch Conjecture: Precise Analysis of an Integer Sequence Arising in Dynamics</a>, Combinatorics, Probability and Computing, 24 (2015), 195-215. Special Issue Honouring the Memory of Philippe Flajolet.
%H A133871 S. Jaidee, S. Stevens and T. Ward, <a href="http://arxiv.org/abs/0801.2082">Mertens' theorem for toral automorphisms</a>, arXiv:0801.2082 [math.DS], 2008-2010.
%H A133871 S. Jaidee, S. Stevens and T. Ward, <a href="https://doi.org/10.1090/S0002-9939-2010-10632-9">Mertens' theorem for toral automorphisms</a>, Proc. Amer. Math. Soc. 139 (2011), no. 5, 1819-1824.
%H A133871 Yasuhiko Kamiyama, <a href="https://www.researchgate.net/profile/Yasuhiko-Kamiyama/publication/388355304_THE_EULER_CHARACTERISTIC_OF_THE_FIBER_PRODUCT_OF_MORSE_FUNCTIONS/">The Euler characteristic of the fiber product of Morse functions</a>, Bull. Korean Math. Soc. (2025) Vol. 62, No. 1, pp. 71-80. See pp. 73, 75.
%H A133871 Mark Daniel Ward, <a href="http://aofa2013.lsi.upc.edu/slides/Ward.pdf">Resolution of T. Ward's Question and the Israel-Finch Conjecture. Precise Asymptotic Analysis of an Integer Sequence Motivated by the Dynamical Mertens' Theorem for Quasihyperbolic Toral Automorphisms</a>, Slides, 2013.
%H A133871 T. Ward, D. W. Cantrell and R. Israel, <a href="http://mathforum.org/kb/message.jspa?messageID=6062821">sci.math.research discussion</a>, 2008.
%F A133871 a(n) = sum of squares of coefficients in Product_{k=1..n} (1-x^k). - _Paul D. Hanna_, Nov 30 2010
%F A133871 a(n) ~ c * d^n / sqrt(n), where d = 1.48770584269062356180051131... and c = 2.40574583936181024... [Ward, 2013]. - _Vaclav Kotesovec_, May 03 2018
%e A133871 a(2) = 4 since Integral_{0..1} sin^2(Pi*x) sin^2(2*Pi*x) dx = 1/4.
%p A133871 a:= n->int(product(4*(sin(Pi*j*x))^2, j=1..n), x=0..1); seq(a(n), n=1..10);
%p A133871 # second Maple program:
%p A133871 A133871:= k -> (-1)^k*coeff(mul((t^j-1)^2,j=1..k),t,k*(k+1)/2);
%p A133871 # _Robert Israel_, Mar 15 2013
%t A133871 p = 1; Table[p = Expand[p*(1 - x^n)^2]; Max[(-1)^n*CoefficientList[p, x]], {n, 1, 100}] (* _Vaclav Kotesovec_, May 03 2018 *)
%t A133871 (* The constant "d" *) Chop[-E^(-I*(Pi^2*(1 + 6*x^2) - 6*PolyLog[2, E^(2*I*Pi*x)]) / (6*Pi*x)) /. x -> (x /. FindRoot[Pi*(Pi*(-1 + 6*x^2) + 12*I*x*Log[1 - E^(2*I*Pi*x)]) + 6*PolyLog[2, E^(2*I*Pi*x)], {x, 4/5}, WorkingPrecision -> 100])] (* _Vaclav Kotesovec_, May 04 2018 *)
%o A133871 (PARI) a(n)=sum(k=0, n*(n+1)/2, polcoeff(prod(m=1, n, 1-x^m+x*O(x^k)), k)^2) \\ _Paul D. Hanna_
%Y A133871 Cf. A005728, A047653.
%K A133871 nonn
%O A133871 1,1
%A A133871 _Thomas Ward_, Jan 07 2008
%E A133871 More terms from _Steven Finch_, Feb 03 2009