zoukankan      html  css  js  c++  java
  • hnu11187

    AC自动机+DP

    #include <cstdio>
    #include <queue>
    #include <cstring>
    using namespace std;
    
    #define D(x) 
    
    const int MAX_D_LEN = 25;
    const int MAX_LEN = 88;
    const int MAX_N = 105;
    const int MAX_CHILD_NUM = 128;
    const int MAX_NODE_NUM = MAX_N * MAX_D_LEN;
    const int INF = 0x3f3f3f3f;
    
    char dna[MAX_LEN];
    
    struct Trie
    {
        int next[MAX_NODE_NUM][MAX_CHILD_NUM];
        int fail[MAX_NODE_NUM];
        bool disease[MAX_NODE_NUM];
        int node_cnt;
        bool vis[MAX_NODE_NUM]; //set it to false
        int root;
        int dp[MAX_LEN][MAX_NODE_NUM];
    
        void init()
        {
            node_cnt = 0;
            root = newnode();
            memset(vis, 0, sizeof(vis));
        }
    
        int newnode()
        {
            for (int i = 0; i < MAX_CHILD_NUM; i++)
                next[node_cnt][i] = -1;
            disease[node_cnt++] = false;
            return node_cnt - 1;
        }
    
        int get_id(char a)
        {
            return (int)a;
        }
    
        void insert(char buf[])
        {
            int now = root;
            for (int i = 0; buf[i]; i++)
            {
                int id = get_id(buf[i]);
                if (next[now][id] == -1)
                    next[now][id] = newnode();
                now = next[now][id];
            }
            disease[now] = true;
        }
    
        void build()
        {
            queue<int>Q;
            fail[root] = root;
            for (int i = 0; i < MAX_CHILD_NUM; i++)
                if (next[root][i] == -1)
                    next[root][i] = root;
                else
                {
                    fail[next[root][i]] = root;
                    Q.push(next[root][i]);
                }
            while (!Q.empty())
            {
                int now = Q.front();
                Q.pop();
                for (int i = 0; i < MAX_CHILD_NUM; i++)
                    if (next[now][i] == -1)
                        next[now][i] = next[fail[now]][i];
                    else
                    {
                        fail[next[now][i]]=next[fail[now]][i];
                        if (disease[fail[next[now][i]]])
                            disease[next[now][i]] = true;
                        Q.push(next[now][i]);
                    }
            }
        }
    
        int work(char* dna)
        {
            int len = strlen(dna);
            for (int i = 0; i < node_cnt; i++)
            {
                if (disease[i])
                    dp[len][i] = INF;
                else
                    dp[len][i] = 0;
            }
            for (int i = len - 1; i >= 0; i--)
            {
                int key = get_id(dna[i]);
                for (int j = 0; j < node_cnt; j++)
                {
                    if (disease[j])
                        continue;
                    dp[i][j] = dp[i + 1][root] + 1;
                    int v = next[j][key];
                    if (disease[v])
                        continue;
                    dp[i][j] = min(dp[i][j], dp[i + 1][v]);
                }
            }
            return dp[0][root];
        }
    
        void debug()
        {
            for(int i = 0;i < node_cnt;i++)
            {
                printf("id = %3d,fail = %3d,end = %3d,chi = [",i,fail[i],disease[i]);
                for(int j = 0;j < MAX_CHILD_NUM;j++)
                    printf("%2d",next[i][j]);
                printf("]
    ");
            }
        }
    }ac;
    
    int n, m;
    char st[MAX_LEN];
    
    int main()
    {
        while (scanf("%d %d", &n, &m), n | m)
        {
            ac.init();
            for (int i = 0; i < n; i++)
            {
                scanf("%s", st);
                getchar();
                ac.insert(st);
            }
            ac.build();
            int ans = 0;
            for (int i = 0; i < m; i++)
            {
                gets(st);
                ans += ac.work(st);
            }
            printf("%d
    ", ans);
        }
        return 0;
    }
    View Code
  • 相关阅读:
    如何编译树莓派内核
    代码导出Reporting Services报表文件
    Bit-Coin收入的一分钱
    如何在树莓派上运行雷神之锤III
    新树莓派入手
    如何通过PowerShell在Visual Studio的Post-build中预热SharePoint站点
    每日一题20201218(389. 找不同)
    每日一题20201217(714. 买卖股票的最佳时机含手续费)
    每日一题20201216(290. 单词规律)
    每日一题20201215(738. 单调递增的数字)
  • 原文地址:https://www.cnblogs.com/rainydays/p/4363967.html
Copyright © 2011-2022 走看看