prog.sf 291 B

12345678910111213141516171819202122
  1. #!/usr/bin/ruby
  2. # a(n) = Sum of the divisors of sigma_n(n).
  3. # https://oeis.org/A347718
  4. include("../../../factordb/auto.sf")
  5. func a(n) {
  6. sigma(sigma(n, n))
  7. }
  8. for n in (1..1000) {
  9. #print(a(n), ", ")
  10. say "#{n} #{a(n)}"
  11. }
  12. __END__
  13. # PARI/GP program:
  14. a(n) = sigma(sigma(n, n));