python求1到100的奇数和
丰色 发自 凹非寺
量子位 | 公众号 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）来模拟图案的生成。这两者之间的关系类似于猎物和捕食者，它们自行移动并有交互：
- 最初，U和V随机分布在一个域上；
- 在每个时间步，它们逐渐扩散到邻近空间；
- 当U和V相遇时，一部分U被V吞噬。因此，V的浓度增加；
- 为了避免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":"cache","pub_time":1662700080000,"pure":"","rawurl":"http://zm.news.so.com/ff7cacee7feffebd5b82b9dcf2d60187","redirect":0,"rptid":"aec5ec457dee7f65","s":"t","src":"量子位","tag":[],"title":"胡渊鸣：import一个“太极”库，让Python代码提速100倍
燕怨家4571求作文——以花木兰的角度(第一人称)叙述自己在战场的经历、感受.300字 -
徐胥承18257521368 ______ 战场上的经历与感受我是花木兰,一个平凡的女子.当我听说父亲被征召上战场时,我决定顶替他去战斗.在战场上,我经历了许多艰辛和考验,也感受到了战争的残酷和无情.战场上的日子充满了紧张和恐惧.我身穿盔甲,手握长,与敌人搏...
燕怨家4571一道数学题:分子与分母之和是100,如果分子加上23,分母加上32,新分数约分后是三分之二.原来分数是多少请各位帮我算算吧!我在线等到一点四十.下... -
徐胥承18257521368 ______[答案] 假设原来分子是x,分母是100-x,那么列方程: 23+x:32+100-x=2:3 264-2x=3x+69 5x=195 x=39 原来分数是39/61
燕怨家4571用1 - 9 九个数字 在不改变他们排列顺序的情况下 要求只在数字中间 添上运算符号 是的数等于100 -
徐胥承18257521368 ______[答案] 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 * 9 = 100 1 - 2 * 3 + 4 * 5 + 6 + 7 + 8 * 9 = 100 1 - 2 * 3 - 4 + 5 * 6 + 7 + 8 * 9 = 100 1 - 2 * 3 - 4 - 5 + 6 * 7 + 8 * 9 = 100 1 * 2 * 3 * 4 + 5 + 6 + 7 * 8 + 9 = 100 1 * 2 * 3 * 4 + 5 + 6 - 7 + 8 * 9 = 100 1 - 2 + 3 * 4 * 5 + 6 * 7 + 8 - 9 = ...
燕怨家4571我国明代珠算家程大位的名著《直指算法统宗》里有一道著名算题:“有100个和尚分100只馒头正好分完.如果大和尚一人分3只小和尚3人分一只,试问大、... -
徐胥承18257521368 ______[选项] A. 25人 B. 30人 C. 50人 D. 75人
燕怨家4571求灵能百分百第二季高清无字幕资源啊!!! -
徐胥承18257521368 ______ 链接: https://pan.baidu.com/s/1WXjImHT5uwUjZYw6rtn3GA 提取码: i677 复制这段内容后打开百度网盘手机App,操作更方便哦 《灵能百分百》是东京电视台2018年1月18日起播出的连续剧,由与田祐希、滨田龙臣主演. 如果资源不正确,或者版本不正确,欢迎追问
燕怨家4571韩国电影《局内人》,求百度云,180分钟.多谢. -
徐胥承18257521368 ______ 链接: https://pan.baidu.com/s/1STbb3XjRvzbEkcQ7kTbsrQ 提取码: vdid 《局内人们》是根据尹胎镐的同名漫画改编的犯罪电影,由吴民浩执导,李秉宪、曹承佑、白允植主演,该片于2015年11月19日在韩国上映 . 该片讲述了可以左右韩国社会的“局内人们”彼此之间义气和背叛的犯罪故事 . 2016年10月7日,该片荣获第25届韩国釜日电影奖男主角奖(李秉宪). 10月24日,由韩国电影评论家协会主办的第36届影评奖,获奖结果公布,《局内人们》入选“影评奖十佳”电影.11月25日,该片荣获第37届韩国青龙电影奖最佳影片、最佳男主角奖(李秉宪).
燕怨家4571我想要一张微信500余额的截图? -
徐胥承18257521368 ______ 你好!如果你想要一张微信500余额的截图可以自己做一个呀,这个别人的肯定不会给你使用的因为这是人家的,而且这样的截图做起来应该也不困难吧.如果你有余额直接截一个就可以如果没有余额你的美图功底好也可以p一个.
燕怨家4571请你运用一种或多种描写方法写一个人物描写的语段不少于100字 -
徐胥承18257521368 ______ 我的同学★瞧她那黝黑的脸蛋,透露着一种不服输的性格,又粗又黑的眉毛下闪着一双玻璃球似的眼睛,小心地注视着对手的一举一动,机警灵活.凹凸的轮廓勾勒出了她对乒乓球的不懈追求.虽然个儿矮,也并不强壮,但她仍然疯狂地热爱运动.正如她那个令人震撼的名字——雷雨点,雷打不断!