A161737 -id:A161737 - cp's OEIS frontend

A090739 Exponent of 2 in 9^n - 1.

3, 4, 3, 5, 3, 4, 3, 6, 3, 4, 3, 5, 3, 4, 3, 7, 3, 4, 3, 5, 3, 4, 3, 6, 3, 4, 3, 5, 3, 4, 3, 8, 3, 4, 3, 5, 3, 4, 3, 6, 3, 4, 3, 5, 3, 4, 3, 7, 3, 4, 3, 5, 3, 4, 3, 6, 3, 4, 3, 5, 3, 4, 3, 9, 3, 4, 3, 5, 3, 4, 3, 6, 3, 4, 3, 5, 3, 4, 3, 7, 3, 4, 3, 5, 3, 4, 3

Offset: 1

Labos Elemer and Ralf Stephan, Jan 19 2004

The exponent of 2 in the factorization of Fibonacci(6n). - T. D. Noe, Mar 14 2014

Records of 3, 4, 5, 6, 7, 8,.. occur at n= 1, 2, 4, 8, 16, 32,... - R. J. Mathar, Jun 28 2025

			For n = 2, we see that -1 + 3^4 = 80 = 2^4 * 5 so a(2) = 4.
For n = 3, we see that -1 + 3^6 = 728 = 2^3 * 7 * 13, so a(3) = 3.

T. D. Noe, Table of n, a(n) for n = 1..1000
T. Lengyel, The order of the Fibonacci and Lucas numbers, Fib. Quart. 33 (1995), 234-239.

Cf. A024101, A069895, A091512, A088660, A090740, A220466.
Cf. A000005, A006519, A120738 (partial sums).
Appears in A161737.

A090739 := proc(n)
    padic[ordp](9^n-1,2) ;
end proc:
seq(A090739(n),n=1..80) ; # R. J. Mathar, Jun 28 2025

Table[Part[Flatten[FactorInteger[ -1+3^(2*n)]], 2], {n, 1, 70}]
Table[IntegerExponent[Fibonacci[n], 2], {n, 6, 600, 6}] (* T. D. Noe, Mar 14 2014 *)

a(n)=valuation(n,2)+3 \\ Charles R Greathouse IV, Mar 14 2014

def A090739(n): return (~n&n-1).bit_length()+3 # Chai Wah Wu, Jul 11 2022

a(n) = A007814(n) + 3.
a((2*n-1)*2^p) = p + 3, p >= 0. - Johannes W. Meijer, Feb 08 2013
a(n) = log_2(A006519(9^n - 1)). - Alonso del Arte, Feb 08 2013
a(n) = 2*tau(4*n)/(tau(4*n) - tau(n)), where tau(n) = A000005(n). - Peter Bala, Jan 06 2021
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 4. - Amiram Eldar, Nov 28 2022

More terms from T. D. Noe, Mar 14 2014

A050605 Column/row 2 of A050602: a(n) = add3c(n,2).

Original entry on oeis.org

0, 0, 1, 1, 0, 0, 2, 2, 0, 0, 1, 1, 0, 0, 3, 3, 0, 0, 1, 1, 0, 0, 2, 2, 0, 0, 1, 1, 0, 0, 4, 4, 0, 0, 1, 1, 0, 0, 2, 2, 0, 0, 1, 1, 0, 0, 3, 3, 0, 0, 1, 1, 0, 0, 2, 2, 0, 0, 1, 1, 0, 0, 5, 5, 0, 0, 1, 1, 0, 0, 2, 2, 0, 0, 1, 1, 0, 0, 3, 3, 0, 0, 1, 1, 0, 0, 2, 2

Offset: 0

Antti Karttunen, Jun 22 1999

nonn

It seems that (n - Sum_{k=1..n} a(k) )/log(n) is bounded. - Benoit Cloitre, Oct 03 2002
2^a(n-1) is the highest power of 2 dividing the triangular number A000217(n) = n*(n+1)/2, for n >= 1. - Benoit Cloitre, Oct 03 2002 [corrected and rewritten by Wolfdieter Lang, Nov 21 2019]
a(n) is the number of trailing 0's in the binary reflected Gray code of n+1 (A014550). - Amiram Eldar, May 15 2021

Vincenzo Librandi, Table of n, a(n) for n = 0..5000
Isabel Cação, Helmuth R. Malonek, Maria Irene Falcão, and Graça Tomaz, Intrinsic Properties of a Non-Symmetric Number Triangle, J. Int. Seq., Vol. 26 (2023), Article 23.4.8.

Bisection gives column/row 1 of A050602: A007814.

Cf. A000217, A001511, A161737, A069834, A014550.

[Valuation(n*(n+1)/2, 2): n in [1..120]]; // Vincenzo Librandi, Aug 11 2017

with(Bits): add3c := proc(a, b) option remember; `if`(0 = And(a, b), 0, 1 + add3c(Xor(a, b), 2*And(a, b))) end: A050605 := n -> add3c(n, 2):
seq(A050605(n), n=0..80); # Johannes W. Meijer, Jun 18 2009; updated by Peter Luschny, Jul 12 2019

Table[IntegerExponent[(n + 1)(n + 2)/2, 2], {n, 0, 100}] (* Jean-François Alcover, Mar 04 2016 *)

PARI
```
a(n)=valuation(n*(n+1)/2,2)
```

a(4*n+2) = A001511(n). - Johannes W. Meijer, Jun 18 2009

a(n) = A007814(n+1) + A007814(n+2) - 1. - Ridouane Oudra, Oct 08 2019

A161736 Denominators of the column sums of the BG2 matrix.

Original entry on oeis.org

1, 9, 75, 1225, 19845, 160083, 1288287, 41409225, 1329696225, 10667118605, 85530896451, 1371086188563, 21972535073125, 176021737014375, 1409850293610375, 90324408810638025, 5786075364399106425, 46326420401234675625, 370882277949065911875, 5938020471163465810125

Offset: 2

Johannes W. Meijer, Jun 18 2009

The BG2 matrix coefficients, see also A008956, are defined by BG2[2m,1] = 2*beta(2m+1) and the recurrence relation BG2[2m,n] = BG2[2m,n-1] - BG2[2m-2,n-1]/(2*n-3)^2 for m = -2, -1, 0, 1, 2, .. and n = 2, 3, .. , with beta(m) = sum((-1)^k/(1+2*k)^m, k=0..infinity). We observe that beta(2m+1) = 0 for m = -1, -2, -3, .. .
A different way to define the matrix coefficients is BG2[2*m,n] = (1/m)*sum(LAMBDA(2*m-2*k,n-1)*BG2[2*k,n], k=0..m-1) with LAMBDA(2*m,n-1) = (1-2^(-2*m))*zeta(2*m)-sum((2*k-1)^(-2*m), k=1..n-1) and BG2[0,n] = Pi/2 for m = 0, 1, 2, .. , and n = 1, 2, 3 .. , with zeta(m) the Riemann zeta function.
The columns sums of the BG2 matrix are defined by sb(n) = sum(BG2[2*m,n], m=0..infinity) for n = 2, 3, .. . For large values of n the value of sb(n) approaches Pi/2.
It is remarkable that if we assume that BG2[2m,1] = 2 for m = 0, 1, .. the columns sums of the modified matrix converge to the original sb(n) values. The first Maple program makes use of this phenomenon and links the sb(n) with the central factorial numbers A008956.
The column sums sb(n) can be linked to other sequences, see the second Maple program.
We observe that the column sums sb(n) of the BG2(n) matrix are related to the column sums sl(n) of the LG2(n) matrix, see A008956, by sb(n) = (-1)^(n+1)*(2*n-1)*sl(n).
a(n+2), for n >= 0, seems to coincide with the numerators belonging to A278145. - Wolfdieter Lang, Nov 16 2016
Suppose that, given values f(x-2*n+1), f(x-2*n+3), ..., f(x-1), f(x+1), ..., f(x+2*n-3), f(x+2*n-1), we approximate f(x) using the first 2*n terms of its Taylor series. Then 1/sb(n+1) is the coefficient of f(x-1) and f(x+1). - Matthew House, Dec 03 2024

			sb(2) = 2; sb(3) = 16/9; sb(4) = 128/75; sb(5) = 2048/1225; etc..

G. C. Greubel, Table of n, a(n) for n = 2..830
R. Arratia, S. Garibaldi, and J. Kilian, Asymptotic distribution for the birthday problem with multiple coincidences, via an embedding of the collision process, arXiv:1310.7055 [math.PR], 2013.

Cf. Numerators A161737.
Cf. A001818, A001902, A008956, A013777, A134372 and A134375.
Cf. A001147, A133221, A161738 and A002474, A278145.

[Denominator((2^(4*n-5)*(Factorial(n-1))^4)/((n-1)*(Factorial(2*n-2))^2)): n in [2..20]]; // G. C. Greubel, Sep 26 2018

nmax := 18; for n from 0 to nmax do A001818(n) := (doublefactorial(2*n-1))^2 od: for n from 0 to nmax do A008956(n, 0):=1 od: for n from 0 to nmax do A008956(n, n) := A001818(n) od: for n from 1 to nmax do for m from 1 to n-1 do A008956(n, m) := (2*n-1)^2*A008956(n-1, m-1) + A008956(n-1, m) od: od: for n from 1 to nmax do for m from 0 to n do s(n, m):=0; s(n, m) := s(n, m)+ sum((-1)^k1*A008956(n, n-k1), k1=0..n-m): od: sb1(n+1) := sum(s(n, k1), k1=1..n) * 2/A001818(n); od: seq(sb1(n), n=2..nmax); # End program 1
nmax1 := nmax; for n from 0 to nmax1 do A001147(n):= doublefactorial(2*n-1) od: for n from 0 to nmax1/2 do A133221(2*n+1) := A001147(n); A133221(2*n) := A001147(n) od: for n from 0 to nmax1 do A002474(n) := 2^(2*n+1)*n!*(n+1)! od: for n from 1 to nmax1 do A161738(n) := ((product((2*n-3-2*k1), k1=0..floor(n/2-1)))) od: for n from 2 to nmax1 do sb2(n) := A002474(n-2) / (A161738(n)*A133221(n-1))^2 od: seq(sb2(n), n=2..nmax1); # End program 2
# Above Maple programs edited by Johannes W. Meijer, Sep 25 2012
r := n -> (1/Pi)*(2*n - 2)*((n - 3/2)!/(n - 1)!)^2: a := n -> numer(simplify(r(n))):
seq(a(n), n = 1..21);  # Peter Luschny, Feb 12 2025

sb[2]=2; sb[n_] := sb[n] = sb[n-1]*4*(n-1)*(n-2)/(2n-3)^2; Table[sb[n] // Denominator, {n, 2, 20}] (* Jean-François Alcover, Aug 14 2017 *)

{a(n) = if( n<2, 0, n--; numerator( binomial( 2*n, n)^2 * n / 2^(n+1) ))}; /* Michael Somos, May 09 2011 */

a(n) = denom(sb(n)) with sb(n) = (2^(4*n-5)*(n-1)!^4)/((n-1)*(2*n-2)!^2) and A161737(n) = numer(sb(n)).
a(n+1) = numerator of C(2*n,n)^2 * n / 2^(n+1). - Michael Somos, May 09 2011
a(n) = A001902(2*n-3). - Mats Granvik, Nov 25 2018
a(n) = numerator((1/Pi)*(2*n - 2)*((n - 3/2)!/(n - 1)!)^2). - Peter Luschny, Feb 13 2025

A380949 a(n) = numerator(r(n)) where r(n) = (n/2)(Pi/2)^cos(Pi(n-1))*((n/2-1/2)!/(n/2)!)^2.

Original entry on oeis.org

0, 1, 1, 4, 9, 64, 75, 256, 1225, 16384, 19845, 65536, 160083, 1048576, 1288287, 4194304, 41409225, 1073741824, 1329696225, 4294967296, 10667118605, 68719476736, 85530896451, 274877906944, 1371086188563, 17592186044416, 21972535073125, 70368744177664, 176021737014375

Offset: 0

Peter Luschny, Feb 11 2025

			r(n) = 0, 1, 1/2, 4/3, 9/16, 64/45, 75/128, 256/175, 1225/2048, ...

Paolo Xausa, Table of n, a(n) for n = 0..1000

Cf. A380950 (denominator), A380910, A380909, A019267 (asymptotic coefficients).

Cf. A161736, A056982, A124399, A161737, A069955, A038534, A278145, A001901.

r := n -> (n/2)*(Pi/2)^cos(Pi*(n-1))*((n/2-1/2)!/(n/2)!)^2:
a := n -> numer(simplify(r(n))): seq(a(n), n = 0..28);
# Alternative:
r := n -> ifelse(n <= 1, n, (n - 1)/(n*r(n - 1))):

Join[{0}, Numerator[FoldList[(#2 - 1)/(#2*#) &, Range[30]]]] (* Paolo Xausa, Feb 14 2025 *)

Product_{k=1..n} a(k) = A380910(n) / A380909(n).
r(n) = (n - 1)/(n*r(n - 1)) for n > 1.
numerator(r(2*n)) = A161736(n).
numerator(r(2*n+1)) = A056982(n).
numerator(r(2*n+1))/4^n = A124399(n).
denominator(r(2*n-2)) = A161737(n).
denominator(r(2*n+1)) = A069955(n).
denominator(r(2*n+1))/(2*n+1) = A038534(n).
denominator(r(2*n+2))/2 = A278145(n).
denominator(r(2*n+2))/2^(2*n+1) = A001901(n).
r(n) ~ (2/Pi)^((-1)^n)*(1 - 1/(2*n) + 1/(8*n^2) + 1/(16*n^3) - 5/(128*n^4) - 23/(256*n^5) ...).

cp's OEIS Frontend

A090739 Exponent of 2 in 9^n - 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A050605 Column/row 2 of A050602: a(n) = add3c(n,2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A161736 Denominators of the column sums of the BG2 matrix.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A380949 a(n) = numerator(r(n)) where r(n) = (n/2)*(Pi/2)^cos(Pi*(n-1))*((n/2-1/2)!/(n/2)!)^2.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A380949 a(n) = numerator(r(n)) where r(n) = (n/2)(Pi/2)^cos(Pi(n-1))*((n/2-1/2)!/(n/2)!)^2.