ihhh题解
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题解的更多相关文章
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
- 2016ACM青岛区域赛题解
A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- poj1399 hoj1037 Direct Visibility 题解 (宽搜)
http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...
- 网络流n题 题解
学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...
- CF100965C题解..
求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...
随机推荐
- Easy UI 面板
驾考园 http://www.jiakaoyuan.com 驾考园信息网 下载(源码)
- ****php:require_once(dirname(__FILE__)."/./config_uc.php");
Q:麻烦清楚地讲解一下这句的意思,具体路径是怎样的,这个文解在 根目录,如果我想放在根目录下的tieba文件夹里,应该怎么修改/./ 这个是表示什么? A: require_once(dirname( ...
- Go的基本示例
有空可以看看, 不知能不能超越JAVA的作法. hello.go package main import "fmt" func main() { s := "hello& ...
- Java程序员要求具备的10项技能
1.语法:必须比较熟悉,在写代码的时候IDE的编辑器对某一行报错应该能够根据报错信息知道是什么样的语法错误并且知道任何修正. 2.命令:必须熟悉JDK带的一些常用命令及其常用选项,命令至少需要熟悉:a ...
- 常用的Java代码汇总
1. 字符串有整型的相互转换 Java 1 2 <strong>Stringa=String.valueOf(2); //integer to numeric ...
- 攻城狮在路上(叁)Linux(十七)--- linux磁盘与文件管理概述
一.复习知识点: 1.扇区是最小的物理存储单位,大小为512bytes. 2.扇区组成一个圆,成为柱面,柱面是分区的最小单位. 3.第一个扇区很重要,因为包含了MBR(446字节)和分区表(64字节) ...
- annotation-config 和 component-scan 的区别
<context:annotation-config> 和 <context:component-scan>是Spring Core里面的两个基础概念,每个使用者都有必要理解怎 ...
- linux使用读写锁pthread_rwlock_t
转自:http://blog.csdn.net/onlyou930/article/details/6755593 使用读写锁 配置读写锁的属性之后,即可初始化读写锁.以下函数用于初始化或销毁读写锁. ...
- PMP 第四章 项目整合管理
1.什么是整合管理,整合什么?如何整合? 项目整合管理包括识别 定义 组合 统一与协调项目管理过组的个过程及项目管理活动二进行的各种过程和活动. 整合兼具统一 合并 连接和一体化的性质,对 ...
- 【SSM】拦截器的原理、实现
一.背景: 走过了双11,我们又迎来了黑色星期五,刚过了黑五,双12又将到来.不管剁手的没有剁手的,估计这次都要剁手了!虽然作为程序猿的我,没有钱但是我们长眼睛了,我们关注到的是我们天猫.淘宝.支付宝 ...