BZOJ 3207 花神的嘲讽计划Ⅰ 哈希+主席树

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不会尴尬



每K个数字分别哈希出一个值,这样子我们把每K个数的哈希值看作是末尾的哈希值,
意思就是全部寄存到末尾。
这样我们就得到了新的K~n的值。

然后我们对这些巨大的哈希值离散化,得到了范围在K~n(或者1~n-K+1)的n-K+1个值。
对于它们的询问,首先当然是求出询问的序列的哈希值;
然后我们只要判断是否存在相应的哈希值,在这段区间里面出现即可。
可以用vector维护直接水水掉,当然我觉得主席树比较ok。
我们只要判断这个位置上的值是不是次数>1即可。
显然建立一个主席树轻松解决了。

主要存在的问题就是如何得出询问时,给出的序列的哈希值后,求出对应的离散化的值。
当然了,我们可以在排好序的哈希值数组里面进行二分。
还是有点麻烦的样子……但是其实不难,
一开始用了map后悔死。

被哈希卡了好久……然后就换成了ull自然溢出。
结果?随便提交了一发竟然AC了还rank44?
问题是……
我刚和网上一个程序拍出来错误样例啊!
哈希这东西,还真的不好说的。



#include<bits/stdc++.h>
#define llt unsigned long long
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0' || ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int 
	N=100005,
	KK=25;  //logN=KK
const llt
	NUM=(llt)107;
int n,m,K,Tcnt,root[N],rank[N];
llt tmp,a[N],t[N];
struct Array{
	llt val;
	int id;
}b[N];
bool cmp(Array x,Array y){
	return x.val<y.val;
}
struct ChairTree{
	int l,r,num;
}ct[N*KK];
int BS(llt x){
	int L=K,R=n,mid;
	while (L<=R){
		mid=(L+R)>>1;
		if (b[mid].val==x) return mid;
		if (b[mid].val<x) L=mid+1;
			else R=mid-1;
	}
	return 0;
}
void insert(int L,int R,int &x,int val){
	ct[Tcnt++]=ct[x];
	x=Tcnt-1;
	ct[x].num++;
	if (L==R) return;
	int mid=(L+R)>>1;
	if (val<=mid) insert(L,mid,ct[x].l,val);
		else insert(mid+1,R,ct[x].r,val);
}
bool query(int L,int R,int ll,int rr,int val){
	if (L==R) return (ct[rr].num-ct[ll].num);
	int mid=(L+R)>>1;
	if (val<=mid) return query(L,mid,ct[ll].l,ct[rr].l,val);
		else return query(mid+1,R,ct[ll].r,ct[rr].r,val);
}
void HASH(int nn){
	t[0]=(llt)0;
	for (int i=1;i<=nn;i++)
		t[i]=t[i-1]*NUM+a[i];	
}
void disc(){
	for (int i=K;i<=n;i++)
		b[i].val=t[i]-t[i-K]*tmp,
		b[i].id=i;
	sort(b+K,b+n+1,cmp);
	int rk=K-1; b[K-1].val=(llt)-1;
	for (int i=K;i<=n;i++)
		rk+=(b[i].val!=b[i-1].val),
		rank[b[i].id]=rk;
}
int main(){
	n=read(),m=read(),K=read();
	for (int i=1;i<=n;i++) a[i]=(llt)read();
	
	tmp=(llt)1; 
	for (int i=1;i<=K;i++) tmp=tmp*NUM;
	HASH(n),disc();
	Tcnt=1,root[K-1]=0;
	for (int i=K;i<=n;i++)
		root[i]=root[i-1],
		insert(K,n,root[i],rank[i]);
		
	int x,y,tt;
	while (m--){
		x=read(),y=read();
		for (int i=1;i<=K;i++) a[i]=(llt)read();
		HASH(K),tt=BS(t[K]);
		if (!tt || !query(K,n,root[x+K-2],root[y],rank[b[tt].id]))
			puts("Yes"); else puts("No");
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值