【题解】魔改版线性基

魔改版线性基解决此类问题。

联系线性空间的性质,我们直接可以构造出这样的基:

\[100000
\\
010000
\\
000010
\\
000001
\]

使得每个基的最高位是唯一的,我们的目的是要能够保证从上往下一直异或一直变大,所以不能使基出现这样的情况:

\[100001
\\
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异或和/魔改版线性基的更多相关文章

  1. BZOJ4671 异或图(容斥+线性基)

    题意 定义两个结点数相同的图 \(G_1\) 与图 \(G_2\) 的异或为一个新的图 \(G\) ,其中如果 \((u, v)\) 在 \(G_1\) 与 \(G_2\) 中的出现次数之和为 \(1 ...

  2. 【题解】 Codeforces 662A Gambling Nim (线性基)

    662A,戳我戳我 Solution: 我们先取\(ans=a[1] \bigoplus a[2] \bigoplus ... \bigoplus a[n]\),然后我们定义\(c[i]=a[i] \ ...

  3. bzoj 2844 albus就是要第一个出场 异或和出现次数 线性基

    题目链接 题意 给定\(n\)个数,将其所有的子集(\(2^n\)个)的异或和按升序排列.给出一个询问\(q\),问\(q\)在该序列中第一次出现位置的下标(下标从\(1\)开始). 题解 结论 记其 ...

  4. BZOJ4671 异或图 斯特林反演+线性基

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4671 题解 半年前刚学计数的时候对这道题怀着深深的景仰,现在终于可以来做这道题了. 类似于一般 ...

  5. P5169 xtq的异或和(FWT+线性基)

    传送门 我咋感觉我学啥都是白学-- 首先可以参考一下这一题,从中我们可以知道只要知道两点间任意一条路径以及整个图里所有环的线性基,就可以得知这两个点之间的所有路径的异或和 然而我好像并不会求线性基能张 ...

  6. 题解——洛谷P3812【模板】线性基

    学了下线性基 使用好像并不复杂 打了板子 但是要注意位运算优先级 #include <cstdio> #include <algorithm> #include <cst ...

  7. 【题解】 luogu 3857 [TJOI2008]彩灯 (线性基)

    luogu3857,懒得复制 Solution: 裸的线性基,往里面添加数,记录添加个数\(sum\),快速幂输出\(2^{sum}\)即可 Code: //It is coded by Ning_M ...

  8. 【题解】 bzoj4004: [JLOI2015]装备购买 (线性基)

    bzoj4004,戳我戳我 Solution: 裸的线性基,这没啥好说的,我们说说有意思的地方(就是我老是wa的地方) Attention: 这题在\(luogu\),上貌似不卡精度,\(bzoj\) ...

  9. hdu3949(线性基,求第k小的异或和

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3949 XOR Time Limit: 2000/1000 MS (Java/Others)    Me ...

随机推荐

  1. ACM的奇计淫巧_输入挂

    什么是输入挂? 众所周知scanf比cin快的多,那么有没有比scanf更快的东西呢?答案就是输入挂,输入挂利用了告诉读取的函数getchar(),然后再人工处理成整数或浮点,比使用scanf快太多. ...

  2. Linux内核网络栈源代码分析

    http://blog.csdn.net/column/details/linux-kernel-net.html

  3. 每天学一点Python(2)

    9月16日(python扩展的安装和使用) 接着上一篇继续.按照之前计划,先分析导出的数据,再做进一步统计. 导出的数据是html类型的,想到的处理方法有: 1.直接readlines然后一行一行找我 ...

  4. 在dedecms后台发表文章显示外部连接栏目

    问题描述:客户的网站,有个顶级栏目,下面包含了几个子栏目,这个顶级栏目不想发布什么内容,点击后进入他的某个子栏目就可以了,这时候把这个顶级栏目设置为“外部连接”就可以了 但是设置顶级栏目为外部连接后, ...

  5. php分页显示文章列表

    <div class="content"> <ul> <?php $querySel = "select * from news where ...

  6. 卸载django

    首先找到django安装路径 Python import sys; sys.path = sys.path[1:]; import django; print(django.__path__) 得到 ...

  7. 24. Spring Boot环境变量读取和属性对象的绑定【从零开始学Spring Boot】

    转:http://blog.csdn.net/linxingliang/article/details/52069509 凡是被spring管理的类,实现接口EnvironmentAware 重写方法 ...

  8. 每天一个 Linux 命令(57):ss命令

    ss是Socket Statistics的缩写.顾名思义,ss命令可以用来获取socket统计信息,它可以显示和netstat类似的内容.但ss的优势在于它能够显示更多更详细的有关TCP和连接状态的信 ...

  9. 浅谈对TDD的看法

    程序猿对TDD这个词一定不陌生.近几年比較火.英文全称Test-Driven Development,測试驱动开发. 它要求在编写某个功能的代码之前先编写測试代码,然后仅仅编写使測试通过的功能代码,通 ...

  10. 可以添加自定义的Select控件

    1.控件dom <select name="WebSiteTarget" id="WebSiteTarget" class="w1" ...