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

A090129 Smallest exponent such that -1 + 3^a(n) is divisible by 2^n.

Original entry on oeis.org

1, 2, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608, 16777216, 33554432, 67108864, 134217728, 268435456, 536870912, 1073741824, 2147483648, 4294967296, 8589934592, 17179869184
Offset: 1

Views

Author

Labos Elemer and Ralf Stephan, Jan 19 2004

Keywords

Comments

A131577 and A011782 are companions, A131577(n) + A011782(n) = 2^n, (and differences each other). - Paul Curtz, Jan 18 2009
A090127 with offset 0: (1, 2, 2, 4, 8, ...) = A(x) / A(x^2), when A(x) = (1 + 2x + 4x^2 + 8x^3 + ...). - Gary W. Adamson, Feb 20 2010
From Wolfdieter Lang, Apr 18 2012: (Start)
a(n) is the order of 3 modulo 2^n. For n=1 and 2 this is obviously 1 and 2, respectively, and for n >= 3 it is 2^(n-2).
For a proof see, e.g., the Graeme McRae link under A068531, the section 'A Different Approach', proposed by Alexander Monnas, the first part, where the result from the expansion of (4-1)^(2^(k-2)) holds only for k >= 3. See also the Charles R Greathouse IV program below where this result has been used.
This means that the cycle generated by 3, taken modulo 2^n, has length a(n), and that 3 is not a primitive root modulo 2^n, if n >= 3 (because Euler's phi(2^n) = 2^(n-1), n >= 1, see A000010).
(End)
Let r(x) = (1 + 2x + 2x^2 + 4x^3 + ...). Then (1 + 2x + 4x^2 + 8x^3 + ...) = (r(x) * r(x^2) * r(x)^4 * r(x^8) * ...). - Gary W. Adamson, Sep 13 2016

Examples

			a(1) = 1 since -1 + 3 = 2 is divisible by 2^1;
a(2) = a(3) = 2 since -1 + 9 = 8 is divisible by 4 = 2^2 and also by 8 = 2^3;
a(5) = 8 since -1 + 6561 = 6560 = 32*205 is divisible by 2^5.
From _Wolfdieter Lang_, Apr 18 2012: (Start)
n=3: the order of 3 (mod 8) is a(3)=2 because the cycle generated by 3 is [3, 3^2==1 (mod 8)].
n=5: a(5) = 2^3 = 8 because the cycle generated by 3 is [3^1=3, 3^2=9, 3^3=27, 17, 19, 25, 11, 1] (mod 32).
  The multiplicative group mod 32 is non-cyclic (see A033949(10)) with the additional four cycles  [5, 25, 29, 17, 21, 9, 13, 1], [7, 17, 23, 1], [15, 1], and [31, 1]. This is the cycle structure of the (Abelian) group Z_8 x Z_2 (see one of the cycle graphs shown in the Wikipedia link 'List of small groups' for the order phi(32)=16, given under A192005).
(End)
		

Crossrefs

Essentially the same as A000079.

Programs

  • Mathematica
    t=Table[Part[Flatten[FactorInteger[ -1+3^(n)]], 2], {n, 1, 130}] Table[Min[Flatten[Position[t, j]]], {j, 1, 10}]
    Join[{1,2},2^Range[30]] (* or *) Join[{1,2},NestList[2#&,2,30]] (* Harvey P. Dale, Nov 08 2012 *)
  • PARI
    a(n)=2^(n+(n<3)-2) \\ Charles R Greathouse IV, Apr 09 2012
    
  • Python
    def A090129(n): return n if n<3 else 1<Chai Wah Wu, Jul 11 2022

Formula

a(n) = 2^(n-2) if n >= 3, 1 for n=1 and 2 for n=2 (see the order comment above).
a(n+2) = A152046(n) + A152046(n+1) = 2*A011782(n). - Paul Curtz, Jan 18 2009

Extensions

a(11) through a(20) from R. J. Mathar, Aug 08 2008
More terms (powers of 2, see a comment above) from Wolfdieter Lang, Apr 18 2012

A090739 Exponent of 2 in 9^n - 1.

Original entry on oeis.org

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

Views

Author

Labos Elemer and Ralf Stephan, Jan 19 2004

Keywords

Comments

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

Examples

			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.
		

Crossrefs

Cf. A000005, A006519, A120738 (partial sums).
Appears in A161737.

Programs

  • Maple
    A090739 := proc(n)
        padic[ordp](9^n-1,2) ;
    end proc:
    seq(A090739(n),n=1..80) ; # R. J. Mathar, Jun 28 2025
  • Mathematica
    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 *)
  • PARI
    a(n)=valuation(n,2)+3 \\ Charles R Greathouse IV, Mar 14 2014
    
  • Python
    def A090739(n): return (~n&n-1).bit_length()+3 # Chai Wah Wu, Jul 11 2022

Formula

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

Extensions

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

A090740 Exponent of 2 in 3^n - 1.

Original entry on oeis.org

1, 3, 1, 4, 1, 3, 1, 5, 1, 3, 1, 4, 1, 3, 1, 6, 1, 3, 1, 4, 1, 3, 1, 5, 1, 3, 1, 4, 1, 3, 1, 7, 1, 3, 1, 4, 1, 3, 1, 5, 1, 3, 1, 4, 1, 3, 1, 6, 1, 3, 1, 4, 1, 3, 1, 5, 1, 3, 1, 4, 1, 3, 1, 8, 1, 3, 1, 4, 1, 3, 1, 5, 1, 3, 1, 4, 1, 3, 1, 6, 1, 3, 1, 4, 1, 3, 1, 5, 1, 3, 1, 4, 1, 3, 1, 7, 1, 3, 1, 4, 1, 3, 1, 5, 1
Offset: 1

Views

Author

Labos Elemer and Ralf Stephan, Jan 19 2004

Keywords

Comments

Also the 2-adic order of Fibonacci(3n) [Lengyel]. - R. J. Mathar, Nov 05 2008

Examples

			n=2: 3^2 - 1 = 8 = 2^3 so a(2)=3.
		

Crossrefs

Programs

  • Maple
    seq(padic:-ordp(3^n-1, 2), n=1..100); # Robert Israel, Dec 28 2015
  • Mathematica
    Table[Part[Flatten[FactorInteger[ -1+3^n]], 2], {n, 1, 70}]
    IntegerExponent[#,2]&/@(3^Range[110]-1) (* Harvey P. Dale, Jan 28 2017 *)
  • PARI
    a(n)=if(n<1,0,if(n%2==0,a(n/2)+1+(n/2)%2,1)) /* Ralf Stephan, Jan 23 2004 */
    
  • PARI
    a(n)=valuation(fibonacci(3*n),2); \\ Joerg Arndt, Oct 28 2012
    
  • PARI
    a(n)=my(t=valuation(n,2)); if(t,t+2,1) \\ Charles R Greathouse IV, Mar 14 2014
    
  • Python
    def A090740(n): return (n&-n).bit_length()+int(not n&1) # Chai Wah Wu, Jul 11 2022

Formula

a(n) = A007814(n) + A059841(n) + 1.
Multiplicative with a(p^e) = e+2 if p = 2; 1 if p > 2. G.f.: A(x) = 1/(1-x^2) + Sum_{k>=0} x^(2^k)/(1-x^(2^k)). - Vladeta Jovovic, Jan 19 2004
G.f.: Sum_{k>=0} t*(1+2*t+t^2+t^3)/(1-t^4) with t=x^2^k. Recurrence: a(2n) = a(n) + 1 + [n odd], a(2n+1) = 1. - Ralf Stephan, Jan 23 2004
a(n) = A337923(3*n). [Lengyel]. - R. J. Mathar, Nov 05 2008
G.f. A(x) satisfies A(x) = A(x^2) + x/(1-x) + x^2/(1-x^4). - Robert Israel, Dec 28 2015
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 5/2. - Amiram Eldar, Nov 28 2022
Dirichlet g.f.: zeta(s)*(2^s+1-1/2^s)/(2^s-1). - Amiram Eldar, Jan 04 2023
Showing 1-3 of 3 results.