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.

Previous Showing 11-20 of 43 results. Next

A094358 Squarefree products of factors of Fermat numbers (A023394).

Original entry on oeis.org

1, 3, 5, 15, 17, 51, 85, 255, 257, 641, 771, 1285, 1923, 3205, 3855, 4369, 9615, 10897, 13107, 21845, 32691, 54485, 65535, 65537, 114689, 163455, 164737, 196611, 274177, 319489, 327685, 344067, 494211, 573445, 822531, 823685, 958467, 974849, 983055
Offset: 1

Views

Author

Robert Munafo, Apr 26 2004

Keywords

Comments

641 is the first member not in sequences A001317, A004729, etc.
Conjectured (by Munafo, see link) to be the same as: numbers n such that 2^^n == 1 mod n, where 2^^n is A014221(n).
It is clear from the observations by Max Alekseyev in A023394 and the Chinese remainder theorem that any squarefree product b of divisors of Fermat numbers satisfies 2^(2^b) == 1 (mod b), hence satisfies Munafo's congruence above. The converse is true iff all Fermat numbers are squarefree. However, if nonsquarefree Fermat numbers exist, the criterion that is equivalent with Munafo's property would be "numbers b such that each prime power that divides b also divides some Fermat number". - Jeppe Stig Nielsen, Mar 05 2014
Also numbers b such that b is (squarefree and) a divisor of A051179(m) for some m. Or odd (squarefree) b where the multiplicative order of 2 mod b is a power of two. - Jeppe Stig Nielsen, Mar 07 2014
From Jianing Song, Nov 11 2023: (Start)
Also squarefree numbers k such that there exists i >= 1 such that k divides 2^^i - 1, where 2^^i = 2^2^...^2 (i times) = A014221(i): 2^^i == 1 (mod k) if and only if ord(2,k) divides 2^^(i-1) (ord(a,k) is the multiplicative order of a modulo k), so such i exists if and only if ord(2,k) is a power of 2. For such k, k divides 2^^i - 1 if and only if 2^^(i-2) >= log_2(ord(2,k)).
Note that 2^^(i-1) divides 2^^i implies that 2^^i - 1 divides 2^^(i+1) - 1, so this sequence is also squarefree numbers k such that k divides 2^^i - 1 for all sufficiently large i. (End)

Examples

			3 is a term because it is in A023394.
51 is a term because it is 3*17 and 17 is also in A023394.
153 = 3*3*17 is not a term because its factorization includes two 3's.
See the Munafo link for examples of the (conjectured) 2^^n == 1 (mod n) property.
		

Crossrefs

Programs

  • Mathematica
    kmax = 10^6;
    A023394 = Select[Prime[Range[kmax]], IntegerQ[Log[2, MultiplicativeOrder[2, #] ] ]&];
    Reap[For[k = 1, k <= kmax, k++, ff = FactorInteger[k]; If[k == 1 || AllTrue[ff, MemberQ[A023394, #[[1]]] && #[[2]] == 1 &], Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, Nov 03 2018 *)
  • PARI
    (  isOK1(n) = n%2==1 && hammingweight(znorder(Mod(2,n)))==1  ) ;  (  isOK2(n) = issquarefree(n) && isOK1(n)  )  \\ isOK1 and isOK2 differ only if n contains a squared prime that divides a Fermat number (none are known) \\ Jeppe Stig Nielsen, Apr 02 2014

Extensions

Edited by T. D. Noe, Feb 02 2009
Example brought in line with name/description by Robert Munafo, May 18 2011

A246674 Run Length Transform of A000225.

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 3, 7, 1, 1, 1, 3, 3, 3, 7, 15, 1, 1, 1, 3, 1, 1, 3, 7, 3, 3, 3, 9, 7, 7, 15, 31, 1, 1, 1, 3, 1, 1, 3, 7, 1, 1, 1, 3, 3, 3, 7, 15, 3, 3, 3, 9, 3, 3, 9, 21, 7, 7, 7, 21, 15, 15, 31, 63, 1, 1, 1, 3, 1, 1, 3, 7, 1, 1, 1, 3, 3, 3, 7, 15, 1, 1, 1, 3, 1, 1, 3, 7, 3, 3, 3, 9, 7, 7, 15, 31, 3, 3, 3, 9, 3, 3, 9, 21, 3, 3, 3, 9, 9, 9, 21, 45, 7, 7, 7, 21, 7, 7, 21, 49, 15, 15, 15, 45, 31, 31, 63, 127, 1
Offset: 0

Views

Author

Antti Karttunen, Sep 08 2014

Keywords

Comments

a(n) can be also computed by replacing all consecutive runs of zeros in the binary expansion of n with * (multiplication sign), and then performing that multiplication, still in binary, after which the result is converted into decimal. See the example below.

Examples

			115 is '1110011' in binary. The run lengths of 1-runs are 2 and 3, thus a(115) = A000225(2) * A000225(3) = ((2^2)-1) * ((2^3)-1) = 3*7 = 21.
The same result can be also obtained more directly, by realizing that '111' and '11' are the binary representations of 7 and 3, and 7*3 = 21.
From _Omar E. Pol_, Feb 15 2015: (Start)
Written as an irregular triangle in which row lengths are the terms of A011782:
1;
1;
1,3;
1,1,3,7;
1,1,1,3,3,3,7,15;
1,1,1,3,1,1,3,7,3,3,3,9,7,7,15,31;
1,1,1,3,1,1,3,7,1,1,1,3,3,3,7,15,3,3,3,9,3,3,9,21,7,7,7,21,15,15,31,63;
...
Right border gives 1 together with the positive terms of A000225.
(End)
		

Crossrefs

Cf. A003714 (gives the positions of ones).
A001316 is obtained when the same transformation is applied to A000079, the powers of two.
Run Length Transforms of other sequences: A071053, A227349, A246588, A246595, A246596, A246660, A246661, A246685, A247282.

Programs

  • Mathematica
    f[n_] := 2^n - 1; Table[Times @@ (f[Length[#]]&) /@ Select[ Split[ IntegerDigits[n, 2]], #[[1]] == 1&], {n, 0, 100}] (* Jean-François Alcover, Jul 11 2017 *)
  • Python
    # uses RLT function in A278159
    def A246674(n): return RLT(n,lambda m: 2**m-1) # Chai Wah Wu, Feb 04 2022

Formula

For all n >= 0, a(A051179(n)) = A247282(A051179(n)) = A051179(n).

A087046 Algebraic order of r_n, the value of r in the logistic map that corresponding to the onset of the period 2^n-cycle.

Original entry on oeis.org

1, 2, 12, 240, 65280, 4294901760, 18446744069414584320, 340282366920938463444927863358058659840, 115792089237316195423570985008687907852929702298719625575994209400481361428480
Offset: 1

Views

Author

Eric W. Weisstein, Aug 04 2003

Keywords

Crossrefs

Cf. A051179 (partial sums).

Programs

  • Mathematica
    Table[If[n <= 1, 1, 2^(2^(n - 1)) - 2^(2^(n - 2))], {n, 1, 10}] (* Cheng Zhang, Apr 02 2012 *)
  • PARI
    a(n) = 1<<(1<<(n-1)) - 1<<(1<<(n-2)); \\ Kevin Ryde, Jan 18 2024

Formula

a(n) = 2^(2^(n - 1)) - 2^(2^(n - 2)) with n>1, a(1)=1. - Cheng Zhang, Apr 02 2012
Sum_{n>=1} 1/a(n) = 1 + A346192. - Amiram Eldar, Jul 18 2021

Extensions

More terms from Cheng Zhang, Apr 03 2012

A190939 Subgroups of nimber addition interpreted as binary numbers.

Original entry on oeis.org

1, 3, 5, 9, 15, 17, 33, 51, 65, 85, 105, 129, 153, 165, 195, 255, 257, 513, 771, 1025, 1285, 1545, 2049, 2313, 2565, 3075, 3855, 4097, 4369, 4641, 5185, 6273, 8193, 8481, 8721, 9345, 10305, 12291, 13107, 15555, 16385, 16705, 17025, 17425, 18465, 20485, 21845
Offset: 0

Views

Author

Tilman Piesk, May 24 2011

Keywords

Comments

Each subgroup {0,a,b,...} of nimber addition can be assigned an integer 1+2^a+2^b+...
These integers ordered by size give this sequence.
Without nimbers the sequence may be defined as follows:
The powerset af a set {0,...,n-1} with the symmetric difference as group operation forms the elementary abelian group (Z_2)^n.
The elements of the group can be numbered lexicographically from 0 to 2^n-1, with 0 representing the neutral element:
{}-->0 , {0}-->2^0=1 , {1}-->2^1=2 , {0,1}-->2^0+2^1=3 , ... , {0,...,n-1}-->2^n-1
So the subgroups of (Z_2)^n can be represented by subsets of {0,...,2^n-1}.
So each subgroup {0,a,b,...} of (Z_2)^n can be assigned an integer 1+2^a+2^b+...
For each (Z_2)^n there is a finite sequence of these numbers ordered by size, and it is the beginning of the finite sequence for (Z_2)^(n+1).
This leads to the infinite sequence:
* 1, (1 until here for (Z_2)^0)
* 3, (2 until here for (Z_2)^1)
* 5, 9, 15, (5 until here for (Z_2)^2)
* 17, 33, 51, 65, 85, 105, 129, 153, 165, 195, 255, (16 until here for (Z_2)^3)
* 257, 513, 771, 1025, 1285, 1545, 2049, 2313, 2565, 3075, 3855, 4097, 4369, 4641, 5185, 6273, 8193, 8481, 8721, 9345, 10305, 12291, 13107, 15555, 16385, 16705, 17025, 17425, 18465, 20485, 21845, 23205, 24585, 26265, 26985, 32769, 33153, 33345, 33825, 34833, 36873, 38505, 39321, 40965, 42405, 43605, 49155, 50115, 52275, 61455, 65535, (67 until here for (Z_2)^4)
* 65537, ...
The number of subgroups of (Z_2)^n is 1, 2, 5, 16, 67, 374, 2825, ... (A006116)
Comment from Tilman Piesk, Aug 27 2013: (Start)
Boolean functions correspond to integers, and belong to small equivalence classes (sec). So a sec can be seen as an infinite set of integers (represented in A227722 by the smallest one). Some secs contain only one odd integer. These unique odd integers, ordered by size, are shown in this sequence. (While the smallest integers from these secs are shown in A227963.)
(End)

Examples

			The 5 subgroups of the Klein four-group (Z_2)^2 and corresponding integers are:
{0      }     -->     2^0                     =   1
{0,1    }     -->     2^0 + 2^1               =   3
{0,  2  }     -->     2^0       + 2^2         =   5
{0,    3}     -->     2^0             + 2^3   =   9
{0,1,2,3}     -->     2^0 + 2^1 + 2^2 + 2^3   =  15
		

Crossrefs

Cf. A227963 (the same small equivalence classes represented by entries of A227722)
Cf. A198260 (number of runs of ones in the binary strings)
Subsequences:
Cf. A051179 (2^2^n-1).
Cf. A083318 (2^n+1).
Cf. A001317 (rows of the Sierpinski triangle read like binary numbers).
Cf. A228540 (rows of negated binary Walsh matrices r.l.b.n.).
Cf. A122569 (negated iterations of the Thue-Morse sequence r.l.b.n.).

Extensions

Offset changed to 0 by Tilman Piesk, Jan 25 2012

A097029 Fixed points when the function f(x) = phi(x) + floor(x/2) is iterated, i.e., solutions to f(x) = x.

Original entry on oeis.org

1, 2, 3, 4, 8, 15, 16, 32, 64, 128, 255, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 65535, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608, 16777216, 33554432, 67108864, 83623935, 134217728, 268435456, 536870912, 1073741824, 2147483648, 4294967295
Offset: 1

Views

Author

Labos Elemer, Aug 27 2004

Keywords

Comments

Trivial fixed points are the powers of 2. How many nontrivial cases exist like 3, 15, 255, 65535: the first 5 terms of A051179. More?
83623935 is the next such term (see also A050474 and A203966). - Michel Marcus, Nov 13 2015

Examples

			For fixed points the cycle lengths are A097026(n=fix)=1, but the reverse is not true because long transients may also lead to 1-cycles.
So, e.g., 1910 is not here because its terminal 1-cycle is prefixed by a long transient: {1910, 1715, 2033, 2924, 2806, 2723, 3689, 4724, 4722, 3933, 4342, 4163, 6041, 8192, 8192}.
		

Crossrefs

Programs

  • PARI
    isok(n) = eulerphi(n) + n\2 == n; \\ Michel Marcus, Nov 13 2015

Extensions

a(30)-a(35) from Michel Marcus, Nov 13 2015
a(36)-a(38) from Jinyuan Wang, Jul 22 2021

A037097 Periodic vertical binary vectors of powers of 3, starting from bit-column 2 (halved).

Original entry on oeis.org

0, 12, 120, 57120, 93321840, 10431955353116229600, 8557304989566294213168677685339060480, 102743047168201563425402150421568484707810385382513037790885688657488312400960
Offset: 2

Views

Author

Antti Karttunen, Jan 29 1999

Keywords

Comments

Conjecture: For n>=3, each term a(n), when considered as a GF(2)[X]-polynomial, is divisible by GF(2)[X] -polynomial (x + 1) ^ A000225(n-1) (= A051179(n-2)). If this holds, then for n>=3, a(n) = A048720bi(A136386(n),A048723bi(3,A000225(n-1))) = A048720bi(A136386(n),A051179(n-2)).

Examples

			When powers of 3 are written in binary (see A004656), under each other as:
000000000001 (1)
000000000011 (3)
000000001001 (9)
000000011011 (27)
000001010001 (81)
000011110011 (243)
001011011001 (729)
100010001011 (2187)
it can be seen that, starting from the column 2 from the right, the bits in the n-th column can be arranged in periods of 2^(n-1): 4, 8, ... This sequence is formed from those bits: 0011, reversed is 11100, which is binary for 12, thus a(3) = 12, 00011110, reversed is 011110000, which is binary for 120, thus a(4) = 120.
		

References

  • S. Wolfram, A New Kind of Science, Wolfram Media Inc., (2002), p. 119.

Crossrefs

a(n) = floor(A037096(n)/(2^(2^(n-1)))). See also A036284, A136386.

Programs

  • Maple
    a(n) := sum( 'bit_n(3^i, n)*(2^i)', 'i'=0..(2^(n-1))-1);
    bit_n := (x, n) -> `mod`(floor(x/(2^n)), 2);

Formula

a(n) = Sum_{k=0..A000225(n-1)} ([A000244(k)/(2^n)] mod 2) * 2^k, where [] stands for floor function.

Extensions

Entry revised Dec 29 2007

A136308 a(n) = (10^2^n - 1)/9.

Original entry on oeis.org

1, 11, 1111, 11111111, 1111111111111111, 11111111111111111111111111111111, 1111111111111111111111111111111111111111111111111111111111111111
Offset: 0

Views

Author

Ctibor O. Zizka, Mar 22 2008

Keywords

Comments

More generally, reading in base B >= 2: a(n) = (B^2^n - 1)/(B-1).
Recurrence: a(n) = a(n-1)*(B^K + 1) and a(0)=1 where K = floor(log_B a(n-1)) + 1.
B = 2 gives A051179; B = 3 gives A059918.

Crossrefs

Cf. A000042 (repunits).
In other bases: A051179, A059918.

Programs

Formula

a(n) = a(n-1)*(10^K + 1) and a(0)=1 where K=floor(log_10 a(n-1)) + 1 = 2^n + 1.
a(n) = A000042(A000079(n)) = A007088(A051179(n)) = A007089(A059918(n)).
A007953(a(n)) = 2^n. - Stefano Spezia, Mar 27 2025

Extensions

Edited by Jason Kimberley, Dec 18 2012

A227722 Smallest Boolean functions from small equivalence classes (counted by A000231).

Original entry on oeis.org

0, 1, 3, 5, 6, 7, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 51, 53, 54, 55, 60, 61, 63, 85, 86, 87, 90, 91, 95, 102, 103, 105, 107, 111, 119, 123, 125, 126, 127, 255, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267
Offset: 0

Views

Author

Tilman Piesk, Jul 22 2013

Keywords

Comments

Two Boolean functions belong to the same small equivalence class (sec) when they can be expressed by each other by negating arguments. E.g., when f(p,~q,r) = g(p,q,r), then f and g belong to the same sec. Geometrically this means that the functions correspond to hypercubes with 2-colored vertices that are equivalent up to reflection (i.e., exchanging opposite hyperfaces).
Boolean functions correspond to integers, so each sec can be denoted by the smallest integer corresponding to one of its functions. There are A000231(n) small equivalence classes of n-ary Boolean functions. Ordered by size they form the finite sequence A_n. It is the beginning of A_(n+1) which leads to this infinite sequence A.

Examples

			The 16 2-ary functions ordered in A000231(2) = 7 small equivalence classes:
a     a(n)    Boolean functions, the left one corresponding to a(n)
0      0      0000
1      1      0001, 0010, 0100, 1000
2      3      0011, 1100
3      5      0101, 1010
4      6      0110, 1001
5      7      0111, 1011, 1101, 1110
6     15      1111
		

Crossrefs

Cf. A227723 (subsequence that does the same thing for big equivalence classes).

Formula

a( A000231 - 1 ) = a(2,6,45,4335...) = 3,15,255,65535... = A051179
a( A000231 ) = a(3,7,46,4336...) = 5,17,257,65537... = A000215

A227723 Smallest Boolean functions from big equivalence classes (counted by A000616).

Original entry on oeis.org

0, 1, 3, 6, 7, 15, 22, 23, 24, 25, 27, 30, 31, 60, 61, 63, 105, 107, 111, 126, 127, 255, 278, 279, 280, 281, 282, 283, 286, 287, 300, 301, 303, 316, 317, 318, 319, 360, 361, 362, 363, 366, 367, 382, 383, 384, 385, 386, 387, 390, 391, 393, 395
Offset: 0

Views

Author

Tilman Piesk, Jul 22 2013

Keywords

Comments

Two Boolean functions belong to the same big equivalence class (bec) when they can be expressed by each other by negating and permuting arguments. E.g., when f(~p,r,q) = g(p,q,r), then f and g belong to the same bec. Geometrically this means that the functions correspond to hypercubes with binarily colored vertices that are equivalent up to rotation and reflection.
Boolean functions correspond to integers, so each bec can be denoted by the smallest integer corresponding to one of its functions. There are A000616(n) big equivalence classes of n-ary Boolean functions. Ordered by size they form the finite sequence A_n. It is the beginning of A_(n+1), which leads to this infinite sequence A.

Examples

			The 16 2-ary functions ordered in A000616(2) = 6 big equivalence classes:
a     a(n)    Boolean functions            hypercube (square)
0      0      0000                         empty
1      1      0001, 0010, 0100, 1000       one in a corner
2      3      0011, 1100, 0101, 1010       ones on a side
3      6      0110, 1001                   ones on a diagonal
4      7      0111, 1011, 1101, 1110       ones in 3 corners
5     15      1111                         full
		

Crossrefs

Cf. A227722 (does the same for small equivalence classes).

Formula

a( A000616 - 1 ) = a(2,5,21,401,...) = 3,15,255,65535,... = A051179

A247282 Run Length Transform of A001317.

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 3, 5, 1, 1, 1, 3, 3, 3, 5, 15, 1, 1, 1, 3, 1, 1, 3, 5, 3, 3, 3, 9, 5, 5, 15, 17, 1, 1, 1, 3, 1, 1, 3, 5, 1, 1, 1, 3, 3, 3, 5, 15, 3, 3, 3, 9, 3, 3, 9, 15, 5, 5, 5, 15, 15, 15, 17, 51, 1, 1, 1, 3, 1, 1, 3, 5, 1, 1, 1, 3, 3, 3, 5, 15, 1, 1, 1, 3, 1, 1, 3, 5, 3, 3, 3, 9, 5, 5, 15, 17, 3, 3, 3, 9, 3, 3, 9, 15, 3, 3, 3, 9, 9, 9, 15, 45
Offset: 0

Views

Author

Antti Karttunen, Sep 22 2014

Keywords

Comments

The Run Length Transform of a sequence {S(n), n>=0} is defined to be the sequence {T(n), n>=0} given by T(n) = Product_i S(i), where i runs through the lengths of runs of 1's in the binary expansion of n. E.g. 19 is 10011 in binary, which has two runs of 1's, of lengths 1 and 2. So T(19) = S(1)*S(2). T(0)=1 (the empty product).
This sequence is obtained by applying Run Length Transform to the right-shifted version of the sequence A001317: 1, 3, 5, 15, 17, 51, 85, 255, 257, ...

Examples

			115 is '1110011' in binary. The run lengths of 1-runs are 2 and 3, thus a(115) = A001317(2-1) * A001317(3-1) = 3*5 = 15.
From _Omar E. Pol_, Feb 15 2015: (Start)
Written as an irregular triangle in which row lengths are the terms of A011782:
1;
1;
1,3;
1,1,3,5;
1,1,1,3,3,3,5,15;
1,1,1,3,1,1,3,5,3,3,3,9,5,5,15,17;
1,1,1,3,1,1,3,5,1,1,1,3,3,3,5,15,3,3,3,9,3,3,9,15,5,5,5,15,15,15,17,51;
...
Right border gives 1 together with A001317.
(End)
		

Crossrefs

Cf. A003714 (gives the positions of ones).
A001316 is obtained when the same transformation is applied to A000079, the powers of two.
Run Length Transforms of other sequences: A071053, A227349, A246588, A246595, A246596, A246660, A246661, A246674, A246685.

Programs

  • Mathematica
    a1317[n_] := FromDigits[ Table[ Mod[Binomial[n-1, k], 2], {k, 0, n-1}], 2];
    Table[ Times @@ (a1317[Length[#]]&) /@ Select[Split[IntegerDigits[n, 2]], #[[1]] == 1&], {n, 0, 100}] (* Jean-François Alcover, Jul 11 2017 *)
  • Python
    # uses RLT function from A278159
    def A247282(n): return RLT(n,lambda m: int(''.join(str(int(not(~(m-1)&k))) for k in range(m)),2)) # Chai Wah Wu, Feb 04 2022

Formula

For all n >= 0, a(A051179(n)) = A246674(A051179(n)) = A051179(n).
Previous Showing 11-20 of 43 results. Next