代码拍卖会 草稿

SDOI2010 代码拍卖会 草稿

题目大意：

给定这样一个定义：对于一个十进制数，如果这个数的数位从高到低非严格单调递增（即形如$111111.....2222...3333...4444....$的数），那么就称这个数叫“同态数”；同态数数位的个数叫做同态数的长度。比如说$1112233456$就是一个长度为$10$的同态数，而$1112233431$就不是。求一个长度为$N$的同态数中，有多少是$P$的倍数。$N\le 10^{18},P\le 500,$答案对$999911659$取余。

当$N$和$P$的值都很小的时候，可以考虑这样一种暴力DP做法：
设$F(i,j,k)$表示填完从右往左的第$i$位，第$i$位的数位为$j$，对$P$取模得$k$的同态数个数。有这样一个方程：
$$F(i,j,k)=\sum _{j\le p\le 9}F(i-1,p,(k-10^{i}j)\mathrm{m}\mathrm{o}\mathrm{d}P)$$
非严格的时间复杂度为$O(9^{2}NP)$。这个是过不了的：其主要原因在于$N$的值过大！
既然如此，我们就只能考虑放弃数位DP，从另一个角度思考这个问题：把同态数进行拆分。

考虑长度为$N$时，最小的同态数为$1111......111(N\text{个}1)$。任意一个同态数，似乎都可以拆分成若干个长度不同的$11111....111$相加！根据这两个性质，我们定义一个概念：同态基。
定义同态基${\beta }_i=\frac{10^i-1}{9}$表示由连续的$i$个$1$组成的同态数。对于任意一个同态数，其均可以表示成${\sum }_{i=1}^N{r}_i{\beta }_i$的形式。注意，选取的同态基个数不能超过$9$个，否则就不满足同态数的要求了！
这种做法的最大优势是取模方便。根据欧拉定理，同态基${\beta }_i$在模$P$意义下是循环的。这意味着当不同的同态基在模$P$意义下相等，那么他们对最终答案的贡献也是相同的。

综上，我们记$M_i$表示满足$\beta \mathrm{m}\mathrm{o}\mathrm{d}P=i$的同态基$\beta$的集合。由于这些同态基对答案的贡献等价，我们可以直接用组合公式来一并计算它们的贡献。
设$F(i,j,k)$表示在集合$M_0\sim M_i$中选取$j$个同态基，模$P$的值为$k$的方案数。由于我们选取的同态基个数不超过$9$个，并且每个同态基长度不相同，我们可以近似使用多重集的组合数公式来计算。
根据乘法原理和组合数公式，我们有：
$$F(i,j,k)=\sum _{0\le q\le j}\left(\genfrac{}{}{0ex}{}{|M_i|+(j-q)-1}{(j-q)}\right)F(i-1,q,(k-i(j-q))\mathrm{m}\mathrm{o}\mathrm{d}P)$$
其中初态有:
$$F(0,j,0)=\left(\genfrac{}{}{0ex}{}{|M_0|+j-1}{j}\right)$$
在最终答案里至少要选取一个长度为$N$的同态基${\beta }_N$。我们可以先扣除这个决策，然后求出最优决策后将这个决策强行赋值给答案。具体来说，我们最多选取8个同态基使得有：
$$ans=\sum _{0\le j\le 8}F(P-1,j,(P-{\beta }_N)\mathrm{m}\mathrm{o}\mathrm{d}P)$$