BZOJ 3207: 花神的嘲讽计划Ⅰ( hash + 可持久化线段树 )

BZOJ 3207: 花神的嘲讽计划Ⅰ( hash + 可持久化线段树 )

O(NK)暴力搞出所有子串的哈希值, 然后就对哈希值离散化建权值线段树, 就是主席树的经典做法了.总时间复杂度O(NK+(N+Q)logN)

--------------------------------------------------------------------------

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cctype>
using namespace std;
typedef unsigned long long ull;
const int maxn = 100009;
const ull P = 1000000007ULL;
inline int read() {
char c = getchar();
int ret = 0;
for(; !isdigit(c); c = getchar());
for(; isdigit(c); c = getchar()) 
ret = ret * 10 + c - '0';
return ret;
}
int seq[maxn], r[maxn];
int N, Q, K, n;
ull H[maxn], Hash[maxn];
void Init() {
N = read(); Q = read(); K = read();
for(int i = 1; i <= N; i++)
seq[i] = read();
n = N - K + 1;
for(int i = 1; i <= n; i++) {
H[i] = 0;
for(int j = 0; j < K; j++)
(H[i] *= P) += seq[i + j];
Hash[i] = H[i];
}
sort(H + 1, H + n + 1);
for(int i = 1; i <= n; i++)
r[i] = lower_bound(H + 1, H + n + 1, Hash[i]) - H;
}
struct Node {
Node *l, *r;
int cnt;
} pool[5000000], *pt, *Root[maxn], *Null;
void Init_sgt() {
pt = pool;
Null = pt++;
Null->l = Null->r = Null;
}
int Val;
Node* Modify(Node* t, int l, int r) {
Node* h = pt++;
if(l == r)
h->cnt = t->cnt + 1;
else {
int m = (l + r) >> 1;
if(Val <= m) {
h->l = Modify(t->l, l, m);
h->r = t->r;
} else {
h->l = t->l;
h->r = Modify(t->r, m + 1, r);
}
}
return h;
}
int Query(Node* t, int l, int r) {
if(l == r)
return t->cnt;
int m = (l + r) >> 1;
if(Val <= m)
return Query(t->l, l, m);
return Query(t->r, m + 1, r);
}
void Solve() {
Init_sgt();
Root[0] = Null;
for(int i = 1; i <= n; i++) {
Val = r[i];
Root[i] = Modify(Root[i - 1], 1, n);
}
while(Q--) {
int l = read(), r = read() - K + 1;
ull hash_val = 0;
for(int i = 0; i < K; i++)
(hash_val *= P) += read();
Val = lower_bound(H + 1, H + n + 1, hash_val) - H;
if(l > r || H[Val] != hash_val) {
puts("Yes");
continue;
}
puts(Query(Root[r], 1, n) > Query(Root[l - 1], 1, n) ? "No" : "Yes");
}
}
int main() {
Init();
Solve();
return 0;
}

--------------------------------------------------------------------------

3207: 花神的嘲讽计划Ⅰ

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1132  Solved: 412
[Submit][Status][Discuss]

Description

背景
花神是神,一大癖好就是嘲讽大J,举例如下:
“哎你傻不傻的!【hqz:大笨J】”
“这道题又被J屎过了!!”
“J这程序怎么跑这么快!J要逆袭了!”
……
描述
这一天DJ在给吾等众蒟蒻讲题,花神在一边做题无聊,就跑到了一边跟吾等众蒟蒻一起听。以下是部分摘录:
1.
“J你在讲什么!”
“我在讲XXX!”
“哎你傻不傻的!这么麻烦,直接XXX再XXX就好了!”
“……”
2.
“J你XXX讲过了没?”
“……”
“那个都不讲你就讲这个了?哎你傻不傻的!”
“……”
DJ对这种情景表示非常无语,每每出现这种情况,DJ都是非常尴尬的。
经过众蒟蒻研究,DJ在讲课之前会有一个长度为N方案,我们可以把它看作一个数列;
同样,花神在听课之前也会有一个嘲讽方案,有M个,每次会在x到y的这段时间开始嘲讽,为了减少题目难度,每次嘲讽方案的长度是一定的,为K。
花神嘲讽DJ让DJ尴尬需要的条件:
在x~y的时间内DJ没有讲到花神的嘲讽方案,即J的讲课方案中的x~y没有花神的嘲讽方案【这样花神会嘲讽J不会所以不讲】。
经过众蒟蒻努力,在一次讲课之前得到了花神嘲讽的各次方案,DJ得知了这个消息以后欣喜不已,DJ想知道花神的每次嘲讽是否会让DJ尴尬【说不出话来】。

Input

第1行3个数N,M,K;
第2行N个数,意义如上;
第3行到第3+M-1行,每行K+2个数,前两个数为x,y,然后K个数,意义如上;

Output

对于每一个嘲讽做出一个回答会尴尬输出‘Yes’,否则输出‘No’

Sample Input

8 5 3
1 2 3 4 5 6 7 8
2 5 2 3 4
1 8 3 2 1
5 7 4 5 6
2 5 1 2 3
1 7 3 4 5

Sample Output

No
Yes
Yes
Yes
No

HINT

题中所有数据不超过2*10^9;保证方案序列的每个数字<=N

2~5中有2 3 4的方案,输出No,表示DJ不会尴尬

1~8中没有3 2 1的方案,输出Yes,表示DJ会尴尬

5~7中没有4 5 6的方案,输出Yes,表示DJ会尴尬

2~5中没有1 2 3的方案,输出Yes,表示DJ会尴尬

1~7中有3 4 5的方案,输出No,表示DJ不会尴尬

Source

上一篇:java 框架 面试


下一篇:delphi TTBXToolBar 增加外部控件