康托展开算法解析与树状数组优化实现

1. 康托展开算法解析与实现

1.1 康托展开原理详解

康托展开是一种将全排列映射到自然数集合的双射算法，它能够计算给定排列在所有可能排列中的字典序排名。这个算法在组合数学和计算机科学中有着重要应用，特别是在需要处理排列相关问题时。

康托展开的核心思想是基于阶乘的权重计算。对于一个长度为N的排列，其康托展开值X可以通过以下公式计算：

X = a₁×(N-1)! + a₂×(N-2)! + ... + a_{N-1}×1! + a_N×0!

其中aᵢ表示在第i位之后比当前数字小的数字个数。这个公式实际上是在计算当前排列之前有多少个"更小"的排列。

举个例子，对于排列[2,1,3]：

• 第一位是2，比2小的未使用数字有1（共1个），贡献1×2! = 2
• 第二位是1，比1小的未使用数字有0个，贡献0×1! = 0
• 第三位是3，固定不变，贡献0
  最终排名为2+0+0+1=3（加1是因为排名从1开始计数）

1.2 树状数组优化实现

直接计算aᵢ（即后面比当前元素小的数字个数）的时间复杂度是O(N²)，对于N=1,000,000的数据规模显然不可行。因此我们需要使用树状数组（Fenwick Tree）来优化这一过程。

树状数组能够在O(logN)时间内完成以下操作：

1. 单点更新：将某个位置的值加1
2. 前缀查询：查询前i个元素的和

在本题中的具体应用步骤：

1. 初始化树状数组所有位置为0
2. 从右向左遍历排列
3. 对于每个元素a[i]，查询树状数组中1到a[i]-1的和，这就是当前元素右侧比它小的元素个数
4. 将a[i]位置的值加1（标记该数字已使用）
5. 将查询结果乘以对应的阶乘，累加到最终结果中

这种优化将算法的时间复杂度降低到了O(N logN)，能够处理大规模数据。

2. 代码实现细节解析

2.1 快速输入输出优化

cpp复制inline void read( rgt int &x ){
    x = 0; while( !isdigit(*p) ) ++p;
    while( isdigit(*p) ) 
        x = x * 10 + ( *p & 15 ), ++p;
}

int main(){
    scanf( "%d", &N ), fac =