Problem C: 重复子串(string)
/*
一个性质? right集合中只有相邻的位置才会有用
那么考虑set启发式合并, 能够求出大概nlogn个有用的对
那么将这些对按照右端点排序, 查询也按照右端点排序就可以离线维护信息
然后需要维护答案的话?? 区间和一次函数取max???
应该有别的办法吧
普通线段树就好了
*/
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
#include<queue>
#include<iostream>
#define ll long long
#define mmp make_pair
#define M 200010
using namespace std;
int read() {
int nm = 0, f = 1;
char c = getchar();
for(; !isdigit(c); c = getchar()) if(c == '-') f = -1;
for(; isdigit(c); c = getchar()) nm = nm * 10 + c - '0';
return nm * f;
}
set<int> st[M];
set<int>::iterator it, tt;
int ch[M][26], fa[M], len[M], sz[M], ans[M], lst = 1, cnt = 1, n, m, tp;
char s[M];
vector<int> to[M];
struct Note {
int l, r, max;
bool operator < (const Note &b) const {
return this->r < b.r;
}
} note[M * 60], que[M];
void insert(int pl, int c) {
int p = ++cnt, f = lst;
lst = p;
len[p] = len[f] + 1;
while(f && !ch[f][c]) ch[f][c] = p, f = fa[f];
if(f == 0) fa[p] = 1;
else {
int q = ch[f][c];
if(len[q] == len[f] + 1) fa[p] = q;
else {
int nq = ++cnt;
len[nq] = len[f] + 1;
fa[nq] = fa[q];
memcpy(ch[nq], ch[q], sizeof(ch[q]));
fa[p] = fa[q] = nq;
while(f && ch[f][c] == q) ch[f][c] = nq, f = fa[f];
}
}
sz[p] = 1;
st[p].insert(pl);
}
int id[M];
bool cmp(int a, int b) {
return sz[a] > sz[b];
}
void dfs(int now) {
id[now] = now;
int flag = sz[now];
for(int i = 0; i < to[now].size(); i++) {
int vj = to[now][i];
dfs(vj);
sz[now] += sz[vj];
}
if(now == 1) return;
sort(to[now].begin(), to[now].end(), cmp);
if(to[now].size() != 0) {
id[now] = id[to[now][0]];
if(flag) {
int v = *st[now].begin();
it = st[id[now]].lower_bound(v);
if(it != st[id[now]].end()) {
note[++tp] = (Note) {
v, *it, len[now]
};
}
if(it != st[id[now]].begin()) {
it--;
note[++tp] = (Note) {
*it, v, len[now]
};
}
st[id[now]].insert(v);
}
for(int i = 1; i < to[now].size(); i++) {
int vj = to[now][i];
it = st[id[vj]].begin();
while(it != st[id[vj]].end()) {
int v = *it++;
tt = st[id[now]].lower_bound(v);
if(tt != st[id[now]].end()) {
note[++tp] = (Note) {
v, *tt, len[now]
};
}
if(tt != st[id[now]].begin()) {
tt--;
note[++tp] = (Note) {
*tt, v, len[now]
};
}
st[id[now]].insert(v);
}
}
}
}
struct Black_ {
#define ls now << 1
#define rs now << 1 | 1
#define lson l, mid, now << 1
#define rson mid + 1, r, now << 1 | 1
int t[M << 2], len[M << 2], ver[M << 2], maxx[M << 2];
void build(int l, int r, int now) {
len[now] = r - l + 1;
ver[now] = -0x3e3e3e3e;
if(l == r) return;
int mid = (l + r) >> 1;
build(lson), build(rson);
}
void add(int now, int v) {
ver[now] = max(ver[now], v);
maxx[now] = max(maxx[now], ver[now] + len[now]);
}
void pushup(int now) {
maxx[now] = max(maxx[now], max(maxx[ls], maxx[rs]));
}
void pushdown(int now) {
add(rs, ver[now]);
add(ls, ver[now] + len[rs]);
}
int que(int l, int r, int now, int ln, int rn) {
if(l > rn || r < ln) return 0;
if(l >= ln && r <= rn) return maxx[now];
int mid = (l + r) >> 1;
pushdown(now);
return max(que(lson, ln, rn), que(rson, ln, rn));
}
void modi(int l, int r, int now, int ln, int rn, int v) {
if(l > rn || r < ln) return;
if(l >= ln && r <= rn) {
add(now, v);
return;
}
int mid = (l + r) >> 1;
pushdown(now);
modi(lson, ln, rn, v + min(max(rn - mid, 0), len[rs]));
modi(rson, ln, rn, v);
pushup(now);
}
void modify(int r, int v) {
modi(1, n, 1, r - v + 1, r, 0);
}
int query(int l, int r) {
return que(1, n, 1, l, r);
}
} Seg;
int main() {
// freopen("string_example_3.in", "r", stdin); freopen("zz.out", "w", stdout);
n = read(), m = read();
scanf("%s", s + 1);
for(int i = 1; i <= n; i++) insert(i, s[i] - 'a');
for(int i = 2; i <= cnt; i++) to[fa[i]].push_back(i);
dfs(1);
for(int i = 1; i <= m; i++) que[i].l = read(), que[i].r = read(), que[i].max = i;
sort(note + 1, note + tp + 1);
// puts(s + 1);
// for(int i = 1; i <= tp; i++) {
// cerr<< note[i].l << " " << note[i].r << " " << note[i].max << "\n";
// }
// cerr << tp << "\n";
sort(que + 1, que + m + 1);
int tp1 = 1, tp2 = 1;
Seg.build(1, n, 1);
for(int i = 1; i <= n; i++) {
while(tp1 <= tp && note[tp1].r <= i) {
Seg.modify(note[tp1].l, note[tp1].max);
tp1++;
}
while(tp2 <= m && que[tp2].r <= i) {
ans[que[tp2].max] = Seg.query(que[tp2].l, que[tp2].r);
tp2++;
}
}
for(int i = 1; i <= m; i++) cout << ans[i] << "\n";
return 0;
}
Problem C: 重复子串(string)的更多相关文章
- POJ 1743 (后缀数组+不重叠最长重复子串)
题目链接: http://poj.org/problem?id=1743 题目大意:楼教主の男人八题orz.一篇钢琴谱,每个旋律的值都在1~88以内.琴谱的某段会变调,也就是说某段的数可以加减一个旋律 ...
- POJ 1743 Musical Theme (后缀数组,求最长不重叠重复子串)(转)
永恒的大牛,kuangbin,膜拜一下,Orz 链接:http://www.cnblogs.com/kuangbin/archive/2013/04/23/3039313.html Musical T ...
- poj 1743 男人八题之后缀数组求最长不可重叠最长重复子串
Musical Theme Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 14874 Accepted: 5118 De ...
- hdu_2668 Daydream O(n)求最长不重复子串
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2668 Daydream Time Limit: 2000/1000 MS (Java/Others) ...
- [Jobdu] 题目1530:最长不重复子串
题目描述: 最长不重复子串就是从一个字符串中找到一个连续子串,该子串中任何两个字符都不能相同,且该子串的长度是最大的. 输入: 输入包含多个测试用例,每组测试用例输入一行由小写英文字符a,b,c... ...
- ACdream 1430——SETI——————【后缀数组,不重叠重复子串个数】
SETI Time Limit: 4000/2000MS (Java/Others) Memory Limit: 128000/64000KB (Java/Others) Submit Statist ...
- POJ3693 Maximum repetition substring —— 后缀数组 重复次数最多的连续重复子串
题目链接:https://vjudge.net/problem/POJ-3693 Maximum repetition substring Time Limit: 1000MS Memory Li ...
- SPOJ - REPEATS —— 后缀数组 重复次数最多的连续重复子串
题目链接:https://vjudge.net/problem/SPOJ-REPEATS REPEATS - Repeats no tags A string s is called an (k,l ...
- [LeetCode] Longest Substring Without Repeating Characters 最长无重复子串
Given a string, find the length of the longest substring without repeating characters. For example, ...
随机推荐
- 第一份offer
11月6日参加的面试,今天签完三方,回头想想,感慨万千. (很多过程没有详细写,只保留了基本的客观事实,避免自吹嫌疑.....) 6号面试,当时来了能有100-200人,以川大和电子科大研究生为主,主 ...
- ngnix——FastCGI 相关参数调优
当 LNMP 组合工作时,首先是用户通过浏览器输入域名请求 Nginx Web 服务,如果请求的是静态资源,则由 Nginx 解析返回给用户:如果是动态请求(如 PHP),那么 Nginx 就会把它通 ...
- 使用RTL-SDR打开车门
在最近几年,入侵汽车在当代社会的黑客圈中成为热点,很多文章表明汽车产业还有很多东西等待完善,在本篇文章中,我会让你熟悉我一直研究的一些概念,以及如何在网状网络中使用一些便宜的部件渗透远程开门系统. 软 ...
- python day18--面向对象,继承
# class Animal: # breath = '呼吸' # # def __init__(self, name, sex, age): # self.name = name # self.se ...
- seajs引入jquery框架问题
seajs引入jquery框架时出现的问题 原因:由于seajs是遵循cmd规范,而jquery是遵循amd规范,所以在seajs调用jquery框架时,需要将amd转换为cmd 转换方法:jquer ...
- IDEA控制台乱码解决
打开Intellij的安装的bin目录(D:\Program Files\JetBrains\IntelliJ IDEA 14.0\bin ),找到上图的两个文件(根据你的系统是32位或64位选择其中 ...
- Spring3(一) 控制反转(IoC)和依赖注入(DI)
几个常用框架之间的关系 1 spring框架概述 1.1 什么是spring Spring是一个开源.轻量级的Java 开发框架.框架的主要优势之一就是其分层架构,分层架构允许使用者选 ...
- C语言笔记1
第一部分 计算机常识 1,人与计算机交流 2,人机交互方式 3,计算机的组成结构 4,计算机的系统组成 5,进制 6,程序和指令 第二部C开发环境 1 人与计算机交流 现实中人与人交流靠语言,那么人 ...
- Linux系统下用find命令查找最近修改过的文件
1.先进入指定目录,比如:cd /var/ftp/ss 2.输入指定命令,相关命令如下介绍: 最近24小时内修改过的文件 find ./ -mtime 0 前48~24小时修改过的文件 find ./ ...
- 第五节《Git基本操作》
我们给原来的数据打一个tag(标签),专业术语叫做“里程碑”,我们先不介绍里程碑的奥秘,只要知道里程碑无非也是一个引用而已. [root@git demo]# pwd/git/my/workspace ...