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 A372214 #13 May 07 2024 07:06:40 %S A372214 1,4,40,487,6376,86629,1203823,16984678,242274280,3484593028, %T A372214 50444222665,734066291974,10728052396111,157349171819155, %U A372214 2314894133906086,34145661019248487,504810905195542504,7478066502444399439,110972913533524676080,1649407167353221551706,24549982881130265421001 %N A372214 a(n) is equal to the n-th order Taylor polynomial (centered at 0) of G(x)^n evaluated at x = 1, where G(x) = (1 - 2*x - sqrt(1 - 8*x + 4*x^2))/(2*x). %C A372214 x*G(x) = (1 - 2*x - sqrt(1 - 8*x + 4*x^2))/2 is the o.g.f. of A047891. %C A372214 The sequence satisfies the Gauss congruences: a(n*p^k) == a(n*p^(k-1)) ( mod p^k ) for all primes p and positive integers n and k. %C A372214 We conjecture that the sequence satisfies the stronger supercongruences a(n*p^k) == a(n*p^(k-1)) ( mod p^(3*k) ) for all primes p >= 5 and positive integers n and k. Examples of these supercongruences are given below. %C A372214 More generally, for each integer m, we conjecture that the sequence {a_m(n) : n >= 0}, defined by setting a_m(n) = the n-th order Taylor polynomial of G(x)^(m*n) evaluated at x = 1, satisfies the same supercongruences. %F A372214 a(n) = [x^n] ( (1 + x)*G(x/(1 + x)) )^n. %F A372214 O.g.f.: ( 1 + x*F'(x)/F(x) )/( 1 - x*F(x) ), where F(x) = (1/x)*Revert( x/G(x) ) = = 1 + 3*x + 21*x^2 + 192*x^3 + 2001*x^4 + ... is the o.g.f. of A219535. %F A372214 Row sums of the Riordan array ( 1 + x*F'(x)/F(x), x*F(x) ) belonging to the Hitting time subgroup of the Riordan group. %F A372214 a(n) ~ sqrt(1 + 17/sqrt(33)) * (59 + 11*sqrt(33))^n / (sqrt(3*Pi*n) * 2^(3*n + 3/2)). - _Vaclav Kotesovec_, May 02 2024 %e A372214 n-th order Taylor polynomial of G(x)^n: %e A372214 n = 0: G(x)^0 = 1 + O(x) %e A372214 n = 1: G(x)^1 = 1 + 3*x + O(x^2) %e A372214 n = 2: G(x)^2 = 1 + 6*x + 33*x^2 + O(x^3) %e A372214 n = 3: G(x)^3 = 1 + 9*x + 63*x^2 + 414*x^3 + O(x^4) %e A372214 n = 4: G(x)^4 = 1 + 12*x + 102*x^2 + 768*x^3 + 5493*x^4 + O(x^5) %e A372214 Setting x = 1 gives a(0) = 1, a(1) = 1 + 3 = 4, a(2) = 1 + 6 + 33 = 40, a(3) = 1 + 9 + 63 + 414 = 487 and a(4) = 1 + 12 + 102 + 768 + 5493 = 6376. %e A372214 The triangle of coefficients of the n-th order Taylor polynomial of G(x)^n, n >= 0, in descending powers of x begins %e A372214 row sums %e A372214 n = 0 | 1 1 %e A372214 n = 1 | 3 1 4 %e A372214 n = 2 | 33 6 1 40 %e A372214 n = 3 | 414 63 9 1 487 %e A372214 n = 4 | 5493 768 102 12 1 6376 %e A372214 ... %e A372214 This is a Riordan array belonging to the Hitting time subgroup of the Riordan group. %e A372214 Examples of supercongruences: %e A372214 a(13) - a(1) = 157349171819155 - 4 = (3^3)*(13^3)*269*9860941 == 0 (mod 13^3). %e A372214 a(2*7) - a(2) = 2314894133906086 - 40 = 2*(3^4)*(7^3)*11*12119*312509 == 0 (mod 7^3). %p A372214 G := x -> (1/2)*(1 - 2*x - sqrt(1 - 8*x + 4*x^2))/x: %p A372214 H := (x, n) -> series(G(x)^n, x, 41): %p A372214 seq(add(coeff(H(x, n), x, k), k = 0..n), n = 0..20); %t A372214 Table[SeriesCoefficient[(2*(1 + x)^2/(1 - x + Sqrt[1 - 6*x - 3*x^2]))^n, {x, 0, n}], {n, 0, 20}] (* _Vaclav Kotesovec_, May 02 2024 *) %o A372214 (PARI) G(x) = (1 - 2*x - sqrt(1 - 8*x + 4*x^2))/(2*x); %o A372214 a(n) = my(x='x+O('x^(n+2))); subst(Pol(Vec(G(x)^n)), 'x, 1); \\ _Michel Marcus_, May 07 2024 %Y A372214 Cf. A047891, A219535, A333090, A333093, A372215. %K A372214 nonn,easy %O A372214 0,2 %A A372214 _Peter Bala_, Apr 23 2024