A090129 Smallest exponent such that -1 + 3^a(n) is divisible by 2^n.
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
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
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).
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
Comments