12
3
2015
0

【bzoj3261】最大异或和

3261: 最大异或和

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 902  Solved: 378
[Submit][Status][Discuss]

Description

     

给定一个非负整数序列 {a},初始长度为 N。       
有   M个操作,有以下两种操作类型:
 
1 、A x:添加操作,表示在序列末尾添加一个数 x,序列的长度 N+1。
2 、Q l r x:询问操作,你需要找到一个位置 p,满足 l<=p<=r,使得:
 
a[p] xor a[p+1] xor ... xor a[N] xor x 最大,输出最大是多少。  

Input

第一行包含两个整数 N  ,M,含义如问题描述所示。   
第二行包含 N个非负整数,表示初始的序列 A 。 
 
接下来 M行,每行描述一个操作,格式如题面所述。   

Output

假设询问操作有 T个,则输出应该有 T行,每行一个整数表示询问的答案。

Sample Input

5 5
2 6 4 3 6
A 1
Q 3 5 4
A 4
Q 5 7 0
Q 3 6 6
对于测试点 1-2,N,M<=5 。
对于测试点 3-7,N,M<=80000 。
对于测试点 8-10,N,M<=300000 。
其中测试点 1, 3, 5, 7, 9保证没有修改操作。
对于 100% 的数据, 0<=a[i]<=10^7。

Sample Output

4
5
6

HINT

对于      100%  的数据,     0<=a[i]<=10^7  

首先我们令b[i]=a[1]^a[2]……^a[i]。因此题目中的a[p]^a[p+1]……^a[n]^x就等于x^b[n]^b[p-1],x^b[n]是定值,我们只要用贪心去得到最大值就可以了。因此我们只要把b数组建一个二进制trie树,然后可持久化,问题就迎刃而解了。

代码:

/**************************************************************
    Problem: 3261
    User: yeweining
    Language: C++
    Result: Accepted
    Time:4148 ms
    Memory:207836 kb
****************************************************************/
 
#include<cstdio>
#include<algorithm>
using namespace std;
int n,m,u,v,a[1000010],b[2000010],cnt,ans,x,y,root[5000010];
char s[10];
struct node{int sum,a[2];}tr[15000010];
void change(int &p,int dep,int x)
{
    tr[++cnt]=tr[p];p=cnt;tr[p].sum++;
    if (dep<0) return;
    int t=(x>>dep)&1;
    change(tr[p].a[t],dep-1,x);
}
void ask(int p1,int p2,int dep,int x)
{
    if (dep<0) return;
    int t=(x>>dep)&1;
    if (tr[tr[p2].a[t^1]].sum!=tr[tr[p1].a[t^1]].sum)
    {
        ans+=(1<<dep);
        ask(tr[p1].a[t^1],tr[p2].a[t^1],dep-1,x);
    }else ask(tr[p1].a[t],tr[p2].a[t],dep-1,x);
}
int main()
{
    scanf("%d%d",&n,&m);
    n++;
    for (int i=2;i<=n;i++)
        scanf("%d",&a[i]);
    for (int i=2;i<=n;i++)
        b[i]=b[i-1]^a[i];
    for (int i=1;i<=n;i++)
    {
        root[i]=root[i-1];
        change(root[i],25,b[i]);
    }
    while (m--)
    {
        scanf("%s",s);
        if (s[0]=='A')
        {
            scanf("%d",&x);
            n++;
            b[n]=(b[n-1]^x);
            root[n]=root[n-1];
            change(root[n],25,b[n]);
        }else
        {
            scanf("%d%d%d",&u,&v,&x);
            ans=0;
            ask(root[u-1],root[v],25,(x^b[n]));
            printf("%d\n",ans);
        }
    }
Category: 可持久化trie树 | Tags:
12
1
2015
0

【JSOI2015】字符串树

这是一道裸的可持久化trie树。这个方法就不多说了,详见“初探主席树”及“cout on a tree”。

代码:

 

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100010
using namespace std;
int n,u,v,cnt,l,q;
int edgenum,vet[2*N],next[2*N],head[2*N],root[N];
int f[N],d[N],num[N],son[N],top[N];
struct node{int sum,son[30];}tr[10*N];
char str[15],s[N][15];
void add(int u,int v)
{
	edgenum++;
	vet[edgenum]=v;
	next[edgenum]=head[u];
	head[u]=edgenum;
}
void dfs(int u,int fa,int dep)
{
	f[u]=fa;
	d[u]=dep;
	num[u]=1;
	int maxnum=0;
	for (int e=head[u];e;e=next[e])
	{
		int v=vet[e];
		if (v!=fa)
		{
			dfs(v,u,dep+1);
			num[u]+=num[v];
			if (num[v]>maxnum){maxnum=num[v];son[u]=v;}
		}
	}
}
void dfs(int u,int number)
{
	top[u]=number;
	if (!son[u]) return;
	dfs(son[u],number);
	for (int e=head[u];e;e=next[e])
	{
		int v=vet[e];
		if (v!=f[u]&&v!=son[u]) dfs(v,v);
	}
}
int getlca(int u,int v)
{
	while (top[u]!=top[v])
	{
		if (d[top[u]]<d[top[v]]) swap(u,v);
		u=f[top[u]];
	}
	if (d[u]<d[v]) return u;else return v;
}
void change(int &p,int k,int dep)
{
	tr[++cnt]=tr[p];p=cnt;tr[p].sum++;
	if (dep==strlen(s[k])) return;
	change(tr[p].son[s[k][dep]-'a'],k,dep+1);
}
int ask(int p,int dep)
{
	if (p==0&&dep) return 0;
	if (dep==l) return tr[p].sum;
	return ask(tr[p].son[str[dep]-'a'],dep+1);
}
void dfs1(int u,int e)
{
	if (u!=1){root[u]=root[f[u]];change(root[u],(e+1)>>1,0);}
	for (int e=head[u];e;e=next[e])
	{
		int v=vet[e];
		if (v!=f[u]) dfs1(v,e);
	}
}
int main()
{
	scanf("%d",&n);
	for (int i=1;i<n;i++)
	{
		scanf("%d%d%s",&u,&v,s[i]);
		add(u,v);add(v,u);
	}
	dfs(1,0,1);
	dfs(1,1);
	dfs1(1,0);
	scanf("%d",&q);
	while (q--)
	{
		scanf("%d%d%s",&u,&v,str);
		l=strlen(str);
		printf("%d\n",ask(root[u],0)+ask(root[v],0)-2*ask(root[getlca(u,v)],0));
	}
}

还有一种我在练习赛上自己yy的一种方法。然而,有些复杂(我是有自知之明的),文字解释起来非常复杂,代码放着,有兴趣的人可以自己研究一下。

代码:

 

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100010
using namespace std;
typedef long long ll;
int n,q,cnt,l;
int edgenum,Edgenum,vet[2*N],next[2*N],head[2*N],Vet[10*N],Next[10*N],Head[10*N],x[N],y[N],a[N],b[N],lca[N];
ll dis[2*N],Ask[10*N],P[10*N],p[10*N],ans[N],mi[20],w[N],c[N],d[N];
int f[N],de[N],num[N],son[N],top[N],tr[10*N];
char s[20];
void add(int u,int v,ll w)
{
	edgenum++;
	vet[edgenum]=v;
	next[edgenum]=head[u];
	head[u]=edgenum;
	dis[edgenum]=w;
}
void dfs(int u,int fa,int dep)
{
	f[u]=fa;
	de[u]=dep;
	num[u]=1;
	int maxnum=0;
	for (int e=head[u];e;e=next[e])
	{
		int v=vet[e];
		if (v!=fa)
		{
			dfs(v,u,dep+1);
			num[u]+=num[v];
			if (num[v]>maxnum){maxnum=num[v];son[u]=v;}
		}
	}
}
void dfs(int u,int number)
{
	top[u]=number;
	if (!son[u]) return;
	dfs(son[u],number);
	for (int e=head[u];e;e=next[e])
	{
		int v=vet[e];
		if (v!=f[u]&&v!=son[u]) dfs(v,v);
	}
}
int getlca(int u,int v)
{
	while (top[u]!=top[v])
	{
		if (de[top[u]]<de[top[v]]) swap(u,v);
		u=f[top[u]];
	}
	if (de[u]<de[v]) return u;else return v;
}
void Add(int u,int v,ll w,ll opt)
{
	Edgenum++;
	Vet[Edgenum]=v;
	Next[Edgenum]=Head[u];
	Head[u]=Edgenum;
	Ask[Edgenum]=w;
	P[Edgenum]=opt;
}
int asksum(ll x)
{
	int ans=0;
	while (x)
	{
		ans+=tr[x];
		x-=((x)&(-x));
	}
	return ans;
}
void tradd(ll x,int val)
{
	while (x<=cnt)
	{
		tr[x]+=val;
		x+=((x)&(-x));
	}
}
void dfs(int u)
{
	for (int e=Head[u];e;e=Next[e])
	{
		int v=Vet[e];
		ans[v]+=(ll)asksum(Ask[e])*P[e];
	}
	for (int e=head[u];e;e=next[e])
	{
		int v=vet[e];
		if (v!=f[u])
		{
			tradd(dis[e],1);
			dfs(v);
			tradd(dis[e],-1);
		}
	}
}
int main()
{
	freopen("strings.in","r",stdin);
	freopen("strings.out","w",stdout);
	scanf("%d",&n);
	mi[0]=1;
	for (int i=1;i<=10;i++) mi[i]=mi[i-1]*31;
	for (int i=1;i<n;i++)
	{
		scanf("%d%d%s",&x[i],&y[i],s);
		l=strlen(s);
		w[i]=0;
		for (int j=0;j<l;j++)
			w[i]=w[i]*31+s[j]-'a'+1;
		w[i]=w[i]*mi[10-l];
		p[++cnt]=(ll)w[i];
		add(x[i],y[i],w[i]);add(y[i],x[i],w[i]);
	}
	dfs(1,0,1);
	dfs(1,1);
	scanf("%d",&q);
	for (int i=1;i<=q;i++)
	{
		scanf("%d%d%s",&a[i],&b[i],s);
		lca[i]=getlca(a[i],b[i]);
		c[i]=0;
		l=strlen(s);
		for (int j=0;j<l;j++)
			c[i]=c[i]*31+s[j]-'a'+1;
		c[i]=c[i]*mi[10-l];
		d[i]=c[i]+(mi[10-l]-1)/30*('z'-'a'+1);
		p[++cnt]=(ll)c[i]-1;
		p[++cnt]=(ll)d[i];
		Add(a[i],i,d[i],1);
		Add(b[i],i,d[i],1);
		Add(lca[i],i,d[i],-2);
		Add(a[i],i,c[i]-1,-1);
		Add(b[i],i,c[i]-1,-1);
		Add(lca[i],i,c[i]-1,2);
	}
	sort(p+1,p+cnt+1);
	cnt=unique(p+1,p+cnt+1)-p-1;
	for (int i=1;i<=edgenum;i++)
		dis[i]=lower_bound(p+1,p+cnt+1,dis[i])-p;
	for (int i=1;i<=Edgenum;i++)
		Ask[i]=lower_bound(p+1,p+cnt+1,Ask[i])-p;
	dfs(1);
	for (int i=1;i<=q;i++)
		printf("%lld\n",ans[i]);
}
Category: 可持久化trie树 | Tags:

Host by is-Programmer.com | Power by Chito 1.3.3 beta | Theme: Aeros 2.0 by TheBuckmaker.com