期望：~~平均 82 ~~‍公测为 80 抽一个当期 up，如果首充全拿完了，一个卡池充 1500 块（648元 × 1 + 198元 × 1 + 6元 × 9，共 120 抽）能保证拿下当期 up 2026/01/16 公测算法更新 ​公测加了个首次 30 抽送一张十连鬼抽（不影响任何水位和保底，也不受水位和保底影响），期望从 82 降到 80，代码变更推上去了，文章细节就不改了，意义不大，DP 矩阵里抽中了就右移没抽中就不动而已 三测的武器池是抽角色池的副产物，只能拿抽角色池得到的代币抽取武器池，暂时没有其他代币来源，因此武器池就先不管了，主要看角色池。 关于角色池抽取，这里有一个详细的机制讲解视频： 【明日方舟:终末地】最清晰抽卡机制讲解 3测迭代再次拉低二游定价？ （多说两句：玩家的总抽卡投入（即所谓“定价”）是要综合福利考虑的，目前三测的福利只能说非常一般，抽卡代价十分高昂。要是海猫早三年坐镇终末地项目组的话抽卡定价可能就不是现在这样了，可惜他在 all in 泡姆泡姆，坐不得） 综合一下，我们可以得到以下信息： 基础六星概率为 0.8 % ，第 66 抽开始每抽概率提升 5 % ，第 80 抽必出六星（不一定是当期 up，也没有大小保底） 如果当前池第 119 抽后依然没有抽到当期 up，则第 120 抽的结果一定为当期 up 当前池每 240 抽会获得一个潜能，但这是独立于抽卡概率的，不会影响第 240 抽的结果 当前池的角色会作为陪跑角色，和其他常驻角色一同出现在后续两个池当中 这里可以看出终末地的抽卡系统是一个高方差+强保底的卡池系统，因为方差高，所以玩家的抽卡过程会更加……刺激，而强保底则让玩家在抽卡过程中不会过于绝望，从而保持一定的抽卡动力（抛开抽卡代币的获取难度不谈的话）。 需要注意在第 80 抽的时候，纸面概率实际上是 0.008 + ( 80 − 65 ) × 0.05 = 0.758 < 1 ，因此在暂时不考虑首个 120 抽保底的情况下，在单个循环（自上一次抽到六星后开始计算，到下一次抽到六星前结束）中第 n 抽出六星的概率 P d r o p ( n ) 计算公式应该分三段考虑，如下： 对于每个六星，还需要计算它是否是当期 up。由于六星中是当期 up 的概率恒定为 50%，所以可以生成一个随机数 r t y p e ∼ U [ 0 , 1 ) ，如果 r t y p e < 0.5 ，则认为该六星是当期 up，否则认为该六星是常驻六星或者陪跑角色（反正不是玩家想要的那个）。 接着引入首个 120 抽保底机制，在每一抽进入通用公式判定之前先判定一次当前积攒的当期卡池抽数是否是 119（即即将计算的这一抽是第 120 抽） 并且没有触发过 120 井，如果是的话，则直接视为抽到了当期 up，否则继续按照上述概率计算公式计算。这个保底是一次性的，一旦通过任意方式（运气好提前抽到，或触发 120 井机制）获得了当期 UP，120 井会立即失效。因为每 240 抽送一个潜能是额外送的（而不是每 240 抽必定出当期六星），所以在第 120 抽之后，概率公式就会恢复为正常的通用公式。 最后我们考虑每 240 抽赠送一个潜能对玩家总抽数的影响。这个潜能赠送机制是独立于抽卡概率的，不会影响第 240 抽的结果，并且如果玩家只想抽一个当前角色的话完全不需要考虑这个机制（因为第 120 抽的时候肯定怎么也得有了），所以它只会在我们计算期望的时候产生影响，在重建抽卡流程的时候完全不用考虑。 至此，我们就得到了完整的抽卡流程，能够精确计算出第 n 抽时出现当期六星的概率了。 一个六星所需要的期望抽数计算 ​ 因为整个抽卡流程是离散的（50% 歪率，120 井，240 赠送），我们无法直接通过纯数学方式计算出抽出 k 个当期 up 角色所需要的准确期望抽数，但是我们可以通过概率密度函数计算在没有井机制的情况下单抽出六星的概率（后续计算新池等歪成功率时需要用到）。 为了直观理解，我们使用期望的等价公式描述这个概率，即“单个出金循环内的期望抽数”等于“第 n 抽还未出金的概率”之和： 由于在单个出金循环中，出金概率会在 66 ≤ n ≤ 79 时每抽提升 5%，因此我们需要分段计算期望抽数。设 r n 为第 n 抽的出金概率（注意这里的出金仅指抽到六星，不代表抽到当期 up），则有 （其实就是上面的 P d r o p ( n ) ，只是为了区分符号重开了一个） 为了方便阅读，我们定义第一阶段（ 1 ≤ n ≤ 65 ）为平稳期，总期望为 Sum 1 ；第二阶段（ 66 ≤ n ≤ 80 ）为上升期，总期望为 Sum 2 ，则出一个六星需要的总期望抽数 E g o l d = Sum 1 + Sum 2 。 其中，平稳期的单抽存活概率 P ( X ≥ n ) 是一个公比为 ( 1 − 0.008 ) = 0.992 的等比数列，因此其期望 Sum 1 为初中数学的等比数列求和公式： 现在再看上升期，由于上升期中每多一抽出金概率就会提升，即 但是你说巧不巧，这个 n 一定是一个 [66，80] 之间的整数，所以我们可以写一个简单的循环求解： 计算得到 Sum 2 ≈ 3.05910 。 因此得到一个六星需要的总期望抽数 向上取整后，抽到一个六星（不考虑歪没歪）需要的期望抽数为 54 抽。 基于以上分析，我们可以定义需要追踪的变量如下： 自上一次出金以后的抽数（可以拆分为“上一个池的水位”和“当前池的水位”） pity ，用于计算当前抽出金的概率 当前池总抽数 bannerPullCount ，用于判定 120 井和 240 赠送 已获得目标角色的数量 targetCopies ，用于追踪目标（需要抽到几个当期六星）完成情况 是否触发了 120 抽保底 oneTwentyUsed ，用于判定当前抽是否保底 额外的，玩家可能还会有当前池子没抽满，希望在后续池子歪出来的想法。但是这个概率 P t a r g e t = 0.5 × 1 N s t d + N e x t r a （ N s t d 是当前池子的常驻六星数量， N e x t r a 是当前池子的陪跑六星数量）实际上是一个静态概率，与抽卡次数无关，所以到时候我们可以直接计算，这里还是专注于考虑当前池子的动态概率计算过程。 基于以上需要追踪的变量，我们可以列出接口结构如下： 要计算这种期望，最简单粗暴但是又很 CS 的方式就是蒙特卡洛模拟，即模拟 k 个玩家，每个玩家进行 i 次抽卡行为，统计每个玩家抽到 j 个当期六星所需的抽数，最后除以 k 得到期望。整个抽卡过程在上文也已经分析过了，这里就直接给出代码，重点关注后面的 DP 实现。 上面说的 MCMC 是一种常用的模拟方式，优点是逻辑简单，并且在模拟次数足够多的情况下，结果会无限接近真实分布。但是缺点也很明显，模拟次数越多的情况下期望计算时间越长，如果把模拟次数提升到 100 万次的话，在我的 Mac mini 上会把浏览器线程硬控将近十秒。并且由于 JS 计算会阻塞页面渲染，所以想给用户展示一个 loading 状态都不好操作。一个折中的解决方案是使用 WebWorker 将计算过程 detach 到后台线程（后面也会给出包含这部分的完整代码），但是终究改变不了 MCMC 是一个计算密集型算法的事实，再怎么优化它在手机上也是跑不快的（更别提我都不知道 iOS Safari 对 WebWorker 能支持到什么情况）。那么有没有一种算法能直接计算概率呢？有的兄弟有的，就是比较掉头发，我们可以祭出 leetcode 禁术之动态规划。 对于动态规划而言，我们需要找出“最小子问题”，然后从最小子问题开始逐步推导出最终问题的解。根据最小子问题的定义不同，一个 DP 问题可以有多种不同的状态定义方式，最终得到不同的状态转移方程。 对于这个抽卡问题而言，每一抽的结果只取决于当前的“卡池状态”（已经抽了多少抽没有出金，抽出了几个当期 up，有没有触发过 120 井），而与之前的抽卡历史无关。即使是完全不同的账号，只要当前的卡池状态一致，下一抽的出率也不会改变。基于这种“将来只取决于现在”的马尔可夫性，我们可以将整个模拟过程简化为对每一抽的概率计算（人话：打表），因此要使用动态规划，我们应当定义清楚什么是所谓的“卡池状态”，以及如何描述它。 就像上一段提到的，在终末地的抽卡模型中，决定下一抽概率和结果的因素有且仅有「当前水位」、「已获得目标数」和「是否触发过 120 井」这三个因素共同决定的，因此这三个状态也会被我们定义成 DP 的三个状态维度：当前水位 dp[p] 、已获得目标角色数 dp[k] 、120 抽保底井 dp[s] 。至于“当前池总抽数”这个变量，它其实是动态规划推进的时间轴 t（或者说外层循环的计数器），不需要作为内部状态存储。 综上所述，我们便可定义动态规划的当前状态 dp[p][k][s] 为：在当前时刻（第 t 抽），账号处于水位为 p、已获得 k 个目标角色、120 井状态为 s 的概率。所有的状态会形成一个名为 dp 的三维空间，其中 p 的取值范围是 0 到 79，k 的取值范围是 0 到（用户设定），s 的取值范围是 1（已使用） 或 0（未使用）。 随着总抽数 t 的增加，我们会不断更新 dp 数组中各个状态的概率值。因为每一抽的概率是上一抽的状态转移过来的，所以我们的最小子问题就是：假设第 t 抽的概率分布为 D P o l d ，推导第 t + 1 抽的概率分布 D P n e w 。 假设当前时刻 t 的状态为 ( p , k , s ) ，则第 t + 1 抽的概率 P c u r r = D P o l d [ p ] [ k ] [ s ] 。我们将第 t + 1 抽（后文称为“当前抽”）的基础出货率记为 R ( p + 1 ) ，则有 而针对 D P n e w 的转移逻辑可以由以下三个顺序执行的步骤计算：120 井判定、常规概率判定，以及 240 抽的赠送判定。 这里有人可能会疑惑为什么还要单独考虑 240 抽的判定（因为它是独立于抽卡概率的，既不重置水位、也不挤占保底）。这个赠送潜能会影响到已获得的目标数 k ，而 k 是在我们 DP 的内部状态中追踪的。当抽数达到 240 的倍数时，玩家手上的角色数量 k 都会强制 +1，相当于对概率分布进行了一次整体右移。这样一来，毕业率曲线在 t = 240 , 480 , ⋯ 处会出现一次垂直跳变。虽然在 MCMC 当中我们不需要关心“已经抽到了几个”，然而在 DP 中，我们是直接计算在 t 抽时正好毕业的概率，因此如果不追踪这一判定的话，我们得到的期望 E 就会明显偏大。

下面我们来推导 D P n e w 的转移逻辑。 正如 流程分析 中提到的，120 井的判定是能够直接改变结果的、优先级最高的规则。其在 DP 中的逻辑为：如果当前总抽数 t + 1 刚好等于 120，并且井未使用（ s = 0 ），则无论 R ( p + 1 ) 的值是多少，都直接视为抽到了当期 up。 其转移目标为：水位 p 转移为 0，目标数转移到 k+1，井标记转移到 1；即状态转移到 ( 0 , k + 1 , 1 ) 。方程可以写作： 并且此时跳过“常规概率判定”和“240 抽赠送判定”（因为不可能到达），直接进入下一抽。 如果没有触发 120 井的话，则会按照上面的基础出货率公式计算 R ( p + 1 ) 的值，并根据其值决定是否抽到当期 up。这里分成三种情况，出金没歪、出金歪了，以及没出金。这里继续分情况讨论。 这一步的概率为 R ( p + 1 ) × 0.5 ，转移目标为：水位 p 转移到 0 ，目标数转移到 k+1 ，井标记转移到 1 ；即状态转移到 ( 0 , k + 1 , 1 ) 。方程可以写作： 出金的两种状态转移概率都是相等的，只是转移目标不同。歪了的概率也为 R ( p + 1 ) × 0.5 ，转移目标为：水位 p 转移到 0 ，目标数转移到 k ，井标记不变；即状态转移到 ( 0 , k , s ) 。方程可以写作： 因为出金和没出金是互斥事件，因此没出金的概率为 1 − R ( p + 1 ) ，转移目标为：水位 p 转移到 p+1 ，目标数不变，井状态不变。方程可以写作： 以上的 3+1 中情况已经涵盖了抽卡过程中所有可能的状态，下面我们可以来考虑每满 240 抽时的潜能 +1 对玩家目标抽的影响。 这一步骤是一个特殊的判定，是在每一抽的常规概率判定后触发的。如果当前的抽数 t + 1 mod 240 = 0 ，则系统会赠送一个六星潜能（相当于目标个数+1），相当于将概率分布中的 k 整体向右平移一位。其转移目标为：对于所有的 p , s ，将 D P n e w [ p ] [ k ] [ s ] 的概率移动到 D P n e w [ p ] [ k + 1 ] [ s ] 。 需要注意的是这里为了防止重复移动（比如 k = 1 移到 2，然后又被重复处理移到 3），最好倒序遍历 k。 至此我们已经考虑到了抽卡过程中的所有情况，接下来我们需要计算整体期望。我们定义 P f i n i s h ( t ) 为“在第 t 抽刚刚达成目标的概率之和”。则可以通过对“在每一步转移时总共有多少概率流向了 k = t a r g e t 的状态”求和来获取最终期望，即 “等歪”状态下的成功率计算 ​ 如果玩家设定的 maxInvestCurrentBanner 用完了，但是 dp 表中依然存在 k < target 的残留概率，那么这部分概率就代表了当前卡池中没抽满的情况。我们称这部分状态为“等歪”状态，因为玩家已经用完了当前池子的所有抽数，只能等待歪出来。在 MCMC 模拟中，我们可以直接继续暴力模拟直到 k >= target 为止；而在 DP 中，由于我们的所有结果都以概率分布的形式存储，所以需要计算这部分概率的期望抽数。 很巧的是我们的状态空间中正好有当前已获得目标角色 k 这个维度，因此假设还需要 N n e e d = Target − k 个目标角色才能达成目标。而因为未来的卡池中针对这个角色没有井也没有送，因此这个概率模型就变成了一个最简单的伯努利试验模型，歪出一个目标角色的期望抽数也是固定的。 我们先不考虑残留概率 P r e s i d u e 。假设常驻六星有 N s t d 个，陪跑六星有 N e x t r a 个。前面我们通过PMF计算得到了获得一个六星需要的期望抽数 E g o l d = 53.89927 抽，并且抽到一个六星是目标角色的概率是 P h i t = 0.5 × 1 N s t d + E e x t r a 。 因此歪出一个目标角色的期望抽数 E w a i t 为： 如果常驻池有 6 个六星，陪跑的有 2 个六星，那么从 0 开始歪出一个目标六星的期望会高达 863 抽（向上取整，因为没有半抽的说法）。 其实算到这一步已经可以不用继续往下算了，但凡一个心智正常的玩家都不会为了抽出一个六星在卡池里豪掷快九百抽。退一步说要是真的有人这么做了，都有这财力了就直接在 up 池抽出来不好吗？ 但是既然纯交（闲）流（的），我们还可以继续最后一步，考虑加入残留抽数 p 的情况。对于每一个未毕业的情况 ( p , k , s ) ，其残留概率为 P r e s i d u e = D P [ p ] [ k ] [ s ] ，其中的 p 就是残留没有出金的抽数，而这个水位的价值 V ( p ) 会让玩家获得下一个六星的成本低于 E g o l d 。因此我们设 E n e x t ( p ) 为在残留水位 p 下获得下一个六星的期望抽数，则有 因此，我们现在可以根据当前状态 ( p , k , s ) 计算出其期望抽数 E s t a t e ( p , k , s ) ： 至此我们就完成了整个抽卡流程的动态规划建模，接下来我们就可以使用这个模型来计算抽卡期望了。 写到这里加上代码也有六千字出头了，本来想用一篇文章的篇幅再介绍一下怎么用 DP 计算最优充值方案的（简单来说就是在有首充的状态下是有限背包，没有首充状态下是无限背包），但是加进来的话恐怕篇幅太长了（vitepress dev 已经开始 OOM 了），就之后再说好了。 附录丨动态规划的数学本质：能不能直接用 PMF 计算期望？ ​ 这时候可能有人要问了，阿婆阿婆，我是个概率论高手但是不会编程，我看你上面推导也写了不少公式‍ （主要是 Gemini 助攻太猛了） ‍，你能不能试试看直接用 PMF 推导出期望？ 答案是能也不能，假如像明日方舟本体那样（没有大小保底和井，每个金有 50% 概率歪）干抽的话，那么设 T t o t a l 为为获得 UP 所需抽数，“求解恰好在第 t 抽毕业的概率 P ( T t o t a l = t ) ”这个问题就变成了一个简单的卷积： 其中 f s i x ( t ) = P ( X ≥ t ) × r t 代表第 t 抽出六星的概率， r t 为第 t 抽的出金概率（注意这里的出金仅指抽到六星，不代表抽到当期 up）。但是由于终末地的 120 保底破坏了这个单纯卷积的独立性，240 抽送一个潜能又破坏了破坏了无记忆性，因此我们就无法使用简单的卷积公式计算期望了。 因此实际上，我们在 DP 求解的过程中构建了一个包含 80(水位) × (Target+1) (目标数) × 2(井状态) 个维度的状态空间的状态机，并且在这个（算上时间轴 t 的话就是）四维空间中推演出了这个状态机随抽数 t 的所有演变过程，通过这种方法得到了 PMF。 具体地说，我们定义状态空间 S 和转移矩阵 M ，并且在任意时刻 t 时，玩家的状态可以由一个状态向量 v t 描述。为了理解简单，我们暂时忽略当前水位 p 的影响，只关注核心目标的状态。 定义 S i , j 表示“已获得 i 个目标角色，且 120 井状态为 j （ j = 0 未激活， j = 1 已激活）”的状态。那么对于单个目标即 k = 1 的情况下，状态向量 v t 包含以下集中状态： S 0 , 0 ：未获得任何目标角色，且 120 井未激活 S 0 , 1 ：未获得任何目标角色，且 120 井已激活 S 1 , 1 ：已获得 1 个目标角色，且 120 井已激活（吸收态） 而每一抽 t → t + 1 的过程，本质上是一个线性变换： 其中 M t 是第 t 抽的转移矩阵，其元素 M i , j 表示从状态 S i , j 转移到状态 S i ′ , j ′ 的概率，而我们的转移函数实际上是在构建这个稀疏矩阵并执行向量乘法。转移分为以下几种情况。 常规转移（非 120 井和 240 抽）： 没出货， P ( S 0 , 0 → S 0 , 0 ) = 1 − r ( p ) ； 歪了， P ( S 0 , 0 → S 0 , 0 ) = r ( p ) × 0.5 ； 以及 P ( S 0 , 0 → S 1 , 1 ) = r ( p ) × 0.5 没出货， P ( S 0 , 0 → S 0 , 0 ) = 1 − r ( p ) ； 歪了， P ( S 0 , 0 → S 0 , 0 ) = r ( p ) × 0.5 ； 以及 P ( S 0 , 0 → S 1 , 1 ) = r ( p ) × 0.5 120 井转移（ t = 119 → 120 ，强制毕业）： P ( S 0 , 0 → S 1 , 1 ) = 1.0 P ( S 0 , 0 → S 1 , 1 ) = 1.0 240 抽转移：上文分析过是向右平移，它会将所有 S i , … 的状态向右平移一个单位到 S i + 1 , … 当 DP 的所有循环执行完毕时，我们的状态空间 S 当中存储的值就正好是概率质量函数 PMF 在 t 处的值，即 而我们的累计分布函数 CDF 可以通过累加 PMF 得到，即 因此，整个抽卡的期望抽数就是 PMF 的加权和。由于我们的抽卡过程被最大允许投入的抽数限制，因此原本的无穷级数 会在 MaxInvest 处被截断，得到最终的期望公式 而 DP 求解的过程本质上是填满这个状态空间，并从后往前推导出 PMF 的过程。