🧩 问题核心与性质分析

解决此题，首先要理解清楚金字塔的构造规则和关键性质：

1. 基本单位：金字塔由 $1 \times 1 \times 1$ 的立方体石块构成。
2. 分层结构：金字塔有 $h$ 层，从下往上编号（第1层是最底层）。每层是一个联通块，且高层石块必须有低层石块作为支撑（不能悬空）。
3. 对称性：每一层都满足严格的左右对称和上下对称，并且所有层的对称轴是重合的。
4. 尺寸要求：每一层的最大长度和最大宽度相等，且为偶数。题目会给出从底层开始的每一层的最大长度 $l_i$，且 $l_i$ 非严格递减。
5. 总石块数：金字塔的总石块数等于 $n$。

数据范围：

• $n \leq 1000$
• $2 \leq l \leq 20$
• $1 \leq h \leq 10$
• 答案对 $1000000007$ 取模。

由于直接计算满足所有条件的金字塔数量比较困难，我们可以考虑利用容斥原理来简化问题。

💡 动态规划结合容斥原理的思路

我们可以通过容斥原理将复杂的连通性条件转化为更容易处理的形式。

🔹 步骤1：忽略连通性条件

首先，我们暂时忽略“每一层必须是联通块”的条件，只考虑对称性、支撑关系和石块总数限制。这样问题就变成了：计算满足以下条件的石块放置方案数：

• 每层在平面上左右、上下对称。
• 高层石块有底层石块支撑。
• 每层的尺寸是给定的偶数且非严格递减。
• 总石块数为 $n$。

这个问题可以通过动态规划来求解。定义 $dp[i][j]$ 表示考虑前 $i$ 层（从下往上），总共使用了 $j$ 个石块的方案数。转移时，我们需要枚举第 $i$ 层放置的石块数 $k$（需满足对称性和支撑条件），然后更新 $dp[i][j] += dp[i-1][j-k]$。

🔹 步骤2：容斥原理应用

在得到了忽略连通性的方案数后，我们需要减去那些不满足连通性的方案。根据容斥原理，我们可以：

1. 计算至少包含 $k$ 个孤立连通块的方案数（赋予权重 $(-1)^{k-1}$ 进行容斥）。
2. 或者，枚举最底层的哪个部分（例如，第一象限的某个特定区域）是不被包含在主要连通块内的，然后减去这些情况。

这通常需要通过枚举断裂点来实现：在每一层的四分之一象限（利用对称性）中，枚举一个分割位置，使得该位置之上的部分与主要部分断开。然后，根据容斥原理，对所有这些不连通的情况进行加权求和。

🔹 步骤3：动态规划状态设计

结合容斥原理，我们的动态规划状态需要包含更多信息。一种可能的状态设计是： $dp[i][j][mask]$ 表示处理到第 $i$ 层，总共使用了 $j$ 个石块，且当前层的连通性状态为 $mask$ 的方案数。 这里的 $mask$ 可以用来表示当前层在四分之一象限中的连通情况，但由于层数较多且需要对称性，实现起来较为复杂。在实际操作中，我们通常需要逐行或逐列递推，并利用括号序列或最小表示法来记录连通性。

🧮 复杂度分析与优化

• 状态数：由于 $l \leq 20$，每层的四分之一象限尺寸不超过 $10 \times 10$，但连通性状态数可能随着宽度指数级增长。因此，我们需要使用插头DP或轮廓线DP等技术来高效处理连通性。
• 优化：利用对称性减少状态数；预处理合法的状态转移；使用滚动数组优化空间。

⚠️ 实现细节提示

1. 对称轴处理：在计算整个层的石块总数时，要小心对称轴上的石块不要重复计算。对于一个 $l \times l$ 的层（$l$ 为偶数），其石块总数可以通过四分之一象限的石块数 $cnt$ 计算为：$total = 4 \times cnt - 2 \times (l/2) + 1$？更稳妥的方法是直接模拟整个层的对称填充。
2. 支撑关系判断：在状态转移时，需要确保第 $i$ 层的每个石块都有第 $i-1$ 层的石块支撑。这可以通过比较相邻两层的形状来判断。
3. 模块化：将问题分解为"生成所有合法的单层四分之一结构"和"计算层与层之间的转移方案数"两个模块，会使代码更清晰。

💎 总结

这道题是一个比较复杂的计数问题，结合了对称性、连通性、拓扑约束（支撑关系）等多个条件，并需要运用动态规划和容斥原理。