利用Python计算1000以内自然数的约数个数，然后通过matplotlib绘制统计图。

下图为约数个数的散点图及其分布情况的条形图。

Python代码：

import collections
import matplotlib.pyplot as plt

def countDivisors(num):
    ans = 1
    x = 2
    while x * x <= num:
        cnt = 1
        while num % x == 0:
            cnt += 1
            num /= x
        ans *= cnt
        x += 1
    return ans * (1 + (num > 1))

MAXNUM = 1000
x = range(1, MAXNUM)
y = map(countDivisors, x)

plt.subplot(2, 1, 1)
plt.title('Divisors Count')
plt.xlim(0, MAXNUM)
plt.ylim(0, max(y) + 1)
plt.scatter(x, y)
plt.grid(True)

plt.subplot(2, 1, 2)
plt.title('Statistics of Divisor Count')
z = collections.Counter(y)
plt.bar(z.keys(), z.values(), align = 'center')
plt.grid(True)

plt.show()

 

本文链接：http://bookshadow.com/weblog/2016/11/28/python-matplotlib-divisor-count-scatter-bar/
请尊重作者的劳动成果，转载请注明出处！书影博客保留对文章的所有权利。

如果您喜欢这篇博文，欢迎您捐赠书影博客： ，查看支付宝二维码