# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a049068
Showing 1-1 of 1

%I A049068 #65 Jul 28 2022 09:07:10
%S A049068 3,5,7,8,10,11,13,14,15,17,18,19,21,22,23,24,26,27,28,29,31,32,33,34,
%T A049068 35,37,38,39,40,41,43,44,45,46,47,48,50,51,52,53,54,55,57,58,59,60,61,
%U A049068 62,63,65,66,67,68,69,70,71,73,74,75,76,77,78,79,80,82,83,84,85,86,87
%N A049068 Complement of quarter-squares (A002620).
%C A049068 Intersection of A000037 and A078358. - _Reinhard Zumkeller_, May 08 2012
%C A049068 Numbers k such that floor(sqrt(k)+1/2) does not divide k. - _Wesley Ivan Hurt_, Dec 01 2020
%H A049068 Reinhard Zumkeller, <a href="/A049068/b049068.txt">Table of n, a(n) for n = 1..10000</a>
%F A049068 a(n) = n + A027434(n).
%F A049068 Other identities and observations. For all n >= 1:
%F A049068 A237347(a(n)) = 2. - _Reinhard Zumkeller_, Mar 18 2014
%F A049068 A240025(a(n)) = 0. - _Reinhard Zumkeller_, Jul 05 2014
%F A049068 a(n) = A080037(n) - 1. - _Peter Kagey_, Dec 08 2015
%F A049068 G.f.: x/(1-x)^2 + Sum_{k>=0} (x^(1+k^2)*(1+x^k))/(1-x)
%F A049068 = (x*Theta3(x)+ x^(3/4)*Theta2(x))/(2-2*x) + (3-x)*x/(2*(1-x)^2) where Theta3 and Theta2 are Jacobi Theta functions. - _Robert Israel_, Dec 09 2015
%F A049068 a(n) = A000037(A000037(n)). - _Gerald Hillier_, Dec 01 2017
%p A049068 A049068:=n->n + ceil(2*sqrt(n)); seq(A049068(n), n=1..100); # _Wesley Ivan Hurt_, Mar 01 2014
%t A049068 max = 100; Complement[Range[0, max], Table[Quotient[n^2, 4], {n, 0, 2*Sqrt[max]}]]  (* _Jean-François Alcover_, Apr 18 2013 *)
%t A049068 Table[n + Ceiling[2 * Sqrt[n]], {n, 100}] (* _Wesley Ivan Hurt_, Mar 01 2014 *)
%o A049068 (PARI) {a(n) = if( n<1, 0, n+1 + sqrtint(4*n - 3))} /* _Michael Somos_, Oct 16 2006 */
%o A049068 (Haskell)
%o A049068 a049068 n = a049068_list !! (n-1)
%o A049068 a049068 = filter ((== 0) . a240025) [0..]
%o A049068 -- _Reinhard Zumkeller_, Jul 05 2014, Mar 18 2014, May 08 2012
%o A049068 (Magma) [n+Ceiling(2*Sqrt(n)): n in [1..70]]; // _Vincenzo Librandi_, Dec 09 2015
%o A049068 (Python)
%o A049068 from math import isqrt
%o A049068 def A049068(n): return n+1+isqrt((n<<2)-1) # _Chai Wah Wu_, Jul 27 2022
%Y A049068 Complement: A002620.
%Y A049068 Cf. A000037, A027434, A078358, A080037, A237347, A240025.
%K A049068 nonn,nice,easy
%O A049068 1,1
%A A049068 _Michael Somos_, Aug 06 1999

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE