login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Complement of A007590, except for initial zeros.
5

%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