login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A166447
a(n) = n*round(sqrt(n)).
3
0, 1, 2, 6, 8, 10, 12, 21, 24, 27, 30, 33, 36, 52, 56, 60, 64, 68, 72, 76, 80, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 186, 192, 198, 204, 210, 216, 222, 228, 234, 240, 246, 252, 301, 308, 315, 322, 329, 336, 343, 350, 357, 364, 371, 378, 385, 392, 456
OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..999 [Offset shifted by Georg Fischer, Aug 30 2021]
FORMULA
a(n) = n*A000194(n). - R. J. Mathar, Nov 18 2009
MAPLE
seq(n*round(sqrt(n)), n=0..60); # Georg Fischer, Aug 28 2021
# second Maple program:
a:= n-> n*isqrt(n):
seq(a(n), n=0..60); # Alois P. Heinz, Aug 29 2021
MATHEMATICA
Table[n Round[Sqrt[n]], {n, 0, 60}] (* Harvey P. Dale, Oct 30 2021 *)
PROG
(Magma) [n*Round(Sqrt(n)): n in [0..60]]; // Vincenzo Librandi, May 15 2016
(Python)
from math import isqrt
def A166447(n): return n*((m:=isqrt(n))+int(n-m*(m+1)>=1)) # Chai Wah Wu, Jun 19 2024
CROSSREFS
Sequence in context: A049637 A284753 A258663 * A075332 A141105 A047395
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Offset corrected by Georg Fischer, Aug 28 2021
STATUS
approved