cp's OEIS Frontend

A352706 G.f. A(x) satisfies: (1 - x*A(x))^7 = 1 - 7*x - x^7*A(x^7).

%I A352706 #9 Sep 01 2022 12:04:34
%S A352706 1,3,13,65,351,1989,11650,69900,427167,2648438,16612947,105215448,
%T A352706 671760933,4318468134,27926126553,181520036178,1185220461867,
%U A352706 7769787812787,51117085998498,337373170647840,2233091755252871,14819626692452231,98582852467595847
%N A352706 G.f. A(x) satisfies: (1 - x*A(x))^7 = 1 - 7*x - x^7*A(x^7).
%C A352706 Essentially an unsigned version of A352705 (after dropping the initial term).
%H A352706 Paul D. Hanna, <a href="/A352706/b352706.txt">Table of n, a(n) for n = 0..1000</a>
%F A352706 G.f. A(x) satisfies:
%F A352706 (1) (1 + x*A(-x))^7 = 1 + 7*x + x^7*A(-x^7).
%F A352706 (2) A(x) = (1 - (1 - 7*x - x^7*A(x^7))^(1/7))/x.
%F A352706 (3) A(x)^7 = A(x^7) (mod 7).
%e A352706 G.f.: A(x) = 1 + 3*x + 13*x^2 + 65*x^3 + 351*x^4 + 1989*x^5 + 11650*x^6 + 69900*x^7 + 427167*x^8 + 2648438*x^9 + 16612947*x^10 + ...
%e A352706 where
%e A352706 (1 - x*A(x))^7 = 1 - 7*x - x^7 - 3*x^14 - 13*x^21 - 65*x^28 - 351*x^35 - 1989*x^42 - 11650*x^49 - 69900*x^56 - 427167*x^63 - 2648438*x^70 + ...
%e A352706 also
%e A352706 (1 - 7*x - x^7*A(x^7))^(1/7) = 1 - x - 3*x^2 - 13*x^3 - 65*x^4 - 351*x^5 - 1989*x^6 - 11650*x^7 - 69900*x^8 - 427167*x^9 - 2648438*x^10 + ...
%e A352706 which equals 1 - x*A(x).
%o A352706 (PARI) {a(n) = my(A=1+3*x); for(i=1,n,
%o A352706 A = (1 - (1 - 7*x - x^7*subst(A,x,x^7) + x*O(x^(n+1)))^(1/7))/x);
%o A352706 polcoeff(A,n)}
%o A352706 for(n=0,30,print1(a(n),", "))
%Y A352706 Cf. A352705, A352702, A352704.
%K A352706 nonn
%O A352706 0,2
%A A352706 _Paul D. Hanna_, Mar 29 2022