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.

A333494 a(1) = 1; a(n) = Sum_{k=1..n-1} ceiling(n/k) * a(k).

This page as a plain text file.
%I A333494 #8 Mar 25 2020 11:35:45
%S A333494 1,2,7,22,69,208,634,1903,5734,17210,51702,155107,465561,1396684,
%T A333494 4190689,12572144,37718360,113155081,339471195,1018413586,3055258062,
%U A333494 9165774828,27497376189,82492128568,247476542954,742429628932,2227289352360,6681868062822,20045605585809
%N A333494 a(1) = 1; a(n) = Sum_{k=1..n-1} ceiling(n/k) * a(k).
%F A333494 G.f. A(x) satisfies: A(x) = x * (1 + (1/(1 - x)) * (A(x) + Sum_{k>=1} A(x^k))).
%F A333494 a(n) ~ c * 3^n, where c = 0.292080665386646518390576592052254840432101999262173908555857806023213143845... - _Vaclav Kotesovec_, Mar 25 2020
%t A333494 a[1] = 1; a[n_] := a[n] = Sum[Ceiling[n/k] a[k], {k, 1, n - 1}]; Table[a[n], {n, 1, 29}]
%t A333494 terms = 29; A[_] = 0; Do[A[x_] = x (1 + (1/(1 - x)) (A[x] + Sum[A[x^k], {k, 1, terms}])) + O[x]^(terms + 1) // Normal, terms + 1]; CoefficientList[A[x], x] // Rest
%Y A333494 Cf. A006590, A014668, A126656, A332846.
%K A333494 nonn
%O A333494 1,2
%A A333494 _Ilya Gutkovskiy_, Mar 24 2020