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()