bzoj 3207 花神的嘲讽计划Ⅰ 主席树+hash

花神的嘲讽计划Ⅰ

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 3112 Solved: 1086
[Submit][Status][Discuss]

Description

背景

花神是神，一大癖好就是嘲讽大J，举例如下：

“哎你傻不傻的！【hqz：大笨J】”

“这道题又被J屎过了！！”

“J这程序怎么跑这么快！J要逆袭了！”

……

描述

这一天DJ在给吾等众蒟蒻讲题，花神在一边做题无聊，就跑到了一边跟吾等众蒟蒻一起听。以下是部分摘录：

“J你在讲什么！”

“我在讲XXX！”

“哎你傻不傻的！这么麻烦，直接XXX再XXX就好了！”

“……”

“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

原创 Memphis

题解：

　　主席树吧，两棵线段树形式相同，满足前缀和的形式是主席树的基本操作

　　这道题k确定，这不就随便hash就可以瞎搞了。

 #include<cmath>

 #include<cstdio>

 #include<cstring>

 #include<iostream>

 #include<algorithm>

 #define ull unsigned ll

 #define ll long long

 #define inf 18446744073709551615UL

 #define N 200007

 #define MX 10000007

 using namespace std;

 int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}

     return x*f;

 }

 int n,m,K,cnt;

 int val[N],b[N],rt[N];

 int sum[MX],a[MX][];

 ull H[N];

 void ins(int x,int &y,unsigned ll val)

 {

     y=++cnt;int tmp=y;

     sum[y]=sum[x]+;

     unsigned ll l=,r=inf,mid;

     while(l!=r)

     {

         mid=l/+r/;

         if((l&)&&(r&))mid++;

         int t=;

         if(val>mid)t=,l=mid+;

         else r=mid-;

         a[y][t]=++cnt;

         a[y][t^]=a[x][t^];

         x=a[x][t];y=a[y][t];sum[y]=sum[x]+;

     }

     y=tmp;

 }

 bool query(int x,int y,unsigned ll val)

 {

     unsigned ll l=,r=inf,mid;

     while(l!=r)

     {

         mid=l/+r/;

         if((l&)&&(r&))mid++;

         int t=;

         if(val>mid) t=,l=mid+;

         else r=mid-;

         x=a[x][t],y=a[y][t];

     }

     return (sum[y]-sum[x])!=;

 }

 int main()

 {

     n=read(),m=read(),K=read();

     for(int i=;i<=n;i++)

         val[i]=read();

     for(int i=;i<=n;i++)

         H[i]=H[i-]*+val[i];

     unsigned ll M=;

     for(int i=;i<=K;i++)M*=;

         for(int i=K;i<=n;i++)

             ins(rt[i-],rt[i],H[i]-H[i-K]*M);

     while(m--)

     {

         int x=read(),y=read();

         for(int i=;i<=K;i++)

             b[i]=read();

         unsigned ll now=;

         for(int i=;i<=K;i++)

             now=now*+b[i];

         if(query(rt[x+K-],rt[y],now))puts("No");

         else puts("Yes");

     }

 }

o(￣ヘ￣o＃)，hzw的代码。