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.
%I A356219 #10 Apr 13 2025 01:45:55
%S A356219 284,287,289,292,294,296,299,301,304,306,309,311,313,316,318,321,323,
%T A356219 325,328,330,333,335,337,340,342,345,347,350,352,354,357,359,362,364,
%U A356219 366,369,371,374,376,379,381,383,386,388,391,393,395,398,400
%N A356219 Intersection of A001952 and A003151.
%C A356219 This is the third of four sequences that partition the positive integers. Starting with a general overview, suppose that u = (u(n)) and v = (v(n)) are increasing sequences of positive integers. Let u' and v' be their complements, and assume that the following four sequences are infinite:
%C A356219 (1) u ^ v = intersection of u and v (in increasing order);
%C A356219 (2) u ^ v';
%C A356219 (3) u' ^ v;
%C A356219 (4) u' ^ v'.
%C A356219 Every positive integer is in exactly one of the four sequences.
%C A356219 Assume that if w is any of the sequences u, v, u', v', then lim_{n->oo} w(n)/n exists and defines the (limiting) density of w. For w = u,v,u',v', denote the densities by r,s,r',s'. Then the densities of sequences (1)-(4) exist, and 1/(r*r') + 1/(r*s') + 1/(s*s') + 1/(s*r') = 1.
%C A356219 For A356219, u, v, u', v', are the Beatty sequences given by u(n) = floor(n*sqrt(2)) and v(n) = floor(((1+sqrt(2))/2)*n), so that r = sqrt(2), s = (1+sqrt(2))/2, r' = (2+sqrt(2))/2, s' = 1 + 1/sqrt(2).
%e A356219 (1) u ^ v = (2, 4, 7, 9, 12, 14, 16, 19, 21, 24, 26, 28, 31, 33, ...) = A003151.
%e A356219 (2) u ^ v' = (1, 5, 8, 11, 15, 18, 22, 25, 29, 32, 35, 39, 42, ...) = A001954.
%e A356219 (3) u' ^ v = (284, 287, 289, 292, 294, 296, 299, 301, 304, 306, ...) = A356219.
%e A356219 (4) u' ^ v' = (3, 6, 10, 13, 17, 20, 23, 27, 30, 34, 37, 40, 44, ...) = A003152.
%t A356219 z = 200;
%t A356219 r = Sqrt[2]; u = Table[Floor[n*r], {n, 1, z}] (* A001951 *)
%t A356219 u1 = Take[Complement[Range[1000], u], z] (* A001952 *)
%t A356219 r1 = 1 + Sqrt[2]; v = Table[Floor[n*r1], {n, 1, z}] (* A003151 *)
%t A356219 v1 = Take[Complement[Range[1000], v], z] (* A003152 *)
%t A356219 t1 = Intersection[u, v] (* A003151 *)
%t A356219 t2 = Intersection[u, v1] (* A001954 *)
%t A356219 t3 = Intersection[u1, v] (* A356219 *)
%t A356219 t4 = Intersection[u1, v1] (* A001952 *)
%Y A356219 Cf. A001951, A001952, A003151, A003152, A001954, A184922 (results of compositions instead of intersections), A341239 (reversed compositions).
%K A356219 nonn,easy
%O A356219 1,1
%A A356219 _Clark Kimberling_, Nov 13 2022