【Trie】[CQOI2016]路由表

题目描述

对于一次查询的一种理解方式是：无视其它所有查询操作，只看添加操作。先清空路由表，然后执行第1到a-1次添加操作。之后再执行第a到b次添加操作过程中，统计匹配改变的次数。

数据范围：

设一条表项的掩码长度为L，数据保证将目的地址转为二进制串后，末尾的32-L位均为0。

分析

在线地，每次在添加表项至路由表的时候，将该地址加入trie，然后在结束的节点标记上添加这个表项的时间。

对于每次询问，在trie上暴力匹配，然后用一个单调栈，使得能够匹配的表项出现的时间随着长度的增加而增加，最后通过二分查找找到在a,b之间改变了多少次即可。

代码

#include<cstdio>
#include<algorithm>
#define MAXN 1000000
#define INF 0xffffffffu
using namespace std;
typedef unsigned int uint;
int n,cnt;
struct query{
uint c;
int len;
inline query(){
}
inline query(uint c,int len):c(c),len(len){
}
}q[MAXN+10];
struct node{
int pos;
int ch[2];
}tree[MAXN*32+10],*root=tree,*tcnt=tree;
template<class T>
void Read(T &x){
char c;
while(c=getchar(),c!=EOF)
if(c>='0'&&c<='9'){
x=c-'0';
while(c=getchar(),c>='0'&&c<='9')
x=x*10+c-'0';
ungetc(c,stdin);
return;
}
}
void read(){
Read(n);
}
void solve(){
char s[20];
uint c,t;
int len,i,a,b,best,ans;
while(n--){
scanf("%s",s);
if(s[0]=='A'){
c=a=0;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(len);
q[++cnt]=query(c,len);
}
else{
c=a=ans=0;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a),Read(b);
best=0;
for(i=1;i<a;i++){
t=INF^((1ll<<(32-q[i].len))-1);
if(q[i].c==(c&t))
if(q[i].len>q[best].len)
best=i;
}
for(;i<=b;i++){
t=INF^((1ll<<(32-q[i].len))-1);
if(q[i].c==(c&t))
if(q[i].len>q[best].len)
best=i,ans++;
}
printf("%d\n",ans);
}
}
}
int s[MAXN+10],tp;
void solve2(){
char ss[20];
int a,b,len,i,x[40],pos,ans;
uint c;
bool d;
node *p;
while(n--){
scanf("%s",ss);
if(ss[0]=='A'){
c=0;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(a);
c=(c<<8)+a;
Read(len);
p=root;
for(i=31;i>=32-len;i--){
d=(c>>i)&1;
if(!p->ch[d])
p->ch[d]=++tcnt-tree;
p=tree+p->ch[d];
}
p->pos=++cnt;
}
else{
len=0;
Read(a);
for(i=7;i>=0;i--)
x[++len]=(a>>i)&1;
Read(a);
for(i=7;i>=0;i--)
x[++len]=(a>>i)&1;
Read(a);
for(i=7;i>=0;i--)
x[++len]=(a>>i)&1;
Read(a);
for(i=7;i>=0;i--)
x[++len]=(a>>i)&1;
Read(a),Read(b);
tp=0;
p=root;
for(i=1;i<=32;i++){
if(p->ch[x[i]])
p=tree+p->ch[x[i]];
else
break;
if(p->pos&&p->pos<=b){
pos=p->pos;
while(tp>0&&s[tp]>pos)
tp--;
s[++tp]=pos;
}
}
ans=tp-(lower_bound(s+1,s+tp+1,a)-s)+1;
printf("%d\n",ans);
}
}
}
int main()
{
read();
if(n<=1000)
solve();
else
solve2();
}