%I #10 Jan 01 2021 11:23:47
%S 1,3,5,6,8,9,13,15,16,19,20,23,26,27,28,29,31,33,34,35,36,38,45,48,50,
%T 51,53,55,56,59,61,63,64,69,71,73,77,78,83,84,85,86,89,91,93,94,96,
%U 100,101,103,104,108,110,115,119,121,124,127,129,131,133,134,135
%N Numbers k such that A129357(8*k) == 2 (mod 4).
%o (PARI) lista(nn) = my(y=prod(k=1, 8*nn, (1-x)^4*sum(j=1, k, binomial(j+2, 3)*x^(j-1)) + x*O(x^(8*nn)))); for(k=1, nn, if(polcoeff(y, 8*k)%4 == 2, print1(k, ", "))); \\ _Jinyuan Wang_, Jan 01 2021
%Y Cf. A129357.
%K nonn
%O 1,2
%A _Paul D. Hanna_, Apr 11 2007
%E More terms from, name edited and offset changed to 1 by _Jinyuan Wang_, Jan 01 2021
|