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


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

<˜—
better upper bound will make it faster
,
F˜Á
using machine language subroutine
5
P˜Ê

V
Z€­À™"Maillet's determinant"
a
dŸÞñ

n¤„PPñ2)„DPñ2)À˜2primes for modulus and DET mod p
¾
x¤‚Detpñô
)À€œ"detp"‚Detp
×
‚§"prime (P>5)"ÁƒP

ŒƒPրуP)Ž™"not a prime."À‡ñ‚
–¨€š
? ƒNÑ(ƒPÞñ
)èñÞñÀƒNNÑ(ƒNÝñ)èñåñ
sª¤‚AƒNÞñ
ƒNNÞñ
)‚WƒNÞñ
ƒNNÞñ
)
…´ˆõCARLITZ.ò
Œ¾™
ªÈ™"Maillet's determinant"
ÍÒ™"=relative class number of"Á
ñÜ™ƒPÁ"th cyclotomic field is"
æ™À™€À…D)À™
𙀚
úª‚W‚A
#
¿
,
˜
=
õCARLITZ
F"
˜ø
Z,
˜make matrix
c6
˜
z@
™"Making matrix"Á
J
‰RÑðŠƒNÞñ

›T
™"."Á
±^
‰SÑðŠƒNÞñ

ûh
‚ARS)Ñ(RÝñ)å(SÝñ)èƒPÞ(RÝñ)å(SÝñ)èƒP
r
Œ
|
ŒÀ™
†
˜©
,
˜Áget upper bound
5š
˜Ê
Q¤
…UboundѳõHadamardcò
Z®
˜ï
w¸
˜ decide primes to use
€Â
˜
ŽÌ
PcÑñŽ
½Ö
ƒPPр΁Pc)À…MуPPÀ„PPð)уPPÀPRÑð
×à
«…MÔñå…Ubound
ê
¡PcÀƒPPр΁Pc)À…MхMåƒPPÀ PRÀ„PPPR)уPP
ô
‘S
Pþ
™"Use"ÁPRÝñ
ÁPRŽ™"primes:"ß™"prime:"
Y˜î
†˜calculate determinant for each prime
˜$
¡&‰IÑðŠPR
Þ0¢‚WðÃNÞñ
ÂðÃNÞñ
)Ñ¢‚AðÃNÞñ
ÂðÃNÞñ
)
:¢‚WðÃNÞñ
ÂðÃNÞñ
)ÞÑñÿ
D¬Içñåñ
<NƒPPфPPI)À™"DET mod"ÁƒPPÁ
aX®‚Detp‚WðÂð)ƒNƒPPƒDP)
“b„DPI)уDPÀ™"is"ÁƒDPÁÀ"IçñÑñŽ™
¬lŒÀPRçñÖñŽ™
µv˜J
Ö€˜kdecide exact determinant
ߊ˜t
ñ”ˆõChineseHò
ž…D҅UboundŽ…DхDޅM
¨•
&²˜
8¼õHadamard
UƺIJ…UboundÑñ

kЉIÑðŠƒNÞñ

wÚ…WÑð
ä‰JÑðŠƒNÞñ

¯î
…WхW݂AIJ)ïñ
¶øŒ
Ò…UboundхUboundå…W
Ù Œ
í•(€Æ…Ubound))
ö ˜
*õChinese
4ºI
/>…MDфPPð)À…DфDPð)
CH‰IÑñ
ŠPR
fRƒPPфPPI)ÀƒD1фDPI)
‚\…WрۅMDƒPP)å…MD
–f…MDхMDåƒPP
Áp…DÑ(ƒD1å…W݅Då(ñ
ޅW))ç…MD
ÈzŒ
Ï„•