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 A331533 #12 Jan 20 2020 12:54:39 %S A331533 0,0,1,0,2,0,1,3,0,4,0,1,3,4,5,0,2,6,0,1,4,7,0,8,0,1,4,8,9,0,2,6,8,10, %T A331533 0,1,3,4,5,7,8,9,11,0,4,12,0,1,3,13,0,2,8,14,0,1,8,15,0,16,0,1,16,17, %U A331533 0,2,8,16,18,0,1,3,8,16,17,19,0,4,12,16,20 %N A331533 Irregular table read by rows; for n >= 0, the n-th row corresponds to the nonnegative integers k such that (n^2) AND (k^2) = k^2, in ascending order (where AND denotes the bitwise AND operator). %C A331533 The n-th row has A331532(n) terms, leading term 0 and last term n. %H A331533 Rémy Sigrist, <a href="/A331533/b331533.txt">Table of n, a(n) for n = 0..7118</a> (rows for n = 0..512) %e A331533 Table begins: %e A331533 0; %e A331533 0, 1; %e A331533 0, 2; %e A331533 0, 1, 3; %e A331533 0, 4; %e A331533 0, 1, 3, 4, 5; %e A331533 0, 2, 6; %e A331533 0, 1, 4, 7; %e A331533 0, 8; %e A331533 ... %o A331533 (PARI) row(n) = select(k -> bitand(n^2, k^2)==k^2, [0..n]) %Y A331533 Cf. A295989, A331532 (row lengths). %K A331533 nonn,tabf,base %O A331533 0,5 %A A331533 _Rémy Sigrist_, Jan 19 2020