poj 1625 dp+大数+ac自动机

题意：

在一个国家当中，给你一个词汇书，你可以用里面所有的字母拼凑出一个单词来，但是在这个国家当中存在一些禁用的单词，那么现在问你在不违反规则的情况下，你能尽可能的多拼凑出相应的单词来。

那么也就是说拼凑的单词当中不能够出现禁用的单词的字段。

后面看到网上教程当中的dp状态转移，傻眼了。。这都可以。

dp[i+1][cnt]=dp[i+1][cnt]+dp[i][j]，我所理解的就是在i+1这个位置上面我放上cnt这个字符，将之前可行的数目累加到长度为i+1的基础上，其实我认为这个与其说是dp，倒不如说是一种递推。

而其中的ac自动机就是用来判断其中新拼凑的单词当中是否出现过类似的字串，出现了就不能够累加到相应的dp值上。

大数当然是因为方案的种类数多了之后难以用int整形存下。

代码：

//大数+dp+ac自动机
#include<iostream>
#include<stdlib.h>
#include<cstring>
#include<cstdio>
#include<queue>

const int Type=55;
const int Max=800;
const int base=100000000;
const int bit=30;

using namespace std;

struct node
{
       node *next[Type];
       node *fail;
       int id,sum;
}tree[Max];

struct BigNum
{
       int num[bit];
       int len;
       BigNum()
       {
           len=1;
           memset(num,0,sizeof(num));
       }
}dp[Type][Max];

int MM(int a,int b)
{
    return a>b?a:b;
}

BigNum operator+(BigNum a,BigNum b)
{
       BigNum c;
       c.len=MM(a.len,b.len);
       int k=0,t;
       for(int i=1;i<=c.len;i++)
       {
           t=a.num[i]+b.num[i]+k;
           k=t/base;
           c.num[i]=t%base;
       }
       while(k)
       {
           c.num[++c.len]=k%base;
           k/=base;
       }
       return c;
}

int n,m,p,tot;
char s[Type],str[Type];
int dx[5*Type];

node *getroot()
{
      tot=0;
      tot++;
      tree[tot].fail=NULL;
      memset(tree[tot].next,NULL,sizeof(tree[tot].next));
      tree[tot].id=tot;
      tree[tot].sum=0;
      return &tree[tot];
}
node *root;

node *getnode()
{
      tot++;
      tree[tot].fail=root;
      memset(tree[tot].next,NULL,sizeof(tree[tot].next));
      tree[tot].id=tot;
      tree[tot].sum=0;
      return &tree[tot];
}

void insert()
{
     node *tmp=root;
     int i=0;
     while(str[i])
     {
           int index=dx[str[i]+128];
           if(tmp->next[index]==NULL)tmp->next[index]=getnode();
           tmp=tmp->next[index];
           i++;
     }
     tmp->sum++;
}

void build_ac()
{
     queue<node*> q;
     while(!q.empty())q.pop();
     q.push(root);
     while(!q.empty())
     {
           node *tmp=q.front();
           q.pop();
           node *ff;
           for(int i=0;i<n;i++)
               if(tmp->next[i]!=NULL)
               {
                  q.push(tmp->next[i]);
                  for(ff=tmp->fail;ff!=NULL;ff=ff->fail)
                  {
                      if(ff->next[i]!=NULL)
                      {
                          tmp->next[i]->fail=ff->next[i];
                          if(ff->next[i]->sum==1)tmp->next[i]->sum=1;
                          break;
                      }
                  }
               }
     }
}

void solution()
{
     dp[0][1].num[1]=1;
     for(int i=0;i<m;i++)
     {
         for(int j=1;j<=tot;j++)
         {
             node *tmp=&tree[j];
             for(int k=0;k<n;k++)
             {
                 if(tmp->next[k]!=NULL)
                 {
                    int index=tmp->next[k]->id;
                    if(tmp->next[k]->sum==0)dp[i+1][index]=dp[i+1][index]+dp[i][j];
                 }
                 else
                 {
                     node *now=tmp->fail;
                     for(;now!=NULL;now=now->fail)
                     {
                         if(now->next[k]!=NULL)
                         {
                             now=now->next[k];
                             break;
                         }
                     }
                     if(now==NULL)now=root;
                     int index=now->id;
                     if(now->sum==0)dp[i+1][index]=dp[i+1][index]+dp[i][j];
                 }
             }
         }
     }
}
void init()
{
     memset(dx,-1,sizeof(dx));
     tot=0;
     gets(s);
     for(int i=0;i<n;i++)
     {
         if(dx[s[i]+128]==-1)
         dx[s[i]+128]=i;
     }
     root=getroot();
     while(p--)
     {
         scanf("%s",str);
         insert();
     }
}
void Print()
{
     BigNum answer;
     for(int i=1;i<=tot;i++)
         answer=answer+dp[m][i];
     int ll=answer.len;
     printf("%d",answer.num[ll]);
     for(int j=ll-1;j>=1;j--)
     {
         if(answer.num[j]<base/10)printf("0");
         if(answer.num[j]<base/100)printf("0");
         if(answer.num[j]<base/1000)printf("0");
         if(answer.num[j]<base/10000)printf("0");
         if(answer.num[j]<base/100000)printf("0");
         if(answer.num[j]<base/1000000)printf("0");
         if(answer.num[j]<base/10000000)printf("0");
         printf("%d",answer.num[j]);
     }
     printf("\n");
}
int main()
{
    scanf("%d %d %d\n",&n,&m,&p);
    init();
    build_ac();
    solution();
    Print();
    return 0;
}