【题解】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 ...
随机推荐
- protobuf3 语法解析
定义一个消息类型 先来看一个非常简单的例子.假设你想定义一个“搜索请求”的消息格式,每一个请求含有一个查询字符串.你感兴趣的查询结果所在的页数,以及每一页多少条查询结果.可以采用如下的方式来定义消息类 ...
- lock与monitor的区别
1.Lock 只能对引用对象加锁 Lock锁定区间内可以对锁定值修改而不发生运行时错误,通常也会采用此种修改方式.这种方式又有点类同于使用Monitor.Wait取得资源,并对这个资源进行操作. 用法 ...
- 线程安全的单例模式还需要对成员变量的set get方法设置锁么
不需要,线程安全的单例模式,在获得对象时已经加锁了,保证每时每刻只有一个线程获得此单例对象.所以不需要再上锁了啊
- iOS 定位坐标不准确的相关整理及解决方案汇总
这几天处理定位相关的代码,彻彻底底的被火星坐标恶心到了. 恶心列表 从 CLLocationManager 取出来的经纬度放到 mapView 上显示,是错的! 从 CLLocationManager ...
- iphone/iOS 访问本地数据库sqlite3
Phone也支持访问本地数据库Sqlite 3.这里简单的介绍一下iPhone上Sqlite 3的使用方法. 首先需要在项目中引用Sqlite 3的开发包,下面是在iPhone SDK 3.0下的目录 ...
- Ajax方式实现注册验证代码
经常用到的经典Ajax代码,记录备用: function CreateAjax() { var XMLHttp; try { XMLHttp = new ActiveXObject("Mic ...
- linux系统故障分析与排查
在处理Linux系统出现的各种故障时,故障的症状是最先发现的,而导致这以故障的原因才是最终排除故障的关键.熟悉Linux系统的日志管理,了解常见故障的分析与解决办法,将有助于管理员快速定位故障点.“对 ...
- [翻译]MySQL 文档: Control Flow Functions(控制流函数)
本文翻译自13.4 Control Flow Functions Table 13.6 Flow Control Operators 名称 描述 CASE Case 运算符 IF() if/else ...
- python:如何判断字符串中的内容是否都为数字并且把字符串转换为数字
使用str.isdigit();有两种使用方法 str.isdigit('12345') =====>True str.isdigit('aaaaa')======>False 或者 '1 ...
- 3DES
3DES是继DESeasy被破解后的DES加密升级版,它属于对称加密. 可指定24位长度的密钥.在java API中也有事实上现,代码例如以下: /** * 3DES 的Java SDK API 实现 ...