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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
|
#!/usr/bin/env perl
use strict;
use warnings;
use Test::More;
use Math::Prime::Util::GMP qw/primorial pn_primorial factorial factorialmod
multifactorial subfactorial factorial_sum
addmod/;
my @small_primes = qw/
2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71
73 79 83 89 97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173
179 181 191 193 197 199 211 223 227 229 233 239 241 251 257 263 269 271 277 281
283 293 307 311 313 317 331 337 347 349 353 359 367 373 379 383 389 397 401 409
419 421 431 433 439 443 449 457 461 463 467 479 487 491 499 503 509 521 523 541
547 557 563 569 571 577 587 593 599 601 607 613 617 619 631 641 643 647 653 659
661 673 677 683 691 701 709 719 727 733 739 743 751 757 761 769 773 787 797 809
811 821 823 827 829 839 853 857 859 863 877 881 883 887 907 911 919 929 937 941
/;
sub nth_prime {
my $n = shift;
return 0 if $n <= 0;
die "Out of range for fake nth_prime: $n" unless defined $small_primes[$n-1];
$small_primes[$n-1];
}
my @pn_primorials = qw/
1
2
6
30
210
2310
30030
510510
9699690
223092870
6469693230
200560490130
7420738134810
304250263527210
13082761331670030
614889782588491410
32589158477190044730
1922760350154212639070
117288381359406970983270
7858321551080267055879090
557940830126698960967415390
40729680599249024150621323470
3217644767340672907899084554130
267064515689275851355624017992790
23768741896345550770650537601358310
2305567963945518424753102147331756070
232862364358497360900063316880507363070
23984823528925228172706521638692258396210
2566376117594999414479597815340071648394470
279734996817854936178276161872067809674997230
31610054640417607788145206291543662493274686990
/;
my @factorials = (qw/
1
1
2
6
24
120
720
5040
40320
362880
3628800
39916800
479001600
6227020800
87178291200
1307674368000
20922789888000
355687428096000
6402373705728000
121645100408832000
2432902008176640000
51090942171709440000
1124000727777607680000
25852016738884976640000
620448401733239439360000
15511210043330985984000000
403291461126605635584000000
10888869450418352160768000000
304888344611713860501504000000
8841761993739701954543616000000
265252859812191058636308480000000
/);
plan tests => 1 # factorial
+ 1 # factorialmod
+ 2 # primorial and pn_primorial
+ 2 # extra primorial tests
+ 1 # subfactorial
+ 1 # factorial_sum
+ 4; # multifactorial
{
my @fact = map { factorial($_) } 0 .. $#factorials;
is_deeply( \@fact, \@factorials, "factorial 0 .. $#factorials" );
}
is_deeply( [map { factorialmod($_>>2,$_) } 10000 .. 10100],
[map { addmod(factorial($_>>2),0,$_) } 10000 .. 10100],
"factorialmod" );
{
my @prim = map { primorial(nth_prime($_)) } 0 .. $#pn_primorials;
my @pnprim = map { pn_primorial($_) } 0 .. $#pn_primorials;
is_deeply(\@prim, \@pn_primorials, "primorial(nth(...)) 0 - $#pn_primorials");
is_deeply(\@pnprim, \@pn_primorials, "pn_primorial(...) 0 - $#pn_primorials");
}
is( primorial(100), '2305567963945518424753102147331756070', "primorial(100)");
is(
primorial(541),
'4711930799906184953162487834760260422020574773409675520188634839616415335845034221205289256705544681972439104097777157991804380284218315038719444943990492579030720635990538452312528339864352999310398481791730017201031090',
"primorial(541)"
);
is_deeply( [ map { subfactorial($_) } 0..23 ],
[qw/1 0 1 2 9 44 265 1854 14833 133496 1334961 14684570 176214841 2290792932 32071101049 481066515734 7697064251745 130850092279664 2355301661033953 44750731559645106 895014631192902121 18795307255050944540 413496759611120779881 9510425471055777937262/],
"subfactoral(n) for 0..23" );
is_deeply( [ map { factorial_sum($_) } 0..22 ],
[qw/0 1 2 4 10 34 154 874 5914 46234 409114 4037914 43954714 522956314 6749977114 93928268314 1401602636314 22324392524314 378011820620314 6780385526348314 128425485935180314 2561327494111820314 53652269665821260314/],
"factorial_sum(n) for 0..22" );
is_deeply( [ map { multifactorial($_,0) } 0..22 ],
[ map { 1 } 0..22 ],
"multifactorial(n,0) for 0..22" );
is_deeply( [ map { multifactorial($_,1) } 0..22 ],
[ map { factorial($_) } 0..22 ],
"multifactorial(n,1) for 0..22" );
is_deeply( [ map { multifactorial($_,2) } 0..26 ],
[qw/1 1 2 3 8 15 48 105 384 945 3840 10395 46080 135135 645120 2027025 10321920 34459425 185794560 654729075 3715891200 13749310575 81749606400 316234143225 1961990553600 7905853580625 51011754393600/],
"multifactorial(n,2) for 0..26" );
is_deeply( [ map { multifactorial($_,3) } 0..29 ],
[qw/1 1 2 3 4 10 18 28 80 162 280 880 1944 3640 12320 29160 58240 209440 524880 1106560 4188800 11022480 24344320 96342400 264539520 608608000 2504902400 7142567040 17041024000 72642169600/],
"multifactorial(n,3) for 0..29" );
|
