【最小异或生成树】密文

3 篇文章 0 订阅
2 篇文章 0 订阅

【题目描述】
有一串长度为n的密文,密文的每一位都可以用一个非负整数来描述,并且每一位都有一个权值ai。你可以进行任意多次操作,每次操作可以选择连续一段密文,花费选择的所有位上权值的异或和的代价获得这段密文每一位的异或和。求至少需要花费多少代价才能将密文的每一位都破解出来。

【输入】
第一行一个正整数n。

第二行n个非负整数,表示ai。

【输出】
输出一个整数,表示答案。

【样例输入】
2
1 3
【样例输出】
3
【提示】
对于10%的数据,n≤10。

对于20%的数据,n≤100。

对于40%的数据,n≤1000。

对于60%的数据,n≤10000。

对于100%的数据, n ≤ 1 0 5 , a i ≤ 1 0 9 n≤10^5,ai≤10^9 n105,ai109

【思路】

这也是一道好题呀。我们需要破解密文的每一位,等价于把密文的前缀异或和 s u m i sum_{i} sumi破解出来。考虑一次询问[l,r]的贡献:如果我们知道了 s u m l − 1 sum_{l-1} suml1,就可以知道 s u m r sum_{r} sumr,反之亦然。形象地说,我们把 s u m l − 1 sum_{l-1} suml1 s u m r sum_{r} sumr联通了,同一个连通块里的前缀可以相互推导。而最初我们只知道 s u m 0 = 0 sum_{0}=0 sum0=0,所以破解所有前缀和就变成了使所有前缀和与 s u m 0 sum_{0} sum0联通。这就变成了一个最小异或生成树问题了。按位贪心,分治即可。

#include<bits/stdc++.h>
#define re register
using namespace std;
const int N=3e6+5;
int n;
inline int red()
{
    int data=0;int w=1; char ch=0;
    ch=getchar();
    while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0' && ch<='9') data=(data<<3)+(data<<1)+ch-'0',ch=getchar();
    return data*w;
}
int val[N];
int ch[N][2],tot=0;
bool v[N];
inline int get()
{
	++tot;
	ch[tot][0]=0;
	ch[tot][1]=0;
	v[tot]=0;
	return tot;
}
inline void add(int id)
{
	int u=0;
	for(int re i=30;i>=0;i--)
	{
		bool nxt=(1ll<<i)&val[id];
		if(!ch[u][nxt])ch[u][nxt]=get();
		u=ch[u][nxt];
	}
	v[u]=1;
}
inline long long check(int id)
{
	long long ret=0;
	int u=0;
	for(int re i=30;i>=0;i--)
	{
		bool nxt=((1ll<<i)&val[id]);
		if(ch[u][nxt])u=ch[u][nxt];
		else u=ch[u][nxt^1],ret+=(1ll<<i);
	}
	return ret;
}
long long ans=0;
void solve(vector<int>g,int pos)
{
	if(g.size()<=1)return;
	if(pos==-1)return;
	vector<int>son[2];
	for(int re i=g.size()-1;i>=0;i--)
		son[(bool)((1ll<<pos)&val[g[i]])].push_back(g[i]);
	if(son[0].size()&&son[1].size())
	{
		tot=0;ch[0][0]=0;ch[0][1]=0;
		long long tmp=1e15;
		for(int re i=son[0].size()-1;i>=0;i--)
			add(son[0][i]);
		for(int re i=son[1].size()-1;i>=0;i--)
			tmp=min(tmp,check(son[1][i]));
		ans+=tmp;
	}
	if(son[0].size())solve(son[0],pos-1);
	if(son[1].size())solve(son[1],pos-1);
}
vector<int>g;
int main()
{
	n=red();val[++n]=0;
	for(int re i=1;i^n;i++)val[i]=red()^val[i-1];
	for(int re i=1;i<=n;i++)g.push_back(i);
	solve(g,30);
	cout<<ans<<"\n";
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值