-
Notifications
You must be signed in to change notification settings - Fork 0
/
from_lcm.pl
30 lines (25 loc) · 1.52 KB
/
from_lcm.pl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
#!/usr/bin/perl
use 5.014;
use Math::GMPz;
use ntheory qw(forsemiprimes forprimes factor forsquarefree random_prime divisors gcd next_prime);
use Math::Prime::Util::GMP qw(mulint is_pseudoprime vecprod divint sqrtint vecprod);
#my $k = "495088126122885";
#my $k = "19976310800932286865"; # record
#my $k = vecprod(3, 5, 17, 23, 29, 43, 53, 89, 113, 127, 157, 257);
#my $k = vecprod(3, 5, 17, 23, 29, 43, 53, 89, 113, 127, 157, 257, 353, 397, 449);
my $k = vecprod(3, 5, 17, 23, 29, 43, 53, 89, 113, 127, 157, 257);
#my $n = "97492561525843375414411984711690826244618218549082363628202417811869034234870968105413397502977";
#my $n = "13757974351340992129905979322929973774675683422410013622969215110517423192617";
#my $n = "492998637615524101988528566254824872883187100795158385773434977471920563361414751";
#my $n = "343007266355769778384701673249549952892082965850160905177012001818058446029240481";
#my $n = "32570940991344830846076116786757514876773522188233729072893528656637375859598588354317";
#my $n = "1783678161425679359147266763387785520358862736924657051419350376424361065071938287397480735353585534563540084656470093";
#my $n = "786769380874976962176734303050188467034887840304275124846721199585748921910091029860031756177";
#my $n = "1611172416406363222443655272298918903139645660167379869658840681191970272250330666161";
my $n = "213259970167095025136189836699535867809435285476249320300587187293474667197712227";
foreach my $d(divisors($n)) {
my $t = mulint($k, $d);
if (is_pseudoprime($t, 2)) {
say $t;
}
}