|
%I #17 Jul 28 2022 12:17:31
%S 1,3,5,6,7,9,10,11,13,14,15,16,17,19,20,21,22,23,25,26,27,28,29,30,31,
%T 33,34,35,36,37,38,39,41,42,43,44,45,46,47,48,49,51,52,53,54,55,56,57,
%U 58,59,61,62,63,64,65,66,67,68,69,70,71,73,74,75,76,77,78,79,80,81,82,83,85,86,87,88,89,90,91,92
%N Complement of A007590, except for initial zeros.
%C A245575(a(n)) mod 2 = 0, or, where even terms occur in A245575; a(A078633(n)) mod 2 = 0. - _Reinhard Zumkeller_, Aug 05 2014
%H Reinhard Zumkeller, <a href="/A182834/b182834.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = n + floor(sqrt(2n-2)).
%t a=2; b=-2;
%t Table[n+Floor[(a*n+b)^(1/2)],{n,80}]
%t Table[n-1+Ceiling[(n*n-b)/a],{n,50}]
%o (Haskell)
%o a182834 n = a000196 (2 * n - 2) + n -- _Reinhard Zumkeller_, Aug 05 2014
%o (Python)
%o from math import isqrt
%o def A182834(n): return n+isqrt(n-1<<1) # _Chai Wah Wu_, Jul 28 2022
%Y Cf. A007590.
%Y Cf. A000196, A078633, A245575.
%K nonn,easy
%O 1,2
%A _Clark Kimberling_, Jan 07 2011
|
