Große Zufallszahlen - Eintrittswahrscheinlichkeit

[OneWayUserAcci]

Hallo,
Ich bin am überlegen sowie auf der Suche nach einer Pseudo-Zufalls-Funktion, die für große Zahlen gedacht ist. Dabei soll man den Range-Wert einstellen können, wie bei (mt_)rand. Die bisherigen Funktionen, die ich gefunden habe, sind ein wenig zu sehr abweichend von der Verteilung der Zahlen in Bezug auf (mt_)rand.

Desweiteren such ich nach einer einfachereren Lösung für folgendes Szenario bzgl. Eintrittswahrscheinlichkeit bei unterschiedlichen Werten. Der Value-Parameter definiert hier die Chance des Eintretens, des jeweiligen Datensatzes.

$users = array(
    'A' => '100',
    'B' => '2000',
    'C' => '30000',
    'D' => '400000',
    'E' => '5000000',
);

$maximum = array_sum_bignum($users);
$number  = mt_rand_bignum('1', $maximum);
$users   = shuffle_keep_keys($users);
$count   = '0';
foreach($users AS $name => $value)
{
    $min   = bcadd($count, '1');
    $count = bcadd($count, $value);
    $max   = $count;

    if(($min <= $number) AND ($max >= $number))
    {
        //echo $min . '' . $max . '' . $value;
        break;
    }
} 

[combie]

Ich mag große Zahlen, auch Zufallszahlen und so wieter, aber von dem was du da möchtest, verstehe ich nix! Auch die Funktionen welche da bei dir zum Einsatz kommen habe ich noch nie gesehen.

[OneWayUserAcci]

Das sind zwei verschiedene Dinge. Beim Ersten gehts nur um den Generator, womit ich größere/kleinere Zahlen rauskriegen möchte als (-)2^31-1. Das Zweite hat mit dem Ersten nur in der Hinsicht was zu tun, weil es um größere Zufallswerte geht. Die Funktionen array_sum_bignum/mt_rand_bignum sind dieselbige wie array_sum/mt_rand, nur mit größeren Werten, wobei die Funktion für die großen Zufallszahlen noch nicht vorhanden ist. Die Funktion shuffle_keep_keys sorgt nur dafür, dass beim Sortieren mit shuffle die Array-Schlüssel erhalten bleiben.

[h3ll]

Das sind zwei verschiedene Dinge. Beim Ersten gehts nur um den Generator, womit ich größere/kleinere Zahlen rauskriegen möchte als (-)2^31-1.

echo mt_rand(0, pow(2, 60));
// 1070457797516197888

echo mt_rand(0, pow(2, 60));
// 426054254793326592

echo mt_rand(0, pow(2, 60));
// 939119261509484544 

[combie]

@h3ll
PHP: mt_getrandmax - Manual

Da hasste wohl ein 64 Bit System
Das hat aber nicht jeder.

echo mt_rand(0, pow(2, 60)).'<br>';
echo mt_rand(0, pow(2, 60)).'<br>';
echo mt_rand(0, pow(2, 60)).'<br>';
echo mt_rand(0, pow(2, 60)).'<br>';
echo mt_rand(0, pow(2, 60)).'<br>'; 

Liefert bei mir IMMER und konsequent 0. Nix anderes.

[h3ll]

@h3ll
PHP: mt_getrandmax - Manual

Da hasste wohl ein 64 Bit System
Das hat aber nicht jeder.

Tja, wer hat, der hat

64-Bit Prozessoren kriegt man inzwischen eh schon nachgeschmissen.

[combie]

Aus meiner Wühlkiste:

function bc_rand($start,$end,$scale=null)
{

    if(is_null($scale))
    {
        if(bccomp($start,$end) >= 0)
          throw new Exception('Left Value must be lower');
    }else
    {
        if(bccomp($start,$end,$scale) >= 0)
          throw new Exception('Left Value must be lower');
    }
    $result = '';
    if(is_null($scale))
    {
        $diff = bcsub($end,$start);
    }else
    {
        $diff = bcsub($end,$start,$scale);
    }
    $begrenzt = true;
    for($i=0;$i<strlen($diff);$i++)
    {
        $ziffer = $diff{$i};
        if($ziffer === '.')
        {
          $result .= '.';
        }else
        {
          if($begrenzt)
          {
           $neue = mt_rand(0,$ziffer);
           if($neue < $ziffer) $begrenzt = false;
          }else
          {
           $neue = mt_rand(0,9);
          }
          $result .= $neue;
        }
    
    }
      
    $result = ltrim($result,'0');
    if(0 === strpos($result,'.')) $result = '0'.$result;
    
    if(is_null($scale))
    {
        return bcadd($start,$result);
    }else
    {
        return bcadd($start,$result,$scale);
    }
}

echo bc_rand('-34534345645635235346332454363453454564','235232389573642345273523451642',6); 

[OneWayUserAcci]

Nun tritt wieder das Problem der Verteilung auf s.u.. Ich nehm mal an, dass im PHP C-Quelltext keine extra Funktion dafür gebaut wurde und einfach eine von C genommen wurde, daher die Begrenzung?! Und sowas wie Zend_Rand gibts glaub auch net, mal schauen, irgendwie kriegt man das schon zamm gewurschtelt.

Array
(
    [10] => 5085
    [7] => 532
    [4] => 524
    [0] => 505
    [8] => 491
    [3] => 485
    [5] => 484
    [9] => 481
    [2] => 480
    [1] => 472
    [6] => 461
) 

[OneWayUserAcci]

Ich hab mich für die folgende Funktion entschieden, welche ich gefunden habe.

function bcrand($min, $max)
{
    if(bccomp($max, $min) != 1)
    {
        return 0;
    }
    $top    = bcsub($max, $min);
    $rand   = bcadd($top, 1);
    $length = strlen($top);
    while(bccomp($rand, $top) == 1)
    {
        unset($rand_part);
        $n = 0;
        while(9 * $n <= $length)
        {
            if($length - 9 * $n >= 9)
            {
                $rand_part[] = mt_rand(0, 999999999);
            }
            else
            {
                $j = 0;
                $foo = '';
                while($j < $length - 9 * $n)
                {
                    $foo .= '9';
                    ++$j;
                }
                $foo += 0;
                $rand_part[] = mt_rand(0, $foo);
            }
            ++$n;
        }
        $i     = 0;
        $rand  = '';
        $count = count($rand_part);
        while($i < $count)
        {
            $rand .= $rand_part[$i];
            ++$i;
        }
    }
    return bcadd($rand, $min, 0);
}