Python程序：求解可产生正确排列的数字之和

假设给定一个数字n，要求编写所有可能的正整数（至n）排列。这些排列按字典序排序，编号从1到n!。其中一个排列被认为是特殊排列。在这个特殊排列中，某些值可能被遗忘，用0代替。我们需要找到可以等于原始排列的排列，然后将它们的编号相加，并将和作为程序的输出。

例如，如果输入的特殊排列 (input_arr) = [0, 2, 0]，n = 3，则输出为7。可能有两种可能的排列：[1, 2, 3]和[3, 2, 1]。这些排列的编号分别为2和5。因此，结果为5 + 2 = 7。

为了解决这个问题，我们将遵循以下步骤：

  • mod := 10^9 + 7
  • i2 := 2 ^ (modulo - 2) mod modulo
  • fact := 一个新的列表，包含值1
  • 对于 x 从 1 到 n+1，执行：
    • 在fact列表末尾插入 (fact的最后一个元素 * x) mod modulo
  • cnt := input_arr中0的出现次数
  • 如果 cnt 等于零，则：
    • res := 0
    • seen_list := 一个新的列表
    • 对于每个索引 i 从 1 开始和 input_arr 中的元素 x，执行：
      • tmp_val := 将x插入seenList保持排序顺序的索引
      • res := res + fact[n-i] *(x - 1 - tmp_val)
      • res := res mod modulo
      • 将x插入seen_list在位置tmp_val
    • 返回 res + 1
  • 否则：
    • ik := (cnt ^ (modulo-2)) mod modulo
    • miss := 一个大小为n的新列表，初始化为True
    • 对于 input_arr 中的每个 x，执行：
      • 如果 x 不等于 0，则：
        • miss[x - 1] := False
    • miss_srtd := 一个新的列表
    • tmp := 0
    • 对于每个索引 i 从 1 开始和 miss 中的元素 x，执行：
      • 如果 x 不为零，则：
        • 在 miss_srtd 末尾插入 i
        • tmp := tmp + i
    • pre := 一个初始化为0的新列表
    • 对于 miss 中的每个 x，执行：
      • 在 pre 末尾插入 pre[-1] + x
    • cnt_cu := 0
    • s := tmp mod modulo * ik mod modulo
    • srtdw := 一个新的列表
    • res := 0
    • z := 0
    • 对于每个索引 i 从 1 开始和 input_arr 中的元素 x，执行：
      • 如果 x 不为零，则：
        • l := 将 x 插入 srtdw 保持排序顺序的位置
        • tmp_val := 将 x 插入 srtdw 保持排序顺序的位置
        • l := l + z * (将 x 插入 miss_srtd 保持排序顺序的位置) mod modulo * ik mod modulo
        • p := x - 1 - l
        • p := p * fact[cnt]
        • p := p mod modulo
        • 将 x 插入 srtdw 在位置 tmp_val
        • cnt_cu := cnt_cu + cnt - pre[x]
      • 否则：
        • l := cnt_cu
        • l := l * ik
        • l := l + z * i2 mod modulo
        • p := s - 1 - l
        • p := p * fact[cnt]
        • p := p mod modulo
        • z := z + 1
        • res := res + p * fact[n-i] mod modulo
        • res := res mod modulo
    • 返回 (res + fact[cnt]) mod modulo

示例

让我们看下面的实现来更好地理解：

import bisect

def solve(input_arr, n):

   modulo = 10 ** 9 + 7
   i2 = pow(2, modulo-2, modulo)
   fact = [1]
   for x in range(1, n+1):
      fact.append(fact[-1] * x % modulo)

   cnt = input_arr.count(0)

   if not cnt:
      res = 0
      seen_list = []
      for i, x in enumerate(input_arr, 1):
         tmp_val = bisect.bisect(seen_list, x)
         res += fact[n-i] * (x - 1 - tmp_val)
         res %= modulo
         seen_list.insert(tmp_val, x)
      return res + 1
   else:
      ik = pow(cnt, modulo-2, modulo)
      miss = [True] * n
      for x in input_arr:
         if x != 0: miss[x-1] = False
      miss_srtd = []
      tmp = 0
      for i, x in enumerate(miss, 1):
         if x:
            miss_srtd.append(i)
            tmp += i
      pre = [0]
      for x in miss:
         pre.append(pre[-1] + x)
      cnt_cu = 0
      s = tmp % modulo * ik % modulo
      srtdw = []
      res = z = 0
      for i, x in enumerate(input_arr, 1):
         if x:
            l = tmp_val = bisect.bisect(srtdw, x)
            l += z * bisect.bisect(miss_srtd, x) % modulo * ik % modulo
            p = x - 1 - l
            p *= fact[cnt]
            p %= modulo
            srtdw.insert(tmp_val, x)
            cnt_cu += cnt - pre[x]
         else:
            l = cnt_cu
            l *= ik
            l += z * i2 % modulo
            p = s - 1 - l
            p *= fact[cnt]
            p %= modulo
            z += 1
         res += p * fact[n-i] % modulo
         res %= modulo
      return (res + fact[cnt])%modulo

print(solve([0, 2, 0], 3))

输入

[0, 2, 0], 3

输出

7

更新于：2021年10月11日

浏览量：135

开启你的职业生涯

完成课程获得认证

开始学习
广告