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