bzoj3530[Sdoi2014]数数 AC自动机+数位dp
Posted GXZlegend
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj3530[Sdoi2014]数数 AC自动机+数位dp相关的知识,希望对你有一定的参考价值。
题目描述
我们称一个正整数N是幸运数,当且仅当它的十进制表示中不包含数字串集合S中任意一个元素作为其子串。例如当S=(22,333,0233)时,233是幸运数,2333、20233、3223不是幸运数。
给定N和S,计算不大于N的幸运数个数。
输入
输入的第一行包含整数N。
接下来一行一个整数M,表示S中元素的数量。
接下来M行,每行一个数字串,表示S中的一个元素。
输出
输出一行一个整数,表示答案模109+7的值。
样例输入
20
3
2
3
14
样例输出
14
题解
AC自动机+数位dp
同学的某道考试题的加强版。。。
由于给定的串多且复杂,要求不能匹配到这些串,所以可以对这些串建立Trie图。
然后设$f[i][j]$表示从$j$开始走$i$个节点,不触碰到危险节点(即得到的是非幸运数)的方案数。
那么如果$j$是危险节点则方案数为0,否则枚举$j$走出去的第一步,用$f[i-1][next[j][k]]$更新$f[i][j]$。
考虑数位dp的过程:
首先处理出位数不满$n$位的数的个数:枚举位数和最高位(非0),然后方案数即为$f[i][next[1][j]]$。
然后考虑位数满$n$位的数的个数:
考虑从高到底的每一位:从0到当前位-1是满的(即后面的数恰好从0到10^{位数}),因此可以直接从$f$中取出。再考虑当前位不满的情况,此时转化为了子问题,按照同样的方法处理即可。
注意最高位是不能包含前导0的,而确定最高位以后其余的位是可以包含前导0的。
时间复杂度$O(10nL)$
细节贼多。。。代码凑合着看吧。。。
#include <queue> #include <cstdio> #include <cstring> #define N 1510 #define mod 1000000007 using namespace std; queue<int> q; int next[N][10] , tot = 1 , fail[N] , tag[N] , f[N][N]; char s[N] , w[N]; void build() { int x , i; for(i = 0 ; i < 10 ; i ++ ) next[0][i] = 1; q.push(1); while(!q.empty()) { x = q.front() , q.pop() , tag[x] |= tag[fail[x]]; for(i = 0 ; i < 10 ; i ++ ) { if(next[x][i]) fail[next[x][i]] = next[fail[x]][i] , q.push(next[x][i]); else next[x][i] = next[fail[x]][i]; } } } int main() { int n , m , i , j , k , t , flag , ans = 0; scanf("%s%d" , s , &m) , n = strlen(s); for(i = 1 ; i <= m ; i ++ ) { scanf("%s" , w); for(j = 0 , t = 1 ; w[j] ; j ++ ) { if(!next[t][w[j] ^ ‘0‘]) next[t][w[j] ^ ‘0‘] = ++tot; t = next[t][w[j] ^ ‘0‘]; } tag[t] = 1; } build(); for(i = 1 ; i <= tot ; i ++ ) f[1][i] = !tag[i]; for(i = 2 ; i <= n ; i ++ ) for(j = 1 ; j <= tot ; j ++ ) if(!tag[j]) for(k = 0 ; k < 10 ; k ++ ) f[i][j] = (f[i][j] + f[i - 1][next[j][k]]) % mod; for(i = 1 ; i < n ; i ++ ) for(j = 1 ; j < 10 ; j ++ ) ans = (ans + f[i][next[1][j]]) % mod; for(i = 0 , t = flag = 1 ; i < n ; i ++ ) { for(j = flag ; j < (s[i] ^ ‘0‘) ; j ++ ) ans = (ans + f[n - i][next[t][j]]) % mod; t = next[t][s[i] ^ ‘0‘]; if(tag[t]) break; flag = 0; } if(!tag[t]) ans = (ans + 1) % mod; printf("%d\n" , ans); return 0; }
以上是关于bzoj3530[Sdoi2014]数数 AC自动机+数位dp的主要内容,如果未能解决你的问题,请参考以下文章
BZOJ.3530.[SDOI2014]数数(AC自动机 数位DP)