zoj 3494 AC自动机+数位DP

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/yqdjl6/article/details/83346947

这道题说实话我到现在也不能完整的把这份代码写出来，大部分都是抄kuangbin博客里的那份代码，其实这份代码让我让我一个人看我差不多能看懂，但让我写我感觉我根本写不了，首先我们对模式串建一个AC自动机，然后构造出来一些不可达的点，再写成一个矩阵bcd[i][j]表示在AC自动机上第i个点加上j后能到达的点，j是十进制，从0到9，然后再在bcd[i][j]上跑数位DP，我还回头写了两道数位DP想了想怎么写，这里面的数位DP注意是可以有前导0的，dp[len][i]表示长度为len在自动机上第J个节点的合法数量。

#include<bits/stdc++.h>
using namespace std;
using LL = int64_t;
const int maxnode=1e6+5;
const int sigma_size=4;
char s[maxnode];
const LL mod=1e9+9;
struct Node {
    int son[sigma_size];
    int val,fail;
}ch[maxnode];
int bcd[2010][10];
struct AC {
    int sz=1;
    queue<int>Q;
    void init(int x) {ch[x].fail=ch[x].val=0;memset(ch[x].son,0,sizeof(ch[x].son));}
    int idx(char c) {return c-'0';}

    void insert(char s[],int v) {
        int u=0,n=strlen(s);
        for(int i=0;i<n;i++) {
            int c=idx(s[i]);
            if(!ch[u].son[c]) {
                init(sz);
                ch[u].son[c]=sz++;
            }
            u=ch[u].son[c];
        }
        ch[u].val=v;
    }

    int change(int pre,int num) {
        if(ch[pre].val) return -1;
        int now=pre;
        for(int i=3;i>=0;i--) {
            if(ch[ch[now].son[(num>>i)&1]].val) return -1;
            now=ch[now].son[(num>>i)&1];
        }
        return now;
    }

    void build() {
        for(int i=0;i<sigma_size;i++) if(ch[0].son[i]) Q.push(ch[0].son[i]);
        while(!Q.empty()) {
            int now=Q.front();Q.pop();
            int fail=ch[now].fail;
            if(ch[fail].val) ch[now].val|=ch[fail].val;
            for(int i=0;i<sigma_size;i++) {
                int nxt=ch[now].son[i];
                if(nxt) {
                    ch[nxt].fail=ch[fail].son[i];
                    Q.push(nxt);
                }
                else ch[now].son[i]=ch[fail].son[i];
            }
        }
        for(int i=0;i<sz;i++) {
            for(int j=0;j<10;j++) {
                bcd[i][j]=change(i,j);//bcd[i][j]表示在自动机上第i个节点，加上j以后到达的节点
            }
        }
    }
};

LL dp[205][2005],cnt[205];

LL dfs(int len,int pos,bool flag, bool zero) {
    if(len==-1) return 1;
    if(flag==false&&dp[len][pos]!=-1) return dp[len][pos];
    LL ans=0;
    if(zero) {
        ans=(ans+dfs(len-1,pos,flag&&cnt[len]==0,true))%mod;
    }
    else {
        if(bcd[pos][0]!=-1) ans=(ans+dfs(len-1,bcd[pos][0],flag&&cnt[len]==0,false))%mod;
    }
    int ends=(flag?cnt[len]:9);
    for(int i=1;i<=ends;i++) {
        if(bcd[pos][i]!=-1) {
            ans=(ans+dfs(len-1,bcd[pos][i],flag&&i==ends,false))%mod;
        }
    }
    if(!flag&&!zero) dp[len][pos]=ans;
    return ans;
}

LL solve(char s[]) {
    memset(cnt,0,sizeof(cnt));
    int len=strlen(s);
    for(int i=0;i<strlen(s);i++) cnt[i]=s[len-1-i]-'0';
    return dfs(len-1,0,1,1);
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    int T;cin>>T;
    while(T--) {
        AC ans;ans.init(0);
        int n;cin>>n;
        for(int i=1;i<=n;i++) {
            cin>>s;
            ans.insert(s,1);
        }
        ans.build();
        memset(dp,-1,sizeof(dp));
        cin>>s;
        for(int i=strlen(s)-1;i>=0;i--) {
            if(s[i]>'0') {
                s[i]--;
                break;
            }
            else s[i]='9';
        }
        LL res=-solve(s);
        cin>>s;
        res=(res%mod+solve(s)%mod+mod)%mod;
        cout<<res <<"\n";
    }
    return 0;
}