BZOJ 3207: 花神的嘲讽计划Ⅰ

3207: 花神的嘲讽计划Ⅰ

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 2145 Solved: 793
[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

20%的数据：1<=N<=1000;1<=M<=1000;
40%的数据：1<=N<=10000;1<=M<=10000;
100%的数据：1<=N<=100000;1<=M<=100000;y-x+1<=N;K<=y-x+1&K<=20;x<=y;

题中所有数据不超过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

原创 Memphis


[Submit][Status][Discuss]

主席树+HASH

因为区间长度固定，先处理出所有长度为K的HASH值，依次插入主席树里（也许你需要离散化一发），然后就是查询区间内是否出现过某个HASH值。

读入很大，快读可以优化不少（不加也能过就是了）。据说莫队也可以水过，然而我的莫队太low（还用STL map），并不能AC。

#include <cstdio>

inline char nextChar(void)
{
static const int siz = 1024;

static char buf[siz];
static char *hd = buf + siz;
static char *tl = buf + siz;

if (hd == tl)
fread(hd = buf, 1, siz, stdin);

return *hd++;
}

inline int nextInt(void)
{
register int ret = 0;
register int neg = false;
register int bit = nextChar();

for (; bit < 48; bit = nextChar())
if (bit == '-')neg ^= true;

for (; bit > 47; bit = nextChar())
ret = ret * 10 + bit - 48;

return neg ? -ret : ret;
}

typedef unsigned long long ull;

const int siz = 200005;
const int tre = 10000005;
const ull inf = 18446744073709551615UL;

int n, m, k;
ull hash[siz];

int tot;
int ls[tre];
int rs[tre];
int sum[tre];
int root[siz];

inline void insert(int &t, int f, ull p)
{
ull lt = 0, rt = inf, mid;

int a = t = ++tot, b = f;

sum[a] = sum[b] + 1;

while (lt != rt)
{
mid = (lt >> 1) + (rt >> 1) + (lt & rt & 1);

if (p <= mid)
rt = mid - 1, ls[a] = ++tot, rs[a] = rs[b], a = ls[a], b = ls[b];
else
lt = mid + 1, rs[a] = ++tot, ls[a] = ls[b], a = rs[a], b = rs[b];

sum[a] = sum[b] + 1;
}
}

inline bool query(int a, int b, ull p)
{
ull lt = 0, rt = inf, mid;

while (lt != rt)
{
mid = (lt >> 1) + (rt >> 1) + (lt & rt & 1);

if (p <= mid)
rt = mid - 1, a = ls[a], b = ls[b];
else
lt = mid + 1, a = rs[a], b = rs[b];
}

return sum[a] != sum[b];
}

signed main(void)
{
//    freopen("taunt.in", "r", stdin);
//    freopen("taunt.out", "w", stdout);

n = nextInt();
m = nextInt();
k = nextInt();

for (int i = 1; i <= n; ++i)
hash[i] = hash[i - 1] * 107UL + nextInt();

ull P = 1UL;

for (int i = 1; i <= k; ++i)
P = P * 107UL;

for (int i = k; i <= n; ++i)
insert(root[i], root[i - 1], hash[i] - hash[i - k] * P);

for (int i = 1, lt, rt; i <= m; ++i)
{
lt = nextInt();
rt = nextInt();

ull h = 0UL;

for (int j = 1; j <= k; ++j)
h = h * 107UL + nextInt();

puts(query(root[lt + k - 2], root[rt], h) ? "No" : "Yes");
}

//    return fclose(stdout), 0;
}

@Author: YouSiki