子串看起来就很SuffixStructures

于是上SAM

本来想着直接LCT

后来发现没法串定位(暴力匹配复杂度不对)

然后就离线吧,先建出来然后链加子树和,树剖就odk。

其实更直接的套路是线段树合并right集合维护。

这个题很像BZOJ3881 那个题更复杂一点,需要树链的并(到时候看看这个题可不可以加强一下扔到校内模拟赛

当然还有正确的ACA写法,直接FAIL树上主席树。(就更像BZOJ3881了

(但感觉都很暴力,我的想法貌似还挺好的)

顺便学习了广义SAM的正确写法。(发现之前的写法死活改不对

//Love and Freedom.
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define ll long long
#define inf 20021225
#define N 200010
#define pb push_back
using namespace std;
int read()
{
    int s=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
    return f*s;
}
struct node{int ch[26],fa,len;}t[N<<1];
struct edge{int to,lt;}e[N<<1];
int lt,poi,cnt,rt,ed[N],in[N<<1];
vector<int> endpos;
void insert(int c,int id)
{
    if(t[lt].ch[c])
    {
        int p=lt,np=lt=t[p].ch[c];
        if(t[np].len==t[p].len+1)    ed[id]=np,endpos.pb(np);
        else
        {
            int nq=++poi; ed[id]=nq; endpos.pb(nq); t[nq].len=t[p].len+1;
            memcpy(t[nq].ch,t[np].ch,sizeof(t[np].ch));
            t[nq].fa=t[np].fa; t[np].fa=nq;
            for(;p&&t[p].ch[c]==np;p=t[p].fa)    t[p].ch[c]=nq;
            lt=nq;
        }
        return;
    }
    int p=lt,np=lt=++poi; endpos.pb(np); ed[id]=np; t[np].len=t[p].len+1;
    for(;p&&!t[p].ch[c];p=t[p].fa) t[p].ch[c]=np;
    if(!p){t[np].fa=rt; return;} int q=t[p].ch[c];
    if(t[q].len==t[p].len+1){t[np].fa=q; return;}
    int nq=++poi; t[nq].fa=t[q].fa; t[q].fa=t[np].fa=nq;
    memcpy(t[nq].ch,t[q].ch,sizeof(t[q].ch)); t[nq].len=t[p].len+1;
    for(;p&&t[p].ch[c]==q;p=t[p].fa)    t[p].ch[c]=nq;
}
#define lowbit(x) (x&-x)
int tr[N<<1];
void add(int x,int v){while(x<=poi) tr[x]+=v,x+=lowbit(x);}
void modify(int l,int r,int v){add(l,v); add(r+1,-v);}
int query(int x){int a=0; while(x) a+=tr[x],x-=lowbit(x); return a;}
int dfn[N<<1],tot,sz[N<<1],son[N<<1],idfn[N<<1],top[N<<1],fa[N<<1];
void addedge(int x,int y){e[++cnt].to=y; e[cnt].lt=in[x]; in[x]=cnt;}
void build(){for(int i=2;i<=poi;i++)    addedge(t[i].fa,i);}
void dfs(int x,int fr)
{
    fa[x]=fr; sz[x]=1;
    for(int i=in[x];i;i=e[i].lt)
    {
        int y=e[i].to; dfs(y,x);
        sz[x]+=sz[y]; if(sz[y]>=sz[son[x]])    son[x]=y;
    }
}
void dfs2(int x,int t)
{
    top[x]=t; dfn[x]=++tot; idfn[tot]=x; if(!son[x])    return;    dfs2(son[x],t);
    for(int i=in[x];i;i=e[i].lt)    if(e[i].to!=son[x])    dfs2(e[i].to,e[i].to);
}
void modify(int x,int v)
{
    while(x)    modify(dfn[top[x]],dfn[x],v), x=fa[top[x]];
}
int ask(int x){return query(dfn[x]);}
#define Q 500010
int n,q; char ch[N]; int ans[Q],pre[N];
struct qry{int id,opt,pos,x;}r[Q<<1]; int qwq;
bool operator<(qry a,qry b){return a.pos<b.pos;}
int main()
{
    n=read(),q=read(); poi=rt=lt=1; int tot=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%s",ch+1); lt=rt; int sn=strlen(ch+1);
        for(int j=1;j<=sn;j++)    ++tot,insert(ch[j]-'a',i);
        pre[i]=tot;
    }
    int L,R,K; build(); dfs(1,0); dfs2(1,1);
    for(int i=1;i<=q;i++)
        L=read(),R=read(),K=read(),
        r[++qwq].id=i,r[qwq].opt=-1,r[qwq].pos=pre[L-1],r[qwq].x=ed[K],
        r[++qwq].id=i,r[qwq].opt=1,r[qwq].pos=pre[R],r[qwq].x=ed[K];
    int p=1; sort(r+1,r+qwq+1);
    while(p<=qwq&&!r[p].pos)    p++;
    for(int i=0;i<endpos.size();i++)
    {
        int pos=endpos[i]; modify(pos,1);
        while(p<=qwq&&r[p].pos==i+1)    ans[r[p].id]+=r[p].opt*ask(r[p].x),p++;
    }
    for(int i=1;i<=q;i++)    printf("%d\n",ans[i]);
    return 0;
}
View Code
12-28 10:40