login
Fourth row of number array A084061.
5

%I #6 Sep 08 2022 08:45:10

%S 1,1,7,54,553,7100,109863,1991752,41426257,972602640,25447064743,

%T 734276705888,23166635069241,793426715543488,29316839407495111,

%U 1162492244159875200,49240280161094287777,2218952409252783579392

%N Fourth row of number array A084061.

%H G. C. Greubel, <a href="/A084065/b084065.txt">Table of n, a(n) for n = 0..300</a>

%F a(n) = ( (n - sqrt(3))^n + (n + sqrt(3))^n )/2.

%p seq( round(((n-sqrt(3))^n + (n+sqrt(3))^n)/2), n=0..20); # _G. C. Greubel_, Jan 11 2020

%t Table[Round[((n+Sqrt[3])^n + (n-Sqrt[3])^n)/2], {n,0,20}] (* _G. C. Greubel_, Jan 11 2020 *)

%o (PARI) vector(21, n, round(((n-1-sqrt(3))^(n-1) + (n-1+sqrt(3))^(n-1))/2) ) \\ _G. C. Greubel_, Jan 11 2020

%o (Magma) [Round(((n-Sqrt(3))^n + (n+Sqrt(3))^n)/2): n in [0..20]]; // _G. C. Greubel_, Jan 11 2020

%o (Sage) [round(((n-sqrt(3))^n + (n+sqrt(3))^n)/2) for n in (0..20)] # _G. C. Greubel_, Jan 11 2020

%o (GAP) List([0..20], n-> ((n-Sqrt(3))^n + (n+Sqrt(3))^n)/2); # _G. C. Greubel_, Jan 11 2020

%Y Cf. A084061, A084062, A084063, A084064, A084095, A084096.

%K easy,nonn

%O 0,3

%A _Paul Barry_, May 11 2003