|
%I
%S 0,1,5,16,53,145,441,1165,3325,8840,24268,63542,172420,446485,1182329,
%T 3055461,7966038,20351233,52507157,133014971,339170500,853849745,
%U 2154442760,5382418592,13481157491,33449252747,83091028844,205013093570,505799029994,1240477207822,3042433058325
%N Row sums of A114994.
%e A114994 begins:
%e 0;
%e 1; => row sum 1;
%e 2,3; => row sum 5;
%e 4,5,7; => row sum 16;
%e 8,9,10,11,15; => row sum 53.
%p isA114994aux := proc(n) local k; for k from 0 do fk1 := 4^(k+1) ; if 2^k mod fk1 = n mod fk1 then return true; elif (2^k mod fk1) > n then return false; end if; end do: end proc:
%p isA114994 := proc(n) option remember ; if n = 0 then return true; elif type(n,'odd') then if n <= 11 then true; else return procname((n-1)/2) ; fi; else return ( isA114994(n/2) and isA114994aux(n/2) ); end if; end proc:
%p A173871 := proc(n) a := 0 ; for k from 2^n to 2^(n+1)-1 do if isA114994(k) then a := a+k; end if; end do; a ; end proc:
%p seq(A173871(n),n=0..19) ; # _R. J. Mathar_, Mar 26 2010
%Y Cf. A114994, A171426.
%K nonn
%O 0,3
%A _Alford Arnold_, Mar 05 2010
%E a(9)-a(19) from _R. J. Mathar_, Mar 26 2010
%E a(0) prepended by and more terms from _Jinyuan Wang_, Feb 23 2020
|
