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.

A372946 Numbers k that divide the k-th NSW number.

This page as a plain text file.
%I A372946 #5 May 18 2024 01:49:25
%S A372946 1,7,217,3937,6727,6847,51943,170671,330337,385687,2484247,2566537,
%T A372946 2904007,3020857,3696967,6465577,9405337,12021439,19384207
%N A372946 Numbers k that divide the k-th NSW number.
%C A372946 Numbers k such that k | A002315(k).
%e A372946 7 is a term since A002315(7) = 275807 = 7 * 39401 is divisible by 7.
%t A372946 seq[kmax_] := Module[{nsw0 = 1, nsw1 = 7, nsw2, s = {1}}, Do[nsw2 = 6*nsw1 - nsw0; If[Divisible[nsw2, k], AppendTo[s, k]]; nsw0 = nsw1; nsw1 = nsw2, {k, 2, kmax}]; s]; seq[52000]
%o A372946 (PARI) lista(kmax) = {my(nsw0 = 1, nsw1 = 7, nsw2); print1("1, "); for(k = 2, kmax, nsw2 = 6*nsw1 - nsw0; if(!(nsw2 % k), print1(k, ", ")); nsw0 = nsw1; nsw1 = nsw2);}
%Y A372946 Cf. A002315, A330276.
%Y A372946 Similar sequences: A014847 (Catalan), A016089 (Lucas), A023172 (Fibonacci), A051177 (partition), A232570 (tribonacci), A246692 (Pell), A266969 (Motzkin).
%K A372946 nonn,more
%O A372946 1,2
%A A372946 _Amiram Eldar_, May 17 2024