利用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/
请尊重作者的劳动成果,转载请注明出处!书影博客保留对文章的所有权利。