experimental-projects/oeis-research/Labos Elemer/a(n) = tau( sigma_n(n) ), where tau is the number of divisors of n/prog.sf

#!/usr/bin/ruby

# a(n) = tau( sigma_n(n) ), where tau is the number of divisors of n.
# https://oeis.org/A064165

include("../../../factordb/auto.sf")

func a(n) {
    sigma(sigma(n,n), 0)
}

for n in (1..1000) {
    #print(a(n), ", ")
    say "#{n} #{a(n)}"
}