zoukankan      html  css  js  c++  java
  • SPOJ LCS Longest Common Substring 和 LG3804 【模板】后缀自动机

    Longest Common Substring

    给两个串A和B,求这两个串的最长公共子串。

    no more than 250000

    分析

    参照OI wiki

    给定两个字符串 S 和 T ,求出最长公共子串,公共子串定义为在 S 和 T 中 都作为子串出现过的字符串 X 。

    我们为字符串 S 构造后缀自动机。

    我们现在处理字符串 T ,对于每一个前缀都在 S 中寻找这个前缀的最长后缀。换句话 说,对于每个字符串 T 中的位置,我们想要找到这个位置结束的 S 和 T 的最长公 共子串的长度。

    为了达到这一目的,我们使用两个变量,当前状态 v 和 当前长度 l 。这两 个变量描述当前匹配的部分:它的长度和它们对应的状态。

    一开始 v=t_0且 l=0 ,即,匹配为空串。

    现在我们来描述如何添加一个字符 T[i] 并为其重新计算答案:

    如果存在一个从 v 到字符 T[i] 的转移,我们只需要转移并让 l 自增一。
    如果不存在这样的转移,我们需要缩短当前匹配的部分,这意味着我们需要按照以下后 缀链接进行转移:

    v=link(v)

    与此同时,需要缩短当前长度。显然我们需要将 l 赋值为 len(v) ,因为经过这个后缀链接后我们到达的状态所对应的最长字符串是一个子串。

    如果仍然没有使用这一字符的转移,我们继续重复经过后缀链接并减小 l ,直到我们 找到一个转移或到达虚拟状态 -1 (这意味着字符 T[i] 根本没有在 S 中出现过, 所以我们设置 v=l=0 )。

    问题的答案就是所有 l 的最大值。

    这一部分的时间复杂度为 O(length(T)) ,因为每次移动我们要么可以使 l 增加一, 要么可以在后缀链接间移动几次,每次都减小 l 的值。

    时间复杂度\(O(|S| + |T|)\)

    co int N=5e5;
    namespace SAM
    {
    	int tot,last;
    	int ch[N][26],fail[N]={-1},len[N];
    	void extend(int k)
    	{
    		int cur=++tot;
    		len[cur]=len[last]+1;
    		int p=last;
    		while(~p&&!ch[p][k])
    		{
    			ch[p][k]=cur;
    			p=fail[p];
    		}
    		if(p==-1)
    			fail[cur]=0;
    		else
    		{
    			int q=ch[p][k];
    			if(len[q]==len[p]+1)
    				fail[cur]=q;
    			else
    			{
    				int clone=++tot;
    				std::copy(ch[q],ch[q]+26,ch[clone]);
    				fail[clone]=fail[q],len[clone]=len[p]+1;
    				while(~p&&ch[p][k]==q)
    				{
    					ch[p][k]=clone;
    					p=fail[p];
    				}
    				fail[cur]=fail[q]=clone;
    			}
    		}
    		last=cur;
    	}
    	void ins(char s[],int n)
    	{
    		for(int i=0;i<n;++i)
    			extend(s[i]-'a');
    	}
    	void solve(char s[],int n)
    	{
    		int ans=0,v=0,l=0;
    		for(int i=0;i<n;++i)
    		{
    			int k=s[i]-'a';
    			if(ch[v][k])
    				v=ch[v][k],++l;
    			else
    			{
    				while(~v&&!ch[v][k])
    					v=fail[v];
    				if(v==-1)
    					v=l=0;
    				else
    					l=len[v]+1,v=ch[v][k];
    			}
    			ans=std::max(ans,l);
    		}
    		printf("%d\n",ans);
    	}
    }
    char buf[N];
    
    int main()
    {
    //	freopen(".in","r",stdin);
    //	freopen(".out","w",stdout);
    	scanf("%s",buf);
    	SAM::ins(buf,strlen(buf));
    	scanf("%s",buf);
    	SAM::solve(buf,strlen(buf));
    	return 0;
    }
    

    【模板】后缀自动机

    给定一个只包含小写字母的字符串S,

    请你求出 S 的所有出现次数不为 1 的子串的出现次数乘上该子串长度的最大值。

    对于100%的数据,|S|<=106

    分析

    学习资料:OI wikiMenci's Blog
    ,张天杨《后缀自动机及其应用》。

    证明后缀自动机的状态数和转移数是线性的:

    不同的状态的Right集只能是不相交或者真包含关系。

    这种性质很容易让人想到树对不对,按照真包含关系做成树的话,大于一个元素的集合必定有大于等于两个儿子(保证了不会出现无用的链),只有一个元素的n个集合就是叶节点,所以总节点数小于等于\(2n-1\)。这样就说明状态数是线性的了。还差一点点,那就是转移数也是线性的。下面就来证明:

    取自动机的一个有向的生成树,使得初始状态出发可以到达所有状态。对于一条不在树上的边,我们可以在树上从初始状态走到它的起点,然后经过它并走到终态,走的这条路径对应的一定是原串的一个后缀(初态→终态,后缀自动机)。所以每条非树边一定可以对应某些后缀,而不同的两条非树边对应的后缀不可能有重复(若重复,跑该后缀的时候状态转移不唯一),所以非树边数小于等于\(n-1\)。算上生成树上的边,我们得到:转移数小于等于\((2n-1-1)+(n-1)=3n-3\)

    我们令叶子节点的size=1.暴力建出parent树然后dfs,求出每个节点的right集合size,然后求len×size的最大值就行了。

    时间复杂度:\(O(|S|)\)

    co int N=2e6;
    namespace SAM
    {
    	int tot,last;
    	int ch[N][26],fail[N],len[N],siz[N];
    	void init()
    	{
    		tot=last=0;
    		fail[0]=-1,len[0]=0;
    	}
    	void extend(int k)
    	{
    		int cur=++tot;
    		len[cur]=len[last]+1,siz[cur]=1;
    		int p=last;
    		while(~p&&!ch[p][k])
    		{
    			ch[p][k]=cur;
    			p=fail[p];
    		}
    		if(p==-1)
    			fail[cur]=0;
    		else
    		{
    			int q=ch[p][k];
    			if(len[p]+1==len[q])
    				fail[cur]=q;
    			else
    			{
    				int clone=++tot;
    				std::copy(ch[q],ch[q]+26,ch[clone]);
    				fail[clone]=fail[q],len[clone]=len[p]+1;
    				while(~p&&ch[p][k]==q)
    				{
    					ch[p][k]=clone;
    					p=fail[p];
    				}
    				fail[q]=fail[cur]=clone;
    			}
    		}
    		last=cur;
    	}
    	void ins(char*s,int n)
    	{
    		for(int i=0;i<n;++i)
    			extend(s[i]-'a');
    	}
    	
    	int nx[N],to[N];
    	ll ans;
    	void build()
    	{
    		for(int i=1;i<=tot;++i)
    			nx[i]=to[fail[i]],to[fail[i]]=i;
    	}
    	void dfs(int x)
    	{
    		for(int i=to[x];i;i=nx[i])
    		{
    			dfs(i);
    			siz[x]+=siz[i];
    		}
    		if(siz[x]>1)
    			ans=std::max(ans,(ll)siz[x]*len[x]);
    	}
    	void solve()
    	{
    		build();
    		dfs(0);
    		printf("%lld\n",ans);
    	}
    }
    char buf[N];
    
    int main()
    {
    //	freopen(".in","r",stdin);
    //	freopen(".out","w",stdout);
    	SAM::init();
    	scanf("%s",buf);
    	SAM::ins(buf,strlen(buf));
    	SAM::solve();
    	return 0;
    }
    
  • 相关阅读:
    Poj 1742 Coins(多重背包)
    Poj 2350 Above Average(精度控制)
    求二进制数中1的个数
    Poj 1659 Distance on Chessboard(国际象棋的走子规则)
    Poj 2411 Mondriaan's Dream(压缩矩阵DP)
    Poj 2136 Vertical Histogram(打印垂直直方图)
    Poj 1401 Factorial(计算N!尾数0的个数——质因数分解)
    poj 2390 Bank Interest(计算本利和)
    Poj 2533 Longest Ordered Subsequence(LIS)
    Poj 1887 Testing the CATCHER(LIS)
  • 原文地址:https://www.cnblogs.com/autoint/p/10345276.html
Copyright © 2011-2022 走看看