优化算法时间复杂度
算法的时间复杂度对程序的执行效率影响***,在 Python 中可以通过选择合适的数据结构来优化时间复杂度,如 list 和 set 查找某一个元素的时间复杂度分别是O(n)和O(1)。不同的场景有不同的优化方式,总得来说,一般有分治,分支界限,贪心,动态规划等思想。
减少冗余数据
如用上三角或下三角的方式去保存一个大的对称矩阵。在0元素占大多数的矩阵里使用稀疏矩阵表示。
合理使用 copy 与 deepcopy
对于 dict 和 list 等数据结构的对象,直接赋值使用的是引用的方式。而有些情况下需要复制整个对象,这时可以使用 copy 包里的 copy 和 deepcopy,这两个函数的不同之处在于后者是递归复制的。效率也不一样:(以下程序在 ipython 中运行)
- import copy
- a = range(100000)
- %timeit -n 10 copy.copy(a) # 运行10次 copy.copy(a)
- %timeit -n 10 copy.deepcopy(a)
- 10 loops, best of 3: 1.55 ms per loop
- 10 loops, best of 3: 151 ms per loop
timeit 后面的-n表示运行的次数,后两行对应的是两个 timeit 的输出,下同。由此可见后者慢一个数量级。
使用 dict 或 set 查找元素
python dict 和 set 都是使用 hash 表来实现(类似c++11标准库中unordered_map),查找元素的时间复杂度是O(1)
- a = range(1000)
- s = set(a)
- d = dict((i,1) for i in a)
- %timeit -n 10000 100 in d
- %timeit -n 10000 100 in s10000 loops, best of 3: 43.5 ns per loop10000 loops, best of 3: 49.6 ns per loop
dict 的效率略高(占用的空间也多一些)。
合理使用生成器(generator)和 yield
- %timeit -n 100 a = (i for i in range(100000))
- %timeit -n 100 b = [i for i in range(100000)]100 loops, best of 3: 1.54 ms per loop100 loops, best of 3: 4.56 ms per loop
使用()得到的是一个 generator 对象,所需要的内存空间与列表的大小无关,所以效率会高一些。在具体应用上,比如 set(i for i in range(100000))会比 set([i for i in range(100000)])快。
但是对于需要循环遍历的情况:
- %timeit -n 100 a = (i for i in range(100000))
- %timeit -n 100 b = [i for i in range(100000)]100 loops, best of 3: 1.54 ms per loop100 loops, best of 3: 4.56 ms per loop
后者的效率反而更高,但是如果循环里有 break,用 generator 的好处是显而易见的。yield 也是用于创建 generator:
- def yield_func(ls):
- for
- i in ls:
- yield i+1
- def not_yield_func(ls):
- return [i+1 for i in ls]
- ls = range(1000000)
- %timeit -n 10 for i in yield_func(ls):pass
- %timeit -n 10 for i in not_yield_func(ls):pass
- 10 loops, best of 3: 63.8 ms per loop
- 10 loops, best of 3: 62.9 ms per loop
对于内存不是非常大的 list,可以直接返回一个 list,但是可读性 yield 更佳(人个喜好)。
python2.x 内置 generator 功能的有 xrange 函数、itertools 包等。
优化循环
循环之外能做的事不要放在循环内,比如下面的优化可以快一倍:
- a = range(10000)
- size_a = len(a)
- %timeit -n 1000 for i in a: k = len(a)
- %timeit -n 1000 for i in a: k = size_a
- 1000 loops, best of 3: 569 µs per loop
- 1000 loops, best of 3: 256 µs per loop
优化包含多个判断表达式的顺序
对于 and,应该把满足条件少的放在前面,对于 or,把满足条件多的放在前面。如:
- a = range(2000)
- %timeit -n 100 [i for i in a if 10 < i < 20 or 1000 < i < 2000]
- %timeit -n 100 [i for i in a if 1000 < i < 2000 or 100 < i < 20]
- %timeit -n 100 [i for i in a if i % 2 == 0 and i > 1900]
- %timeit -n 100 [i for i in a if i > 1900 and i % 2 == 0]
- 100 loops, best of 3: 287 µs per loop
- 100 loops, best of 3: 214 µs per loop
- 100 loops, best of 3: 128 µs per loop
- 100 loops, best of 3: 56.1 µs per loop
使用 join 合并迭代器中的字符串
- In [1]: %%timeit
- ...: s = ''
- ...: for i in a:
- ...: s += i
- ...:10000 loops, best of 3: 59.8 µs per loopIn [2]: %%timeit
- s = ''.join(a)
- ...:100000 loops, best of 3: 11.8 µs per loop
join 对于累加的方式,有大约5倍的提升。
选择合适的格式化字符方式
- s1, s2 = 'ax', 'bx'
- %timeit -n 100000 'abc%s%s' % (s1, s2)
- %timeit -n 100000 'abc{0}{1}'.format(s1, s2)
- %timeit -n 100000 'abc' + s1 + s2
- 100000 loops, best of 3: 183 ns per loop
- 100000 loops, best of 3: 169 ns per loop
- 100000 loops, best of 3: 103 ns per loop
不借助中间变量交换两个变量的值
- In [3]: %%timeit -n 10000
- a,b=1,2
- ....: c=a;a=b;b=c;
- ....:10000 loops, best of 3: 172 ns per loop
- In [4]: %%timeit -n 10000
- a,b=1,2
- a,b=b,a
- ....:
- 10000 loops, best of 3: 86 ns per loop
使用a,b=b,a而不是c=a;a=b;b=c;来交换a,b的值,可以快1倍以上。
使用 if is
- a = range(10000)
- %timeit -n 100 [i for i in a if i == True]
- %timeit -n 100 [i for i in a if i is True]
- 100 loops, best of 3: 531 µs per loop
- 100 loops, best of 3: 362 µs per loop
使用 if is True 比 if == True 将近快一倍。
使用级联比较x < y < z
- x, y, z = 1,2,3
- %timeit -n 1000000 if x < y < z:pass
- %timeit -n 1000000 if x < y and y < z:pass
- 1000000 loops, best of 3: 101 ns per loop
- 1000000 loops, best of 3: 121 ns per loop
x < y < z效率略高,而且可读性更好。
while 1 比 while True 更快
- def while_1():
- n = 100000
- while 1:
- n -= 1
- if n <= 0: break
- def while_true():
- n = 100000
- while True:
- n -= 1
- if n <= 0: break
- m, n = 1000000, 1000000
- %timeit -n 100 while_1()
- %timeit -n 100 while_true()
- 100 loops, best of 3: 3.69 ms per loop
- 100 loops, best of 3: 5.61 ms per loop
while 1 比 while true 快很多,原因是在 python2.x 中,True 是一个全局变量,而非关键字。
使用**而不是 pow
- %timeit -n 10000 c = pow(2,20)
- %timeit -n 10000 c = 2**2010000 loops, best of 3: 284 ns per loop10000 loops, best of 3: 16.9 ns per loop
**就是快10倍以上!
使用 cProfile, cStringIO 和 cPickle 等用c实现相同功能(分别对应profile, StringIO, pickle)的包
- import cPickle
- import pickle
- a = range(10000)
- %timeit -n 100 x = cPickle.dumps(a)
- %timeit -n 100 x = pickle.dumps(a)
- 100 loops, best of 3: 1.58 ms per loop
- 100 loops, best of 3: 17 ms per loop
由c实现的包,速度快10倍以上!
使用***的反序列化方式
下面比较了 eval, cPickle, json 方式三种对相应字符串反序列化的效率:
- import json
- import cPickle
- a = range(10000)
- s1 = str(a)
- s2 = cPickle.dumps(a)
- s3 = json.dumps(a)
- %timeit -n 100 x = eval(s1)
- %timeit -n 100 x = cPickle.loads(s2)
- %timeit -n 100 x = json.loads(s3)
- 100 loops, best of 3: 16.8 ms per loop
- 100 loops, best of 3: 2.02 ms per loop
- 100 loops, best of 3: 798 µs per loop
可见 json 比 cPickle 快近3倍,比 eval 快20多倍。
使用C扩展(Extension)
目前主要有 CPython(python最常见的实现的方式)原生API, ctypes,Cython,cffi三种方式,它们的作用是使得 Python 程序可以调用由C编译成的动态链接库,其特点分别是:
CPython 原生 API: 通过引入 Python.h 头文件,对应的C程序中可以直接使用Python 的数据结构。实现过程相对繁琐,但是有比较大的适用范围。
ctypes: 通常用于封装(wrap)C程序,让纯 Python 程序调用动态链接库(Windows 中的 dll 或 Unix 中的 so 文件)中的函数。如果想要在 python 中使用已经有C类库,使用 ctypes 是很好的选择,有一些基准测试下,python2+ctypes 是性能***的方式。
Cython: Cython 是 CPython 的超集,用于简化编写C扩展的过程。Cython 的优点是语法简洁,可以很好地兼容 numpy 等包含大量C扩展的库。Cython 的使得场景一般是针对项目中某个算法或过程的优化。在某些测试中,可以有几百倍的性能提升。
cffi: cffi 的就是 ctypes 在 pypy(详见下文)中的实现,同进也兼容 CPython。cffi提供了在 python 使用C类库的方式,可以直接在 python 代码中编写C代码,同时支持链接到已有的C类库。
使用这些优化方式一般是针对已有项目性能瓶颈模块的优化,可以在少量改动原有项目的情况下大幅度地提高整个程序的运行效率。
并行编程
因为 GIL 的存在,Python 很难充分利用多核 CPU 的优势。但是,可以通过内置的模块 multiprocessing 实现下面几种并行模式:
多进程:对于 CPU 密集型的程序,可以使用 multiprocessing 的 Process,Pool 等封装好的类,通过多进程的方式实现并行计算。但是因为进程中的通信成本比较大,对于进程之间需要大量数据交互的程序效率未必有大的提高。
多线程:对于 IO 密集型的程序,multiprocessing.dummy 模块使用 multiprocessing 的接口封装 threading,使得多线程编程也变得非常轻松(比如可以使用 Pool 的 map 接口,简洁高效)。
分布式:multiprocessing 中的 Managers 类提供了可以在不同进程之共享数据的方式,可以在此基础上开发出分布式的程序。
不同的业务场景可以选择其中的一种或几种的组合实现程序性能的优化。
终级大杀器:PyPy
PyPy 是用 RPython(CPython 的子集)实现的 Python,根据官网的基准测试数据,它比 CPython 实现的 Python 要快6倍以上。快的原因是使用了 Just-in-Time(JIT)编译器,即动态编译器,与静态编译器(如gcc,javac等)不同,它是利用程序运行的过程的数据进行优化。由于历史原因,目前 pypy 中还保留着 GIL,不过正在进行的 STM 项目试图将 PyPy 变成没有 GIL 的 Python。
如果 python 程序中含有C扩展(非cffi的方式),JIT 的优化效果会大打折扣,甚至比 CPython 慢(比 Numpy)。所以在 PyPy 中***用纯 Python 或使用 cffi 扩展。
随着STM,Numpy 等项目的完善,相信 PyPy 将会替代 CPython。
使用性能分析工具
除了上面在 ipython 使用到的 timeit 模块,还有 cProfile。cProfile 的使用方式也非常简单: python -m cProfile filename.py,filename.py 是要运行程序的文件名,可以在标准输出中看到每一个函数被调用的次数和运行的时间,从而找到程序的性能瓶颈,然后可以有针对性地优化。
参考
[1] http://www.ibm.com/developerworks/cn/linux/l-cn-python-optim/
[2] http://maxburstein.com/blog/speeding-up-your-python-code/