Original entry on oeis.org
1, 6, 9, 15, 19, 24, 27, 31, 36, 42, 45, 51, 55, 60, 66, 69, 73, 78, 81, 87, 91, 96, 99, 103, 108, 114, 117, 121, 126, 129, 135, 139, 144, 150, 153, 159, 163, 168, 171, 175, 180, 186, 189, 195, 199, 204, 210, 213, 217, 222, 225, 231, 235, 240, 246, 249, 255
Offset: 1
(1) u ^ v = (3, 5, 8, 10, 12, 14, 16, 18, 21, 23, 26, 28, 30, 33, ...) = A360394
(2) u ^ v' = (1, 6, 9, 15, 19, 24, 27, 31, 36, 42, 45, 51, 55, 60, ...) = A360395
(3) u' ^ v = (7, 11, 17, 20, 25, 29, 32, 38, 43, 47, 53, 56, 62, ...) = A360396
(4) u' ^ v' = (2, 4, 13, 22, 34, 40, 49, 58, 64, 76, 85, 94, 106, ...) = A360397
-
z = 400;
u = Accumulate[1 + ThueMorse /@ Range[0, z]]; (* A026430 *)
u1 = Complement[Range[Max[u]], u]; (* A356133 *)
v = u + 2 ; (* A360392 *)
v1 = Complement[Range[Max[v]], v]; (* A360393 *)
Intersection[u, v] (* A360394 *)
Intersection[u, v1] (* A360395 *)
Intersection[u1, v] (* A360396 *)
Intersection[u1, v1] (* A360397 *)
Original entry on oeis.org
1, 3, 6, 9, 14, 19, 23, 28, 33, 36, 41, 46, 51, 54, 60, 63, 68, 73, 77, 82, 87, 90, 96, 99, 103, 109, 114, 117, 121, 128, 130, 136, 141, 144, 149, 154, 159, 162, 168, 171, 175, 181, 186, 189, 194, 199, 203, 209, 213, 216, 222, 225, 230, 235, 239, 245, 249
Offset: 1
(1) u o v = (5, 8, 10, 12, 15, 16, 18, 21, 24, 26, 27, 30, 31, 35, 37, 39, ...) = A360398
(2) u o v' = (1, 3, 6, 9, 14, 19, 23, 28, 33, 36, 41, 46, 51, 54, 60, 63, 68, ...) = A360399
(3) u' o v = (7, 13, 20, 22, 29, 32, 34, 40, 47, 49, 53, 58, 62, 67, 74, 76, ...) = A360400
(4) u' o v' = (2, 4, 11, 17, 25, 38, 43, 56, 64, 71, 79, 92, 101, 106, 119, ...) = A360401
-
z = 2000;
u = Accumulate[1 + ThueMorse /@ Range[0, 600]]; (* A026430 *)
u1 = Complement[Range[Max[u]], u]; (* A356133 *)
v = u + 2; (* A360392 *)
v1 = Complement[Range[Max[v]], v]; (* A360393 *)
zz = 100;
Table[u[[v[[n]]]], {n, 1, zz}] (* A360398 *)
Table[u[[v1[[n]]]], {n, 1, zz}] (* A360399 *)
Table[u1[[v[[n]]]], {n, 1, zz}] (* A360400 *)
Table[u1[[v1[[n]]]], {n, 1, zz}] (* A360401 *)
Original entry on oeis.org
1, 4, 9, 13, 19, 22, 24, 31, 36, 40, 42, 49, 51, 58, 64, 66, 73, 76, 78, 85, 87, 94, 99, 103, 106, 112, 117, 121, 126, 129, 133, 139, 144, 148, 150, 157, 159, 166, 171, 175, 178, 184, 189, 193, 199, 202, 204, 210, 213, 220, 225, 229, 235, 238, 240, 246, 253
Offset: 1
(1) v o u = (3, 7, 10, 11, 14, 16, 17, 20, 23, 25, 26, 29, 30, 33, 37, ...) = A360402
(2) v' o u = (1, 4, 9, 13, 19, 22, 24, 31, 36, 40, 42, 49, 51, 58, 64, ...) = A360403
(3) v o u' = (5, 8, 12, 18, 21, 28, 32, 35, 39, 46, 50, 53, 59, 62, 67, ...) = A360404
(4) v' o u' = (2, 6, 15, 27, 34, 45, 55, 60, 69, 81, 91, 96, 108, 114, ...) = A360405
-
z = 2000; zz = 100;
u = Accumulate[1 + ThueMorse /@ Range[0, 600]]; (* A026430 *)
u1 = Complement[Range[Max[u]], u]; (* A356133 *)
v = u + 2; (* A360392 *)
v1 = Complement[Range[Max[v]], v]; (* A360393 *)
Table[v[[u[[n]]]], {n, 1, zz}] (* A360402 *)
Table[v1[[u[[n]]]], {n, 1, zz}] (* A360403 *)
Table[v[[u1[[n]]]], {n, 1, zz}] (* A360404 *)
Table[v1[[u1[[n]]]], {n, 1, zz}] (* A360405 *)
-
def A360393(n):
if n < 3: return [0, 1, 2][n]
return 3*n - 5 - (n-3).bit_count() % 2
def A026430(n): return n+(n-1>>1)+(n-1&1|(n.bit_count()&1^1))
def A360403(n): return A360393(A026430(n)) # Winston de Greef, Mar 24 2023
Comments