zoukankan      html  css  js  c++  java
  • hdu 2825 Wireless Password AC自动机+状态DP

    时间卡得紧,写成递推可以做一些优化

    //#pragma comment(linker, "/STACK:1024000000,1024000000")
    #include<cstdio>
    #include<cstring>
    #include<cstdlib>
    #include<algorithm>
    #include<iostream>
    #include<sstream>
    #include<cmath>
    #include<climits>
    #include<string>
    #include<map>
    #include<queue>
    #include<vector>
    #include<stack>
    #include<set>
    using namespace std;
    typedef long long ll;
    typedef unsigned long long ull;
    typedef pair<int,int> pii;
    #define pb(a) push_back(a)
    #define INF 0x1f1f1f1f
    #define lson idx<<1,l,mid
    #define rson idx<<1|1,mid+1,r
    #define PI  3.1415926535898
    template<class T> T min(const T& a,const T& b,const T& c) {
        return min(min(a,b),min(a,c));
    }
    template<class T> T max(const T& a,const T& b,const T& c) {
        return max(max(a,b),max(a,c));
    }
    void debug() {
    #ifdef ONLINE_JUDGE
    #else
    
        freopen("d:\in.txt","r",stdin);
       // freopen("d:\out1.txt","w",stdout);
    #endif
    }
    int getch() {
        int ch;
        while((ch=getchar())!=EOF) {
            if(ch!=' '&&ch!='
    ')return ch;
        }
        return EOF;
    }
    
    const int MAX_NODE=100;
    const int SIGMA_SIZE=26;
    const int mod = 20090717;
    
    int ch[MAX_NODE][SIGMA_SIZE];
    int val[MAX_NODE];
    int fail[MAX_NODE];
    int sz;
    
    void init()
    {
        memset(ch[0],0,sizeof(ch[0]));
        sz=1;
        val[0]=0;
    }
    int idx(char c){return c-'a';}
    void insert(const char *s,int k)
    {
        int u=0;
        for(int i=0;s[i]!='';i++)
        {
            int v=idx(s[i]);
            if(!ch[u][v])
            {
                memset(ch[sz],0,sizeof(ch[sz]));
                val[sz]=0;
                ch[u][v]=sz++;
            }
            u=ch[u][v];
    
        }
        val[u]=1<<k;
    }
    
    void construct ()
    {
        fail[0]=0;
        queue<int> q;
        for(int c=0;c<SIGMA_SIZE;c++)
        {
            int u=ch[0][c];
            if(u){fail[u]=0;q.push(u);}
        }
    
        while(!q.empty())
        {
            int r=q.front();q.pop();
    
            for(int c=0;c<SIGMA_SIZE;c++)
            {
                int u=ch[r][c];
                if(!u)
                {
                    ch[r][c]=ch[fail[r]][c];
                    continue;
                }
                q.push(u);
    
                int v=fail[r];
                while(v&&!ch[v][c])v=fail[v];
                fail[u]=ch[v][c];
                val[u]|=val[fail[u]];
            }
        }
    }
    
    int check[1<<10];
    void prework()
    {
        for(int st=0;st<(1<<10);st++)
        {
            check[st]=0;
            int nst=st;
            while(nst)
            {
                check[st]+=nst%2;
                nst/=2;
            }
        }
    }
    int dp[MAX_NODE][26][1<<10];
    
    int DP(int n,int m,int k)
    {
        memset(dp,0,sizeof(dp));
        dp[0][0][0]=1;
        for(int i=0,k=0;i<n;i++,k^=1)
        {
            for(int u = 0;u<sz;u++)
            {
                for(int st=0;st<(1<<m);st++)if(dp[u][i][st]>0)
                {
                    for(int c=0;c<SIGMA_SIZE;c++)
                    {
                        dp[ch[u][c]][i+1][st|val[ch[u][c]]]=(dp[ch[u][c]][i+1][st|val[ch[u][c]]]+dp[u][i][st])%mod;
                    }
                }
            }
        }
        int res=0;
        for(int u=0;u<sz;u++)
        {
            for(int st=(1<<k)-1;st<=(1<<m);st++)
                if(check[st]>=k)
                    res=(res+dp[u][n][st])%mod;
        }
        return res;
    }
    int main()
    {
        prework();
        int n,m,k;
        while(scanf("%d%d%d",&n,&m,&k)!=EOF&&(n+m+k))
        {
            init();
            for(int i=1;i<=m;i++)
            {
                char buf[15];
                scanf("%s",buf);
                insert(buf,i-1);
            }
            construct();
            int res=DP(n,m,k);
            printf("%d
    ",res);
        }
        return 0;
    }
    View Code
  • 相关阅读:
    DP——背包问题(三)
    堆——练习题
    DP——背包问题(二)
    二叉树的后序遍历(暴力版) 小白菜oj 1034
    树状数组2模板 Luogu 3368
    树状数组1模板 Luogu 3374
    DP——最长上升子序列(n^2与n log n)
    线段树(区间修改)模板题 Luogu 2357 守墓人
    c语言学习摘录
    python 学习摘录
  • 原文地址:https://www.cnblogs.com/BMan/p/3415152.html
Copyright © 2011-2022 走看看