Bzoj3207 花神的嘲讽计划Ⅰ

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2498  Solved: 903

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

看错tag啦!以为是可持久化trie树才点进来,想了半天觉得这和trie树没关系啊……用hash就行了嘛

百度一下果然是hash

主席树+HASH

discuss里说K<=20,那么$O(nk)$暴力找出所有区间的hash值,丢到主席树里,挨个回答询问即可。

WA了一串,以为是hash有问题,越发起劲地换模数

突然冷静一想,是不是主席树写错了?

仔细一看果然,下标错了一位。

好气呀

 1 /*by SilverN*/
 2 #include<iostream>
 3 #include<algorithm>
 4 #include<cstring>
 5 #include<cstdio>
 6 #include<cmath>
 7 #include<vector>
 8 #define LL unsigned long long
 9 using namespace std;
10 //const LL MX=998244353LL*157+1;
11 const LL MX=998244353LL*157*23333+7;
12 const int mxn=100010;
13 LL read(){
14     LL x=0;char ch=getchar();
15     while(ch<'0' || ch>'9'){ch=getchar();}
16     while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
17     return x;
18 }
19 struct node{
20     int l,r;
21     int val;
22 }t[mxn*90];
23 int rot[mxn],sz=0;
24 void insert(LL p,int v,LL l,LL r,int y,int &rt){
25     rt=++sz;
26     t[rt]=t[y];
27     t[rt].val+=v;
28     if(l==r)return;
29     LL mid=(l>>1)+(r>>1);
30     if((l&1)&&(r&1))mid++;
31     if(p<=mid)insert(p,v,l,mid,t[y].l,t[rt].l);
32     else insert(p,v,mid+1,r,t[y].r,t[rt].r);
33     return;
34 }
35 int query(LL val,LL l,LL r,int y,int rt){
36     if(!rt)return 0;
37     if(l==r){
38         if(t[rt].val-t[y].val>0)return 1;
39         else return 0;
40     }
41     LL mid=(l>>1)+(r>>1);
42     if((l&1)&&(r&1))mid++;
43     if(val<=mid)return query(val,l,mid,t[y].l,t[rt].l);
44     else return query(val,mid+1,r,t[y].r,t[rt].r);
45 }
46 int n,m,K;
47 LL a[mxn];
48 int main(){
49     int i,j;
50     n=read();m=read();K=read();
51     for(i=1;i<=n;i++)a[i]=read();
52     for(i=1;i<=n-K+1;i++){
53         LL res=0;
54         for(j=0;j<K;j++)
55             res=(res*113+a[i+j])%MX;
56         insert(res,1,0,MX,rot[i-1],rot[i]);
57     }
58     int x,y;
59     while(m--){
60         x=read();y=read();LL res=0;
61         for(j=1;j<=K;j++)res=(res*113+read())%MX;
62         y-=K-1;
63         if(x<=y && query(res,0,MX,rot[x-1],rot[y]))
64             printf("No
");
65         else printf("Yes
");
66     }
67     return 0;
68 }
原文地址:https://www.cnblogs.com/SilverNebula/p/6771242.html