¿Cuál sería una manera elegante y pitónica de implementar cumsum? Alternativamente
- si there'a ya una forma integrada de hacerlo, que sería aún mejor por supuesto ...Elegante cumsum pitónico
Respuesta
Está disponible en Numpy:
>>> import numpy as np
>>> np.cumsum([1,2,3,4,5])
array([ 1, 3, 6, 10, 15])
O utilice itertools.accumulate
desde Python 3.2 :
>>> from itertools import accumulate
>>> list(accumulate([1,2,3,4,5]))
[ 1, 3, 6, 10, 15]
Si Numpy no es una opción, un circuito generador sería la solución más elegante que se me ocurre:
def cumsum(it):
total = 0
for x in it:
total += x
yield total
Ej.
>>> list(cumsum([1,2,3,4,5]))
>>> [1, 3, 6, 10, 15]
nice and simple with the yield total :) –
para los bucles son Pythonic
def cumsum(vec):
r = [vec[0]]
for val in vec[1:]:
r.append(r[-1] + val)
return r
@larsmans: ¿qué sería 'r [-1]' en ese caso? – WolframH
a=[1,2,3,4,5]
def cumsum(a):
a=iter(a)
cc=[next(a)]
for i in a:
cc.append(cc[-1]+i)
return cc
print cumsum(a)
"[1, 3, 6, 10, 15]"
en su lugar:
a=[1,2,3,4,5]
def cumsum(a):
for i in range(1,len(a)):
a[i]+=a[i-1]
cumsum(a)
print a
"[1, 3, 6, 10, 15]"
a=[1,2,3,4,5]
def cumsum(a):
b=a[:]
for i in range(1,len(a)):
b[i]+=b[i-1]
return b
print cumsum(a)
"[1, 3, 6, 10, 15]"
a = [1, 2, 3 ,4, 5]
# Using list comprehention
cumsum = [sum(a[:i+1]) for i in range(len(a))] # [1, 3, 6, 10, 15]
# Using map()
cumsum = map(lambda i: sum(a[:i+1]), range(len(a))) # [1, 3, 6, 10, 15]
Si bien esto es agradable y fácil, ¡adviértase que en una lista grande esto será muy lento! :) –
Mi idea era utilizar reducir de manera funcional:
from operator import iadd
reduce(lambda acc, itm: iadd(acc, [acc[-1] + itm]), [1, 2, 3, 4, 5], [0])[1:]
>>> [1, 3, 6, 10, 15]
iadd del módulo del operador tiene la propiedad exclusiva de realizar una adición in situ y devolver el destino como resultado.
Si eso [1:] copia hace que usted está nervioso, que asimismo podría hacer:
from operator import iadd
reduce(lambda acc, itm: (iadd(acc[0], [acc[1] + itm]), acc[1] + itm),
[1, 2, 3, 4, 5], ([], 0))[0]
>>> [1, 3, 6, 10, 15]
Pero descubrí que localmente el primer ejemplo es mucho más rápido y generadores de la OMI son más Pythonic de programación funcional como 'reducir ':
reduce(lambda acc, itm: (iadd(acc[0], [acc[1] + itm]), acc[1] + itm), values_ten, ([], 0))[0]
Average: 6.4593828736e-06
reduce(lambda acc, itm: (iadd(acc[0], [acc[1] + itm]), acc[1] + itm), values_mil, ([], 0))[0]
Average: 0.727404361961
reduce(lambda acc, itm: iadd(acc, [acc[-1] + itm]), values_ten, [0])[1:]
Average: 5.16271911336e-06
reduce(lambda acc, itm: iadd(acc, [acc[-1] + itm]), values_mil, [0])[1:]
Average: 0.524223491301
cumsum_rking(values_ten)
Average: 1.9828751369e-06
cumsum_rking(values_mil)
Average: 0.234241141632
list(cumsum_larsmans(values_ten))
Average: 2.02786211569e-06
list(cumsum_larsmans(values_mil))
Average: 0.201473119335
Aquí está la secuencia de comandos de referencia, tu caso es distinto:
from timeit import timeit
def bmark(prog, setup, number):
duration = timeit(prog, setup=setup, number=number)
print prog
print 'Average:', duration/number
values_ten = list(xrange(10))
values_mil = list(xrange(1000000))
from operator import iadd
bmark('reduce(lambda acc, itm: (iadd(acc[0], [acc[1] + itm]), acc[1] + itm), \
values_ten, ([], 0))[0]',
setup='from __main__ import iadd, values_ten', number=1000000)
bmark('reduce(lambda acc, itm: (iadd(acc[0], [acc[1] + itm]), acc[1] + itm), \
values_mil, ([], 0))[0]',
setup='from __main__ import iadd, values_mil', number=10)
bmark('reduce(lambda acc, itm: iadd(acc, [acc[-1] + itm]), \
values_ten, [0])[1:]',
setup='from __main__ import iadd, values_ten', number=1000000)
bmark('reduce(lambda acc, itm: iadd(acc, [acc[-1] + itm]), \
values_mil, [0])[1:]',
setup='from __main__ import iadd, values_mil', number=10)
def cumsum_rking(iterable):
values = list(iterable)
for pos in xrange(1, len(values)):
values[pos] += values[pos - 1]
return values
bmark('cumsum_rking(values_ten)',
setup='from __main__ import cumsum_rking, values_ten', number=1000000)
bmark('cumsum_rking(values_mil)',
setup='from __main__ import cumsum_rking, values_mil', number=10)
def cumsum_larsmans(iterable):
total = 0
for value in iterable:
total += value
yield total
bmark('list(cumsum_larsmans(values_ten))',
setup='from __main__ import cumsum_larsmans, values_ten', number=1000000)
bmark('list(cumsum_larsmans(values_mil))',
setup='from __main__ import cumsum_larsmans, values_mil', number=10)
Y aquí está mi versión de Python cadena:
Python 2.7 (r27:82525, Jul 4 2010, 09:01:59) [MSC v.1500 32 bit (Intel)] on win32
+1 para el benchmarking minucioso – Jonathan
def cumsum(a):
return map(lambda x: sum(a[0:x+1]), range(0, len(a)))
cumsum([1,2,3])
> [1, 3, 6]
- 1. Diseño pitónico sin enums
- 2. ¿Es pitónico usar bools como ints?
- 3. ¿Es pitónico utilizar interfaces/clases base abstractas?
- 4. ¿Es pitónico imitar la sobrecarga de métodos?
- 5. ¿Por qué csv.reader no es pitónico?
- 6. setattr con kwargs, pitónico o no?
- 7. equivalente pitónico de unshift o rehacer?
- 8. funcionalidad "isnotnan" en numpy, ¿puede ser esto más pitónico?
- 9. ¿Cómo diseñarías un marco de interfaz de usuario muy "pitónico"?
- 10. Convirtiendo a mayúsculas. ¿Qué camino es más pitónico?
- 11. ¿Cómo implementar un equivalente pitónico de cola -F?
- 12. Elegante modelo AST
- 13. Comparación de objeto elegante
- 14. Parallel.ForEach - Elegante cancelación
- 15. Regex elegante coincidencia de patrón
- 16. Manera elegante de comparar secuencias
- 17. PowerShell: una manera elegante de crear cierres
- 18. La construcción de lazo más elegante?
- 19. Buscando el código más elegante despachador
- 20. Rieles: ¿Manera elegante de manejar la navegación?
- 21. manera elegante al sesgo booleano aleatoria
- 22. Manera elegante de validar los valores
- 23. Elegante foreach - else construir en Razor
- 24. Rieles: ¿alguna forma elegante de manejar 404?
- 25. Cómo deshabilitar un botón más elegante
- 26. Sobrecarga elegante del operador en D
- 27. Solución isNumeric() más elegante para java
- 28. Salida elegante de un programa en Python?
- 29. Solución elegante para saltos de línea (PHP)
- 30. ¿Cómo diseñar correctamente una ProgressBar vertical elegante?
cumSum como en http://docs.scipy.org/doc/numpy/reference/generated/numpy.cumsum.html? –
http://docs.python.org/dev/library/itertools.html#itertools.accumulate – ChessMaster