|
%I #11 Jul 13 2022 14:45:42
%S 3,4,6,7,9,10,11,13,14,15,16,18,19,20,21,23,24,25,26,27,29,30,31,32,
%T 33,34,36,37,38,39,40,41,43,44,45,46,47,48,49,51,52,53,54,55,56,57,58,
%U 60,61,62,63,64,65,66,67,69,70,71,72,73,74,75,76,77,79,80,81,82,83,84,85,86,87,88,90,91,92,93,94,95,96,97,98,99,101,102,103,104,105,106
%N a(n) = n+floor(sqrt(3n+1)); complement of A183861.
%F a(n)=n+floor(sqrt(3n+1)).
%t a=3;b=1;
%t Table[n+Floor[(a*n+b)^(1/2)],{n,90}] (* A183860 *)
%t Table[n-1+Ceiling[(n*n-b)/a],{n,70}] (* A183861 *)
%o (PARI) a(n) = n + sqrtint(3*n+1); \\ _Michel Marcus_, Jul 13 2022
%Y Cf. A183861.
%K nonn,easy
%O 1,1
%A _Clark Kimberling_, Jan 07 2011
|
