• HDU 4641


    动态更新后缀自动机,每次不断依据当前添加的节点不断往前寻找父节点上字符串最多可出现的次数

    这里为了减少运算,当父节点已经达到k次就不在往前寻找,因为之前的必然达到k次,也已经统计在内

     1 #include <cstdio>
     2 #include <cstring>
     3 #include <iostream>
     4 
     5 using namespace std;
     6 #define N 500010
     7 #define M 26
     8 
     9 struct SamNode{
    10     SamNode *son[26] , *f;
    11     int l , sc;
    12     void init(){
    13         f = NULL;
    14         for(int i=0 ; i<26 ; i++) son[i] = NULL;
    15         l = sc = 0;
    16     }
    17 }*root , *last , sam[N];
    18 
    19 int cnt , n , m , k , ret;
    20 char s[N];
    21 
    22 void init(){
    23     sam[0].init();
    24     root = last = &sam[cnt=0];
    25     ret = 0;
    26 }
    27 
    28 void add(int x)
    29 {
    30     SamNode *p = &sam[++cnt] , *jp=last;
    31     p->init();
    32     p->l = jp->l+1;
    33     last = p;
    34     for( ; jp&&!jp->son[x] ; jp=jp->f) jp->son[x]=p;
    35     if(!jp) p->f = root;
    36     else{
    37         if(jp->l+1 == jp->son[x]->l) p->f = jp->son[x];
    38         else{
    39             SamNode *r = &sam[++cnt] , *q = jp->son[x];
    40             r->init();
    41             *r = *q;
    42             r->l = jp->l+1;
    43             q->f = p->f = r;
    44             for( ; jp && jp->son[x]==q ; jp=jp->f) jp->son[x]=r;
    45         }
    46     }
    47     while(p!=root && p->sc<k){
    48         p->sc++;
    49         if(p->sc == k) ret+=(p->l-p->f->l);
    50         p = p->f;
    51     }
    52 }
    53 
    54 void solve()
    55 {
    56     for(int i=0 ; i<n ; i++) add(s[i]-'a');
    57  //   for(int i=0 ; i<=5 ; i++) cout<<i<<" :"<<sam[i].sc<<endl;
    58     int op;
    59     char c[3];
    60     for(int i=0 ; i<m ; i++){
    61         scanf("%d" , &op);
    62         if(op == 2) printf("%d
    " , ret);
    63         else{
    64             scanf("%s" , c);
    65             add(c[0]-'a');
    66         }
    67         //debug
    68     //    for(int i=0 ; i<=5 ; i++) cout<<i<<" :"<<sam[i].sc<<endl;
    69     }
    70 }
    71 
    72 int main()
    73 {
    74   //  freopen("a.in" , "r" , stdin);
    75     while(~scanf("%d%d%d" , &n , &m , &k))
    76     {
    77         scanf("%s" , s);
    78         init();
    79         solve();
    80     }
    81     return 0;
    82 }
  • 相关阅读:
    75分以下是文盲
    罗永浩最近怎么了
    北京奥运会赛事项目竞赛日程表
    从今天开始我的blog增加计数器
    07工作总结
    五道脑筋测试题,全答对的是天才
    转:国家名字的含义
    我的新装备双狂
    上海海鲜自助:喜多屋vs古象大酒店
    因为有了爱
  • 原文地址:https://www.cnblogs.com/CSU3901130321/p/4590851.html
Copyright © 2020-2023  润新知