• BZOJ 2946 [Poi2000]公共串 (二分+Hash/二分+后缀数组/后缀自动机)


    求多串的最长公共字串.
    法1: 二分长度+hash 传送门
    法2: 二分+后缀数组 传送门
    法3: 后缀自动机
    拿第一个串建自动机,然后用其他串在上面匹配.每次求出SAM上每个节点的最长匹配长度后,再在全局取最小值(因为是所有串的公共串)就行了.

    CODE

    #include<bits/stdc++.h>
    using namespace std;
    char cb[1<<15],*cs=cb,*ct=cb;
    #define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
    template<class T>inline void read(T &res) {
    	char ch; int flg = 1; while(!isdigit(ch=getchar()))if(ch=='-')flg=-flg;
    	for(res=ch-'0';isdigit(ch=getchar());res=res*10+ch-'0'); res*=flg;
    }
    const int MAXN = 2005;
    const int C = 26;
    int m, lk[MAXN<<1], ch[MAXN<<1][C], len[MAXN<<1], mx[MAXN<<1], mn[MAXN<<1], bin[MAXN<<1], seq[MAXN<<1], sz, last;
    inline void init() {
        last = sz = 0; ++sz;
        lk[0] = -1; len[0] = 0;
    }
    inline void Copy(int A, int B) {
        lk[A] = lk[B];
        memcpy(ch[A], ch[B], sizeof ch[B]);
    }
    inline int insert(int p, int c) {
        int cur = sz++; last = cur;
        len[cur] = len[p] + 1;
        for(; ~p && !ch[p][c]; p = lk[p]) ch[p][c] = cur;
        if(p == -1) lk[cur] = 0;
        else {
            int q = ch[p][c];
            if(len[p] + 1 == len[q]) lk[cur] = q;
            else {
                int x = sz++;
                len[x] = len[p] + 1;
                Copy(x, q);
                lk[cur] = lk[q] = x;
                for(; ~p && ch[p][c] == q; p = lk[p]) ch[p][c] = x;
            }
        }
        return last;
    }
    inline void Match(char *s) {
        int p = 0, now = 0, j = 0;
        for(int c; s[j]; ++j) {
            if(ch[p][c=s[j]-'a'])
                ++now, p = ch[p][c];
            else {
                for(; ~p && !ch[p][c]; p = lk[p]);
                if(p == -1) now = p = 0;
                else now = len[p] + 1, p = ch[p][c];
            }
            mx[p] = max(mx[p], now);
        }
        for(int i = sz-1; i; --i) {
            mn[seq[i]] = min(mn[seq[i]], mx[seq[i]]);
            mx[lk[seq[i]]] = max(mx[lk[seq[i]]], mx[seq[i]]);
            mx[seq[i]] = 0;
        }
    }
    char s[MAXN];
    int main() {
    	read(m); scanf("%s", s);
    	if(--m == 0) return printf("%d
    ", strlen(s)), 0;
    	init(); int p = 0;
    	for(int i = 0; s[i]; ++i) p = insert(p, s[i]-'a');
        for(int i = 0; i < sz; ++i) ++bin[mn[i] = len[i]];
        for(int i = 1; i < sz; ++i) bin[i] += bin[i-1];
        for(int i = sz-1; ~i; --i) seq[--bin[len[i]]] = i;
    	while(m--) scanf("%s", s), Match(s);
        int ans = 0;
    	for(int i = 1; i < sz; ++i)
            ans = max(ans, mn[i]);
        printf("%d
    ", ans);
    }
    
    
    
  • 相关阅读:
    Mac之雕虫小技
    Python 合并两个字典(Dictionary)中相同key的value的方法
    linux只保留最新n个文件
    自动化测试笔记
    python实现一个判断时间的装饰,并抛出自定义异常
    python在字典列表中计算具有特定键的元素
    canvas游戏
    教为学:Oracle 11g OCP之路(七):数据结构
    教为学:Oracle 11g OCP之路(八):用户权限管理
    ADO.NET入门教程(六) 谈谈Command对象与数据检索
  • 原文地址:https://www.cnblogs.com/Orz-IE/p/12039314.html
Copyright © 2020-2023  润新知