zoukankan      html  css  js  c++  java
  • Substrings(SPOJ8222) (sam(后缀自动机))

    You are given a string (S) which consists of 250000 lowercase latin letters at most. We define (F(x)) as the maximal number of times that some string with length (x) appears in (S). For example for string 'ababa' (F(3)) will be 2 because there is a string 'aba' that occurs twice. Your task is to output (F(i)) for every (i) so that (1<=i<=|S|).

    Input

    String (S) consists of at most 250000 lowercase latin letters.

    Output

    Output (|S|) lines. On the (i-th) line output (F(i)).

    Example

    Input:

    ababa
    

    Output:

    3
    2
    2
    1
    1
    

    题意:

    给定一个字符串(S),令(F(x))表示(S)中所有长度为(x)的子串出现最多的次数,

    题解:

    把串塞进一个后缀自动机,在(parent)树上反向拓扑,求出每个长度时的最大值。

    #include<bits/stdc++.h>
    using namespace std;
    const int N=2000010;
    char s[N];
    int a[N],c[N],ans[N];
    struct SAM{
    	int last,cnt;
    	int size[N],ch[N][26],fa[N<<1],l[N<<1];
    	void ins(int c){
    		int p=last,np=++cnt;last=np;l[np]=l[p]+1;
    		for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;
    		if(!p)fa[np]=1;
    		else{
    			int q=ch[p][c];
    			if(l[p]+1==l[q])fa[np]=q;
    			else{
    				int nq=++cnt;l[nq]=l[p]+1;
    				memcpy(ch[nq],ch[q],sizeof ch[q]);
    				fa[nq]=fa[q];fa[q]=fa[np]=nq;
    				for(;ch[p][c]==q;p=fa[p])ch[p][c]=nq;
    			}
    		}
    		size[np]=1;
    	}
    	void build(char s[]){
    		int len=strlen(s+1);
    		last=cnt=1;
    		for(int i=1;i<=len;++i)ins(s[i]-'a');
    	}
    	void calc(int len){
    		for(int i=1;i<=cnt;++i)c[l[i]]++;
    		for(int i=1;i<=cnt;++i)c[i]+=c[i-1];
    		for(int i=1;i<=cnt;++i)a[c[l[i]]--]=i;
    		for(int i=cnt;i;--i){
    			int p=a[i];
    			size[fa[p]]+=size[p];
    			ans[l[p]]=max(ans[l[p]],size[p]);
    		}
    		for(int i=1;i<=len;++i){
    			printf("%d
    ",ans[i]);
    		}
    	}
    }sam;
    int main(){
    	cin>>s+1;
    	sam.build(s);
    	sam.calc(strlen(s+1));
    }
    
  • 相关阅读:
    11.关于django的content_type表
    6.re正则表达式
    1.关于python 的hmac加密
    4.瀑布流js
    10.django的一些方法理解
    7.一些比较有用的网站
    准备辞职了,走之前想解决的问题ptr 为空
    Oracle 代码生成小工具免费下载
    Js 中一个判断类型的方法
    jQuery 加上最后自己的验证
  • 原文地址:https://www.cnblogs.com/zhenglier/p/10097776.html
Copyright © 2011-2022 走看看