Creo que para solucionar este problema no es necesario ningún bucle.
Esta es la solución más simple que he podido encontrar:
from __future__ import division
def sumar_multiplos_k(a1, an,d):
n = an/d
S = ((a1+an)/2)*n
return S
print(sumar_multiplos_k(3,999,3)+sumar_multiplos_k(5,995,
Claro están siempre las ecuaciones de sumatorias, el mejor método para este
caso, de hecho en euler lo comentan hasta el cansancio, en mi caso prefiero
buscar una soluciòn más del mundo del algoritmo que de las matematicas para
jugar con el lenguaje. Seguiré la idea de cython ya que me parece
inter