|
%I #19 Mar 06 2023 20:40:23
%S 0,3,8,26,32,46,54,92,115,199,268,306,326,368,390,436,796,828,861,963,
%T 1071,1146,1304,1345,1429,1472,1560,2193,2246,2300,2465,2636,2873,
%U 2995,3120,3312,3377,3443,3713,3782,3922,3993,4137,4738,4816,5297,6062,6239
%N Numbers k such that floor(Pi*k) is a square.
%H Harry J. Smith, <a href="/A061812/b061812.txt">Table of n, a(n) for n = 0..1000</a>
%e Floor(Pi*861) = 2704 = 52^2, so 861 is a term;
%e floor(Pi*1146) = 3600 = 60^2, so 1146 is a term.
%t Select[Range[0,7000],IntegerQ[Sqrt[Floor[Pi #]]]&] (* _Harvey P. Dale_, May 14 2011 *)
%o (PARI) for(n=0,1200, if(issquare(floor(Pi*n)),print(n)))
%o (PARI) { n=-1; for (a=0, 3714376, if (issquare(floor(Pi*a)), write("b061812.txt", n++, " ", a)) ) } \\ _Harry J. Smith_, Jul 28 2009
%Y Cf. A022844.
%K easy,nice,nonn
%O 0,2
%A _Jason Earls_, Jul 01 2001
%E More terms from _Naohiro Nomoto_, Oct 04 2001
|
