bzoj 3473 后缀自动机多字符串的子串处理方法
后缀自动机处理多字符串字串相关问题。
首先,和后缀数组一样,用分割符连接各字符串,然后建一个后缀自动机。
我们定义一个节点代表的字符串为它原本代表的所有串去除包含分割符后的串。每个节点代表的字符串的数量可以用DP来计算(不能用right集合来算了)。
对于原来n个串中的一个串,其所有前缀可以通过将该串放到自动机上跑来获得,对于某个前缀,其所有后缀包括在该前缀本身的节点以及parent树的祖先节点中。这样我们就获得访问某个串所有子串的技能了。
对于这道题,我们可以先建出后缀自动机,然后对于n个串中的每个串,找到包含其子串的所有节点(可以保证所有子串一定且唯一出现在某个节点中)。然后将它们的计数器+1。弄完后,对于每个节点,我们就可以知道其代表的串是n个串中多少个串的子串。
最后再对于每个串,找出所有字串(不同位置要区分),统计答案。
如果不同位置不区分,那么我们得到的节点不能重复。如果要区分,对于每个前缀,其parent树上的节点都要计算,即使以前被计算过(因为他们的结束位置不同,所以肯定要计算)。
/**************************************************************
Problem: 3473
User: idy002
Language: C++
Result: Accepted
Time:728 ms
Memory:120928 kb
****************************************************************/ #include <cstdio>
#include <cstring>
#include <cassert>
#include <algorithm>
#define N 200010
#define S 500010
#define P 18
using namespace std; typedef long long dnt; int n, k;
char buf[N], *shead[N];
int son[S][], pnt[S], val[S], ntot, last;
int head[S], dest[S], next[S], etot;
int dfn[S], dep[S], anc[S][P+], idgr[S], stk[S], qu[S], top, bg, ed, idc;
dnt dp[S], eff[S];
int log[S]; void init() {
ntot = last = ;
pnt[] = -;
}
void append( int c ) {
int p = last;
int np = ++ntot;
val[np] = val[p]+;
while( p!=- && !son[p][c] )
son[p][c]=np, p=pnt[p];
if( p==- ) {
pnt[np] = ;
} else {
int q=son[p][c];
if( val[q]==val[p]+ ) {
pnt[np] = q;
} else {
int nq = ++ntot;
memcpy( son[nq], son[q], sizeof(son[nq]) );
val[nq] = val[p]+;
pnt[nq] = pnt[q];
pnt[q] = pnt[np] = nq;
while( p!=- && son[p][c]==q )
son[p][c]=nq, p=pnt[p];
}
}
last = np;
}
void make_topo() {
for( int u=; u<=ntot; u++ ) {
for( int c=; c<=; c++ ) {
int v=son[u][c];
if( !v ) continue;
idgr[v]++;
}
}
qu[bg=ed=] = ;
while( bg<=ed ) {
int u=qu[bg++];
for( int c=; c<=; c++ ) {
int v=son[u][c];
if( !v ) continue;
idgr[v]--;
if( idgr[v]== )
qu[++ed] = v;
}
}
}
void dodp() {
make_topo();
dp[] = ;
for( int i=; i<=ed; i++ ) {
int u=qu[i];
for( int c=; c<=; c++ ) {
int v=son[u][c];
if( !v ) continue;
dp[v] += dp[u];
}
}
}
void adde( int u, int v ) {
etot++;
dest[etot] = v;
next[etot] = head[u];
head[u] = etot;
}
void build() {
for( int u=; u<=ntot; u++ )
adde( pnt[u], u );
}
void dfs( int u ) {
dfn[u] = ++idc;
for( int p=; p<=P && anc[u][p-]; p++ )
anc[u][p] = anc[anc[u][p-]][p-];
for( int t=head[u]; t; t=next[t] ) {
int v=dest[t];
dep[v] = dep[u]+;
anc[v][] = u;
dfs(v);
}
}
int lca( int u, int v ) {
if( dep[u]<dep[v] ) swap(u,v);
int t=dep[u]-dep[v];
for( int p=; t; p++,t>>= )
if( t& ) u=anc[u][p];
if( u==v ) return u;
for( int p=log[dep[u]]; anc[u][]!=anc[v][]; p-- )
if( anc[u][p]!=anc[v][p] ) u=anc[u][p],v=anc[v][p];
return anc[u][];
}
void fetch( char *s ) {
top = ;
int u = ;
for( int i=; s[i]; i++ ) {
int c=s[i]-'a'+;
u = son[u][c];
assert(u!=);
stk[++top] = u;
}
}
bool cmp( int u, int v ) {
return dfn[u]<dfn[v];
}
void effort( char *s ) {
fetch(s);
sort( stk+, stk++top, cmp );
eff[stk[]] += ;
for( int i=; i<=top; i++ ) {
int u=stk[i];
int ca=lca(u,stk[i-]);
eff[u] += ;
eff[ca] -= ;
}
}
void bfs() {
qu[bg=ed=] = ;
while( bg<=ed ) {
int u=qu[bg++];
for( int t=head[u]; t; t=next[t] ) {
int v=dest[t];
qu[++ed] = v;
}
}
for( int i=ed; i>=; i-- ) {
int u=qu[i];
for( int t=head[u]; t; t=next[t] ) {
int v=dest[t];
eff[u] += eff[v];
}
}
dp[] = ;
for( int i=; i<=ed; i++ ) {
int u=qu[i];
if( eff[u]>=k ) {
dp[u] = dp[pnt[u]]+dp[u];
} else {
dp[u] = dp[pnt[u]];
}
}
}
void query( char *s ) {
fetch(s);
sort( stk+, stk++top, cmp );
dnt rt = ;
for( int i=; i<=top; i++ ) {
int u=stk[i];
rt += dp[u];
}
printf( "%lld ", rt );
}
int main() {
// input and build sam
scanf( "%d%d", &n, &k );
init();
char *buf_cur = buf;
for( int i=; i<=n; i++ ) {
shead[i] = buf_cur;
scanf( "%s", shead[i] );
for( int j=; shead[i][j]; j++ )
append( shead[i][j]-'a'+ );
append( );
buf_cur += strlen(shead[i]) + ;
}
log[] = -;
for( int i=; i<=ntot; i++ ) log[i] = log[i>>]+;
// dodp to calc the number of real substring
dodp();
// build parent tree and its dfs order
build();
dfs();
// calc echo string's effort
for( int i=; i<=n; i++ )
effort( shead[i] );
// bfs to calc the sum of subans
bfs();
// query the answer of eacho string
for( int i=; i<=n; i++ )
query( shead[i] );
printf( "\n" );
}
bzoj 3473 后缀自动机多字符串的子串处理方法的更多相关文章
- POJ 1509 Glass Beads 后缀自动机 模板 字符串的最小表示
http://poj.org/problem?id=1509 后缀自动机其实就是一个压缩储存空间时间(对节点重复利用)的储存所有一个字符串所有子串的trie树,如果想不起来长什么样子可以百度一下找个图 ...
- Lexicographical Substring Search (spoj7259) (sam(后缀自动机)+第k小子串)
Little Daniel loves to play with strings! He always finds different ways to have fun with strings! K ...
- SPOJ LCS 后缀自动机找最大公共子串
这里用第一个字符串构建完成后缀自动机以后 不断用第二个字符串从左往右沿着后缀自动机往前走,如能找到,那么当前匹配配数加1 如果找不到,那么就不断沿着后缀树不断往前找到所能匹配到当前字符的最大长度,然后 ...
- [SPOJ1811]Longest Common Substring 后缀自动机 最长公共子串
题目链接:http://www.spoj.com/problems/LCS/ 题意如题目,求两个串的最大公共子串LCS. 首先对其中一个字符串A建立SAM,然后用另一个字符串B在上面跑. 用一个变量L ...
- SPOJ - SUBST1 New Distinct Substrings —— 后缀数组 单个字符串的子串个数
题目链接:https://vjudge.net/problem/SPOJ-SUBST1 SUBST1 - New Distinct Substrings #suffix-array-8 Given a ...
- bzoj 3676 后缀自动机+马拉车+树上倍增
思路:用马拉车把一个串中的回文串个数降到O(n)级别,然后每个串在后缀自动机上倍增找个数. #include<bits/stdc++.h> #define LL long long #de ...
- 不在B中的A的子串数量 HDU - 4416 (后缀自动机模板题目)
题目: 给定一个字符串a,又给定一系列b字符串,求字符串a的子串不在b中出现的个数. 题解: 先将所有的查询串放入后缀自动机(每次将sam.last=1)(算出所有子串个数) 然后将母串放入后缀自动机 ...
- BZOJ 4327 JSOI2012 玄武密码(后缀自动机)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=4327 [题目大意] 求每个子串在母串中的最长匹配 [题解] 对母串建立后缀自动机,用每 ...
- 【算法专题】后缀自动机SAM
后缀自动机是用于识别子串的自动机. 学习推荐:陈立杰讲稿,本文记录重点部分和感性理解(论文语言比较严格). 刷题推荐:[后缀自动机初探],题目都来自BZOJ. [Right集合] 后缀自动机真正优于后 ...
随机推荐
- 微信小程序实战--开发一个简单的快递单号查询
功能如图: 虽然工作中只负责小程序后台开发,但是还是小程序开发产生了浓厚的兴趣,官方文档也是超级详细了 这里就简单做一个快递单号的查询: 新建一个page: 接着就可以写wxml了.这里用一个简单的i ...
- vue2.0组件之间的传值
1.父子组件--props props:需要注意的是,虽然的是单向的数据流,但是如果传递的是数组或是对象这样的引用类型,子组件数据变化,父组件的数据通也会变化 子组件代码 <template&g ...
- Python [SSL: CERTIFICATE_VERIFY_FAILED] certificate verify failed 解决方法
一个搭建在SAE上的Django应用,使用新浪微博提供的Python SDK已经稳定运行一年有余,但最近开始持续出现微博认证失败的状况. 摘录微博python SDK的错误提示如下所示: ERROR: ...
- http://code52.org/DownmarkerWPF/
http://code52.org/DownmarkerWPF/ http://kb.cnblogs.com/page/132209/
- IntelliJ IDEA 修改IDE字体、代码字体。
IntelliJ IDEA 默认的 IDE 菜单字体太小,看着不舒服 ,我们调节下: ==============以上修改 仅仅针对的IDE字体,对代码的字体不生效. 所以如果代码 你觉得小 还得修改 ...
- Mysql建表语句
create table messages( id int unsigned primary key auto_increment, sender varchar(64) not null, gett ...
- Linux下fastbin利用小结——fd覆盖与任意地址free(House of Spirit)
linux下的fastbin是ctf中pwn题的重点出题点.去年(2015)中,XCTF就有两站是使用fastbin的利用作为pwn400的压轴题来出现,这也是我刚开始接触fastbin的利用,参考了 ...
- Html5和Css3扁平化风格网页
前言 扁平化概念的核心意义 去除冗余.厚重和繁杂的装饰效果.而具体表现在去掉了多余的透视.纹理.渐变以及能做出3D效果的元素,这样可以让“信息”本身重新作为核心被凸显出来.同时在设计元素上,则强调了抽 ...
- js各浏览器兼容取的元素的位置X坐标 Y坐标
JS code: function getElementPos(elementId) { var ua = navigator.userAgent.toLowerCase(); var isOpera ...
- #define const typedef
#define用法 1. 定义简单的常数:定义常量,便于修改 #define N 1000 2. 定义简单的函数:注意多使用括号 define可以像函数那样接受一些参数,如下: #define max ...