异或序列 [set优化DP]

也许更好的阅读体验

$\mathcal{D}escription$
有一个长度为 $n$的自然数序列 $a$，要求将这个序列分成至少 $m$ 个连续子段
每个子段的价值为该子段的所有数的按位异或
要使所有子段的价值按位与的结果最大，输出这个最大值

$T$组询问
$T\le 10,n,m\le 1000,a_i\le 2^{30}$
$\mathcal{S}olution$
实际上数据范围可开大很多

我们贪心的一位一位的确定最终答案，即看当前考虑的位能否为$1$
记$s_i$表示前$i$个数的异或和，$⨁$表示异或
设当前考虑到了第$b$位
令$res=ans|(1<<b)$
一段区间$[j+1,i]$如果是一个合法的区间，可以得到
$\left(s_i⨁s_j\right)\&res=res$
于是我们得到了一个$n^{2}log$的$DP$方程
$f_i=max{f}_i,f_j+1$其中$\left(s_i⨁s_j\right)=res$
枚举位是$log$的，这样就可以$AC$此题了

实际这个$DP$可以进一步优化
$\left(s_i⨁s_j\right)\&res=res$可以推出
$\left(s_i\&res\right)⨁\left(s_j\&res\right)=res$
$\Rightarrow s_i\&res=\left(s_j\&res\right)⨁res$
即要将$s_i$到$s_j$这段作为一个子段必须满足上面的条件
因为题目是至少$m$段，所以分的越多越好
则我们可以考虑完$s_i$的最优答案后将$s_i⨁res$作为第一关键字存进$set$
$f_i=find(s_i⨁res)$
这样一次转移就是$log$的
复杂度为$nlo{g}^2$

$\mathcal{C}ode$

/*******************************
Author:Morning_Glory
LANG:C++
Created Time:2019年10月26日 星期六 09时18分19秒
*******************************/
#include <cstdio>
#include <fstream>
#include <cstring>
#include <set>
#define mp make_pair
using namespace std;
const int maxn = 2003;
//{{{cin
struct IO{
	template<typename T>
	IO & operator>>(T&res){
		res=0;
		bool flag=false;
		char ch;
		while((ch=getchar())>'9'||ch<'0')	flag|=ch=='-';
		while(ch>='0'&&ch<='9')	res=(res<<1)+(res<<3)+(ch^'0'),ch=getchar();
		if (flag)	res=~res+1;
		return *this;
	}
}cin;
//}}}
int n,m,T,ans;
int a[maxn],s[maxn];
set < pair<int,int> > v;
set < pair<int,int> > :: iterator it,nx;
//{{{solve
void solve (int x)
{
	int res=ans|(1<<x);
	bool flag;
	v.clear();
	for (int i=1;i<=n;++i){
		int val=s[i]&res;
		v.insert(mp(val,0));
		nx=it=v.lower_bound(mp(val,0));
		++nx;
		while (nx!=v.end()&&nx->first==val){
			v.erase(it);
			it=nx,++nx;
		}
		if (it->second==0){
			if (val==res){
				v.insert(mp(val^res,1));
				if (i==n)	flag=it->second+1>=m;
			}
		}
		else{
			v.insert(mp(val^res,(it->second)+1));
			if (i==n)	flag=it->second+1>=m;
		}
	}
	if (flag)	ans=res;
}
//}}}
int main()
{
	cin>>T;
	while (T--){
		cin>>n>>m;
		ans=0;
		for (int i=1;i<=n;++i){
			cin>>a[i];
			s[i]=s[i-1]^a[i];
		}

		for (int i=29;~i;--i)	solve(i);
		printf("%d\n",ans);
	}
	return 0;
}

如有哪里讲得不是很明白或是有错误，欢迎指正
如您喜欢的话不妨点个赞收藏一下吧