【BZOJ】【P3578】【GTY的人类基因组计划2】【题解】【线段树套set or 线段树+set+hash】

传送门：www.lydsy.com/JudgeOnline/problem.php?id=3578

做为出题人当然要写题解啦……

数据结构题大水题……，题意不用我说了吧……注意1,2,3做了实验，1,2做还会产生2点点数，1,2,3,4做还会产生4点点数

判重是重点，蒟蒻已经在下方注明了：善用STL……，于是map<set<int>,bool> 即可……呃……

于是线段树维护区间，套平衡树，显然只有叶节点才需要用平衡树,区间做实验既是查询又是修改，需要打标记，然后传标记就好了……呃，真是大水题一个

数据全为随机生成，只有两个极限数据……不过暴力当然是不能过的了……haha

第一次在bzoj上出题，如果数据错了请联系我，我会认真改正

数据+标程：pan.baidu.com/s/1hqj3ZYO

===========================================

刚上bzoj 5min就被发现数据弱了……呃……在n=100000,m=2,q=100000的数据后标程T成翔，一看内存都吃了一个G了！！！，看来蒟蒻还是得好好想想了……

===========================================

作为出题人蒟蒻没想到会有这么多种做法，无限YM各位神犇，再列出两种做法:

1.

hash集合原来可以这样……OTZ

2.(这是目前Rank1 wwx神犇的做法)：

目前还没读懂，200+ms跪烂

===========================================

PS:hash选种子真TM蛋疼，过不了的可以尝试改改srand()的种子,或者……对着数据改……

Code：

/*
	ID:zky
*/
#include<map>
#include<set>
#include<cstdio>
#include<iostream>
#include<algorithm>
#define lson i<<1,l,mid
#define rson i<<1|1,mid+1,r
#define L i<<1
#define R i<<1|1
using namespace std;
const int maxn=1e5+10;
const int maxm=1e5+10;
const int maxq=1e5+10;
int n,m,q;
map<set<int>,int>M;
int home[maxn];
struct seg_tree{
	struct node{
		set<int>s;
		int size,lazy;
		node(){
			size=lazy=0;
		}
	};
	node t[maxn<<2];
	void deb(int i){
//		printf("#%d size:%d lazy:%d\n",i,t[i].size,t[i].lazy);
//		for(set<int>::iterator it=t[i].s.begin();it!=t[i].s.end();it++)
//			cout<<*it<<" ";
//		cout<<endl;
	}
	void pushdown(int i,int l,int r){
		if(t[i].lazy==0)return;
		int mid=l+r>>1;
		if(l==mid)
			M[t[L].s]=1;
		if(r==mid+1)
			M[t[R].s]=1;
		t[L].lazy=1;
		t[L].size=0;
		t[R].lazy=1;
		t[R].size=0;
		t[i].lazy=0;
	}
	void rz(int i){
		t[i].size=0;
//		if(!M.count(t[L].s))
		t[i].size+=t[L].size;
//		1if(!M.count(t[R].s))		
		t[i].size+=t[R].size;
	}
	void Change(int i,int l,int r,int pos,int val){
		deb(i);
		if(l==r){
			deb(i);
			t[i].s.insert(val);
			if(!M.count(t[i].s))
				t[i].size=t[i].s.size();
			else
				t[i].size=0;
			deb(i);
			return ;
		}
		int mid=l+r>>1;
		pushdown(i,l,r);
		if(pos<=mid)Change(lson,pos,val);
		else Change(rson,pos,val);
		rz(i);
		deb(i);
	}
	void Erase(int i,int l,int r,int pos,int val){
		deb(i);
		if(l==r){
			deb(i);
			t[i].s.erase(val);
			if(!M.count(t[i].s))
				t[i].size=t[i].s.size();
			else
				t[i].size=0;
			deb(i);
			return ;
		}
		int mid=l+r>>1;
		pushdown(i,l,r);
		if(pos<=mid)Erase(lson,pos,val);
		else Erase(rson,pos,val);
		rz(i);
		deb(i);
	}
	int qsum(int i,int l,int r,int l0,int r0){
		deb(i);
		if(l0<=l&&r0>=r){
			deb(i);
			M[t[i].s]=1;
			int s=t[i].size;
			t[i].size=0;
			t[i].lazy=1;
			deb(i);
			return s;
		}
		pushdown(i,l,r);
		int mid=l+r>>1;

		int ans=0;
		if(l0<=mid)ans+=qsum(lson,l0,r0);
		if(r0>mid) ans+=qsum(rson,l0,r0);
		rz(i);
		deb(i);
		return ans;
	}
}T;
void Change(int i,int j){
	if(home[i])
	T.Erase(1,1,n,home[i],i);
	T.Change(1,1,n,(home[i]=j),i);
}
int Qsum(int l,int r){
	return T.qsum(1,1,n,l,r);
}
int getint(){
	int res=0,ok=0;char ch;
	while(1){
		ch=getchar();
		if(ch<='9'&&ch>='0'){
			res*=10;res+=ch-'0';ok=1;
		}else if(ok)break;
	}return res;
}
int main(){
	n=getint();m=getint();q=getint();
	for(int i=1;i<=n;i++)Change(i,1);
	while(q--){
		char opt[3];scanf("%s",opt);
		if(opt[0]=='C'){
			int i=getint(),j=getint();
			Change(i,j);
		}else{
			int l=getint(),r=getint();
			printf("%d\n",Qsum(l,r));
		}
//		for(map<set<int>,int>::iterator it=M.begin();it!=M.end();it++)
//			for(set<int>::iterator i=it->first.begin();i!=it->first.end();i++)
//				cout<<*i<<" ";
//		cout<<endl;
	}
	return 0;
}

hash:

/*
	ID:zky
*/
#include<set>
#include<cctype>
#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#define lson i<<1,l,mid
#define rson i<<1|1,mid+1,r
#define L i<<1
#define R i<<1|1
using namespace std;
const int maxn=1e5+10;
typedef unsigned long long lld;
lld MOD=100000000000007LL;
set<lld>hash;
lld mp[maxn];
int room[maxn];
int n,m,q;
lld t[maxn<<2];
int size[maxn<<2];
int rsize[maxn<<2];
int lazy[maxn<<2];
struct seg_tree{

	seg_tree(){
		memset(t,0,sizeof(t));
		memset(size,0,sizeof(size));
		memset(lazy,0,sizeof(lazy));
		memset(rsize,0,sizeof(rsize));
	}
	void pushdown(int i,int l,int r){
		if(!lazy[i])return;
		int mid=l+r>>1;
		if(l==mid)
			hash.insert(t[L]);
		if(r==mid+1)
			hash.insert(t[R]);
		size[L]=0;
		size[R]=0;
		lazy[L]=1;
		lazy[R]=1;		
		lazy[i]=0;
	}
	void rz(int i){
		size[i]=size[L]+size[R];
	}
	void Erase(int i,int l,int r,int pos,int val){
		if(l==r){
			t[i]^=mp[val];
			rsize[i]--;
			if(!hash.count(t[i]))
				size[i]=rsize[i];
			else
				size[i]=0;
			return;
		}
		pushdown(i,l,r);
		int mid=l+r>>1;
		if(pos<=mid)Erase(lson,pos,val);
		else Erase(rson,pos,val);
		rz(i);
	}
	void Change(int i,int l,int r,int pos,int val){
		if(l==r){
			t[i]^=mp[val];
			rsize[i]++;
			if(!hash.count(t[i]))
				size[i]=rsize[i];
			else
				size[i]=0;
			return;
		}
		pushdown(i,l,r);
		int mid=l+r>>1;
		if(pos<=mid)Change(lson,pos,val);
		else Change(rson,pos,val);
		rz(i);
	}
	int qsum(int i,int l,int r,int l0,int r0){
		if(l0<=l&&r0>=r){
			hash.insert(t[i]);
			int s=size[i];
			size[i]=0;
			lazy[i]=1;
			return s;			
		}
		pushdown(i,l,r);
		int mid=l+r>>1;
		int ans=0;
		if(l0<=mid)ans+=qsum(lson,l0,r0);
		if(r0>mid) ans+=qsum(rson,l0,r0);
		rz(i);
		return ans;		
	}
}T;
void Change(int i,int j){
	if(room[i])
	T.Erase(1,1,n,room[i],i);
	T.Change(1,1,n,(room[i]=j),i);
}
int Qsum(int l,int r){
	return T.qsum(1,1,n,l,r);
}
int getint(){
	int res=0,ok=0;char ch;
	while(1){
		ch=getchar();
		if(ch<='9'&&ch>='0'){
			res*=10;res+=ch-'0';ok=1;
		}else if(ok)break;
	}return res;
}
int main(){
	srand(112);
//	freopen("input10.txt","r",stdin);
	n=getint();m=getint();q=getint();
	for(int i=1;i<=n;i++)mp[i]=(rand()*13131+21313*(rand()+756753))*10007;
	for(int i=1;i<=n;i++)Change(i,1);
	while(q--){
		char opt[3];scanf("%s",opt);
		if(opt[0]=='C'){
			int i=getint(),j=getint();
			Change(i,j);
		}else{
			int l=getint(),r=getint();
			printf("%d\n",Qsum(l,r));
		}
	}
	return 0;
}

3.Orz wwx

/*
	ID:zky
*/
#include<set>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1e5+10;
typedef unsigned long long lld;
set<lld>hash;
set<int>s;
lld mp[maxn];
lld room[maxn];
int size[maxn];
int bel[maxn];
int n,m,q;
int getint(){
	int res=0,ok=0;char ch;
	while(1){
		ch=getchar();
		if(ch<='9'&&ch>='0'){
			res*=10;res+=ch-'0';ok=1;
		}else if(ok)break;
	}return res;
}
void Change(int i,int j){
	room[bel[i]]^=mp[i];
	size[bel[i]]--;
	if(!hash.count(room[bel[i]]))
		s.insert(bel[i]);
	else
		s.erase(bel[i]);
	room[(bel[i]=j)]^=mp[i];
	size[bel[i]]++;
	if(!hash.count(room[bel[i]]))
		s.insert(bel[i]);
	else
		s.erase(bel[i]);	
}
int Qsum(int l,int r){
	int ans=0;
	
	for(;;){
		set<int>::iterator it=s.lower_bound(l);
		if(it==s.end()||*it>r)break;
		int deb=*it;
		ans+=size[*it];
		hash.insert(room[*it]);
		s.erase(it);
	}
	
	return ans;
}
int main(){
	srand(112);
	n=getint();m=getint();q=getint();
	for(int i=1;i<=n;i++)mp[i]=(rand()*13131+21313*(rand()+756753))*10007;
	size[1]=n;
	for(int i=1;i<=n;i++){
		room[1]^=mp[i];
		bel[i]=1;
	}
	s.insert(1);
	while(q--){
		char opt[3];scanf("%s",opt);
		if(opt[0]=='C'){
			int i=getint(),j=getint();
			Change(i,j);
		}else{
			int l=getint(),r=getint();
			printf("%d\n",Qsum(l,r));
		}
	}		
	return 0;
}