|
%I #19 Dec 07 2021 18:16:54
%S 0,0,1,1,2,3,6,10,18,31,56,96,172,299,530,929,1646,2893,5126,9044,
%T 16028,28362,50328,89249,158598,281830,501538,892857,1591282,2837467,
%U 5064334,9044023,16163946,28906213,51729844,92628401,165967884,297541263,533731692,957921314
%N Inverse Euler transform of the tribonacci numbers A000073.
%p read transforms; # https://oeis.org/transforms.txt
%p arow := len -> EULERi([seq(A000073(n), n = 0..len)]): arow(39);
%p # second Maple program:
%p t:= n-> (<<0|1|0>, <0|0|1>, <1|1|1>>^n)[1, 3]:
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
%p add(binomial(a(i)+j-1, j)*b(n-i*j, i-1), j=0..n/i)))
%p end:
%p a:= proc(n) option remember; t(n-1)-b(n, n-1) end:
%p seq(a(n), n=1..40); # _Alois P. Heinz_, Dec 05 2021
%t (* EulerInvTransform is defined in A022562. *)
%t EulerInvTransform[LinearRecurrence[{1, 1, 1}, {0, 0, 1}, 40]]
%o (SageMath)
%o def euler_invtrans(A) :
%o L = []; M = []
%o for i in range(len(A)) :
%o s = (i+1)*A[i] - sum(L[j-1]*A[i-j] for j in (1..i))
%o L.append(s)
%o s = sum(moebius((i+1)/d)*L[d-1] for d in divisors(i+1))
%o M.append(s/(i + 1))
%o return M
%o @cached_function
%o def a(n): return a(n-1) + a(n-2) + a(n-3) if n > 2 else [0,0,1][n]
%o print(euler_invtrans([a(n) for n in range(40)]))
%o (Python) # After the Maple program of _Alois P. Heinz_.
%o from functools import cache
%o from math import comb
%o def binomial(n, k):
%o if n == -1: return 1
%o return comb(n, k)
%o @cache
%o def A000073(n):
%o if n <= 1: return 0
%o if n == 2: return 1
%o return A000073(n-1) + A000073(n-2) + A000073(n-3)
%o @cache
%o def b(n, i):
%o if n == 0: return 1
%o if i < 1: return 0
%o return sum(binomial(a(i) + j - 1, j) *
%o b(n - i * j, i - 1) for j in range(1 + n // i))
%o @cache
%o def a(n): return (A000073(n - 1) - b(n, n - 1))
%o print([a(n) for n in range(1, 41)])
%o (PARI)
%o InvEulerT(v)={my(p=log(1+x*Ser(v))); dirdiv(vector(#v,n,polcoef(p,n)), vector(#v,n,1/n))}
%o seq(n) = InvEulerT(Vec(x^2/(1 - x - x^2 - x^3) + O(x^n), -n)) \\ _Andrew Howroyd_, Dec 05 2021
%Y Column k=2 of A349802.
%Y Cf. A000073, A057597 (tribonacci numbers for n <= 0), A006206 and A060280.
%Y Cf. A349903, A349977.
%K nonn
%O 1,5
%A _Peter Luschny_, Dec 05 2021
|
