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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
|
#!/Users/colincrain/perl5/perlbrew/perls/perl-5.32.0/bin/perl
#
# cubans-in-their-prime.pl
#
# First Series Cuban Primes
# Submitted by: Mohammad S Anwar
# Write a script to compute first series Cuban Primes <= 1000.
# Please refer wikipedia page for more informations.
#
# Output:
# 7, 19, 37, 61, 127, 271, 331, 397, 547, 631, 919.
#
# © 2022 colin crain
## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
use warnings;
use strict;
use utf8;
use feature ":5.26";
use feature qw(signatures);
no warnings 'experimental::signatures';
sub nth_series_cubic ( $dif ) {
my @c;
my $x;
for my $y ( 1..100 ) {
$x = $y + $dif;
push @c, ($x**3 - $y**3) / ($x - $y) ;
}
return @c;
}
sub is_prime ($num) {
for ( 2..sqrt $num ) {
return 0 unless $num % $_;
}
return 1;
}
## output section
my $out;
for my $diff ( 1..10 ) {
$out .= sprintf " %3s | %s\n",
$diff,
join ' ', grep { is_prime($_) && $_ <= 1000 } nth_series_cubic( $diff ) ;
}
say<<~"END";
series | sequence
-------+------------------------------------------
$out
END
|
