File: lucasuv.pl

package info (click to toggle)
libmath-prime-util-perl 0.73-2
  • links: PTS, VCS
  • area: main
  • in suites: bookworm, forky, sid, trixie
  • size: 2,796 kB
  • sloc: perl: 24,676; ansic: 11,471; makefile: 26; python: 24
file content (15 lines) | stat: -rw-r--r-- 370 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
#!/usr/bin/env perl
use warnings;
use strict;
use ntheory ":all";
use Math::GMP;

for my $n (8539783 .. 8539788) {
  for (1..2000) {
    my($u,$v) = lucas_sequence($n, 1, -1, $_);
    my $u1 = Math::GMP->new(lucasu(1,-1,$_)) % $n;
    my $v1 = Math::GMP->new(lucasv(1,-1,$_)) % $n;
    die "U $_ $n $u $u1" unless $u==$u1;
    die "V $_ $n $v $v1" unless $v==$v1;
  }
}