10分做法:

由于空间卡得紧,所以给了10分暴力分0.0
所以大家很容易就知道暴力就是线段树套ac自动机辣
时间:$O((\sum |qSi| + \sum |nSi|)*log Q)$
空间:$O((\sum |nSi|)*log n)$

40-50分做法:

发现我们可以用莫队做,可是有删除操作不好搞。
发现其实删除操作就是在trie上把一个串的trie节点的size统统删去1,然后用跑的时候这个不更新size=0的终止节点就行辣,由于fail的性质,可以证明这样做是正确的。于是维护size我们用lct。
虽然理论复杂度和满分做法一样,可是区间修改的splay十分蛋疼,你懂的。
我比较良心因此数据给了一些分给此做法~因此期望分就是40~50分辣。
时间:$O(N^{1.5}log(max\{|nSi|\})+Qlog(max\{|qSi|\}))$ 
空间:$O(\sum |nSi| + \sum |qSi|)$

满分做法:

其实发现我们可以在终止节点打个标记= =所以lct从维护trie变成维护fail树,然后就能通过全部数据辣...
时间:$O(N^{1.5}log(max\{|nSi|\})+Qlog(max\{|qSi|\}))$
空间:$O(\sum |nSi| + \sum |qSi|)$
维护树还能用树链剖分+线段树,这种做法虽然在修改上多了一个log,但是由于常熟十分小,期望分数还是100分的,只要你没写残...
时间:$O(N^{1.5}log(max\{|nSi|\})+Qlog^2(max\{|qSi|\}))$
空间:$O(\sum |nSi| log(\sum |nSi|) + \sum |qSi|)$
妈呀你们做题才发现这题能被分块水掉............
完了完了我的人生啊.............

完了完了saffah大爷说了一种秒标程的做法:

同样是建线段树,然后我们查询一个区间的时候先标记,总共会标记$O(QlogN)$个节点,然后我们对每个标记的节点再建ac自动机,查询完一个就rebuild自动机....这样空间就不会爆了......orz
这样每个串最多被build到$O(log(N))$个自动机内,因此查询复杂度为$O(log(N) \sum |nSi|)$

乱搞?

1、暴力修改终止节点:然而fail树高度最坏$O(\max\{S\})$。(咦好像数据并没有出全a的数据,没水过的估计被卡常了?) 、、
2、分块大发:似乎这是正解...然而还是被卡常了.....

PS:

复杂度和最长串长度有关?然而好像我造的随机数据并没有太长QAQ

满分第一种做法(ac自动机+莫队+lct):

#include <bits/stdc++.h>
using namespace std;
const int N=1e5+10, nQ=1e5+10, nSz=26, S=5*1e5+10;
struct node *null, *root;
struct node {
node *c[2], *f, *ch[nSz], *fail;
int n, mx, tag;
void init() { c[0]=c[1]=f=null; fail=root; for(int i=0; i<nSz; ++i) ch[i]=root; mx=n=tag=0; }
bool d() { return f->c[1]==this; }
bool isson() { return f->c[0]==this || f->c[1]==this; }
void setc(node *x, bool d) { c[d]=x; x->f=this; }
void up() { mx=max(c[0]->mx, c[1]->mx); if(tag) mx=max(mx, n); }
}Po[S], *iT=Po, *pos[N], *s[S];
node *newnode() { iT->init(); return iT++; }
void rot(node *x) {
node *f=x->f; bool d=x->d();
if(f->isson()) f->f->setc(x, f->d());
else x->f=f->f;
f->setc(x->c[!d], d);
x->setc(f, !d);
f->up();
}
void splay(node *x) {
while(x->isson())
if(!x->f->isson()) rot(x);
else x->d()==x->f->d()?(rot(x->f), rot(x)):(rot(x), rot(x));
x->up();
}
node *access(node *x) { node *y=null; for(; x!=null; y=x, x=x->f) splay(x), x->c[1]=y; return y; }
void add(char *s, int id) {
int n=strlen(s); node *now=root;
for(int i=0; i<n; ++i) {
int x=s[i]-'a';
if(now->ch[x]==root) now->ch[x]=newnode();
now=now->ch[x];
}
now->n=n; pos[id]=now;
}
void bfs() {
int fr=0, ta=0;
s[ta++]=root; node *x;
while(fr!=ta) {
x=s[fr++];
for(int i=0; i<nSz; ++i)
if(x->ch[i]!=root) {
s[ta++]=x->ch[i];
if(x==root) continue;
x->ch[i]->fail=x->fail->ch[i];
x->ch[i]->f=x->ch[i]->fail;
}
else x->ch[i]=x->fail->ch[i];
}
}
int getans(char *s, int n) {
int ans=0; node *now=root;
for(int i=0; i<n; ++i) {
int x=s[i]-'a';
now=now->ch[x];
access(now); splay(now); ans=max(ans, now->mx);
}
return ans;
}
void update(int x, int f) { splay(pos[x]); pos[x]->tag+=f; pos[x]->up(); }
int blc[nQ], n, Q, nowid=0;
char a[S];
struct qq { int l, r, id, n, pos; }q[nQ];
bool cmp1(const qq &a, const qq &b) { return blc[a.l]==blc[b.l]?a.r<b.r:a.l<b.l; }
bool cmp2(const qq &a, const qq &b) { return a.id<b.id; }
void init() {
null=iT++; root=iT++; null->init(); root->init();
for(int i=0, sq=sqrt(n+0.5); i<n; ++i) blc[i]=i/sq;
}
int main() {
scanf("%d%d", &n, &Q);
init();
for(int i=0; i<n; ++i) { scanf("%s", a); add(a, i); }
for(int i=0; i<Q; ++i) { int l, r; scanf("%d%d%s", &l, &r, a+nowid); --l; --r; q[i]=(qq){l, r, i, (int)strlen(a+nowid), nowid}; nowid+=q[i].n; }
bfs();
sort(q, q+Q, cmp1);
int l=0, r=-1;
for(int i=0; i<Q; ++i) {
int nl=q[i].l, nr=q[i].r;
while(nl<l) update(--l, 1);
while(nl>l) update(l++, -1);
while(nr<r) update(r--, -1);
while(nr>r) update(++r, 1);
blc[q[i].id]=getans(a+q[i].pos, q[i].n);
}
for(int i=0; i<Q; ++i) printf("%d\n", blc[i]);
return 0;
}

  

满分第二种做法(ac自动机+莫队+树剖+线段树):

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm> using namespace std; struct edge {
int t;
edge *next;
};
struct qry {
char* str;
int len, l, r, bl, i;
};
struct seg {
int l, r, s, c, v;
seg *ls, *rs;
};
inline bool cmpQry(const qry& a, const qry& b) {
if (a. bl == b. bl)
return a. r > b. r;
else
return a. bl < b. bl;
} const int maxn = 100009;
const int maxlen = 500009; int n, m, bsz, ans[maxn];
int tr[maxlen][27], fl[maxlen], ep[maxn], tn, dep[maxlen];
int d[maxlen], sz[maxlen], tc, fc[maxlen], ch[maxlen], cl[maxlen];
char g[maxlen], cbuf_arr[maxlen], *cbuf(cbuf_arr);
qry q[maxn];
edge ebuf_arr[maxlen << 1], *ebuf(ebuf_arr), *head[maxlen];
seg sbuf_arr[maxn << 2], *sbuf(sbuf_arr), *rt[maxlen]; int trieIns(char* a) {
int p(1);
for (; *a; ++ a) {
int d(*a - 97);
if (!tr[p][d])
dep[tr[p][d] = ++ tn] = dep[p] + 1;
p = tr[p][d];
}
return p;
}
inline void addEdge(int u, int v) {
ebuf-> t = v;
ebuf-> next = head[u];
head[u] = ebuf ++;
} #define cpos(_p_) (d[_p_]-d[ch[fc[_p_]]])
#define midp(_p_) ((_p_->l+_p_->r)>>1)
inline seg* sgtBuild(int l, int r) {
seg* p(sbuf ++);
p-> l = l, p-> r = r;
p-> s = 0, p-> c = 0;
if (l + 1 < r) {
p-> ls = sgtBuild(l, midp(p));
p-> rs = sgtBuild(midp(p), r);
}
return p;
}
inline void sgtChg(seg* p, int po, int v) {
if (p-> l + 1 == p-> r) {
p-> c += v;
if (p-> c)
p-> s = p-> v;
else
p-> s = 0;
}
else {
if (po < midp(p))
sgtChg(p-> ls, po, v);
else
sgtChg(p-> rs, po, v);
p-> s = max(p-> ls-> s, p-> rs-> s);
}
}
inline void sgtSet(seg* p, int po, int vo) {
if (p-> l + 1 == p-> r)
p-> v = vo;
else if (po < midp(p))
sgtSet(p-> ls, po, vo);
else
sgtSet(p-> rs, po, vo);
}
inline int sgtQry(seg* p, int po) {
if (p-> r == po)
return p-> s;
else if (po <= midp(p))
return sgtQry(p-> ls, po);
else
return max(p-> ls-> s, sgtQry(p-> rs, po));
} void acaBuild() {
static int q[maxlen];
int hd(0), tl(1);
fl[q[0] = 1] = 0;
d[1] = 1;
while (hd < tl) {
int u(q[hd ++]), v;
for (int i = 0; i < 26; ++ i)
if ((v = tr[u][i])) {
int w;
for (w = fl[u]; w && !tr[w][i]; w = fl[w]);
if (!w)
fl[v] = 1;
else
fl[v] = tr[w][i];
d[v] = d[fl[v]] + 1;
addEdge(fl[v], v);
q[tl ++] = v;
}
}
tc = 0;
for (int i = tl - 1; i >= 0; -- i) {
int p(q[i]), z(-1);
sz[p] = 1;
for (edge* e = head[p]; e; e = e-> next) {
sz[p] += sz[e-> t];
if (z == -1 || sz[e-> t] > sz[z])
z = e-> t;
}
if (z == -1)
cl[fc[p] = ++ tc] = 1;
else
++ cl[fc[p] = fc[z]];
ch[fc[p]] = p;
}
for (int i = 1; i <= tc; ++ i)
rt[i] = sgtBuild(0, cl[i]);
for (int i = 1; i <= tn; ++ i)
sgtSet(rt[fc[i]], cpos(i), dep[i]);
} int getTrans(int p, int w) {
for (; p && !tr[p][w]; p = fl[p]);
if (p)
return tr[p][w];
else
return 1;
}
void cptMo() {
sort(q, q + m, cmpQry);
int l(1), r(1);
sgtChg(rt[fc[ep[1]]], cpos(ep[1]), 1);
for (int i = 0; i < m; ++ i) {
while (r < q[i]. r) {
++ r;
sgtChg(rt[fc[ep[r]]], cpos(ep[r]), 1);
}
while (l > q[i]. l) {
-- l;
sgtChg(rt[fc[ep[l]]], cpos(ep[l]), 1);
}
while (r > q[i]. r) {
sgtChg(rt[fc[ep[r]]], cpos(ep[r]), -1);
-- r;
}
while (l < q[i]. l) {
sgtChg(rt[fc[ep[l]]], cpos(ep[l]), -1);
++ l;
}
int g(0);
for (int p = 1, j = 0; j < q[i]. len; ++ j) {
p = getTrans(p, q[i]. str[j] - 97);
for (int u = p; u; u = fl[ch[fc[u]]])
g = max(g, sgtQry(rt[fc[u]], cpos(u) + 1));
}
ans[q[i]. i] = g;
}
} int main() {
scanf("%d%d", &n, &m);
bsz = (int)pow(n, 0.5) + 1;
dep[tn = 1] = 0;
for (int i = 1; i <= n; ++ i) {
scanf("%s", g);
ep[i] = trieIns(g);
}
acaBuild();
for (int i = 0; i < m; ++ i) {
scanf("%d%d%s", &q[i]. l, &q[i]. r, cbuf);
q[i]. len = strlen(q[i]. str = cbuf);
cbuf += q[i]. len;
q[i]. bl = q[i]. l / bsz;
q[i]. i = i;
}
cptMo();
for (int i = 0; i < m; ++ i)
printf("%d\n", ans[i]);
}

  

满分第三种做法(还没写...):

ihhh题解的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

随机推荐

  1. hdu 1754:I Hate It(线段树,入门题,RMQ问题)

    I Hate It Time Limit: 9000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total S ...

  2. 蛋疼的Fedora17

    在公司给同事要安装ASM oracle 11g,说要在forder17上安装,于是乎我先在自己的虚拟机上安装了一个forder17,遇到了几个很蛋疼的问题:      1. 安装的时候没有创建普通用户 ...

  3. [Java] 使用Java Visual VM寻找PermGen Space的解决办法

    在Eclipse使用tomcat运行3个项目时,老是报这个错误,以下为错误详情: 2014-5-28 13:47:41 org.apache.catalina.core.StandardWrapper ...

  4. WebStorm中将Project分享到GitHub时报“Error Running Git”错误的解决办法

    错误信息 Cannot run program "git.exe":CreateProcess error=2,系统找不到指定的文件. 解决办法 从错误信息就可以知道,WebSto ...

  5. 智能车学习(十四)——K60单片机GPIO学习

    一.头文件: #ifndef __MK60_GPIO_H__ #define __MK60_GPIO_H__ #include "MK60_gpio_cfg.h" /* * 定义管 ...

  6. python特殊函数 __len__(self):

    __len__ 如果一个类表现得像一个list,要获取有多少个元素,就得用 len() 函数. 要让 len() 函数工作正常,类必须提供一个特殊方法__len__(),它返回元素的个数. 例如,我们 ...

  7. Liferay 6.2 改造系列之十一:默认关闭CDN动态资源

    在行业客户中,一般无法提供CDN服务,因此默认关闭CDN动态资源功能: 在/portal-master/portal-impl/src/portal.properties文件中,有如下配置: # # ...

  8. 在C#程序中实现插件架构

    阅读提示:这篇文章将讲述如何利用C#奇妙的特性,实现插件架构,用插件(plug-ins)机制建立可扩展的解决方案. 在.NET框架下的C#语言,和其他.NET语言一样提供了很多强大的特性和机制.其中一 ...

  9. java的重载、覆盖和隐藏的区别

    重载:方法名相同,但参数不同的多个同名函数 注意:1.参数不同的意思是参数类型.参数个数.参数顺序至少有一个不同 2.返回值和异常以及访问修饰符,不能作为重载的条件(因为对于匿名调用,会出现歧义,eg ...

  10. Color the ball(线段树)

    Time Limit: 9000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission( ...