Divisores

x1 = x1 / m ‘no hay que usar un IF porque es divisible con seguridad
Wend
smar2 = n
End Function
Con esta nueva implementación podemos calcular S(x) para valores
mayores. Lo hemos intentado para comprobar que S(200001)=409 y se
ha conseguido de forma prácticamente instantánea. Coincide con el
resultado obtenido con PARI.
El problema está en que necesitas la función MCD. Aquí tienes una:
Public Function mcd(a1, b1)
Dim a, b, r
'Halla el MCD de a1 y b1
r = 1
a = a1
b = b1
If b = 0 Then b = 1
If a = 0 Then a = 1
While r 0
r = a - b * Int(a / b)
If r 0 Then a = b: b = r
Wend
mcd = b
End Function
Puedes estudiar esta versión muy sintética en PARI:
a(n)={local(m=1,x=n);while(x>1,m++;x=x/gcd(x,m));m}
72