110 Diophantine reciprocals II -- oeis.pl 473 B

1234567891011121314151617181920212223242526
  1. #!/usr/bin/perl
  2. # Daniel "Trizen" Șuteu
  3. # Date: 17 February 2017
  4. # Website: https://github.com/trizen
  5. # https://projecteuler.net/problem=110
  6. # Runtime: 1.807s
  7. use 5.010;
  8. use strict;
  9. use LWP::Simple qw(get);
  10. use Math::AnyNum qw(:overload);
  11. use ntheory qw(divisors);
  12. for my $line (split(/\R/, get('https://oeis.org/A002093/b002093.txt'))) {
  13. my $n = (split(' ', $line))[-1] >> 1;
  14. if ((divisors($n * $n) + 1) >> 1 > 4e6) {
  15. say $n;
  16. last;
  17. }
  18. }