「CodePlus 2017 11 月赛」可做题

题目描述
qmqmqm 希望给 sublinekelzrip 出一道可做题。于是他想到了这么一道题目：给一个长度为n的非负整数序列ai​​，你需要计算其异或前缀和bi，满足条件b1=a1​​，bi=bi−1 xor ai(i≥2)。
但是由于数据生成器出现了问题，他生成的序列a的长度特别长，并且由于内存空间不足，一部分ai已经丢失了，只剩余m个位置的元素已知。现在 qmqmqm 找到你，希望你根据剩余的ai​​，计算出所有可能的a序列对应的b序列中∑ni=1bi的最小值。

输入格式
输入第一行两个非负整数n、m，分别表示原始序列a的长度及剩余元素的个数。
之后m行，每行2个数i、ai，表示一个剩余元素的位置和数值。

输出格式
输出一个整数表示可能的最小值。

样例输入
5 3
4 0
3 7
5 0

样例输出
7

样例解释
已知的a序列为：X,X,7,0,0，其中X表示这个位置丢失了。一种可能的a序列为0,7,7,0,0，对应的b序列为0,7,0,0,0，和最小为7。可以证明不存在和更小的情况。

数据范围与提示
1≤n≤109​​，0≤m≤min{n,105}，0≤ai≤109
注意未知的ai可以超过已知ai的范围。
保证输入中所有的i不同，且满足1≤i≤n。
来自 CodePlus 2017 11 月赛，清华大学计算机科学与技术系学生算法与竞赛协会 荣誉出品。
Credit：idea/卢政荣　命题/卢政荣　验题/何昊天
Git Repo：https://git.thusaac.org/publish/CodePlus201711
感谢腾讯公司对此次比赛的支持。

思路
首先想想序列中只有0和1的情况，可以很容易想到dp的做法：fi,j表示第i位填j的代价。转移自己想一想就好了。时间复杂度O(n)。
但是i很大时，数组是存不下的，这个时候可以发现：在连续的一段长度大于2的空白区间内，最优方案一定是除了头和尾以外其他的位都为0，对答案造成的贡献也是0，这样可以把所有的长度大于2的区间都变成长度为2的区间。时间复杂度O(m)。
再来想想序列中不只有0和1的情况：容易看出，异或的每一二进制位都是独立的，那么可以把数字的二进制位拆开，分开做一遍上述算法。时间复杂度O(mlog(max{ai}))。

代码

#include <cstdio> #include <algorithm>const int maxn=100000; const int inf=0x3f3f3f3f;struct data {int x,y;bool operator <(const data &other) const{return x<other.x;} };int n,m,p[maxn<<3],f[maxn<<3][2]; long long ans; data d[maxn+10];inline int read() {int x=0,f=1;char ch=getchar();while((ch<'0')||(ch>'9')){if(ch=='-'){f=-f;}ch=getchar();}while((ch>='0')&&(ch<='9')){x=x*10+ch-'0';ch=getchar();}return x*f; }int main() {n=read();m=read();for(register int i=1; i<=m; ++i){d[i].x=read();d[i].y=read();}std::sort(d+1,d+m+1);d[0].x=0;d[0].y=0;d[m+1].x=n+2;d[m+1].y=0;n=0;for(register int i=1; i<=m+1; ++i)//将每个大于2的间隔都压成2{if(d[i].x-d[i-1].x==1){++n;p[n]=d[i].y;}else if(d[i].x-d[i-1].x==2){++n;p[n]=-1;++n;p[n]=d[i].y;}else{++n;p[n]=-1;++n;p[n]=-1;++n;p[n]=d[i].y;}}for(register int k=30; k>=0; --k)//枚举每个二进制位{f[0][0]=0;f[0][1]=inf;for(register int i=1; i<=n; ++i)//dp{if(p[i]==-1){f[i][0]=std::min(f[i-1][0],f[i-1][1]);f[i][1]=std::min(f[i-1][1],f[i-1][0])+1;}else if(p[i]&(1<<k)){f[i][0]=f[i-1][1];f[i][1]=f[i-1][0]+1;}else{f[i][0]=f[i-1][0];f[i][1]=f[i-1][1]+1;}}ans=(ans<<1)+std::min(f[n][0],f[n][1]);}printf("%I64d\n",ans);return 0; }

转载于:https://www.cnblogs.com/Canopus-wym/p/10376281.html

总结

以上是生活随笔为你收集整理的「CodePlus 2017 11 月赛」可做题的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。