首页 >>  正文

python打出1到100所有数字

来源:baiyundou.net   日期:2024-09-21

丰色 发自 凹非寺

量子位 | 公众号 QbitAI

众所周知,Python的简单和易读性是靠牺牲性能为代价的——

尤其是在计算密集的情况下,比如多重for循环。

不过现在,大佬胡渊鸣说了:

只需import 一个叫做“Taichi”的库,就可以把代码速度提升100倍

不信?

来看三个例子。

计算素数的个数,速度x120

第一个例子非常非常简单,求所有小于给定正整数N的素数。

标准答案如下:

我们将上面的代码保存,运行。

当N为100万时,需要2.235s得到结果:

现在,我们开始施魔法。

不用更改任何函数体,import“taichi”库,然后再加两个装饰器:

Bingo!同样的结果只要0.363s,快了将近6倍。

如果N=1000万,则只要0.8s;要知道,不加它可是55s,一下子又快了70倍

不止如此,我们还可以在ti.init()中加个参数变为ti.init(arch=ti.gpu) ,让taich在GPU上进行计算。

那么此时,计算所有小于1000万的素数就只耗时0.45s了,与原来的Python代码相比速度就提高了120倍

厉不厉害?

什么?你觉得这个例子太简单了,说服力不够?我们再来看一个稍微复杂一点的。

动态规划,速度x500

动态规划不用多说,作为一种优化算法,通过动态存储中间计算结果来减少计算时间。

我们以经典教材《算法导论》中的经典动态规划案例“最长公共子序列问题(LCS)”为例。

比如对于序列a = [0, 1, 0, 2, 4, 3, 1, 2, 1]和序列b = [4, 0, 1, 4, 5, 3, 1, 2],它们的LCS就是:

LCS(a, b) = [0, 1, 4, 3, 1, 2]。

用动态规划的思路计算LCS,就是先求解序列a的前i个元素和序列b的前j个元素的最长公共子序列的长度,然后逐步增加i或j的值,重复过程,得到结果。

我们用f[i, j]来指代这个子序列的长度,即LCS((prefix(a, i), prefix(b, j)。其中prefix(a, i) 表示序列a的前i个元素,即a[0], a[1], …, a[i - 1],得到如下递归关系:

完整代码如下:

现在,我们用Taichi来加速:

结果如下:

胡渊鸣电脑上的程序最快做到了0.9秒内完成,而换成用NumPy来实现,则需要476秒,差异达到了超500倍!

最后,我们再来一个不一样的例子。

反应 - 扩散方程,效果惊人

自然界中,总有一些动物身上长着一些看起来无序但实则并非完全随机的花纹。

图灵机的发明者艾伦·图灵是第一个提出模型来描述这种现象的人。

在该模型中,两种化学物质(U和V)来模拟图案的生成。这两者之间的关系类似于猎物和捕食者,它们自行移动并有交互:

  1. 最初,U和V随机分布在一个域上;
  2. 在每个时间步,它们逐渐扩散到邻近空间;
  3. 当U和V相遇时,一部分U被V吞噬。因此,V的浓度增加;
  4. 为了避免U被V根除,我们在每个时间步添加一定百分比 (f) 的U并删除一定百分比 (k) 的V。

上面这个过程被概述为“反应-扩散方程”:

其中有四个关键参数:Du(U的扩散速度),Dv(V的扩散速度),f(feed的缩写,控制U的加入)和k(kill的缩写,控制V的去除)

如果Taichi中实现这个方程,首先创建网格来表示域,用vec2表示每个网格中U, V的浓度值。

拉普拉斯算子数值的计算需要访问相邻网格。为了避免在同一循环中更新和读取数据,我们应该创建两个形状相同的网格W×H×2。

每次从一个网格访问数据时,我们将更新的数据写入另一个网格,然后切换下一个网格。那么数据结构设计就是这样:

一开始,我们将U在网格中的浓度设置为 1,并将V放置在50个随机选择的位置:

那么实际计算就可以用不到10行代码完成:

","force_purephv":"0","gnid":"937e299d4c2c2a176","img_data":[{"flag":2,"img":[{"desc":"","height":"920","title":"","url":"https://p0.ssl.img.360kuai.com/t01e23794c94540b7d3.jpg","width":"920"},{"desc":"","height":"765","title":"","url":"https://p0.ssl.img.360kuai.com/t0144fb37a5bf3f53de.jpg","width":"1080"},{"desc":"","height":"227","title":"","url":"https://p0.ssl.img.360kuai.com/t0150b2fef7620e3dbc.jpg","width":"1080"},{"desc":"","height":"839","title":"","url":"https://p0.ssl.img.360kuai.com/t016b17cfbc73cc0857.jpg","width":"1080"},{"desc":"","height":"218","title":"","url":"https://p0.ssl.img.360kuai.com/t0179edf8404f9ae37d.jpg","width":"1080"},{"desc":"","height":"106","title":"","url":"https://p0.ssl.img.360kuai.com/t01e32851aaad11b65f.jpg","width":"1080"},{"desc":"","height":"179","title":"","url":"https://p0.ssl.img.360kuai.com/t01bec382de36100f5c.jpg","width":"1080"},{"desc":"","height":"1048","title":"","url":"https://p0.ssl.img.360kuai.com/t01d29f206b72759e5c.jpg","width":"952"},{"desc":"","height":"220","title":"","url":"https://p0.ssl.img.360kuai.com/t013e42423d8b222dd8.jpg","width":"1080"},{"desc":"","height":"200","title":"","url":"https://p0.ssl.img.360kuai.com/t01d2b2fb1751d67f82.jpg","width":"1080"},{"desc":"","height":"290","title":"","url":"https://p0.ssl.img.360kuai.com/t01a6e79e099eafb7c4.jpg","width":"732"},{"desc":"","height":"103","title":"","url":"https://p0.ssl.img.360kuai.com/t0120a647a3f436ce6b.jpg","width":"1080"},{"desc":"","height":"321","title":"","url":"https://p0.ssl.img.360kuai.com/t01758422991b7876fd.jpg","width":"1080"}]}],"original":0,"pat":"art_src_3,fts0,sts0","powerby":"hbase","pub_time":1662700080000,"pure":"","rawurl":"http://zm.news.so.com/ff7cacee7feffebd5b82b9dcf2d60187","redirect":0,"rptid":"aec5ec457dee7f65","s":"t","src":"量子位","tag":[],"title":"胡渊鸣:import一个“太极”库,让Python代码提速100倍

戚矩界4531怎么用python做一个输出10个文字排列的文本文件 -
边钟肥19830303507 ______ #coding =utf-8 reader = open('test.txt').readlines() n=0 for len in reader: writer=open(str(reader.index(len)+1)+'.txt','w') if (len.find('\n')>0): for n in range(0,10): writer.write(len[:-1]+reader[n]) else: for n in range(0,10): writer.write(len+reader[n]) writer.close() 这个基本能满足你要的功能,你把输入和输出的文件名改一下就行了

戚矩界4531怎么用Python 的for循环和while循环算平均数分别用for和while 在1到10里任意输入五个数然后算出他们的平均值.求大神 -
边钟肥19830303507 ______[答案] m = 0for i in range(5): s = raw_input() m += float(s)/5.0 print 'the average is: ',m i = 0m = 0while i解析看不懂?免费查看同类题视频解析查看解答

戚矩界4531python编程 输入一组整数数组,求出两两之差的最小绝对值.只需得出最小值 如:输入:[10,3,12,9] 输入:1 -
边钟肥19830303507 ______ 你好, 代码如下: ------ a = [10, 8, 2, 45, 69, 38, 11, 15] #假设该列表为需要输入的一组数 a.sort(reverse = True) #首先对这组数进行从大到小的排序 print a #输出排序结果 min = a[0] #令min变量记录该列表中最大的值 for i in range( len(a) -1 ): #i用...

戚矩界4531python输入十个数,打印每个数对应的立方数 -
边钟肥19830303507 ______ 以下是示例代码: for i in range(10): n=input('输入一个数:\n') fn=float(n) print(fn**3)也可以一次输入10个数,中间用空格分开: nnum=input('输入10个数,中间用空格分开:\n') nnums=nnum.split() ns=[float(n) for n in nnums] for n in ns: print(n**3)不明白可追问

戚矩界4531键盘输入一个数(5到10)如果不是之间的数,允许用户重新输入编写一个... -
边钟肥19830303507 ______ #!/usr/bin/python# Filename : factorial.pyfrom functools import reduceresult = 0numb = int(input('请输入一个(5-10)...

戚矩界4531怎样在python中生成一串0到n的数字名相加 -
边钟肥19830303507 ______ def print_plus(n): t = n+1 for i in range(t): ostr = str(i) if i == 0: strs = ostr + '+' elif i strs += ostr + '+' else: strs += ostr else: return strs print print_plus(6)代码如上.结果如图:

戚矩界4531怎样用python编程将二进制转成10进制 -
边钟肥19830303507 ______ 十进制整数转换为二进制整数采用"除2取余,逆序排列"法. #include <stdio.h> int main() { int n,a[100],i=0,j; scanf("%d",&n); while(n) 扩展资料: 基本符号是0到9十个数字.要表示这十个数的10倍,就将这些数字左移一位,用0补上空位,即10,20,30,...,90;要表示这十个数的10倍,就继续左移数字的位置,即100,200,300,....要表示一个数的1/10,就右移这个数的位置,需要时就0补上空位:1/10位0.1,1/100为0.01,1/1000为0.001. 参考资料来源:百度百科-十进制

戚矩界4531编写一个程序,求1到10之间所有偶数的和及其所有奇数的和 -
边钟肥19830303507 ______ 以python语言为例: #!/usr/bin/env python # -*- coding: utf-8 -*- #编写一个程序,求1到10之间所有偶数的和及其所有奇数的和 sum_1 = 0 sum_2 = 0 for i in range(1,11): if (i % 2==0): sum_2=sum_2+i else: sum_1=sum_1+i print "奇数和是%d" ...

戚矩界4531python 3 语言写一个计算输入数每位数字和的.例如输入128 就算出1+2+8.要求见补充 -
边钟肥19830303507 ______ print('Sum of digits :',sum(map(int,input('Enter a number :'))))

(编辑:自媒体)
关于我们 | 客户服务 | 服务条款 | 联系我们 | 免责声明 | 网站地图 @ 白云都 2024