danog/PrimeModule
1
0
Fork 0
mirror of https://github.com/danog/PrimeModule.git synced 2024-11-26 20:34:37 +01:00
Code Issues Projects Releases Wiki Activity

Implement GNU factor factorization

Browse Source
This commit is contained in:
Daniil Gentili 2023-04-12 16:19:09 +02:00
parent 5e2718b969
commit 1fcfddc253
2 changed files with 38 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
lib/danog/PrimeModule.php
View File

@ -112,6 +112,33 @@ class PrimeModule
return false;
}
public static function factor_single($what)
{
if (function_exists('shell_exec')) {
$res = trim(shell_exec('timeout 10 factor '.$what.' 2>&1') ?? '');
$res = explode(' ', $res);
if (count($res) !== 3) {
return false;
}
return (int) $res[1];
}
return false;
}
public static function factor($what)
{
$res = [self::factor_single($what)];
if ($res[0] === false) {
return false;
}
while (array_product($res) !== $what) {
$res[] = self::factor_single($what / array_product($res));
}
return $res;
}
public static function python($what)
{
$res = [self::python_single($what)];

14
tests/testing.php
View File

@ -28,6 +28,10 @@ function test($n)
$GLOBALS['medium']['native_cpp'] += $time;
echo '| '.str_pad($result, 6, ' ', STR_PAD_RIGHT).' | native cpp | '.str_pad($time, 20, ' ', STR_PAD_RIGHT).' |'.PHP_EOL;
list($time, $result) = get_time(['\danog\PrimeModule', 'factor'], $n, true);
$GLOBALS['medium']['factor'] += $time;
echo '| '.str_pad($result, 6, ' ', STR_PAD_RIGHT).' | factor | '.str_pad($time, 20, ' ', STR_PAD_RIGHT).' |'.PHP_EOL;
list($time, $result) = get_time(['\danog\PrimeModule', 'python_alt'], $n, true);
$GLOBALS['medium']['python_alt'] += $time;
echo '| '.str_pad($result, 6, ' ', STR_PAD_RIGHT).' | python alt | '.str_pad($time, 20, ' ', STR_PAD_RIGHT).' |'.PHP_EOL;
@ -58,6 +62,10 @@ function test_single($n, $messy = false)
$GLOBALS['medium']['native_cpp'] += $time;
echo '| '.str_pad($result, strlen($n), ' ', STR_PAD_RIGHT).' | native cpp | '.str_pad($time, 20, ' ', STR_PAD_RIGHT).' |'.PHP_EOL;
list($time, $result) = get_time(['\danog\PrimeModule', 'factor_single'], $n);
$GLOBALS['medium']['factor'] += $time;
echo '| '.str_pad($result, strlen($n), ' ', STR_PAD_RIGHT).' | factor | '.str_pad($time, 20, ' ', STR_PAD_RIGHT).' |'.PHP_EOL;
list($time, $result) = get_time(['\danog\PrimeModule', 'python_single_alt'], $n);
$GLOBALS['medium']['python_alt'] += $time;
echo '| '.str_pad($result, strlen($n), ' ', STR_PAD_RIGHT).' | python alt | '.str_pad($time, 20, ' ', STR_PAD_RIGHT).' |'.PHP_EOL;
@ -95,7 +103,7 @@ function gen_payload()
}
echo PHP_EOL.'----------- HUGE SEMIPRIME TESTS (100 semiprimes) ----------'.PHP_EOL;
$GLOBALS['medium'] = ['python' => 0, 'python_alt' => 0, 'wolfram' => 0, 'native' => 0, 'native_cpp' => 0];
$GLOBALS['medium'] = ['python' => 0, 'python_alt' => 0, 'wolfram' => 0, 'native' => 0, 'native_cpp' => 0, 'factor' => 0];
$tg = fsockopen('tcp://149.154.167.51:443');
fwrite($tg, chr(239));
stream_set_timeout($tg, 1);
@ -112,7 +120,7 @@ foreach ($medium as $type => $total) {
}
echo PHP_EOL.'------------------- SMALL MULTIPLE FACTOR TESTS -------------------'.PHP_EOL;
$GLOBALS['medium'] = ['python' => 0, 'python_alt' => 0, 'wolfram' => 0, 'native' => 0, 'native_cpp' => 0];
$GLOBALS['medium'] = ['python' => 0, 'python_alt' => 0, 'wolfram' => 0, 'native' => 0, 'native_cpp' => 0, 'factor' => 0];
foreach ([200, 327, 35, 13589] as $multiple) {
test($multiple);
@ -123,7 +131,7 @@ foreach ($medium as $type => $total) {
}
echo PHP_EOL.'------------------- HUGE SEMIPRIME TESTS (MESSY) ------------------'.PHP_EOL;
$GLOBALS['medium'] = ['python' => 0, 'python_alt' => 0, 'wolfram' => 0, 'native' => 0];
$GLOBALS['medium'] = ['python' => 0, 'python_alt' => 0, 'wolfram' => 0, 'native' => 0, 'factor' => 0];
$m = [1724114033281923457, 2189285106422392999, 3510535493073971807, 1756377470921216651, 1767867620107504387, 2149465210997855797];
foreach ($m as $messy) {
test_single($messy, true);