【数位DP】题集

1.【HDOJ2089

题意:求区间内不出现4和62的数的个数

解法:模板题

2.【HDOJ3555

题意:求区间内不出现49的数的个数

解法:模板题

3.【HDOJ5179

题意:对于一个十进制数,把每一位拆出来变成一个数组,要求数组右边的数都小于等于左边的数,且左数能整除右数。求区间内满足条件的个数

解法:dp[i = 20][j = 10]: j为前一位上的数
按要求转移即可,要考虑前导零,注意防止%0的情况出现

4.【HDOJ3652

题意:求区间内含13或者是13倍数的数的个数

解法:dp[i = 20][j = 13][k = 3]: j为对13的余数,k为状态(0: 前一位不是1 | 1: 前一位是1 | 2: 前面出现过13)
最后当j为0,k为2时加入计数,计算不含13并且不是13倍数的数
输出时为N + 1 - solve(N) (0也在统计内)

5.【HDOJ3709

题意:定义平衡数,以某一数位作为原点,左边数位上数字 × 距离之和与右边数位上的相等。求区间内平衡数的个数

解法:dp[i = 20][j = 5000][k = 20]: j为左边数位权值之和,k为枚举的原点
需要枚举原点,然后统计左边权值的减去右边的权值最后是否为零,由于 0 对于每个枚举位置都会被统计,最后要再减去枚举次数 - 1个重复的0

6.【HDOJ6148

题意:对于一个十进制数,要求从左到右的数位没有出现先递增接着递减的情况,求区间内满足条件的个数

解法:dp[i = 105][j = 10][k = 2]: j为前一个数,k表示前面是否已出现递增
按要求转移即可,要考虑前导零,这题的0是不合法的,额外减1就行了,注意数据范围10^100

7.【HDOJ5898

题意:对于一个十进制数,要求数位中奇数段长度都为偶数,偶数段长度都为奇数,求区间内满足条件的个数

解法:dp[i = 20][j = 10][k = 20]: j为前一个数,k为前一个数所属段的长度
要考虑前导零,转移时,当前数位若与前一位奇偶相反,表示前一段已结束,判断长度是否满足要求剪枝,并重置k为1继续向下枚举,最后统计时还要再判断一次j与k,奇偶相反则加入统计

8.【HDOJ4734

题意:定义:F(x) = An * 2^n-1 + An-1 * 2^n-2 + ... + A2 * 2 + A1 * 1,Ai是十进制数位,求区间内满足F(i)<=F(a)的i的个数

解法:dp[i = 20][j = 20000]: j表示枚举到当前位时,在满足条件情况下后面的数位还能有多少权值(初始值为F(a))
最后当j >= 0加入计数,中间状态若j已经小于0则可剪枝

9.【HDOJ4352

题意:求区间内数位从左到右满足LIS = m的数的个数

解法:dp[i = 20][j = 1 << 10][k = 11]: j状态压缩表示当前LIS状态,k表示所求LIS长度(不加k会导致每组数组要重新初始化,超时)
要考虑前导零,在数位上跑nlogn的LIS

int updateState(int state, int digit) {
    // i=digit严格递增,状态压缩无法计算不减递增
    for (int i = digit; i <= 9; i++) {
        if (state & (1 << i)) {
            state ^= (1 << i);
            break;
        }
    }
    return state | (1 << digit);
}
int getLIS(int state) {
    int cnt = 0;
    while (state) {
        if (state & 1) {
            cnt++;
        }
        state >>= 1;
    }
    return cnt;
}

10.【HDOJ4507

题意:求区间内和7无关的数字的平方和,结果对1e9+7取模

解法:dp[i = 20][j = 7][k = 7]: j表示整个数%7的余数,k表示数位和%7的余数
要求的是所有满足条件的数的平方和,要对一个数按照平方和公式拆分。
假设满足条件的数有234, 245, 266,那么

234^2 + 245^2 + 266^2 =
(200 + 34)^2 + (200 + 45)^2 + (200 + 66)^2 =
3 * 200^2 + 2 * 200 * (34 + 45 + 66) + (34^2 + 35^2 + 66^2)

因此在枚举到2的时候,表达式里未知的只有
3, (34 + 45 + 66), (34^2 + 35^2 + 66^2)

结构体数组dp中记录三个值,cnt代表满足条件的数的个数,sum表示满足条件的数的求和,sum2表示满足条件的数的平方和

不带取模的转移方程:

ans.cnt += tmp.cnt  
ans.sum += (tmp.sum + tmp.cnt * (i * p[pos]))
ans.sum2 += (tmp.cnt * (i * p[pos])^2 + 2 * (i * p[i]) * tmp.sum + tmp.sum2) 

为了防止出现负数,输出要(solve(r) - solve(l - 1) + mod) % mod

11.【POJ3252

题意:求区间内满足二进制下0的数位 >= 1的数位的数的个数

解法:dp[i = 32][j = 32][k = 32]: j表示0的个数,k表示1的个数
要考虑前导零,按照题意转移即可,注意j + len - pos < k时可剪枝

12.【HDOJ3886

题意:给出一个字符串,只含/, -, ,要求一个数上的数位按相应字符上升、不变或下降(如12345554321符合/-),求区间内满足条件的数的个数,输出时保留最后八位

解法:dp[i = 105][j = 10][k = 105]: j表示前一个数,k表示匹配到字符串哪个位置
要注意的点比较多:
① 输入数据可能有前导零,范围10^100
② 字符串数据的大整数减法(L - 1
③ 输出时保留最后八位,即对1e8取模,同时按照指定格式输出
printf("%08lld ",(solve(r) - solve(l) + mod) % mod);
④ 对每组数据都要初始化dp,因为要匹配的字符串变了
⑤ 枚举时要考虑前导零,只能在无前导零状态时进行匹配
⑥ 当能匹配下一个字符时,尽量先匹配,如果不能再判断是否符合之前的字符

原文地址:https://www.cnblogs.com/stolf/p/9477411.html