(MAGMAMagma) k:=30; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 25 2019
(MAGMAMagma) k:=30; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 25 2019
editing
approved
nonn,easy
approved
editing
editing
approved
Join[{1}, NestList[29#&, 30, 20]] (* Harvey P. Dale, Aug 27 2020 *)
approved
editing
proposed
approved
editing
proposed
editing
proposed
GExpansion of g.f.: (1+x)/(1-29*x).
a(n) = Sum_{k=0..n} A097805(n,k)*(-1)^(n-k)*30^k. [_- _Philippe Deléham_, Dec 04 2009]
a(0) = 1; for n > 0, a(n) = 30*29^(n-1). [_- _Vincenzo Librandi_, Dec 05 2009]
E.g.f.: (30*exp(29*x) -1)/29. - G. C. Greubel, Sep 25 2019
k:=30; seq(`if`(n=0, 1, k*(k-1)^(n-1)), n = 0..25); # G. C. Greubel, Sep 25 2019
With[{k=30}, Table[If[n==0, 1, k*(k-1)^(n-1)], {n, 0, 25}]] (* G. C. Greubel, Sep 25 2019 *)
(Python) for i in range(100131):print(i, 30*29**(i-1) if i>0 else 1) # Kenny Lau, Aug 03 2017
(PARI) vector(26, n, k=30; if(n==1, 1, k*(k-1)^(n-2))) \\ G. C. Greubel, Sep 25 2019
(MAGMA) k:=30; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 25 2019
(Sage) k=30; [1]+[k*(k-1)^(n-1) for n in (1..25)] # G. C. Greubel, Sep 25 2019
(GAP) k:=30;; Concatenation([1], List([1..25], n-> k*(k-1)^(n-1) )); # G. C. Greubel, Sep 25 2019
Cf. A003945.
approved
editing