按字典顺序n k

时间:2019-05-21 08:39:56

标签: php math permutation combinatorics

我想编写以字典顺序n选择k返回下一个排列的函数。

到目前为止,我发现以下算法


class Permutation
{
    public function swap(& $a, $first, $second)
    {
        $tmp = $a[$first];
        $a[$first] = $a[$second];
        $a[$second] = $tmp;
    }


    function nextPermutation(& $a, $n, $k)
    {
        do {
            $first = $n - 2;
            while ($first != -1 && $a[$first] >= $a[$first + 1]) $first--;
            if ($first == -1)
                return false;
            $second = $n - 1;
            while ($a[$first] >= $a[$second]) $second--;
            $this->swap($a, $first, $second);
            $left = $first + 1;
            $right = $n - 1;
            while ($left < $right) {
                $this->swap($a, $left++, $right--);
            }

        } while ($first > $k - 1);
        return true;

    }

    public function run()
    {
        $n=10;
        $k=4;
        $a = array_merge(range(ord(0), ord(9)), range(ord('A'), ord('Z')), range(ord('a'), ord('z')));
        $i = 0;
        while ($this->nextPermutation($a, $n, $k)) {
            $i++;
            echo($i . ")  ");
            for ($j = 0; $j < $k; $j++) {
                echo $a[$j] . " ";
            }
            echo("\n");
        }
    }

}

但是它仅在n到13时有效,然后变慢。 我需要使其适用于n = 62和k = 6。有可能做到吗?

1 个答案:

答案 0 :(得分:1)

此组合对象在英语组合语中没有特定名称,在俄语和法语组合语中称为“安排” A(n,k)

有44 261 653 680个排列A(62,6)-因此它们的生成和输出(〜500 GB)花费的时间太长。在像C / C ++ / Delphi这样的编译语言中(无输出)生成可能要花费一些时间(大约每秒估计10 ^ 6-10 ^ 7个结果),但是在PHP(解释语言)中,此过程应该较慢。

您的代码使用常规置换算法,根据所需的n!排列生成n!/(n-k)!结果。

用于生成安排的简单递归Python代码。对于n = 25,k = 5,工作约3秒。不是最佳选择,在书中会更好地检查。

n = 4
k = 2
ar = [0] * k
used = set()

def genArr(idx):
    for i in range(n):
        if not i in used:
            used.add(i)
            ar[idx] = i
            if idx == k - 1:
                print(ar) #comment for large n
            else:
                genArr(idx + 1)
            used.remove(i)

genArr(0)
print('OK')

[0, 1, 2]
[0, 1, 3]
[0, 2, 1]
[0, 2, 3]
......
[2, 3, 1]
[3, 0, 1]
[3, 0, 2]
[3, 1, 0]
[3, 1, 2]
[3, 2, 0]
[3, 2, 1]