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

最新推荐文章于 2020-09-19 22:44:34 发布
最新推荐文章于 2020-09-19 22:44:34 发布
阅读量2.5k 收藏
点赞数
分类专栏： OI 一入OI深似海 文章标签： 省选 bzoj
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/iamzky/article/details/25807361
版权

传送门：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;
}


iamzky
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
bzoj 3730. 震波（动态点分治 + vector树状数组）
01-06
因为子树树形被破坏，在做减法时，子节点子树对父节点的贡献不能用子树维护的信息 + 连向父节点的边的贡献得到，考虑在每个节点再维护一个线段树，按距离建树用来统计子树内的节点到其父节点距离为 p 的点权和，这样...
数据结构 树套树 线段树套set
fdxgcw的博客
07-20 129
原题 #include <iostream> #include <queue> #include <cstring> #include <set> #include <algorithm> using namespace std; int n, m; const int N = 50010, M = N * 4, INF = 1e9; struct Tree { int l, r; multiset<int>s; } tr[M];
树套树——线段树套set
行走天涯的豆沙包的博客
09-19 315
请你写出一种数据结构，来维护一个长度为 n 的序列，其中需要提供以下操作： 1 pos x，将 pos 位置的数修改为 x。 2 l r x，查询整数 x 在区间 [l,r] 内的前驱(前驱定义为小于 x，且最大的数)。 数列中的位置从左到右依次标号为 1∼n。 区间 [l,r] 表示从位置 l 到位置 r 之间（包括两端点）的所有数字。 区间内排名为 k 的值指区间内从小到大排在第 k 位的数值。（位次从 1 开始） 输入格式 第一行包含两个整数 n,m，表示数列长度以及操作次数。 第二行包含 n 个整数
Codeforces Round #397 F. Souvenir线段树套set
now_ing的博客
10-28 231
博客：http://www.cnblogs.com/qscqesze/p/6405272.html这算是一个时间换空间的做法。 对于每一个树上的点弄一个set。然后离线从右到左， 每次更新的时候如果 if((p==t.end()||(p-val)>=mi)&&(p==t.begin()||(val-(–p))>=mi)) 这样就不用单个点地查询了。。
BZOJ3578 GTY的人类基因组计划2
玫葵之蝶
11-15 221
Link Difficulty 算法难度4，思维难度6，代码难度4 Description 一开始给定nnn个元素，都位于第一个集合中 有mmm个集合，编号1~m qqq次操作，分两种类型： 将一个元素移动到另一个集合里 询问编号在[l,r][l,r][l,r]之间的以前没有询问过的集合的大小之和 1≤n,m,q≤1051\le n,m,q\le 10^51≤n,m,q≤105 Solutio...
BZOJ3578:GTY的人类基因组计划2(集合hash,STL)
weixin_33828101的博客
02-17 93
Description GTY召唤了n个人来做实验，GTY家的房子很大，有m个房间一开始所有人都在1号房间里，GTY会命令某人去某个房间等待做实验，或者命令一段区间的房间开始实验，实验会获得一些实验信息点数，点数为房间里的人数，如果一个房间里的一群人已经做过实验了那么这些人将不会增加实验信息点数(不会增加是针对这一群人的，不是对这群人中的每个人，即1,2,3做了实验，1,2再做实验还会增...
BZOJ泡泡堂题解.doc
09-29
题解 , 文档 , 资料 BZOJ 泡泡堂
wangjunrui666#bzoj-problem##5165. 树上倍增1
07-25
您需要写一个树上倍增算法，以实现如下操作：A x 新建一个节点，将它作为x节点的儿子，编号为当前节点总数+1。保证询问不超过1000次输出格式对于每一个Q操作，
数论题解：「BZOJ1053」反素数/「Violet5」樱花
01-22
「BZOJ1053」反素数/「Violet5」樱花 详细题解
CTSC 2011 无穷图的桥(BZOJ 2307) 题解.ppt
03-17
CTSC 2011 无穷图的桥(BZOJ 2307) 题解.ppt
BZOJ 3578: GTY的人类基因组计划2 STL
BlackJack
08-29 350
3578: GTY的人类基因组计划2 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 485  Solved: 206 [Submit][Status][Discuss] Description   GTY召唤了n个人来做实验，GTY家的房子很大，有m个房间一开始所有人都在1号房间里，GTY会命令某人去某个房间等待做实验，或者命令一段
BZOJ 3578 GTY的人类基因组计划2 set
(╯°口°)╯(┴—┴
01-13 1190
题目大意：有一些人一开始都在一号房间。有一些指令使他们区别的房间。每次选定一段区间作试验，如果有的房间的组合已经做过试验，那么就不计入最后的试验点数，否则计入试验点数为试验的人数。问最后会获得多少试验点数。 思路：对于每一个人随即一个long long，一个集合所代表的数字就是这些long long的异或和。然后用set来暴力判重。 CODE： #include #
P2087 GTY的人类基因组计划2 set map
weixin_30823227的博客
05-23 69
　　 题目描述 GTY召唤了n个人来做实验，GTY家的房子很大，有m个房间一开始所有人都在1号房间里，GTY会命令某人去某个房间等待做实验，或者命令一段区间的房间开始实验，实验会获得一些实验信息点数，点数为房间里的人数，如果一个房间里的一群人已经做过实验了那么这些人将不会增加实验信息点数(不会增加是针对这一群人的，不是对这群人中的每个人，即1,2,3做了实验，1,2再做实验还会增加2点实验点数...
bzoj 3578: GTY的人类基因组计划2 (STL+hash)
clover_hxy的博客
07-13 406
3578: GTY的人类基因组计划2 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 312  Solved: 139 [Submit][Status][Discuss] Description   GTY召唤了n个人来做实验，GTY家的房子很大，有m个房间一开始所有人都在1号房间里，GTY会命令某人去某个房间等待做实验，或者命令一段
CF 19D, 线段树套set/KD-tree
weixin_30556959的博客
07-28 104
题目大意：有n个操作，要么在平面上增删一个点，要么询问某个坐标右上侧（严格大于等于）第一个点是谁（x相同取y最小）。 解：非常智障看错题了，没看到要在右上，所以写了个树套set，发现看错题嘀咕了一下艹这询问set搞不来，想了一下似乎kd-tree可以做，就写了棵，常数卡得很死过了。后来看题解还真是树套set。离散化后每个线段都是一个set，然后插入删除的时候沿途把set搞一下，询问的时候用个小技...
BZOJ3578 GTY的人类基因组计划2 HASH
a18700013354的博客
02-28 103
题意：有N个人M个房间，初始时所有人都在1号房间，维护：1、让某个人去某个房间 2、假如[l,r]中所有的人还未一起统计过，Ans+=人数。输出Ans 题解： 集合HASH get 根据异或的性质，a^x^x=a，因此一个集合的HASH值就是其所有元素的异或起来的值 接下来就是STL大法好 将每个房间看成一个集合，每个人赋一个初始rand，set存储当前所有HASH值没有出现...
bzoj4829 [TJOI2017]dna（后缀数组）
Icefox的博客
07-23 547
luogu3763有题面可以去看的说。。大牛分站开o2优化了（没错我就是这样才过的。）后缀数组处理一下lcp，然后枚举三次就好了。bzoj上我的总是T。嗯，一定是姿势不对。。不过我看神犇们都是用SAM写的，贼快。我也要去学习一下。。 #include #include #define N 200005 char s[N]; int tst,n,len1,len2,m,st[N][19],lo
（线段树+加map）1287 加农炮
black_horse2018的博客
05-11 152
1287 加农炮 1 秒 131,072 KB 40 分 4 级题 一个长度为M的正整数数组A，表示从左向右的地形高度。测试一种加农炮，炮弹平行于地面从左向右飞行，高度为H，如果某处地形的高度大于等于炮弹飞行的高度H（A[i] >= H），炮弹会被挡住并落在i - 1处，则A[i - 1] + 1。如果H <= A[0]，则这个炮弹无效，如果H > 所有...
【BZOJ 4892】DNA
C20181503csy的博客
10-07 216
又一个DNA 【题意】 加里敦大学的生物研究所，发现了决定人喜不喜欢吃藕的基因序列S，有这个序列的碱基序列就会表现出喜欢吃藕的性状，但是研究人员发现对碱基序列S，任意修改其中不超过3个碱基，依然能够表现出吃藕的性状。现在研究人员想知道这个基因在DNA链S0上的位置。所以你需要统计在一个表现出吃藕性状的人的DNA序列S0上，有多少个连续子 串可能是该基因，即有多少个S0的连续子串修改小于等于三个字母...
bzoj 1040 zjoi2008 骑士 树形dp
最新发布
08-08
我们可以使用树形dp来解决这个问题，具体来说，我们可以将这个图看作是一棵树，每个点的父节点是它的前驱或者后继，然后我们从根节点开始，依次向下遍历，对于每个节点，我们可以考虑它的两个儿子，如果它的两个儿子...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值