Perhaps some one can find useful this function to compute the modular inverse of a integer (extended euclidean algorithm):
function invmod($a,$b) {
$n=$b;
$x=0; $lx=1; $y=1; $ly=0;
while ($b) {
$t=$b;
$q=bcdiv($a,$b,0);
$b=bcmod($a,$b);
$a=$t;
$t=$x; $x=bcsub($lx,bcmod(bcmul($q,$x),$n)); $lx=$t;
$t=$y; $y=bcsub($ly,bcmod(bcmul($q,$y),$n)); $ly=$t;
}
if (bccomp($lx,0) == -1)
$lx=bcadd($lx,$n);
return $lx;
}
// verify
$n="2447995268898324993537772139997802321";
$t="64941057316178801556773346239351236811";
$m="123456789";
$i=invmod($t,$n);
// (t*m)*inv(t) is m
echo bcmod(bcmul(bcmod(bcmul($t,$m),$n),$i),$n) == $m;
bcdiv
(PHP 4, PHP 5)
bcdiv — Division zweier Zahlen beliebiger Genauigkeit
Beschreibung
string bcdiv
( string $left_operand
, string $right_operand
[, int $scale
] )
Dividiert left_operand durch right_operand .
Parameter-Liste
- left_operand
-
Der linke Operand in Stringform.
- right_operand
-
Der rechte Operand in Stringform.
- scale
-
Dieser optionale Parameter kann benutzt werden um die Anzahl der Nachkommastellen im Ergebnis festzulegen. Sie können mit Hilfe der bcscale() Funktion auch einen globalen Vorgabewert für alle Funktionen festlegen.
Rückgabewerte
Gibt das Ergebnis der Division als String zurück. Wenn der Parameter right_operand 0 ist, wird NULL zurückgegeben.
Beispiele
Beispiel #1 bcdiv()-Beispiel
<?php
echo bcdiv('105', '6.55957', 3); // 16.007
?>
bcdiv
MM
18-Jan-2008 03:29
18-Jan-2008 03:29
Gautam
29-Aug-2007 07:56
29-Aug-2007 07:56
<?php
//converting in to required precision of decimal points
$result= bcdiv(89.99999999997,2.57865741235478,2);
echo "$result"; // 34.90 result with 2 decimal points
?>
cristianDOTzuddas]NOSPAM[gmailDOTcom
24-Jul-2005 05:10
24-Jul-2005 05:10
Decimal to binary conversion, using BC Math.
Note: this function is VERY slow if the decimal number is too big!
<?
function bc_decbin($dec_str) {
if (strlen($dec_str)>0) {
$bin_str = '';
do {
if (((int)$dec_str[strlen($dec_str)-1] % 2) === 0)
$bin_str .= '0';
else
$bin_str .= '1';
$dec_str = bcdiv($dec_str, '2');
} while ($dec_str!='0');
return strrev($bin_str);
}
else
return null;
}
?>
-----
Cristian
www.CodeFlower.com
