OFFSET
0,3
LINKS
Chai Wah Wu, Table of n, a(n) for n = 0..14
EXAMPLE
A038003(4) = 9694845 which has 96 divisors.
PROG
(Python)
from operator import mul
from functools import reduce
from sympy import factorint
A255132_list, c, s = [1, 1], {}, 3
for n in range(2, 2**15):
for p, e in factorint(4*n-2).items():
if p in c:
c[p] += e
else:
c[p] = e
for p, e in factorint(n+1).items():
if c[p] == e:
del c[p]
else:
c[p] -= e
if n == s:
c2 = reduce(mul, [e+1 for e in c.values()])
A255132_list.append(c2)
s = 2*s+1
CROSSREFS
KEYWORD
nonn
AUTHOR
Chai Wah Wu, Feb 15 2015
STATUS
approved
