空罐Cans (ac自动机 dp)

原创 2013年12月03日 12:26:57

忘记了 这一句  val[u] |= val[f[u]];///!!!

参考:http://blog.csdn.net/no__stop/article/details/8943271

另有:

hdu 3341 Lost's revenge关于压缩的技巧

lightoj 1427 Substring Frequency (II) (ac自动机)在fail树上的拓扑

。。。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include <ctime>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <string>
#include <set>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;

#define FE(i, a, b) for(int i = (a); i <= (b); ++i)
#define FD(i, b, a) for(int i = (b); i >= (a); --i)
#define REP(i, N) for(int i = 0; i < (N); ++i)
#define CLR(a, v) memset(a, v, sizeof(a))
#define PB push_back
#define MP make_pair
typedef long long LL;
const int INF = 0x3f3f3f3f;
const int MOD = 10007;

const int maxn = 1555;
const int SIG = 4;

int dp[2][111][maxn];
int ans1, ans2;

int n, p;

struct AC{
    int ch[maxn][SIG];
    int f[maxn];
    int dep[maxn];
    int val[maxn];
    int sz;
    int newnode(int deep)///!!!
    {
        CLR(ch[sz], 0);
        dep[sz] = deep;
        val[sz] = 0;
        return sz++;
    }
    void init()
    {
        sz = 0;
        newnode(0);
    }
    int idx(char c) { return c-'a'; }
    void insert(char *s)
    {
        int u=0,n=strlen(s);
        int d = 0;
        for(int i=0;i<n;i++)
        {
            d++;
            int c=idx(s[i]);
            if(!ch[u][c]) ch[u][c] = newnode(d);
            u = ch[u][c];
        }
        val[u] = 1;
    }
    void getFail()
    {
        queue<int> q;
        f[0] = 0;
        for(int c = 0; c < SIG; c++)
        {
            int u = ch[0][c];
            if(u) { f[u] = 0; q.push(u); }
        }

        while(!q.empty())
        {
            int r = q.front(); q.pop();
            for(int c = 0; c < SIG; c++)
            {
                int  u = ch[r][c];
                if(!u) { ch[r][c] = ch[f[r]][c]; continue; }
                q.push(u);

                int v = f[r];
                while(v && !ch[v][c]) v = f[v];
                f[u] = ch[v][c];
                val[u] |= val[f[u]];///!!!
            }
        }
    }
    void solve(char *s)
    {
        int len = strlen(s);
        CLR(dp, 0);
        int fla = 0;
        int uu = 0;
        for (int i = 0; i < len; i++)
        {
            int c = idx(s[i]);
            uu = ch[uu][c];
            if (val[uu]) fla = 1;
        }
        if (fla) { ans1 = 0; ans2 = 1; return; }
        dp[0][len][uu] = 1;

        int now, next;
        now = 0;

        for (int i = 0; i <= p; i++)
        {
            next = now ^ 1;
            for (int r = 0; r < sz; r++)
            {
                ans1 += dp[now][0][r];
                ans1 %= MOD;///!!!
                for (int j = 1; j <= 100; j++)
                {
                    int tmp = dp[now][j][r];
                    if (!tmp) continue;
                    if (val[r]) { ans2 += tmp; ans2 %= MOD; continue; }///!!!
                    for (int g = 0; g < 4; g++)
                    {
                        int u = ch[r][g];
                        dp[next][j + 1][u] += tmp;
                        dp[next][j + 1][u] %= MOD;
                    }
                    if (dep[r] <= j - 1) dp[next][j - 1][r] += tmp, dp[next][j - 1][r] %= MOD;
                    else dp[next][j - 1][f[r]] += tmp, dp[next][j - 1][f[r]] %= MOD;
                }
            }
            CLR(dp[now], 0);
            now ^= 1;
        }
    }

}ac;

char s[maxn];
char ca[maxn];
int main()
{
    while (~scanf("%s", s))
   {
        scanf("%d%d", &p, &n);
        ac.init();
        for (int i = 1; i <= n; i++)
        {
            scanf("%s", ca);
            ac.insert(ca);
        }
        ac.getFail();
        ans1 = ans2 = 0;
        ac.solve(s);
        printf("%d %d\n", ans1, ans2);
   }
    return 0;
}


相关文章推荐

空罐 AC自动机 + DP 好题

链接:http://zerojudge.tw/ShowProblem?problemid=b179 做这道题首先要对ac自动机理解透彻!! 蛮好的一道题目,每一天字符串除了会在尾部扩展一个,开头的...

hdu 5955 Guessing the Dice Roll 概率DP,ac自动机,高斯消元 ★★

掷骰子,n个人,每人预测一个长度为L的序列,直至筛子序列的最后L个数与某个人预测的一致为止游戏结束(每个人预测的序列不一样,且长度均为L)。...

POJ2778 DNA Sequence AC自动机,DP,矩阵加速

题目大意,找长度为n的不包含某些zi fu chang d

HDU 4511 小明系列故事——女友的考验 (AC自动机+DP)

题意:   终于放寒假了,小明要和女朋友一起去看电影。这天,女朋友想给小明一个考验,在小明正准备出发的时候,女朋友告诉他,她在电影院等他,小明过来的路线必须满足给定的规则:   1、假设小明在的位置...
  • Ezereal
  • Ezereal
  • 2016年08月08日 09:19
  • 159

poj 2778 (AC自动机+dp+矩阵快速幂)

断断续续看了一周AC自动机和Trie图,总算明白一些了,A的不易啊 /* * ============================================================...
  • dlut_ju
  • dlut_ju
  • 2013年08月25日 01:27
  • 702

AC自动机+DP小结 (一)

好久没有更新博客了,最近真是懒到家了,南京赛前重点复习了下AC自动机+DP方面的题,写下来总结一下。 HDU 2457  DNA repair http://acm.hdu.edu.cn/showpr...

hdu2457 DNA repair(AC自动机+dp)

题目链接:点击打开链接 题意描述:给定一些带有病毒的DNA序列,然后给出一串DNA,问最少修改多少个字符使DNA中不含有带病的序列? 解题思路:AC自动机+dp 分析: 首先这是一个多字符串问...

hdu4758(ac自动机,状态压缩dp)

Walk Through Squares Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65535/65535 K (Java/Oth...

HDU - 2457 DNA repair AC自动机+dp

题意:给出n个模式串,给出一个主串,要求在主创上改动尽可能少的字符,使任一模式串不在主串中出现。 AC自动机+dp 按照一般的 AC自动机+dp转移 的模式(套路)来做,与1月底做的一个经典例题题差不...

NOJ1222-English Game(AC自动机+dp)

问题描述This English game is a simple English words connection game. The rules are as follows: there are...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:空罐Cans (ac自动机 dp)
举报原因:
原因补充:

(最多只允许输入30个字)