https://oeis.org/A153823

factdivisors.py

from sympy import *
for n in range(2,14):
    r=1
    for p in primerange(n+1):
        q=1
        x=0
        while p**q<=n:
            x+=n//(p**q)
            q+=1
        r*=(x+1)
        r%=m
    print(r-1,end=', ')