for i := 2 to n do
    begin
      min:=maxlongint;
      for j := 1 to trunc(sqrt(i)) do
        if min>f[i-j*j] then
          min:=f[i-j*j];
      f[i]:=min+1;
    end;