Rosetta Code

Zsigmondy numbers: Difference between revisions

Page Discussion

Zsigmondy numbers (view source)

Revision as of 22:46, 25 November 2022

213 bytes removed ,  1 year ago
m
→‎{{header|Raku}}: minor simplification
No edit summary
m (→‎{{header|Raku}}: minor simplification)
Line 469:
 
sub Zsigmondy ($a, $b) {
my @aexp = 1, $a, * × $a … *;
my @bexp = 1, $b, * × $b … *;
(1..∞).map: -> $n {
my @divisors = divisors(@aexp[$n] - @bexp[$n]).&divisors.sort(-*).first: -*;> $d {
@divisors.first: -> $d { all (1..^$n).map: -> $m { (@aexp[$m_] - @bexp[$m_]) gcd $d == 1 } }
}
}
}
 
for '064078', (2,1), '064079', (3,1), '064080', (4,1), '064081', (5,1), '064082', (6,1),
for 'A064078: Zsigmondy(n,2,1)', (2,1),
'064083', (7,1), '109325', (3,2), '109347', (5,3), '109348', (7,3), '109349', (7,5)
'A064079: Zsigmondy(n,3,1)', (3,1),
-> $oeis, $seq {
'A064080: Zsigmondy(n,4,1)', (4,1),
'A064081say "\nA$oeis: Zsigmondy(n,5$seq[0],$seq[1])',:\n" ~ Zsigmondy(5,1|$seq),[^20]
}</syntaxhighlight>
'A064082: Zsigmondy(n,6,1)', (6,1),
'A064083: Zsigmondy(n,7,1)', (7,1),
'A109325: Zsigmondy(n,3,2)', (3,2),
'A109347: Zsigmondy(n,5,3)', (5,3),
'A109348: Zsigmondy(n,7,3)', (7,3),
'A109349: Zsigmondy(n,7,5)', (7,5)
-> $name, $seq { say "\n$name:\n" ~ Zsigmondy(|$seq)[^20] }</syntaxhighlight>
{{out}}
<pre>A064078: Zsigmondy(n,2,1):
Thundergnat
10,333

edits

Retrieved from "https://rosettacode.org/wiki/Special:MobileDiff/333875"