前 k 大异或值问题

最新推荐文章于 2024-05-21 23:16:53 发布
最新推荐文章于 2024-05-21 23:16:53 发布
阅读量825 收藏 1
点赞数 2
分类专栏： 字典树 文章标签： 算法 数据结构 c++ 贪心算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_51547083/article/details/122747322
版权

多区间第 k 大异或值(k较小)

链接：[P5283 十二省联考 2019] 异或粽子 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

题意：给定序列 a，求 n 2 n^2 n2 个子区间中前 k 大子区间异或和的和。

题解：首先区间异或和可以用前缀异或和 O ( 1 ) O(1) O(1) 查询，那么题目转化为求前 k 大的点对异或的和。用堆来维护最大值，遍历序列求出所有固定右端点的第 1 大异或值，丢进堆里。每次取最大值，假设当前为 第5位的第 t 大值，则往堆内放入第 5 第 t+1 位最大值。注意同一个右端点加入次数有限制，比如第 1 位只有一个最大值。

//#pragma GCC optimize("O3")
//#pragma GCC optimize("unroll-loops")
#include<iostream>
#include<algorithm>
#include<queue>

using namespace std;
typedef long long ll;
typedef pair<ll,ll> P;
const int maxn=5e5+5;

ll va[maxn<<6],kt[maxn],a[maxn],n,k,ans;
int p[maxn<<6][2],ct[maxn<<6],rt[maxn],now;
priority_queue<P,vector<P>,less<P> >pq;

void insert(int u,int v,ll x)
{
	rt[v]=++now;
	int l=rt[u],r=rt[v];
	for(int i=32;i>=0;i--)
	{
		int k=x>>i&1;
		p[r][k]=++now;
		p[r][k^1]=p[l][k^1];
		ct[r]=ct[l]+1;
		l=p[l][k],r=p[r][k];
	}
	va[r]=x; ct[r]=ct[l]+1;
	return;
}

ll query(int u,ll x,int y)
{
	int r=rt[u];
	for(int i=32;i>=0;i--)
	{
		int q=x>>i&1;
		if(ct[p[r][q^1]]>=y)r=p[r][q^1];
		else y-=ct[p[r][q^1]],r=p[r][q];
	}
	return va[r]^x;
}


int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin>>n>>k; insert(0,1,0);
	for(int i=2;i<=n+1;i++)
	{
		cin>>a[i],a[i]^=a[i-1],insert(i-1,i,a[i]);
	}
	for(int i=2;i<=n+1;i++)pq.push(P(query(i-1,a[i],++kt[i]),i));
	while(k--)
	{
		P x=pq.top(); pq.pop();
		ans+=x.first;
		if(kt[x.second]<x.second-1)
		{
			pq.push(P(query(x.second-1,a[x.second],++kt[x.second]),x.second)); 
		}
	}
	cout<<ans<<"\n";
	return 0;
}

多区间第 k 大异或值加强版

链接：CF241B Friends - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

与上面相比主要是 k 值扩大，无法枚举 k，采取二分的形式，先确定 k 值位置，然后将大于 k 值的和累加。每个点考虑层贡献。求取大于 k 值的所有数与一个数异或和的和可以先维护每个节点下边所有层带的 1 和 0 的个数，然后就可以直接遍历这个节点的 log 层考虑每层的贡献即可。复杂度 O ( n l o g 2 w ) O(nlog^2w) O(nlog2w)。

//#pragma GCC optimize("O3")
//#pragma GCC optimize("unroll-loops")
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>

using namespace std;
typedef long long ll;
const int maxn=5e4+5;
const int mod=1e9+7;

ll a[maxn],n,m,sum,res;
int p[maxn<<5][2],ct[maxn<<5],rt[maxn],l[maxn],r[maxn],now;
unsigned short int sm[maxn<<5][31][2];

inline void insert(int u,int v,ll x)
{
	rt[v]=++now;
	int l=rt[u],r=rt[v];
	for(int i=30;i>=0;i--)
	{
		int q=x>>i&1;
		p[r][q]=++now;
		ct[r]=ct[l]+1;
		p[r][q^1]=p[l][q^1];
		for(int j=i;j>=0;j--)
		{
			sm[r][j][x>>j&1]=sm[l][j][x>>j&1]+1;
			sm[r][j][(x>>j&1)^1]=sm[l][j][(x>>j&1)^1];
		}
		l=p[l][q],r=p[r][q];
	}
	ct[r]=ct[l]+1;
	return;
}

inline int query(int u,ll x,ll mid)
{
	int r=rt[u],cnt=0; 
	for(int i=30;i>=0;i--)
	{
		int q=x>>i&1;
		if(mid>>i&1)r=p[r][q^1];
		else
		{
			cnt+=ct[p[r][q^1]];
			r=p[r][q];
		}
	}
	cnt+=ct[r]; 
	return cnt;
}

inline bool check(ll mid)
{
	ll cnt=0;
	for(int i=2;i<=n;i++)cnt+=query(i-1,a[i],mid);
	return cnt>=m;
}

void Query(int u,ll x,ll y)
{
	int r=rt[u],dep[31]={0};
	for(int i=30;i>=0;i--)
	{
		int q=x>>i&1;
		if(y>>i&1)r=p[r][q^1];
		else
		{
			res+=ct[p[r][q^1]];
			for(int j=30;j>i;j--)dep[j]+=ct[p[r][q^1]]*(y>>j&1);
			dep[i]+=ct[p[r][q^1]];
			for(int j=i-1;j>=0;j--)dep[j]+=sm[p[r][q^1]][j][(x>>j&1)^1];
			r=p[r][q];
		}
	}
	res+=ct[r]; 
	for(int i=30;i>=0;i--)dep[i]+=ct[r]*(y>>i&1);
	for(int i=30;i>=0;i--)sum+=dep[i]*(1ll<<i); 
	return;
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>a[i],insert(i-1,i,a[i]);
	}
	ll l=0,r=1<<30;
	while(l<r)
	{
		ll mid=l+r+1>>1;
		if(check(mid))l=mid;
		else r=mid-1;
	}
	for(int i=2;i<=n;i++)Query(i-1,a[i],l);
	cout<<((sum-(res-m)*l)%mod+mod)%mod<<"\n";
	return 0;
}

更优解法

因为二分出来的是一个数，所有点都跑这一个数，所以直接边把所有点的数往同一个方向移动即可。注意可持久化要左右都要继承或修改，不要忘了继承这一项。

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>

using namespace std;
typedef long long ll;
const int maxn=5e4+5;
const int mod=1e9+7;

ll a[maxn],n,k,cnt,ans,res;
int p[maxn<<5][2],ct[maxn<<5],rt[maxn],r[maxn],now;
short unsigned int sm[maxn<<5][31][2];

inline void insert(int u,int v,int x)
{
	rt[v]=++now;
	int l=rt[u],r=rt[v];
	for(int i=30;i>=0;i--)
	{
		int q=x>>i&1;
		p[r][q]=++now;
		ct[r]=ct[l]+1;
		for(int j=i;j>=0;j--)
		{
			sm[r][j][x>>j&1]=sm[l][j][x>>j&1]+1;
			sm[r][j][(x>>j&1)^1]=sm[l][j][(x>>j&1)^1];
		}
		p[r][q^1]=p[l][q^1];
		l=p[l][q],r=p[r][q];
	}
	ct[r]=ct[l]+1;
	return; 
} 

inline void query()
{
	for(int i=2;i<=n;i++)r[i]=rt[i-1];
	ll dep[31]={0},y=0;
	for(int i=30;i>=0;i--)
	{
		int cp=0;
		for(int j=2;j<=n;j++)
		{
			cp+=ct[p[r[j]][(a[j]>>i&1)^1]];
		}
		if(cnt+cp>=k)
		{
			for(int j=2;j<=n;j++)
			{
				r[j]=p[r[j]][(a[j]>>i&1)^1];
			}
			y+=(1ll<<i);
		}
		else
		{
			cnt+=cp;
			for(int j=2;j<=n;j++)
			{
				int nr=p[r[j]][(a[j]>>i&1)^1];
				for(int o=30;o>i;o--)dep[o]+=ct[nr]*(y>>o&1);
				dep[i]+=ct[nr];
				for(int o=i-1;o>=0;o--)dep[o]+=sm[nr][o][(a[j]>>o&1)^1];
				r[j]=p[r[j]][a[j]>>i&1];
			}
		}
	}
	for(int i=2;i<=n;i++)
	{
		cnt+=ct[r[i]];
		for(int j=30;j>=0;j--)
		{
			dep[j]+=ct[r[i]]*(y>>j&1);
		}
	}
	for(int i=30;i>=0;i--)ans+=(1ll<<i)*dep[i];
	res=y;
	return;
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin>>n>>k;
	for(int i=1;i<=n;i++)cin>>a[i],insert(i-1,i,a[i]);
	query();
	cout<<((ans-(cnt-k)*res)%mod+mod)%mod<<"\n";
	return 0;
}

多区间第 k 大异或值更优解法

链接：[P5795 THUSC2015]异或运算 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

上面通过枚二分来寻找第 k 大异或值。可以发现每次对 mid 都是使所有数异或出来为 mid，则可以直接将所有数进行同时转移，判定是该使其异或 x i x_{i} xi​ 为哪一位去同时走。复杂度 O ( q n l o g w ) O(qnlogw) O(qnlogw)。

//#pragma GCC optimize("O3")
//#pragma GCC optimize("unroll-loops")
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>

using namespace std;
typedef long long ll;
const int maxn=3e5+5;

int x[maxn],l[maxn],r[maxn],n,m,t;
int p[maxn<<5][2],ct[maxn<<5],rt[maxn],now;

inline void insert(int u,int v,int x)
{
	rt[v]=++now;
	int l=rt[u],r=rt[v];
	for(int i=30;i>=0;i--)
	{
		int q=x>>i&1;
		p[r][q]=++now;
		ct[r]=ct[l]+1;
		p[r][q^1]=p[l][q^1];
		l=p[l][q],r=p[r][q];
	}
	ct[r]=ct[l]+1;
	return;
}

inline int query(int u,int v,int ql,int qr,int k)
{
	int ans=0;
	for(int i=ql;i<=qr;i++)l[i]=rt[u],r[i]=rt[v];
	for(int i=30;i>=0;i--)
	{
		int res=0;
		for(int j=ql;j<=qr;j++)
		{
			int q=x[j]>>i&1;
			res+=ct[p[r[j]][q^1]]-ct[p[l[j]][q^1]];
		}
		if(res>=k)
		{
			ans+=1<<i;
			for(int j=ql;j<=qr;j++)
			{
				int q=x[j]>>i&1;
				l[j]=p[l[j]][q^1],r[j]=p[r[j]][q^1];	
			}
		}
		else
		{
			k-=res;
			for(int j=ql;j<=qr;j++)
			{
				int q=x[j]>>i&1;
				l[j]=p[l[j]][q],r[j]=p[r[j]][q];
			}
		} 
	}
	return ans;
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)cin>>x[i];
	for(int i=1;i<=m;i++)
	{
		int u; cin>>u,insert(i-1,i,u);
	}
	cin>>t;
	for(int i=1;i<=t;i++)
	{
		int u,d,l,r,k;
		cin>>u>>d>>l>>r>>k;
		cout<<query(l-1,r,u,d,k)<<"\n";
	}
	return 0;
}
__LazyCat__
关注
  • 2
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
php使用异或实现的加密解密实例
12-18
则b=a ^ c (^是异或的意思)，php在处理异或的字符时先把字符转化为二进制的ascii值，对这些值进行异或，获取结果后在将ascii值转化为字符，原理说晚了直接贴实现的代码：复制代码 代码如下:echo ‘<meta charset=...
#114. k 大异或和
有顶天
11-25 773
师从 zeleqian ZlycerQan前辈，学习k大异或和 链接： https://www.cnblogs.com/ZlycerQan/p/7223893.html 题目描述 这是一道模板题。 给由 n n n 个数组成的一个可重集 S S S，每次给定一个数 k k k，求一个集合 T⊆S T \subseteq S T⊆S，使得集合 T T T 在 S S S 的所有非空子集的不...
【LeetCode】﹝异或ி﹞解码异或、区间查询、第k大异或值
热门推荐
像写诗一样优雅地写代码~
12-04 1万+
【LeetCode】﹝异或ி﹞解码异或、区间查询、第k大异或值 文章目录【LeetCode】﹝异或ி﹞解码异或、区间查询、第k大异或值解码异或后的数组★解码异或后的排列★★找出所有子集的异或总和再求和★形成两个异或相等数组的三元组数目★★子数组异或查询★★数组中两个数的最大异或值★★找出第 K 大的异或坐标值★★ 解码异或后的数组★ 1720. 解码异或后的数组 【题目】未知 整数数组 arr 由 n 个非负整数组成。 经编码后变为长度为 n - 1 的另一个整数数组 encoded ，其中 encoded
【十二省联考2019】异或粽子【01Trie】【堆】【前k大套路】
Lstdo的博客
04-15 187
题意：给定长度为nnn的序列aaa，求前kkk大的区间异或和之和。 n≤5×105,k≤min⁡(n(n−1)2,2×105),ai<232n\leq 5\times 10^5,k\leq \min(\frac{n(n-1)}2,2\times10^5),a_i<2^{32}n≤5×105,k≤min(2n(n−1)​,2×105),ai​<232 原来省选题这么裸，i了i了 显...
各大高校操作系统考试题（相当全面）
06-26
含全国各大高校的操作系统试题，很精彩，如果需要就下来看看，不容错过
Excel公式与函数大辞典.宋翔(带书签高清文字版).pdf
03-02
3.1.4 两位数年份问题 121 3.2 返回当前的日期、时间和指定的日期、时间 121 3.2.1 NOW——返回当前日期和时间 121 3.2.2 TODAY——返回当前日期 123 3.2.3 DATE——返回指定日期的序列号 124 3.2.4 TIME——...
leetcode中文版-LeetCode-Solution-PPT:我给LeetCode中文版制作的题解PPT
06-29
题：数组中两个数的最大异或值 题解地址： LeetCode 第 1079 题：活字印刷 题解地址： LeetCode 第 1080 题：根到叶路径上的不足节点 题解地址： LeetCode 第 1081 题：不同字符的最小子序列 题解地址：
leetcode题库-work:C++后端开发面经整理
06-29
出现前缀字符串查询或者异或值直接用 高精度 ★ 大数乘法或大数加法 注意 数值较大的结果考虑用long long 动态规划写不来，先写暴力过部分测试样例 特殊输出先提交一下，骗点分了再说 练习处理诡异输入的能力 面试常...
链接表存储图(C++注释详解): 构建表 & 深度优先遍历 (DFS)
H13420972436的博客
05-15 435
代码实现: 构建表 & 深度优先遍历 (DFS)
【算法】二分算法——x的平方根
2302_79031646的博客
05-17 799
用二分边界算法来求解——x的平方根，从暴力求解到二分的优化细说
day09
2202_76101487的博客
05-17 738
正值少年的日常打卡！！！
BERT系列算法解读
qq_52053775的博客
05-19 59
从而提高模型在SQuAD、MNLI和SST-2数据集上的性能。总的来看，这些优化措施显著增强了模型的表现。
基于地理坐标的高阶几何编辑工具算法（8）——整形面
HouGISer的博客
05-19 610
整形面算法常用于快速修改一个几何面的局部形状。
【算法】前缀和——二维前缀和模板题
最新发布
2302_79031646的博客
05-21 219
用“二维前缀和模板题”来介绍二位前缀和算法模板
算法练习第22天|39. 组合总和、40.组合总和II
Brillian123的博客
05-15 825
39. 组合总和39. 组合总和 - 力扣（LeetCode）
代码随想录算法训练营第五十五天|LeetCode583.两个字符串的删除操作 、LeetCode72.编辑距离
m0_64514622的博客
05-16 880
dp[i][j]:表示以下标i-1为结尾的word1,和以j-1为结尾的字符串word2，想要达到相等，所需要删除元素的最少次数。dp[i][j]:表示以下标i-1为结尾的word1,和以j-1为结尾的字符串word2，最近编辑距离为dp[i][j]word2为空字符串，以i-1为结尾的字符串word1要要删除i个元素才能和word2相同。word1为空字符串，以j-1为结尾的字符串word2要要删除j个元素才能和word1相同。操作三：word1替换word1[i-1]，使其与word2[j-1]相同。
力扣-219. 存在重复元素 II
ls3614140的博客
05-19 103
给你一个整数数组 nums 和一个整数 k ，判断数组中是否存在两个 不同的索引 i 和 j ，满足 nums[i] == nums[j] 且 abs(i - j) <= k。如果存在，返回 true；否则，返回 false。输入：nums = [1,2,3,1,2,3], k = 2。输入：nums = [1,2,3,1], k = 3。输入：nums = [1,0,1,1], k = 1。
算法001：移动零
m0_62319039的博客
05-17 390
moveZeroes
向量间相互异或verilog代码
06-02
每个异或模块中的`assign`语句可以使用`^`运算符将所有位上的值进行异或运算，得到一个10位二进制数，作为输出数组中对应位置的值。 需要注意的是，由于在生成异或模块时使用了循环，因此这个模块只能处理100个向量...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值