【题解】kth异或和/魔改版线性基
【题解】魔改版线性基
魔改版线性基解决此类问题。
联系线性空间的性质,我们直接可以构造出这样的基:
\\
010000
\\
000010
\\
000001
\]
使得每个基的最高位是唯一的,我们的目的是要能够保证从上往下一直异或一直变大,所以不能使基出现这样的情况:
\\
000001
\]
一个不能从上往下一直异或一直变大的例子。
考虑如何构造\(kth\) 大(小),考虑这样的性质,我们记\(a_i\)表示从下往上第\(i\)个基,显然从\(0\)开始,如果我们异或了\(a_x\),那么我们可以保证我们会比这个线性空间内\(2^{x-1}\)个值都要大,直接考虑从上往下查询,类似树上倍增。
考虑无解的情况,显然当我们发现\(k \ge 2^{|E|}\)时就无解了。
然而我们有一个问题,那就是我们默认\(0\)是可以被表示出来的,然而我们知道,当给定元素全部互为不相关时,这个时候就不存在一个\(0\)会被表示出来了。那么我们可以特判一下即可 。
好,那么如何构造从上往下一直异或一直变大的基呢?魔改一下就好了。具体看代码。
分析复杂度:我们把魔改操作的\(log^2\)看做常数吧...那就是\(O(nlogn)\)
当然你如果不承认魔改操作是常数,那就继续魔改\(upd\)函数,也可以做到\(nlogn\)。具体就是加入的时候从小往大加。
下面的第\(k\)小,转换为第\(k\)大用\(2^{|E|}\)减去即可。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
template <class ccf>
inline ccf qr(ccf ret) {
ret = 0;
register char c = getchar();
while (c < 48 || c > 57) c = getchar();
while (c >= 48 && c <= 57) ret = ret * 10 + c - 48, c = getchar();
return ret;
}
namespace bs {
ll num[51], sav[51];
struct BASE {
ll data[51];
int cnt;
ll p;
BASE() {
memset(data, 0, sizeof data);
cnt = 0;
p = 1;
}
inline ll size() { return 1LL << cnt; }
inline void set() {
for (register int t = 1; t <= 50; ++t)
if (data[t])
for (register int i = t + 1; i <= 50; ++i)
if (data[i] & num[t])
data[i] ^= data[t];
}
inline void upd(ll x) {
for (register int t = 50; t >= 1; --t) {
if (x & num[t]) {
if (not data[t]) {
data[t] = x;
++cnt;
break;
} else
x ^= data[t];
}
if(not x) {
p=0;
break;
}
}
}
inline ll que(ll k) {
register ll ret = 0, sav = 0;
if(cnt==50) --k;
if (k > this->size())
return -1;
for (register int t = 50, f = 0; t >= 1; --t) {
if (data[t]) {
++f;
if (sav + (1LL << (cnt - f)) <= k)
sav += 1LL << (cnt - f), ret = ret ^ data[t];
}
}
return ret;
}
inline void clear() {
memset(data, 0, sizeof data);
cnt = 0;
}
inline void debug() {
for (register int t = 50; t >= 1; --t) {
if (data[t])
cout << ((bitset<10>)data[t]) << ' ' << data[t] << endl;
}
cout << endl;
}
};
inline void init() {
for (num[0] = 2, num[1] = 1; num[0] <= 50; ++num[0]) num[num[0]] = num[num[0] - 1] << 1;
}
}
bs::BASE qaq;
int main() {
bs::init();
int n = qr(1);
for (register int t = 1; t <= n; ++t) qaq.upd(qr(1ll));
int m = qr(1);
qaq.set();
for (register int t = 1; t <= m; ++t) printf("%lld\n", qaq.que(qr(1ll)));
return 0;
}
【题解】kth异或和/魔改版线性基的更多相关文章
- BZOJ4671 异或图(容斥+线性基)
题意 定义两个结点数相同的图 \(G_1\) 与图 \(G_2\) 的异或为一个新的图 \(G\) ,其中如果 \((u, v)\) 在 \(G_1\) 与 \(G_2\) 中的出现次数之和为 \(1 ...
- 【题解】 Codeforces 662A Gambling Nim (线性基)
662A,戳我戳我 Solution: 我们先取\(ans=a[1] \bigoplus a[2] \bigoplus ... \bigoplus a[n]\),然后我们定义\(c[i]=a[i] \ ...
- bzoj 2844 albus就是要第一个出场 异或和出现次数 线性基
题目链接 题意 给定\(n\)个数,将其所有的子集(\(2^n\)个)的异或和按升序排列.给出一个询问\(q\),问\(q\)在该序列中第一次出现位置的下标(下标从\(1\)开始). 题解 结论 记其 ...
- BZOJ4671 异或图 斯特林反演+线性基
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4671 题解 半年前刚学计数的时候对这道题怀着深深的景仰,现在终于可以来做这道题了. 类似于一般 ...
- P5169 xtq的异或和(FWT+线性基)
传送门 我咋感觉我学啥都是白学-- 首先可以参考一下这一题,从中我们可以知道只要知道两点间任意一条路径以及整个图里所有环的线性基,就可以得知这两个点之间的所有路径的异或和 然而我好像并不会求线性基能张 ...
- 题解——洛谷P3812【模板】线性基
学了下线性基 使用好像并不复杂 打了板子 但是要注意位运算优先级 #include <cstdio> #include <algorithm> #include <cst ...
- 【题解】 luogu 3857 [TJOI2008]彩灯 (线性基)
luogu3857,懒得复制 Solution: 裸的线性基,往里面添加数,记录添加个数\(sum\),快速幂输出\(2^{sum}\)即可 Code: //It is coded by Ning_M ...
- 【题解】 bzoj4004: [JLOI2015]装备购买 (线性基)
bzoj4004,戳我戳我 Solution: 裸的线性基,这没啥好说的,我们说说有意思的地方(就是我老是wa的地方) Attention: 这题在\(luogu\),上貌似不卡精度,\(bzoj\) ...
- hdu3949(线性基,求第k小的异或和
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3949 XOR Time Limit: 2000/1000 MS (Java/Others) Me ...
随机推荐
- bzoj 5125: [Lydsy1712月赛]小Q的书架
新学了一波 决策单调性 dp 套路.... 这种dp一般是长这样的 => f[i][j] = max/min { f[i-1][k] + cost(k+1,j)} ,其中cost函数满足四边形 ...
- Factory Method 和AbstractFactory
对应慕课视频的连接:https://www.imooc.com/video/5316 1,工厂模式的应用场景 有一组类似的对象需要被创建 在编码时不能预见需要被创建哪种类的实例 在系统需要考虑扩展性的 ...
- Delphi中@,^,#,$分别表示什么?
@:取址运算符; var int:integer; p:^integer; new(P); int:=24; p:=@int; dispose(P); ^:指针的引用解析操作符; var pint:^ ...
- PCA原理
http://blog.csdn.net/shizhixin/article/details/51181379
- Codis的安装
其他环境准备: 安装JDK,安装Zookeeper 1.创建codis账户 useradd codis passwd codis 2.解压codis3.1.3-go1.7.4-linux.tar.gz ...
- 【前端阅读】——《程序员思维修炼》摘记&读后感&思维导图
前言:这是一本介绍如何用脑的书,并从思维的角度(以程序员为例),介绍如何从新手成为专家.作者带领着读者(我)共同经历一次有关认知科学.神经学.学习和行为理论的旅程,探索人类大脑令人 惊奇的工作的机制, ...
- POJ 3254 Corn Fields 状态压缩
这题对我真的非常难.实在做不出来,就去百度了,搜到了一种状压DP的方法.这是第一种 详细见凝视 #include <cstdio> #include <cstring> #in ...
- JAVA Eclipse开发Android如何设置滚动条最大值最小值
最小值默认为0,你最好在实现逻辑中修改 最大值为max 初始值为progress <SeekBar android:id="@+id/seekBarSpeedMov ...
- hdu1394 Minimum Inversion Number(最小逆序数)
Minimum Inversion Number Time Limit : 2000/1000ms (Java/Other) Memory Limit : 65536/32768K (Java/O ...
- Ant Design 3.0 使用案例
代码地址如下:http://www.demodashi.com/demo/12309.html 本文适合对象 有过React使用经验. 有过webpack使用经验. 了解node. DEMO使用方式 ...