f[x_] := Sum[k, {k, 1, x}] x = 0; While[Length[Divisors[f[x]]] < 501, x = x + 1]; x 12375 f[12375] 76576500