Category : BASIC Source Code
Archive   : UBAS830.ZIP
Filename : MAILLET3.UB


˜MAILLET3
D˜Carlitz's version of Maillet's determinant
˜cf. Rev Numb Th LeVEQUE, C20-18 by Carlitz & Olson
µ(˜sample of calculation of integer determinants
Ó2˜by using multi-modulo

<˜better upper bound will make it faster

F˜
,
P€­À™"Maillet's determinant"
7
ZŸÞñ
s
d¤„PPñ2)„DPñ2)À˜primes for modulus and DET mod p
Ž
n§"prime (P>5)"ÁƒPÀ™
¸
xƒPրуP)Ž™"not a prime."À‡ñn
Á
‚¨€š
Þ
ŒƒNÑ(ƒPÞñ
)èñÞñ
–¤‚AƒNÞñ
ƒNÞñ
)„WƒNÞñ
ƒNÞñ
)
$ ˆõCARLITZXò
+ª™
I´™"Maillet's determinant"
k¾™"=relative class number of"
È™ƒPÁ"th cyclotomic field is"
Ò™€À…D)
¨Ü™À™€š
¸æª„W‚A
¿ð¿
Èú˜
Ù
õCARLITZ
â
˜
ö
˜make matrix
ÿ"
˜
,
‰RÑðŠƒNÞñ

+6
‰SÑðŠƒNÞñ

u@
‚ARS)Ñ(RÝñ)å(SÝñ)èƒPÞ(RÝñ)å(SÝñ)èƒP
’J
™€ìñ)‚ARS)Á
›T
ŒÀ™
¤^
ŒÀ™
­h
˜
År
˜get upper bound
Î|
˜
ê†
…UboundѳõHadamardò
ó
˜
š
˜decide primes to use
¤
˜
K®
ƒPPрÞñïñ)À…MуPPÀ„PPð)уPPÀPRÑð
e¸
…MÔñå…Ubound
›Â
ƒPPрރPP)À…MхMåƒPPÀ PRÀ„PPPR)уPP
¤Ì
‘
ØÖ
™"Use"ÁPRÝñ
ÁPRŽ™"primes:"™"prime:"
áà
˜
ê
˜calculate determinant for each prime
ô
˜
)þ
‰IÑðŠPR
f¢„WðÃNÞñ
ÂðÃNÞñ
)Ñ¢‚AðÃNÞñ
ÂðÃNÞñ
)
‹ƒPPфPPI)À™"DET mod"ÁƒPPÁ
¥ƒDPѳõDETP;ò(ƒPP)
Ä&„DPI)уDPÀ™"is"ÁƒDP
Ë0Œ
Ô:˜
õD˜decide exact determinant
þN˜
XˆõChineserò
5b…D҅UboundŽ…DхDޅM
<l•
Ev˜
ˆ€õDETP(ƒP)À˜ assume every elements >=0 and §ŠºIJKƒDPƒA
DPÑñ
À˜ determinant mod p
ž‰JÑðŠƒNÞñ
G¨‰IсJŠƒNÞñ
À˜ search non_zero
f² „WIJ)Ž”‰À‡ñÐ
m¼Œ
Æ”‰ÀƒDPÑðÀ‡ñH
™ÐIցJŽ´ƒDP
ÑÚÀ¥¢„WIJÃNÞñ
)„WJJÃNÞñ
)
êäƒAфWJJ)
îƒDPуDPåƒAçƒP
%øƒInvрۄWJJ)ƒP)
B‰IсJÝñ
ŠƒNÞñ

[  ƒCфWIJ)
} ƒCŽƒCуCåƒInvçƒP
Ö À‰KсJÝñ
ŠƒNÞñ
À„WIK)Ñ(„WIK)ރCå„WJK))çƒPÀŒ
Ý*Œ
ä4Œ
 >ƒDPуDPå„WƒNÞñ
ƒNÞñ
)çƒP
 H•(ƒDP)
% R˜
7 \õHadamard
T fºIJ…UboundÑñ

j p‰IÑðŠƒNÞñ

v z…WÑð
Œ „‰JÑðŠƒNÞñ

® Ž
…WхW݂AIJ)ïñ
µ ˜Œ
Ñ ¢…UboundхUboundå…W
Ø ¬Œ
ì ¶•(€Æ…Ubound))
õ À˜

ÊõChinese

ÔºI
.
Þ…MDфPPð)À…DфDPð)
B
艁IÑñ
ŠPR
e
òƒPPфPPI)ÀƒD1фDPI)

ü…WрۅMDƒPP)å…MD
•
…MDхMDåƒPP
À
…DÑ(ƒD1å…W݅Då(ñ
ޅW))ç…MD
Ç
Œ
Î
$•