BZOJ.5092.[Lydsy1711月赛]分割序列(高维前缀和)
\(Description\)
\(Solution\)
首先处理\(a_i\)的前缀异或和\(s_i\)。那么在对于序列\(a_1,...,a_n\),在\(i\)位置处分开的价值为:\(s_i+s_i\ ^{\wedge}s_n\)。
虽然有个加,但依旧可以考虑按位计算。如果\(s_n\)的第\(k\)位为\(1\),那\(s_i\)的第\(k\)位为\(0\)或是\(1\)贡献都是\(2^k\)(贡献即\(s_i+s_i\ ^{\wedge}s_n\)在第\(k\)位上是否为\(1\));如果\(s_n\)的第\(k\)位为\(0\),那么\(s_i\)第\(k\)位为\(0\)则贡献为\(0\),为\(1\)则贡献为\(2*2^k\)。
\(n\)就是指我们当前处理的前缀是\(a[1...n]\)。然后从高到低枚举每一位\(k\),如果\(s_n\)在这一位为\(1\),显然答案一定可以得到\(2^k\)的贡献;否则\(s_n\)在这一位为\(0\),我们应尽量让分割位置\(i\)满足\(s_i\)在第\(k\)位为\(1\),也就是找一个\(n\)前面的位置\(i\)满足\(s_i\)在第\(k\)位为\(1\),如果找得到,答案就可以得到\(2^{k+1}\)的贡献,并限制了\(s_i\)的第\(k\)位为\(1\)。
继续枚举更低位\(k'\)时,在第二种情况\(s_i\)不仅要满足\(s_i\)在第\(k'\)位为\(1\),还要满足之前的第\(k\)位为\(1\),也就是找是否存在\(s_i\)第\(k,k'\)位同时为\(1\)的\(n\)前面的位置\(i\)。
之后同理。
也就是说我们要求是否存在\(i\leq n\),\(s_i\)的第\(k\)位为\(1\)且前\(k\)位都满足之前的限制(某些位必须为\(1\))。
不妨去求,第\(k\)位为\(1\)且满足限制的最靠前的位置\(i\),判断是否有\(i\leq n\)。
因为限制就是某些位必须为\(1\),其它位任意,也就是超集。所以用高维前缀和维护满足某种限制的集合中,最靠前的位置就可以了。
复杂度\(O(2^kk)\)。
//6576kb 1748ms
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 500000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=(1<<20)+5;
int s[300005],f[N];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
int main()
{
static int pw[30];
pw[0]=1;
for(int i=1; i<=21; ++i) pw[i]=pw[i-1]<<1;
int n=read(),mx=0;
memset(f,0x3f,sizeof f);
for(int i=1,t; i<=n; ++i)
s[i]=s[i-1]^read(), mx=std::max(mx,s[i]), f[s[i]]=std::min(f[s[i]],i);
int bit=1;
while(pw[bit]<=mx) ++bit;
for(int i=0,lim=1<<bit; i<bit; ++i)
for(int s=0; s<lim; ++s)
if(!(s&pw[i])) f[s]=std::min(f[s],f[s|pw[i]]);
for(int i=1; i<=n; ++i)
{
int ans=0,t=0;
for(int j=bit-1; ~j; --j)
if(s[i]&pw[j]) ans|=pw[j];
else if(f[t|pw[j]]<=i) t|=pw[j], ans+=pw[j+1];//+= not |=...
printf("%d\n",ans);
}
return 0;
}
BZOJ.5092.[Lydsy1711月赛]分割序列(高维前缀和)的更多相关文章
- bzoj 5092 [Lydsy1711月赛]分割序列——高维前缀和
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=5092 套路地弄一个前缀异或和,就变成 f[ i ]=max_{j=0}^{i} { s[ j ...
- bzoj 5092: [Lydsy1711月赛]分割序列
5092: [Lydsy1711月赛]分割序列 Time Limit: 5 Sec Memory Limit: 256 MBSubmit: 219 Solved: 100[Submit][Stat ...
- bzoj 5092 [Lydsy1711月赛]分割序列 贪心高维前缀和
[Lydsy1711月赛]分割序列 Time Limit: 5 Sec Memory Limit: 256 MBSubmit: 213 Solved: 97[Submit][Status][Dis ...
- BZOJ:5092 [Lydsy1711月赛]分割序列(贪心&高维前缀和)
Description 对于一个长度为n的非负整数序列b_1,b_2,...,b_n,定义这个序列的能量为:f(b)=max{i=0,1,...,n}((b_1 xor b _2 xor...xor ...
- BZOJ5092:[Lydsy1711月赛]分割序列(贪心,高维前缀和)
Description 对于一个长度为n的非负整数序列b_1,b_2,...,b_n,定义这个序列的能量为:f(b)=max{i=0,1,...,n}((b_1 xor b_2 xor...xor b ...
- bzoj 5092 分割序列 —— 高维前缀和
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=5092 首先,处理出异或前缀和 s[i],i 位置的答案就是 s[j] + s[j]^s[i] ...
- bzoj 5094 [Lydsy1711月赛]硬盘检测 概率dp
[Lydsy1711月赛]硬盘检测 Time Limit: 1 Sec Memory Limit: 256 MBSubmit: 273 Solved: 75[Submit][Status][Dis ...
- bzoj 5093 [Lydsy1711月赛]图的价值 NTT+第二类斯特林数
[Lydsy1711月赛]图的价值 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 245 Solved: 128[Submit][Status][D ...
- BZOJ 5093[Lydsy1711月赛]图的价值 线性做法
博主曾更过一篇复杂度为$O( k· \log k)$的多项式做法在这里 惊闻本题有$ O(k)$的神仙做法,说起神仙我就想起了于是就去学习了一波 幂与第二类斯特林数 推导看这里 $$ x^k=\sum ...
随机推荐
- MySQL 导出数据库,出现 “mysqldump: Got error: 1146”
出现场景 在 cmd 导出数据库时: mysqldump -hlocalhost -uroot -p student_db > C:\student_db.sql 出现: mysqldump: ...
- (BFS) leetcode 279. Perfect Squares
Given a positive integer n, find the least number of perfect square numbers (for example, 1, 4, 9, 1 ...
- Docker:企业级私有仓库harbor[十六]
一.安装配置 1.下载安装包 链接:https://pan.baidu.com/s/1Z9I7zYXSt-8ve3lFT2YCeg 提取码:iuqj 2.安装docker和docker-compose ...
- Entity Framework入门教程(3)---EF中的上下文简介
1.DbContext(上下文类) 在DbFirst模式中,我们添加一个EDM(Entity Data Model)后会自动生成一个.edmx文件,这个文件中包含一个继承DbContext类的上下文实 ...
- Quartz.net 3.x使用总结(二)——Db持久化和集群
上一篇简单介绍了Quartz.net的概念和基本用法,这一篇记录一下Quartz.net通过数据库持久化Trigger和Jobs等数据,并简单配置Quartz.net的集群. 1.JobStore介绍 ...
- 2018-2019-2 《Java程序设计》第3周学习总结
20175319 2018-2019-2 <Java程序设计>第3周学习总结 教材学习内容总结 第三周通过课本与蓝墨云上的视频学习了<Java2实用教程>第四章类与对象 成功激 ...
- docker学习------docker私有仓库的搭建
192.168.138.102:23451.私有仓库的搭建(docker pull registry),拉取最新的镜像 2.查看拉取的仓库镜像(docker images) 3.启用registry镜 ...
- Mysql的跨服务器 关联查询--Federated引擎
1.确认开启Federated引擎 查询FEDERATED功能是否开启: show ENGINES; 2.如果状态为NO则需修改my.ini文件,增加一行federated配置: ...
- 用juery的ajax方法调用aspx.cs页面中的webmethod方法
首先在 aspx.cs文件里建一个公开的静态方法,然后加上WebMethod属性. 如: [WebMethod] public static string GetUserName() { //.... ...
- Python中的 @staticmethod@classmethod方法
python类中有三种方法,常见的是实例方法,另外两种是staticmethod装饰的静态方法,和classmethod装饰的类方法. 1.对比 流畅的python里,用一个例子进行了对比: (1)两 ...