|
%I #16 Sep 08 2022 08:45:28
%S 22,471260364628084305,6457022669043550542502557676,
%T 105149403852520725445003265581519105,
%U 41911381174488637014293971538580334000626
%N a(n) = Sum_{k=1..21} n^A001358(k).
%H G. C. Greubel, <a href="/A123665/b123665.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = 1 +n^4 +n^6 +n^9 +n^10 +n^14 +n^15 +n^21 +n^22 +n^25 +n^26 + n^33 +n^34 +n^35 +n^38 +n^39 +n^46 +n^49 +n^51 +n^55 +n^57 +n^58.
%t Table[1 + n^4 + n^6 + n^9 + n^10 + n^14 + n^15 + n^21 + n^22 + n^25 +
%t n^26 + n^33 + n^34 + n^35 + n^38 + n^39 + n^46 + n^49 + n^51 +
%t n^55 + n^57 + n^58, {n, 1, 50}] (* _G. C. Greubel_, Oct 26 2017 *)
%o (PARI) for(n=1,50, print1(1 + n^4 + n^6 + n^9 + n^10 + n^14 + n^15 + n^21 + n^22 + n^25 + n^26 + n^33 + n^34 + n^35 + n^38 + n^39 + n^46 + n^49 + n^51 + n^55 + n^57 + n^58, ", ")) \\ _G. C. Greubel_, Oct 26 2017
%o (Magma) [1 + n^4 + n^6 + n^9 + n^10 + n^14 + n^15 + n^21 + n^22 + n^25 +
%o n^26 + n^33 + n^34 + n^35 + n^38 + n^39 + n^46 + n^49 + n^51 + n^55 + n^57 + n^58: n in [1..50]]; // _G. C. Greubel_, Oct 26 2017
%Y Cf. A001358, A002523, A091940, A123177, A123656-A123665.
%K easy,nonn,less
%O 1,1
%A _Jonathan Vos Post_, Oct 04 2006
%E Better name from _Joerg Arndt_, May 23 2021
| 