题意

题目链接

Sol

神仙题Orz

后缀自动机 + 线段树合并

首先对所有的\(t_i\)建个广义后缀自动机,这样可以得到所有子串信息。

考虑把询问离线,然后把\(S\)拿到自动机上跑,同时维护一下最长能匹配的位置,对于每个以\(i\)位置为右端点的询问我们需要找到\(len\)最小的状态满足\(len[sta] >= pr - pl + 1\)(这部分把每个以\(i\)为端点的询问排序后暴力跳即可,复杂度\(O(n \sqrt{n})\))。那么现在的问题就是对于每个状态,如何知道他在每个\(T_i\)中的出现次数。

直接线段树合并一下就好啦

总复杂度:\(O(大常数的nlogn + 小常数的n\sqrt{n})\)

疯狂wa的原因居然是pair<int, int>不支持自定义小于号

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 2e6 + 10;
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, M, Q;
struct Pair {
int fi, se;
};
bool operator < (const Pair &a, const Pair &b) {
return (a.fi < b.fi) || (a.fi == b.fi && a.se > b.se);
}
Pair operator + (const Pair &a, const Pair &b) {
return {a.fi + b.fi, a.se};
}
Pair ans[MAXN];
char S[MAXN];
string T[MAXN];
struct Query {
int pl, l, r, id;
bool operator < (const Query &rhs) const {
return pl < rhs.pl;
}
};
vector<Query> qry[MAXN];
namespace Seg {
int root[MAXN], ls[MAXN], rs[MAXN], cnt;
Pair mx[MAXN]; void update(int k) {
mx[k] = max(mx[ls[k]], mx[rs[k]]);
}
void IntAdd(int &k, int l, int r, int p, int v) {
if(!k) k = ++cnt;
if(l == r) {mx[k].fi++; mx[k].se = l; return ;}
int mid = l + r >> 1;
if(p <= mid) IntAdd(ls[k], l, mid, p, v);
else IntAdd(rs[k], mid + 1, r, p, v);
update(k);
}
int Merge(int x, int y) {
if(!x || !y) return x ^ y;
int nw = ++cnt; mx[nw] = mx[x];
if(!ls[x] && !rs[x]) {mx[nw].fi += mx[y].fi; return nw;}
ls[nw] = Merge(ls[x], ls[y]);
rs[nw] = Merge(rs[x], rs[y]);
update(nw);
return nw;
}
Pair Query(int k, int l, int r, int ql, int qr) {
if(!k) return {0, 0};
if(ql <= l && r <= qr) return mx[k];
int mid = l + r >> 1;
if(ql > mid) return Query(rs[k], mid + 1, r, ql, qr);
else if(qr <= mid) return Query(ls[k], l, mid, ql, qr);
else {
Pair al = Query(ls[k], l, mid, ql, qr);
Pair ar = Query(rs[k], mid + 1, r, ql, qr);
if((al.fi > ar.fi) || (al.fi == ar.fi && al.se < ar.se)) return al;
else return ar;
}
}
}
namespace SAM { int fa[MAXN], ch[MAXN][26], len[MAXN], root = 1, las = 1, tot = 1;
vector<int> par[MAXN];
void insert(int x, int id) {
int now = ++tot, pre = las; las = now; len[now] = len[pre] + 1;
if(id) Seg::IntAdd(Seg::root[now], 1, M, id, 1);
for(; pre && !ch[pre][x]; pre = fa[pre]) ch[pre][x] = now;
if(!pre) fa[now] = root;
else {
int q = ch[pre][x];
if(len[pre] + 1 == len[q]) fa[now] = q;
else {
int nq = ++tot; fa[nq] = fa[q]; len[nq] = len[pre] + 1;
memcpy(ch[nq], ch[q], sizeof(ch[q]));
for(; pre && ch[pre][x] == q; pre = fa[pre]) ch[pre][x] = nq;
fa[q] = fa[now] = nq;
}
} }
void Build() {
for(int i = 1; i <= tot; i++) par[fa[i]].push_back(i);
}
void dfs(int x) {
for(auto &to : par[x]) {
dfs(to);
Seg::root[x] = Seg::Merge(Seg::root[x], Seg::root[to]);
}
}
void work() {
int now = root, dl = 0;
for(int i = 1; i <= N; i++) {
int nxt = S[i] - 'a';
while(!ch[now][nxt] && now) now = fa[now], dl = len[now];
if(!now) {
now = 1; dl = 0;
for(auto &q : qry[i]) ans[q.id].se = q.l;
continue;
}
now = ch[now][nxt]; dl++;
int t = now;
for(auto &q : qry[i]) {
if(dl < i - q.pl + 1) {ans[q.id].se = q.l; continue;}
while(len[fa[t]] >= i - q.pl + 1) t = fa[t];
ans[q.id] = Seg::Query(Seg::root[t], 1, M, q.l, q.r);
if(!ans[q.id].fi) ans[q.id].se = q.l;
}
}
} }
int main() {
// freopen("a.in", "r", stdin);
scanf("%s", S + 1);
N = strlen(S + 1);
cin >> M;
for(int i = 1; i <= M; i++) {
cin >> T[i];
string &ns = T[i];
for(int j = 0; j < ns.length(); j++) SAM::insert(ns[j] - 'a', i);
SAM::las = 1;
}
SAM::Build();
SAM::dfs(1);
// for(int i = 1; i <= Seg::cnt; i++) printf("%d ", Seg::mx[i]);
cin >> Q;
for(int i = 1; i <= Q; i++) {
int l = read(), r = read(), pl = read(), pr = read();
qry[pr].push_back({pl, l, r, i});
}
for(int i = 1; i <= N; i++) stable_sort(qry[i].begin(), qry[i].end());
SAM::work();
for(int i = 1; i <= Q; i++) {
printf("%d %d\n", ans[i].se, ans[i].fi);
}
return 0;
}

cf666E. Forensic Examination(广义后缀自动机 线段树合并)的更多相关文章

  1. 【CF666E】Forensic Examination - 广义后缀自动机+线段树合并

    广义SAM专题的最后一题了……呼 题意: 给出一个长度为$n$的串$S$和$m$个串$T_{1\cdots m}$,给出$q$个询问$l,r,pl,pr$,询问$S[pl\cdots pr]$在$T_ ...

  2. CF 666E Forensic Examination——广义后缀自动机+线段树合并

    题目:http://codeforces.com/contest/666/problem/E 对模式串建广义后缀自动机,询问的时候把询问子串对应到广义后缀自动机的节点上,就处理了“区间”询问. 还要处 ...

  3. [CF666E]Forensic Examination:后缀自动机+线段树合并

    分析 用到了两个小套路: 使用线段树合并维护广义后缀自动机的\(right\)集合. 查询\(S[L,R]\)在\(T\)中的出现次数:给\(T\)建SAM,在上面跑\(S\),跑到\(R\)的时候先 ...

  4. CF666E Forensic Examination(后缀自动机+线段树合并)

    给你一个串S以及一个字符串数组T[1..m],q次询问,每次问S的子串S[pl..pr]在T[l..r]中的哪个串里的出现次数最多,并输出出现次数. 如有多解输出最靠前的那一个. 我们首先对m个字符串 ...

  5. Codeforces.666E.Forensic Examination(广义后缀自动机 线段树合并)

    题目链接 \(Description\) 给定串\(S\)和\(m\)个串\(T_i\).\(Q\)次询问,每次询问\(l,r,p_l,p_r\),求\(S[p_l\sim p_r]\)在\(T_l\ ...

  6. CF666E Forensic Examination 广义SAM、线段树合并、倍增、扫描线

    传送门 朴素想法:对\(M\)个匹配串\(T_1,...,T_M\)建立广义SAM,对于每一次询问,找到这个SAM上\(S[pl...pr]\)对应的状态,然后计算出对于每一个\(i \in [l,r ...

  7. CF666E Forensic Examination 广义后缀自动机_线段树合并_树上倍增

    题意: 给定一个串 $S$ 和若干个串 $T_{i}$每次询问 $S[pl..pr]$ 在 $Tl..Tr$ 中出现的最多次数,以及出现次数最多的那个串的编号. 数据范围: 需要离线 题解:首先,很常 ...

  8. Codeforces 666E Forensic Examination(广义后缀自动机+线段树合并)

    将所有串(包括S)放一块建SAM.对于询问,倍增定位出该子串所在节点,然后要查询的就是该子串在区间内的哪个字符串出现最多.可以线段树合并求出该节点在每个字符串中的出现次数. #include<b ...

  9. BZOJ3413: 匹配(后缀自动机 线段树合并)

    题意 题目链接 Sol 神仙题Orz 后缀自动机 + 线段树合并... 首先可以转化一下模型(想不到qwq):问题可以转化为统计\(B\)中每个前缀在\(A\)中出现的次数.(画一画就出来了) 然后直 ...

随机推荐

  1. Jquery 图片走马灯效果原理

    本篇只讲解水平走马灯效果,垂直向上走马灯效果不讲解,原理一样,但是水平走马灯效果有一个小坑.待会讲解 照例先上代码: HTML: <div class="box"> & ...

  2. iOS编码规范(简版)

    1. 总体指导原则 [规则1-1]首先是为人编写程序,其次才是计算机. 说明:这是软件开发的基本要点,软件的生命周期贯穿产品的开发.测试.生产.用户使用.版本升级和后期维护等长期过程,只有易读.易维护 ...

  3. springboot之定时任务

    定时线程 说到定时任务,通常会想到JDK自带的定时线程来执行,定时任务. 回顾一下定时线程池. public static ScheduledExecutorService newScheduledT ...

  4. MyBatis 源码分析 - 配置文件解析过程

    * 本文速览 由于本篇文章篇幅比较大,所以这里拿出一节对本文进行快速概括.本篇文章对 MyBatis 配置文件中常用配置的解析过程进行了较为详细的介绍和分析,包括但不限于settings,typeAl ...

  5. 通俗理解N-gram语言模型。(转)

    从NLP的最基础开始吧..不过自己看到这里,还没做总结,这里有一篇很不错的解析,可以分享一下. N-gram语言模型 考虑一个语音识别系统,假设用户说了这么一句话:“I have a gun”,因为发 ...

  6. RSA实现JS前端加密,PHP后端解密

    web前端,用户注册与登录,不能直接以明文形式提交用户密码,容易被截获,这时就引入RSA. 前端加密 需引入4个JS扩展文件,jsbn.js.prng4.js.rng.js和rsa.js. <h ...

  7. postgres Date/Time 学习笔记

    一.Date/Time Types 参考文档:https://www.postgresql.org/docs/9.2/static/datatype-datetime.html Types 别名 ti ...

  8. LeetCode:104_Maximum Depth of Binary Tree | 二叉树的最大深度 | Easy

    要求:求二叉树的深度(二叉树的深度为最远叶子节点到根节点的距离,即根节点到最远叶子节点的距离) Given a binary tree, find its maximum depth. The max ...

  9. .NET 构造DataTable返回多个json值

    有时候我们使用Ajax链接一般处理程序需要返回多个值,然而这些数据并非在一个查询表内,此时便想到构造一个虚拟的DataTable,这样就可以返回多个值了(当然有很多办法,这是其中一种 ). 首先我们需 ...

  10. 使用JDBC连接各种数据库

    项目需要连多钟数据库,把链接的关键部分在这里记录下来. 1.Oracle8/8i/9i数据库(thin模式) Class.forName("oracle.jdbc.driver.Oracle ...