皮卡丘的梦想（线段树+位运算）

皮卡丘的梦想2

1000 ms 65536 KiB

Submit Status My Status Origin

Description

一天，一只住在 501 实验室的皮卡丘决定发奋学习，成为像 LeiQ 一样的巨巨，于是他向镇上的贤者金桔请教如何才能进化成一只雷丘。

金桔告诉他需要进化石才能进化，并给了他一个地图，地图上有 n 个小镇，他需要从这些小镇中收集进化石。

接下来他会进行 q 次操作，可能是打听进化石的信息，也可能是向你询问第 l 个小镇到第 r 个小镇之间的进化石种类。

如果是打听信息，则皮卡丘会得到一个小镇的进化石变化信息，可能是引入了新的进化石，也可能是失去了全部的某种进化石。

如果是向你询问，你需要回答他第 l 个小镇到第 r 个小镇之间的进化石种类。

Input

首先输入一个整数 T (1 <= T <= 10)，代表有 T 组数据。

每组数据的第一行输入一个整数 n (1 <= n <= 100000) 和一个整数 q (1 <= q <= 100000)，分别代表有 n 个小镇，表皮卡丘有 q 次操作。

接下来输入 q 行，对于每次操作，先输入操作类型，然后根据操作类型读入：

1: 紧接着输入 2 个整数 a (1 <= a <= n), b (1 <= b <= 60)，表示第 a 个小镇引入了第 b 种进化石
2: 紧接着输入 2 个整数 a (1 <= a <= n), b (1 <= b <= 60)，表示第 a 个小镇失去了全部第 b 种进化石
3: 紧接着输入 2 个整数 l, r (1 <= l <= r <= n)，表示他想询问从第 l 个到第 r 个小镇上可收集的进化石有哪几种

Output

对于每组输入，首先输出一行 "Case T:"，表示当前是第几组数据。

对于每组数据中的每次 3 操作，在一行中按编号升序输出所有可收集的进化石。如果没有进化石可收集，则输出一个 MeiK 的百分号 "%"（不包括引号）。

Sample

Input

Output

CopyCase 1:
%
50
50
20 50
%
%

Source

【“师创杯”山东理工大学第九届ACM程序设计竞赛正式赛】bLue

题解：题目咋一看感觉是线段树的模板题，但是麻烦在需要你输出石头的种类，这就不是单纯的线段树了，还需要加入一点二进制的思想（好像还是第一次接触到这样的二进制思想题目），用二进制从右向左0代表没有1代表有，如01001代表这个村庄有1号和4号石头，如此便可快速的输出石头的种类，同时线段树的区间合并需要用或运算来维护，区间删除则且上需要删除种类的二进制取反（想一下不难理解）。

#include<stdio.h>
#include <algorithm>
#include<iostream>
#include<string.h>
#include<vector>
#include<stdlib.h>
#include<math.h>
#include<queue>
#include<deque>
#include<ctype.h>
#include<map>
#include<set>
#include<stack>
#include<string>
#include<algorithm>
#define INF 0x3f3f3f3f
#define FAST_IO ios::sync_with_stdio(false)
const double PI = acos(-1.0);
const double eps = 1e-6;
const int MAX=1e5+10;
const int mod=1e9+7;
typedef long long ll;
using namespace std;
#define gcd(a,b) __gcd(a,b)
inline ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
inline ll qpow(ll a,ll b){ll r=1,t=a; while(b){if(b&1)r=(r*t)%mod;b>>=1;t=(t*t)%mod;}return r;}
inline ll inv1(ll b){return qpow(b,mod-2);}
inline ll exgcd(ll a,ll b,ll &x,ll &y){if(!b){x=1;y=0;return a;}ll r=exgcd(b,a%b,y,x);y-=(a/b)*x;return r;}
inline ll read(){ll x=0,f=1;char c=getchar();for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;for(;isdigit(c);c=getchar()) x=x*10+c-'0';return x*f;}const int maxn=100000;
ll a[maxn+5],sum[maxn<<2];
void pushup(int st){sum[st]=sum[st<<1] | sum[st<<1|1];return;}void update(int p,int c,int l,int r,int st)
{if(l==r){sum[st]|=1ll<<(c-1);return;}int m=(l+r)>>1;if(p<=m) update(p,c,l,m,st<<1);else update(p,c,m+1,r,st<<1|1);pushup(st);
}void del(int p,int c,int l,int r,int st)
{if(l==r){sum[st] &= ~(1ll<<(c-1));return;}int m=(l+r)>>1;if(p<=m) del(p,c,l,m,st<<1);else del(p,c,m+1,r,st<<1|1);pushup(st);
}ll query(int L,int R,int l,int r,int st)
{if(l>=L && r<=R) return sum[st];int m=(l+r)>>1;ll ans=0;if(L<=m) ans|=query(L,R,l,m,st<<1);if(R>m) ans|=query(L,R,m+1,r,st<<1|1);return ans;
}int main()
{int t;t=read();for(int k=1;k<=t;k++){memset(sum,0,sizeof(sum));printf("Case %d:\n",k);int n=read();int q=read();while(q--){int node=read();int pos=read();int val=read();if(node==1){update(pos,val,1,n,1);}else if(node==2){del(pos,val,1,n,1);}else{ll ans=query(pos,val,1,n,1);ll cnt=0;int flag=1,ot=1;while(ans){if(ans&1){if(flag) flag=0;else putchar(' ');cout<<ot;}ot++;ans>>=1;}if(flag){printf("%%\n");}elseprintf("\n");}}}return 0;
}