/*
一个性质? 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)的更多相关文章

  1. POJ 1743 (后缀数组+不重叠最长重复子串)

    题目链接: http://poj.org/problem?id=1743 题目大意:楼教主の男人八题orz.一篇钢琴谱,每个旋律的值都在1~88以内.琴谱的某段会变调,也就是说某段的数可以加减一个旋律 ...

  2. POJ 1743 Musical Theme (后缀数组,求最长不重叠重复子串)(转)

    永恒的大牛,kuangbin,膜拜一下,Orz 链接:http://www.cnblogs.com/kuangbin/archive/2013/04/23/3039313.html Musical T ...

  3. poj 1743 男人八题之后缀数组求最长不可重叠最长重复子串

    Musical Theme Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 14874   Accepted: 5118 De ...

  4. hdu_2668 Daydream O(n)求最长不重复子串

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2668 Daydream Time Limit: 2000/1000 MS (Java/Others)  ...

  5. [Jobdu] 题目1530:最长不重复子串

    题目描述: 最长不重复子串就是从一个字符串中找到一个连续子串,该子串中任何两个字符都不能相同,且该子串的长度是最大的. 输入: 输入包含多个测试用例,每组测试用例输入一行由小写英文字符a,b,c... ...

  6. ACdream 1430——SETI——————【后缀数组,不重叠重复子串个数】

    SETI Time Limit: 4000/2000MS (Java/Others) Memory Limit: 128000/64000KB (Java/Others) Submit Statist ...

  7. POJ3693 Maximum repetition substring —— 后缀数组 重复次数最多的连续重复子串

    题目链接:https://vjudge.net/problem/POJ-3693 Maximum repetition substring Time Limit: 1000MS   Memory Li ...

  8. SPOJ - REPEATS —— 后缀数组 重复次数最多的连续重复子串

    题目链接:https://vjudge.net/problem/SPOJ-REPEATS REPEATS - Repeats no tags  A string s is called an (k,l ...

  9. [LeetCode] Longest Substring Without Repeating Characters 最长无重复子串

    Given a string, find the length of the longest substring without repeating characters. For example, ...

随机推荐

  1. Android 音视频深入 七 学习之路的总结和资料分享

    说个实话一开始我对基于Android如何开发音视频很迷茫,甚至对音视频开发都不是很明白,我看了Android 音视频开发入门指南 http://blog.51cto.com/ticktick/1956 ...

  2. Java启动工程时,加载固定数据到Map中(不用每次访问数据库)

    如下:启动工程时,加载配置文件的用户名和密码到map中, 通过getKeyByValue()获取即可 import java.util.HashMap; import java.util.Map; i ...

  3. firewall-cmd.man

    FIREWALL-CMD(1) firewall-cmd FIREWALL-CMD(1) NAME firewall-cmd - firewalld command line client SYNOP ...

  4. bootstrap表格鼠标悬停与状态类

    今天在学习的过程中遇到在表格一章中 鼠标悬停 与 状态类 无效的问题,是因为在css文件中默认颜色都为白色造成的,解决方式如下: 在bootstrap文件夹中找到bootstrap.css文件(我的版 ...

  5. vue 关于solt得用法

    solt 第一种用法 父组件 <templateSolt></templateSolt> <templateSolt> <p>slot分发了内容< ...

  6. C#磁性窗体设计

    using System;using System.Collections.Generic;using System.ComponentModel;using System.Data;using Sy ...

  7. 入门项目 A5-3 interface-user 第三方接口3

    ''' 用户接口层 ''' # 导入数据库包下面的处理数据模,为了使用其内部名称空间 from db import db_handler # 注册接口函数,接收名字与密码两个参数 def regist ...

  8. webView 获取内容高度不准确的原因是因为你设置了某个属性

    不管是UIWebView 还是 WKWebView 这里 获取js属性 获取高度的方法 我就不一一细说了 ,本文最主要不说这个 ,网上有太多的方法 我最不摘几个 CGFloat webViewHeig ...

  9. 让MySQL数据库跑的更快的7个优化建议!

    随着容量和负载的增加,MySQL 的性能会日趋缓慢.这里有七点建议能够保证 MySQL 的平稳运行. 性能是我们衡量应用的一种方式,而应用性能的一项指标就是用户体验,也就是平时我们常说的:“用户需要等 ...

  10. PTA——洗牌

    PTA 7-43 Shuffling Machine #include<stdio.h> int main() { int i,n,*result; scanf("%d" ...