STATUS
proposed
approved
The OEIS is supported by the many generous donors to the OEIS Foundation.
proposed
approved
editing
proposed
Jairo Bochi and Piotr Laskawiec, <a href="https://arxiv.org/abs/2301.12574">Spectrum maximizing products are not generically unique</a>, arXiv:2301.12574 [math.OC], 2023.
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
return 1 if n<1 else ((2**(n//2+1) if n%2 + else 3)/4*2**(n//2-1)) + sum(totient(n//d)*2**d for d in divisors(n))/(2*/n)//2
proposed
editing
editing
proposed
def a000029(n):
def a000029(n): return 1 if n<1 else (n%2 + 3)/4*2**int(n//2) + sum([totient(n//d)*2**d for d in divisors(n)])/(2*n)
def a(n):
def a(n): return 1 if n<1 else sum([mobius(d)*a000029(n//d) for d in divisors(n)])
print ([a(n) for n in range(51)] ) # Indranil Ghosh, Apr 23 2017
approved
editing
print [a(n) for n in xrangerange(51)] # Indranil Ghosh, Apr 23 2017