OFFSET
1,2
COMMENTS
If [M/k^2] is changed to [M/k] we get A100617.
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
10 -> 10 - [10/4] = 8 -> 8 - [8/9] = 8, which is now fixed, so a(10) = 8.
MAPLE
f:=proc(n) local i, j, k; k:=n; for i from 2 to 10000 do j := floor(k/(i^2)); if j < 1 then break; fi; k := k-j; od; k; end;
MATHEMATICA
a[n_] := (k = 2; FixedPoint[# - Floor[#/(k++)^2]&, n]); Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Feb 10 2018 *)
PROG
(Haskell)
a100618 n = f 2 n where
f k n | n' == 0 = n
| otherwise = f (k+1) (n-n') where n' = div n (k^2)
-- Reinhard Zumkeller, Sep 15 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Dec 03 2004
STATUS
approved