cf666E. Forensic Examination(广义后缀自动机 线段树合并)
时间:2022-06-19
本文章向大家介绍cf666E. Forensic Examination(广义后缀自动机 线段树合并),主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
题意
Sol
神仙题Orz
后缀自动机 + 线段树合并
首先对所有的(t_i)建个广义后缀自动机,这样可以得到所有子串信息。
考虑把询问离线,然后把(S)拿到自动机上跑,同时维护一下最长能匹配的位置,对于每个以(i)位置为右端点的询问我们需要找到(len)最小的状态满足(len[sta] >= pr - pl + 1)(这部分把每个以(i)为端点的询问排序后暴力跳即可,复杂度(O(n sqrt{n})))。那么现在的问题就是对于每个状态,如何知道他在每个(T_i)中的出现次数。
直接线段树合并一下就好啦
总复杂度:(O(大常数的nlogn + 小常数的nsqrt{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 %dn", ans[i].se, ans[i].fi);
}
return 0;
}
- 第九章:Shiro的Web——深入浅出学Shiro细粒度权限开发框架
- 第十章:Shiro的Cache——深入浅出学Shiro细粒度权限开发框架
- Appboy基于MongoDB的数据密集型实践
- 微信企业号登录授权Java实现获取员工userid根据userid换openid
- 微信支付-微信红包Java版本
- Universe入门
- 分享一款值得分享的写作工具
- 微信二次开发Java自定义菜单事件实现
- 微信OAuth授权获取用户OpenId-JAVA(个人经验)
- 【手写文字识别】-JavaAPI示例代码
- 【Python3-API】情感倾向分析示例代码
- SpringMVC+Hibernate +MySql+ EasyUI实现CRUD(一)
- 【Python3-API】通用文字识别示例代码
- Python入门教程之安装MyEclipse插件和安装Python环境
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 详解Numpy中的数组拼接、合并操作
- 批量更改图像尺寸到统一大小
- 张东升,我知道是你!如何使用GAN做一个秃头生产器
- Roslyn 在 NuGet 包中放注释 xml 文件的方法
- OpenCV3 图像的加载、修改、显示与保存
- 打卡群刷题总结0727——搜索旋转排序数组 II
- dotnet 的 TaskCompletionSource 的 TrySetResult 是线程安全
- OpenCV 摄像头与视频读取并显示
- OpenCV 边缘检测
- Spring Boot 集成 WebSocket 实现服务端推送消息到客户端
- OpenCV 绘制图形
- OpenCV 图片缩放
- AI树莓派——构建树莓派大脑(NCNN环境搭建)
- PC人脸识别登录,出乎意料的简单
- TensorLayer学习之简介与安装