native.pl 1.1 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546
  1. #!/usr/bin/perl
  2. # Floor of log of product of all primes between n!+1 and (n+1)!.
  3. # https://oeis.org/A294196
  4. use 5.014;
  5. use Math::GMPz qw();
  6. use Math::MPFR qw(MPFR_RNDZ);
  7. use ntheory qw(forprimes factorial);
  8. for my $n (2..100) {
  9. my $f = 0;
  10. forprimes {
  11. $f += log($_);
  12. } factorial($n)+1, factorial($n+1);
  13. say "floor($f) = ", int($f);
  14. }
  15. __END__
  16. floor(2.70805020110221) = 2
  17. floor(15.821901318467) = 15
  18. floor(87.846705739409) = 87
  19. floor(579.092934941327) = 579
  20. floor(4276.65691765426) = 4276
  21. floor(35103.7037849562) = 35103
  22. floor(322168.829509598) = 322168
  23. floor(3264471.9003149) = 3264471
  24. floor(36285842.9959193) = 36285842
  25. floor(439070392.951992) = 439070392
  26. floor(5747983086.45772) = 5747983086
  27. floor(2.70805020110221007) = 2
  28. floor(15.8219013184670268) = 15
  29. floor(87.8467057394090264) = 87
  30. floor(579.092934941326837) = 579
  31. floor(4276.65691765426939) = 4276
  32. floor(35103.7037849562028) = 35103
  33. floor(322168.82950959723) = 322168
  34. floor(3264471.90031498688) = 3264471
  35. floor(36285842.9959180593) = 36285842
  36. floor(439070392.95186203) = 439070392
  37. floor(5747983086.4605573) = 5747983086