prog.sf 434 B

12345678910111213141516171819202122232425262728
  1. #!/usr/bin/ruby
  2. # Perfect powers that are equal to the sum of the first k perfect powers > 1 for some k.
  3. # https://oeis.org/A364948
  4. # Known terms:
  5. # 4, 121, 2548735225
  6. # Next term, if it exists, is > 979846576384.
  7. var sum = 0
  8. var sums = Set()
  9. var n = 1
  10. for k in (2..1e9) {
  11. n = n.next_perfect_power
  12. sum += n
  13. sums << sum
  14. if (sums.has(n)) {
  15. say [n, k]
  16. }
  17. }
  18. __END__
  19. [4, 2]
  20. [121, 14]
  21. [2548735225, 51912]