Bzoj3998 弦论

物理题目传送门

求第k大的子串？SAM模板题啊

CLJ的论文都讲了怎么做啊，把自动机看成一个后缀Trie求出size让后像多叉平衡树那样乱搞就好了~

比前两个哈希的题好多了~ （顺便，hdu高亮好好看啊）

#pragma GCC opitmize("O3")
#pragma G++ opitmize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 1000010
using namespace std;
char c[N];
int s[N][26],mx[N],sz[N],f[N];
int n,m,T,k,cnt=1,lst=1,v[N],r[N],w[N];
inline int extend(int c){
    int p=lst,np=lst=++cnt,q,nq;
    mx[np]=mx[p]+1; sz[np]=1;
    for(;p&&!s[p][c];p=f[p]) s[p][c]=np;
    if(!p) return f[np]=1;
    q=s[p][c];
    if(mx[q]==mx[p]+1) f[np]=q;
    else{
        nq=++cnt;
        mx[nq]=mx[p]+1;
        f[nq]=f[q]; f[q]=f[np]=nq;
        memcpy(s[nq],s[q],26<<2);
        for(;p&&s[p][c]==q;p=f[p]) s[p][c]=nq;
    }
}
inline void dfs(int x){
    if(k<=sz[x]) return; else k-=sz[x];
    for(int j=0;j<26;++j)
        if(k>w[s[x][j]])k-=w[s[x][j]];
        else { putchar(j+'a'); dfs(s[x][j]); break; }
}
int main(){
    scanf("%s%d%d",c+1,&T,&k); n=strlen(c+1);
    for(int i=1;i<=n;++i) extend(c[i]-'a');
    for(int i=1;i<=cnt;++i) ++v[mx[i]];
    for(int i=1;i<=n;++i) v[i]+=v[i-1];
    for(int i=cnt;i;--i) r[v[mx[i]]--]=i;
    if(T) for(int i=cnt;i;--i) sz[f[r[i]]]+=sz[r[i]];
    else  for(int i=cnt;i;--i) sz[i]=1; sz[1]=*sz=0; memcpy(w,sz,sizeof w);
    for(int i=cnt;i;--i) for(int j=0;j<26;++j) w[r[i]]+=w[s[r[i]][j]];
    if(k>w[1]) puts("-1"); else dfs(1);
}