怎么全是广义后缀自动机,我\(AC\)自动机不服

这道题可以使用的算法很多,\(SA\)或者\(SAM\)应该都可以

但是我都不会

但是这毕竟是一个多串匹配问题,\(AC\)自动机还是可以刚一刚的

我们先考虑一下暴力做法

先将操作离线下来,之后对于所有的询问串建立\(AC\)自动机,之后我们把\(n\)个给出串放到\(AC\)自动机上跑一遍,统计一下有多少个询问串可以是其子串,给这些询问串打上标记,最后输出每一个询问串被打了几次标记就好了

但是这个做法需要我们暴力去跳\(fail\)指针,于是很容易被卡成平方级别

思考一下跳\(fail\)的实质就是在\(fail\)树上暴力从当前点访问到根,而最后的查询也只是查询一个子树内部不同的标记数

于是我们熟练地想到了把\(fail\)树抽离出来,现在的问题转化为一个子树内部有多少个不同颜色的标记

是不是非常眼熟,这不就是HH的项链多了个\(dfs\)序吗

于是现在莫队或者树状数组都可以上了,不过可以树状数组谁写莫队啊

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<vector>
#define re register
#define maxn 360005
#define lowbit(x) ((x)&(-x))
int n,m,num,__,cnt;
struct E {int v,nxt;}e[maxn<<1];
struct Ask{int x,y,rk;}a[maxn];
char S[maxn];
std::string s[10005];
std::vector<int> v[maxn],pre[maxn];
int son[360005][26],fail[maxn],to[maxn];
int deep[maxn],dfn[maxn],head[maxn],sum[maxn],lst[maxn];
int c[maxn],f[maxn],E[100005],Ans[maxn];
inline void add_edge(int x,int y){e[++num].v=y,e[num].nxt=head[x],head[x]=num;}
inline int add(int x,int val){for(re int i=x;i<=__;i+=lowbit(i)) c[i]+=val;}
inline int ask(int x){int ans=0;for(re int i=x;i;i-=lowbit(i)) ans+=c[i];return ans;}
inline void ins(int x)
{
int now=0;
scanf("%s",S+1);
int len=strlen(S+1);
for(re int i=1;i<=len;i++)
{
if(!son[now][S[i]-'a']) son[now][S[i]-'a']=++cnt;
now=son[now][S[i]-'a'];
}
E[x]=now;
}
inline void Build()
{
std::queue<int> q;
for(re int i=0;i<26;i++) if(son[0][i]) q.push(son[0][i]);
while(!q.empty())
{
int k=q.front();
q.pop();
if(k) add_edge(fail[k],k);
for(re int i=0;i<26;i++)
if(son[k][i]) fail[son[k][i]]=son[fail[k]][i],q.push(son[k][i]);
else son[k][i]=son[fail[k]][i];
}
}
void dfs(int x)
{
dfn[++__]=x,to[x]=__;
sum[x]=1;
for(re int i=head[x];i;i=e[i].nxt)
if(!deep[e[i].v])
{
deep[e[i].v]=deep[x]+1;
dfs(e[i].v);
sum[x]+=sum[e[i].v];
}
}
inline void query(int x)
{
int len=s[x].size();
int now=0;
for(re int i=0;i<len;i++)
{
now=son[now][s[x][i]-'a'];
if(!f[now]) v[now].push_back(x);
if(!f[fail[now]])
v[fail[now]].push_back(x);
f[fail[now]]=f[now]=1;
}
now=0;
for(re int i=0;i<len;i++)
{
now=son[now][s[x][i]-'a'];
f[now]=f[fail[now]]=0;
}
}
inline int cmp(Ask A,Ask B){return A.y<B.y;}
int main()
{
scanf("%d%d",&n,&m);
for(re int i=1;i<=n;i++) std::cin>>s[i];
for(re int i=1;i<=m;i++) ins(i);
Build(),deep[0]=1,dfs(0);
for(re int i=1;i<=n;i++) query(i);
for(re int i=1;i<=m;i++) a[i].x=to[E[i]],a[i].y=to[E[i]]+sum[E[i]]-1,a[i].rk=i;
std::sort(a+1,a+m+1,cmp);
for(re int i=1;i<=__;i++)
for(re int j=0;j<v[dfn[i]].size();j++)
pre[i].push_back(lst[v[dfn[i]][j]]),lst[v[dfn[i]][j]]=i;
int top=1;
for(re int i=1;i<=__;i++)
{
for(re int j=0;j<pre[i].size();j++)
{
add(i,1);
if(pre[i][j])
add(pre[i][j],-1);
}
while(a[top].y==i)
Ans[a[top].rk]=ask(a[top].y)-ask(a[top].x-1),top++;
}
for(re int i=1;i<=m;i++) printf("%d\n",Ans[i]);
return 0;
}

SPOJ8093【JZPGYZ - Sevenk Love Oimaster】的更多相关文章

  1. SP8093 JZPGYZ - Sevenk Love Oimaster 解题报告

    SP8093 JZPGYZ - Sevenk Love Oimaster 题目大意 给定\(n(n\le 10000)\)个模板串,以及\(m(m\le 60000)\)个查询串(模板串总长\(\le ...

  2. 【洛谷 SP8093】 JZPGYZ - Sevenk Love Oimaster(后缀自动机)

    题目链接 广义sam.. #include <cstdio> #include <cstring> #include <algorithm> using names ...

  3. SPOJ 8093 JZPGYZ - Sevenk Love Oimaster

    思路 可以用复杂度不对的做法水过去 相当于求parent树子树中的颜色种数,可以离线后树状数组(HH的项链,询问右端点排序之后维护last),dsu on tree,莫队都可以 但是也可以记录每个点上 ...

  4. SP8093 JZPGYZ - Sevenk Love Oimaster(广义后缀自动机)

    题意 题目链接 Sol 广义后缀自动机板子题..和BZOJ串那个题很像 首先建出询问串的SAM,然后统计一下每个节点被多少个串包含 最后直接拿询问串上去跑就行了 #include<bits/st ...

  5. SP8093 JZPGYZ - Sevenk Love Oimaster(SAM)

    /* 打模板题啊 每个串影响到的集合直接枚举跳parent处理即可 */ #include<cstdio> #include<algorithm> #include<cs ...

  6. SP8093 JZPGYZ - Sevenk Love Oimaster

    传送门 广义后缀自动机-- 其实也不是很难理解,就是每次SAM插入一个串之后,插入新的串的时候,要把last重新调到1的位置,共用一些节点. 这个题我们首先要预处理出来每个状态被多少个串共用.挺暴力的 ...

  7. BZOJ 2780: [Spoj]8093 Sevenk Love Oimaster [广义后缀自动机]

    JZPGYZ - Sevenk Love Oimaster     Oimaster and sevenk love each other.       But recently,sevenk hea ...

  8. 【BZOJ2780】【SPOJ】Sevenk Love Oimaster(后缀自动机)

    [BZOJ2780][SPOJ]Sevenk Love Oimaster(后缀自动机) 题面 BZOJ 洛谷 题解 裸的广义后缀自动机??? 建立广义后缀自动机建立出来之后算一下每个节点被几个串给包括 ...

  9. 【BZOJ2780】[Spoj]8093 Sevenk Love Oimaster 广义后缀自动机

    [BZOJ2780][Spoj]8093 Sevenk Love Oimaster Description Oimaster and sevenk love each other.     But r ...

随机推荐

  1. 【OpenCV-python】CV2和PIL按box信息实现图像裁剪

    # 用cv2实现裁剪 import cv2 import os img = cv2.imread("./test_and_verification/1406170100001.jpg&quo ...

  2. Python 字符串 (str)

    作者博文地址:https://www.cnblogs.com/liu-shuai/ Python字符串的常用操作包括以下但不限于以下操作: 1 字符串的替换.删除.切片.复制.连接.比较.查找.分割等 ...

  3. elasticsearch清空type下的所有数据

    DELETE /twitter/tweet/_query{"query": { "match_all": {} }} 注:twitter为index,tweet ...

  4. Vue2.0以后,有哪些变化

    最近移动端项目版本升级,Vue由之前的1.0升级到2.3,那么,Vue2.0之后,有哪些细节的变化呢,现在总结如下: 1.在每个组件模板,不再支持片段代码 组件中模板: 之前: <templat ...

  5. vs2017取消起始页(设定起始页)/(.ashx文件的添加)

    1.设定起始页:在要设定起始页的视图(如Index视图)上直接右击,然后点击“设为起始页” 2.取消起始页:点击菜单栏“调试”,然后点击最后一行“项目名+属性”,打开后选中左侧栏中的“web”,选中“ ...

  6. 如何查询日志文件中的所有ip,正则表达式

    IPV4必须满足以下四条规则: 1.任何一个1位或2位数字,即0-99: 2.任何一个以1开头的3位数字,即100-199: 3.任何一个以2开头.第2位数字是0-4之间的3位数字,即200-249: ...

  7. Android开发ListView嵌套ImageView实现单选按钮

    做Android开发两年的时间,技术稍稍有一些提升,刚好把自己实现的功能写出来,记录一下,如果能帮助到同行的其他人,我也算是做了件好事,哈哈!!废话不多说,先上个图. 先上一段代码: if (last ...

  8. POJ1659 Frogs' Neighborhood(Havel–Hakimi定理)

    题意 题目链接 \(T\)组数据,给出\(n\)个点的度数,问是否可以构造出一个简单图 Sol Havel–Hakimi定理: 给定一串有限多个非负整数组成的序列,是否存在一个简单图使得其度数列恰为这 ...

  9. for计算位置坐标

    代码1效果图: 代码2效果图: <!DOCTYPE html> <html> <head> <meta charset="UTF-8"&g ...

  10. bootstrap模态框实现相对定位拖拽

    1.正常的拖拽是用绝对定位absolute来实现的,可是bootstrap的模态框是用relative,为了统一更改方便,就照着相对定位来实现拖拽效果. $(".modal .modal-h ...