本文实例讲述了Python利用matplotlib绘制约数个数统计图。分享给大家供大家参考,具体如下
本文实例讲述了Python利用matplotlib绘制约数个数统计图。分享给大家供大家参考,具体如下:
利用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()
Python
matplotlib
绘制
约数个数
统计图