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.

A353199 a(1)=a(2)=1. a(n) = the sum of the indices of the terms, from among the first (n-2) terms of the sequence, which divide a(n-1).

This page as a plain text file.
%I A353199 #16 Jun 05 2022 08:27:00
%S A353199 1,1,1,3,6,10,6,15,10,12,22,6,22,17,6,34,20,21,10,21,28,6,49,6,71,6,
%T A353199 95,6,121,6,149,6,179,6,211,6,245,29,6,281,6,320,57,10,40,101,6,361,6,
%U A353199 408,497,31,6,457,6,510,681,10,84,634,6,565,6,626,6,689,6,754,44,30,965
%N A353199 a(1)=a(2)=1. a(n) = the sum of the indices of the terms, from among the first (n-2) terms of the sequence, which divide a(n-1).
%e A353199 a(10) = 12. Of the terms before that, a(1), a(2), a(3), a(4), a(5), and a(7) divide 12. Hence, a(11) = 1+2+3+4+5+7 = 22.
%t A353199 a[1] = a[2] = 1; a[n_] := a[n] = Total[Select[Range[n - 2], Divisible[a[n - 1], a[#]] &]]; Array[a, 100] (* _Amiram Eldar_, Apr 30 2022 *)
%o A353199 (PARI) first(n) = { n = max(n, 2); res = List([1,1]); for(i = 3, n, v = Vec(select(x->res[#res]%x == 0, res, 1))[^-1]; listput(res, vecsum(v)) ); res } \\ _David A. Corneth_, Apr 30 2022
%Y A353199 Cf. A124063.
%K A353199 nonn
%O A353199 1,4
%A A353199 _Eric Fox_, Apr 30 2022