%I #20 Sep 18 2024 16:33:03
%S 0,1,2,4,6,8,9,10,12,14,16,18,20,22,24,25,26,28,30,32,34,36,38,40,42,
%T 44,46,48,49,50,52,54,56,58,60,62,64,66,68,70,72,74,76,78,80,81,82,84,
%U 86,88,90,92,94,96,98,100,102,104,106,108,110,112,114,116,118,120,121,122
%N Union of squares and the even numbers.
%C Numbers k such that k^k is a square. - _Chai Wah Wu_, Sep 18 2024
%t With[{upto=150},Union[Join[Range[0,upto,2],Range[Sqrt[upto]]^2]]] (* _Harvey P. Dale_, Apr 29 2015 *)
%o (PARI) is_A176693(n)=!bittest(n, 0)||issquare(n) \\ _M. F. Hasler_, Apr 19 2015
%o (Python)
%o from math import isqrt
%o def A176693(n):
%o def f(x): return n-1+(x>>1)+(x&1)-(isqrt(x)+1>>1)
%o m, k = n-1, f(n-1)
%o while m != k: m, k = k, f(k)
%o return m # _Chai Wah Wu_, Sep 18 2024
%Y Cf. A128201 (Union of squares and the odd numbers).
%K nonn,easy
%O 1,3
%A _Vincenzo Librandi_, Apr 24 2010
%E Initial a(1)=0 added by _M. F. Hasler_, Apr 20 2015