[POJ1625]Censored!(AC自动机+DP+高精度)

Censored!

Time Limit: 5000MS		Memory Limit: 10000K
Total Submissions: 10824		Accepted: 2966

Description

The alphabet of Freeland consists of exactly N letters. Each sentence of Freeland language (also known as Freish) consists of exactly M letters without word breaks. So, there exist exactly N^M different Freish sentences.

But after recent election of Mr. Grass Jr. as Freeland president
some words offending him were declared unprintable and all sentences
containing at least one of them were forbidden. The sentence S contains a
word W if W is a substring of S i.e. exists such k >= 1 that S[k] =
W[1], S[k+1] = W[2], ...,S[k+len(W)-1] = W[len(W)], where k+len(W)-1
<= M and len(W) denotes length of W. Everyone who uses a forbidden
sentence is to be put to jail for 10 years.

Find out how many different sentences can be used now by freelanders without risk to be put to jail for using it.

Input

The
first line of the input file contains three integer numbers: N -- the
number of letters in Freish alphabet, M -- the length of all Freish
sentences and P -- the number of forbidden words (1 <= N <= 50, 1
<= M <= 50, 0 <= P <= 10).

The second line contains exactly N different characters -- the
letters of the Freish alphabet (all with ASCII code greater than 32).

The following P lines contain forbidden words, each not longer than
min(M, 10) characters, all containing only letters of Freish alphabet.

Output

Output the only integer number -- the number of different sentences freelanders can safely use.

Sample Input

2 3 1
ab
bb

Sample Output

5

Source

Northeastern Europe 2001, Northern Subregion

同HDU2222，只是需要高精度

感觉可能有个问题，就是这题题目没有规定屏蔽词包含了所有字母，但是好像对解题没有什么影响。

代码用时：1.5h 高精度输出写炸了

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define rep(i,l,r) for (int i=l; i<=r; i++)
 using namespace std;

 const int N=;
 int n,m,p,nd,fail[N],b[N],q[N],c[N][N];
 char s[N],S[N];

 int get(char ch){ rep(i,,n) if (S[i]==ch) return i; return -; }

 void ins(char S[]){
     int x=,len=strlen(s+);
     rep(i,,len){
         int k=get(s[i]);
         if (!c[x][k]) ++nd,c[x][k]=nd;
         x=c[x][k];
     }
     b[x]=;
 }

 void getfail(){
     int st=,ed=;
     rep(i,,n) if (c[][i]) q[++ed]=c[][i];
     while (st!=ed){
         int x=q[++st];
         rep(i,,n)
             if (!c[x][i]) c[x][i]=c[fail[x]][i];
                 else q[++ed]=c[x][i],fail[c[x][i]]=c[fail[x]][i];
         b[x]|=b[fail[x]];
     }
 }

 struct D{
     int v[N],len;
     D(int x=){
         memset(v,,sizeof(v));
         for (len=; x>; x/=) v[len++]=x%;
         len--;
     }
     D operator +(const D &a){
         D ans;
         ans.len=max(len,a.len);
         for (int i=; i<=ans.len; i++){
             ans.v[i]+=v[i]+a.v[i]; ans.v[i+]+=ans.v[i]/; ans.v[i]%=;
         }
         while (ans.v[ans.len+]) ans.len++;
         return ans;
     }
     void print(){
         if (len==-) { printf("%d\n",); return; }
         for (int i=len; ~i; i--) printf("%d",v[i]);
         printf("\n");
     }
 }dp[][N];

 int check(int k,int j){ int ans=; rep(i,,n) if (c[k][i]==j) ans++; return ans; }

 int main(){
     freopen("poj1625.in","r",stdin);
     freopen("poj1625.out","w",stdout);
     scanf("%d%d%d",&n,&m,&p); scanf("%s",S+);
     rep(i,,p) scanf("%s",s+),ins(s);
     getfail(); dp[][]=D();
     rep(i,,m) rep(j,,nd) if (!b[j])
         rep(k,,nd) if (!b[k])
             for (int ti=check(k,j); ti--; ) dp[i][j]=dp[i][j]+dp[i-][k];
     D ans; rep(i,,nd) ans=ans+dp[m][i];
     ans.print();
     return ;
 }