|
|
|
|
#20 by Joerg Arndt at Tue Aug 13 01:55:58 EDT 2024
|
|
|
|
#19 by Michel Marcus at Mon Aug 12 23:41:02 EDT 2024
|
|
|
|
#18 by Chai Wah Wu at Mon Aug 12 22:26:39 EDT 2024
|
|
|
|
#17 by Chai Wah Wu at Mon Aug 12 22:26:32 EDT 2024
|
| PROG
|
return -(1<<len(g:=factorint(m).values()))+prod(e+1 for e in g) # Chai Wah Wu, Aug 12 2024
|
|
|
|
#16 by Chai Wah Wu at Mon Aug 12 22:26:23 EDT 2024
|
| PROG
|
(Python)
from math import prod, isqrt
from sympy import mobius, factorint
def A174961(n):
def f(x): return n+sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
m, k = n, f(n)
while m != k:
m, k = k, f(k)
return -(1<<len(g:=factorint(m).values()))+prod(e+1 for e in g) # Chai Wah Wu, Aug 12 2024
|
| STATUS
|
approved
editing
|
|
|
|
#15 by OEIS Server at Sat Dec 09 07:07:11 EST 2023
|
| LINKS
|
Amiram Eldar, <a href="/A174961/b174961_1.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from G. C. Greubel)
|
|
|
|
#14 by Peter Luschny at Sat Dec 09 07:07:11 EST 2023
|
|
|
Discussion
|
Sat Dec 09
| 07:07
| OEIS Server: Installed new b-file as b174961.txt. Old b-file is now b174961_1.txt.
|
|
|
|
#13 by Joerg Arndt at Sat Dec 09 04:23:26 EST 2023
|
|
|
|
#12 by Amiram Eldar at Sat Dec 09 03:41:41 EST 2023
|
|
|
|
#11 by Amiram Eldar at Sat Dec 09 03:38:11 EST 2023
|
| EXAMPLE
|
For n = 4, the fourth nonsquarefree number is A013929(4) = 12 which has 2 non-unitary divisors, 42 and 126. Therefore a(4) = 2.
|
|
|
|