[bzoj5337][loj2576][TJOI2018]str【字符串哈希】【后缀自动机】【dp】

【题目链接】
  https://www.lydsy.com/JudgeOnline/problem.php?id=5337
  https://loj.ac/problem/2576
【题解】
  随便怎么做都可以。
  比较简单的做法是:记 f[i][j] f [ i ] [ j ] 表示考虑了前 i i 个串,当前的末尾匹配到了第j位。那么可以用字符串哈希判断一段是否能匹配做到 O(1) O ( 1 ) 转移。
  时间复杂度 O(LA) O ( L ∗ ∑ A )
  但我做的时候脑回路比较清奇,把状态记到了后缀自动机上。。转移变成在自动机上游走。
  奇怪的复杂度 O(L) O ( L ∗ ∑ 子 串 长 度 )
【代码】

# include <bits/stdc++.h>
# define    ll      long long
# define    inf     0x3f3f3f3f
# define    N       30010
# define    K       110
using namespace std;
int read(){
    int tmp = 0, fh = 1; char ch = getchar();
    while (ch < '0' || ch > '9') {if (ch == '-') fh = -1; ch = getchar(); }
    while (ch >= '0' && ch <= '9'){tmp = tmp * 10 + ch - '0'; ch = getchar(); }
    return tmp * fh;
}
struct Trie{
    int son[26], step, pre;
}T[N];
const int  P = 1e9 + 7;
int f[K][N], n, l, num[K], place, use[N], sum[N], tag[N], las;
char s[N], tmp[N];
int extend(char x, int l, int las){
    x -= 'A';
    int now = ++place;
    T[now].step = l;
    while (T[las].son[x] == 0){
        T[las].son[x] = now;
        if (las == 0){
            T[now].pre = 0;
            return now;
        }
        else las = T[las].pre;
    }
    int q = T[las].son[x], p = las;
    if (T[p].step + 1 == T[q].step)
        T[now].pre = q;
        else {
            int np = ++place;
            T[np].step = T[p].step + 1;
            T[np].pre = T[q].pre;
            for (int i = 0; i < 26; i++)
                T[np].son[i] = T[q].son[i];
            T[q].pre = np, T[now].pre = np;
            for (; T[p].son[x] == q; p = T[p].pre)
                T[p].son[x] = np;
        }
    return now;
}
void solve(int x){
    if (use[x]) return;
    use[x] = true;
    if (tag[x]) sum[x] = 1;
    for (int i = 0; i < 26; i++){
        if (T[x].son[i] == 0) continue;
        solve(T[x].son[i]);
        sum[x] = (sum[x] + sum[T[x].son[i]]) % P;
    }
}
int main(){
//  freopen("A.in", "r", stdin);
//  freopen(".out", "w", stdout);
    n = read();
    scanf("\n%s", s + 1);
    l = strlen(s + 1), las = 0;
    for (int i = 1; i <= l; i++)
        las = extend(s[i], i, las);
    while (las != 0) tag[las] = true, las = T[las].pre;
    f[0][0] = 1;
    for (int i = 1; i <= n; i++){
        num[i] = read();
        for (int j = 1; j <= num[i]; j++){
            scanf("\n%s", tmp + 1);
            int len = strlen(tmp + 1);
            for (int k = 1; k <= len; k++) tmp[k] = tmp[k] - 'A';
            for (int k = 0; k <= place; k++){
                if (f[i - 1][k] == 0) continue;
                int now = k, flag = true;
                for (int t = 1; t <= len; t++){
                    if (T[now].son[tmp[t]] == 0){
                        flag = false;
                        break;
                    }
                    now = T[now].son[tmp[t]];
                }
                if (flag == true){
                    f[i][now] = f[i][now] + f[i - 1][k];
                    f[i][now] = (f[i][now] >= P) ? (f[i][now] - P) : f[i][now];
                }
            }
        }
    }
    int ans = 0;
    solve(0);
    for (int i = 0; i <= place; i++)
        ans = (ans + 1ll * f[n][i] * sum[i]) % P;
    printf("%d\n", ans);
    return 0;
}
### BZOJ1461 字符串匹配 题解 针对BZOJ1461字符串匹配问题解决方法涉及到了KMP算法以及树状数组的应用。对于此类问题,朴素的算法无法满足时间效率的要求,因为其复杂度可能高达O(ML²),其中M代表模式串的数量,L为平均长度[^2]。 为了提高效率,在这个问题中采用了更先进的技术组合——即利用KMP算法来预处理模式串,并通过构建失配树(也称为失败指针),使得可以在主串上高效地滑动窗口并检测多个模式串的存在情况。具体来说: - **前缀函数与KMP准备阶段**:先对每一个给定的模式串执行一次KMP算法中的pre_kmp操作,得到各个模式串对应的next数组。 - **建立失配树结构**:基于所有模式串共同构成的一棵Trie树基础上进一步扩展成带有失配链接指向的AC自动机形式;当遇到某个节点不存在对应字符转移路径时,则沿用该处失配链路直至找到合适的目标或者回到根部重新开始尝试其他分支。 - **查询过程**:遍历整个待查文本序列的同时维护当前状态处于哪一层级下的哪个子结点之中,每当成功匹配到完整的单词就更新计数值至相应位置上的f_i变量里去记录下这一事实。 下面是简化版Python代码片段用于说明上述逻辑框架: ```python from collections import defaultdict def build_ac_automaton(patterns): trie = {} fail = [None]*len(patterns) # 构建 Trie 树 for i,pattern in enumerate(patterns): node = trie for char in pattern: if char not in node: node[char]={} node=node[char] node['#']=i queue=[trie] while queue: current=queue.pop() for key,value in list(current.items()): if isinstance(value,int):continue if key=='#': continue parent=current[key] p=fail[current is trie and 0 or id(current)] while True: next_p=p and p.get(key,None) if next_p:break elif p==0: value['fail']=trie break else:p=fail[id(p)] if 'fail'not in value:value['fail']=next_p queue.append(parent) return trie,fail def solve(text, patterns): n=len(text) m=len(patterns) f=[defaultdict(int)for _in range(n)] ac_trie,_=build_ac_automaton(patterns) state=ac_trie for idx,char in enumerate(text+'$',start=-1): while True: trans=state.get(char,state.get('#',{}).get('fail')) if trans!=None: state=trans break elif '#'in state: state[state['#']['fail']] else: state=ac_trie cur_state=state while cur_state!={}and'#'in cur_state: matched_pattern_idx=cur_state['#'] f[idx][matched_pattern_idx]+=1 cur_state=cur_state['fail'] result=[] for i in range(len(f)-1): row=list(f[i].values()) if any(row): result.extend([sum((row[:j+1]))for j,x in enumerate(row[::-1])if x>0]) return sum(result) patterns=["ab","bc"] text="abc" print(solve(text,text)) #[^4] ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值