POJ 1625 自动机

来源：互联网发布：linux有ping命令吗编辑：程序博客网时间：2024/05/01 02:23
给出包含n个可见字符的字符集，以下所提字符串均由该字符集中的字符构成。给出p个长度不超过10的字符串，求长为m且不包含上述p个字符串的字符串有多少个。
g++提交
int    mat[108][108] ;int    matn ;int    N  ;map<char ,int> to ;//ACconst  int  maxm = 108 ;const  int  kind = 256 ;int    next[maxm][kind] ;int    fail[maxm] ;int    id[maxm] ;struct AC{       int root , n ;       int newnode(){           id[n] = 0 ;           for(int i = 0 ; i < kind ; i++) next[n][i] = -1 ;           return n++ ;       }       void Init(){            n = 0 ;            root = newnode() ;       }       void Insert(char *s){            int now = root  , k ;            while(*s){                 k = to[*s] ;                 if(next[now][k] == -1)                    next[now][k] = newnode() ;                 now = next[now][k] ;                 s++ ;            }            id[now] = 1 ;       }       void  makeAC(){             queue<int> q ;             fail[root] = root ;             int now , i ;             for(i = 0 ; i < kind ; i++){                 if(next[root][i] == -1)                     next[root][i] = root ;                 else{                     fail[next[root][i]] = root ;                     q.push(next[root][i]) ;                 }             }             while(! q.empty()){                  now = q.front() ; q.pop() ;                  if(id[fail[now]]) id[now] = 1 ;                  for(i = 0  ; i < kind ; i++){                      if(next[now][i] == -1)                          next[now][i] = next[fail[now]][i] ;                      else{                          fail[next[now][i]] = next[fail[now]][i] ;                          q.push(next[now][i]) ;                      }                  }             }       }       void  getMat(){             matn = n ;             memset(mat , 0 , sizeof(mat)) ;             for(int i = 0 ; i < n ; i++){                 for(int j = 0 ; j < N ; j++){                     if(id[next[i][j]] == 0)                          mat[i][next[i][j]]++ ;                 }             }       }}ac ;//EndACconst int maxn = 9999 ; // 每位储存4位const int dig = 4 ;class BigNum{ public  :       int  a[7150] ;       int  len ;       BigNum(){len = 1 ; memset(a , 0 , sizeof(a)) ;}       BigNum(const int) ;       BigNum(const char*) ;       BigNum operator *(const BigNum &) const ;       BigNum operator +(const BigNum &) const ;       friend ostream& operator <<(ostream &out , const BigNum&T) ;};BigNum::BigNum(const int T){      int t = T ;      len = 0 ;  memset(a , 0 , sizeof(a)) ;      while(t > maxn){            a[len++] = t - t/(maxn+1)*(maxn+1) ;            t = t/(maxn+1) ;      }      a[len++] = t ;}BigNum::BigNum(const char* s){      memset(a , 0 , sizeof(a)) ;      int t , i , j , k , id = 0 , l = strlen(s) ;      len = l/dig ;      if(l%dig) len++ ;      for(i = l-1 ; i >= 0 ; i -= dig){          t = 0 ;          k = max(0  , i - dig + 1);          for(j = k ; j <= i ; j++)  t = t*10 + s[j] - '0' ;          a[id++] = t ;      }}ostream& operator <<(ostream &out , const BigNum&T){      printf("%d" , T.a[T.len-1]) ;      for(int i = T.len-2 ; i >= 0 ; i--)  printf("%04d" ,T.a[i]) ;      return out ;}BigNum BigNum::operator*(const BigNum &T) const{       BigNum s ;       int i , j , up  , t ;       for(i = 0 ; i < len ; i++){           up = 0 ;           for(j = 0 ; j < T.len ; j++){                t = a[i]*T.a[j] + s.a[i+j] + up ;                if(t > maxn){                     s.a[i+j] = t - t/(maxn+1)*(maxn+1) ;                     up = t/(maxn+1) ;                }                else{                    up = 0 ;                    s.a[i+j] = t ;                }           }           if(up != 0) s.a[i+j] = up ;       }       s.len = len + T.len ;       while(s.a[s.len-1] == 0 && s.len>1) s.len-- ;       return s ;}BigNum BigNum::operator + (const BigNum &T) const{       BigNum s ;       s.len = max(len , T.len) ;       for(int i = 0 ; i < s.len ; i++){           s.a[i] += a[i] + T.a[i] ;           if(s.a[i] > maxn){                s.a[i+1]++ ;                s.a[i] %= (maxn+1) ;           }       }       if(s.a[s.len]!=0) s.len++ ;       return s ;}char  str[58] ;BigNum dp[2][108] ;int  main(){     int i , m , p , j  , k ;     while(scanf("%d%d%d",&N ,&m ,&p) != EOF){          scanf("%s" , str) ;          to.clear() ;          for(i = 0 ; i < N ; i++) to[str[i]] = i ;          ac.Init() ;          while(p--){              scanf("%s" , str) ;              ac.Insert(str) ;          }          ac.makeAC() ;          ac.getMat() ;          int now = 0 ;          dp[now][0] = 1 ;          for(i = 1 ; i < matn ;i++)              dp[now][i] = BigNum(0) ;          for(i = 0 ; i < m ; i++){              now ^= 1 ;              for(j = 0 ; j < matn ; j++)                 dp[now][j] = BigNum(0) ;              for(j = 0 ; j < matn ; j++){                  for(k = 0 ; k < matn ; k++){                      if(mat[j][k] > 0)                         dp[now][k] = dp[now][k] + dp[now^1][j] * BigNum(mat[j][k]) ;                  }              }          }          BigNum sum = BigNum(0) ;          for(i = 0 ; i < matn ; i++)              sum = sum + dp[now][i]  ;          cout<< sum <<endl ;     }     return 0 ;}
0 0