线性基学习

拖了好久,还是看掉吧QAQ

果然还是这篇讲的最好了:https://www.cnblogs.com/ljh2000-jump/p/5869991.html

以下是转载于上面这篇神犇的博客

1、线性基:

  若干数的线性基是一组数a1,a2,...ana1,a2,...an,其中axax的最高位的11在第xx位。

  通过线性基中元素xorxor出的数的值域与原来的数xorxor出数的值域相同。

 

2、线性基的构造法:

  对每一个数pp从高位到低位扫,扫到第xx位为11时,若axax不存在,则ax=pax=p并结束此数的扫描,否则令p=pp=p   xorxor ax。ax。

 

3、查询:

  用线性基求这组数xorxor出的最大值:从高往低扫axax,若异或上axax使答案变大,则异或。

 

4、判断:

  用线性基求一个数能否被xorxor出:从高到低,对该数每个是11的位置xx,将这个数异或上axax(注意异或后这个数为1的位置和原数就不一样了),若最终变为00,则可被异或出。当然需要特判00(在构造过程中看是否有p变为0即可)。例子:(11111,10001)(11111,10001)的线性基是a5=11111a5=11111,a4=01110a4=01110,要判断1111111111能否被xorxor出,1111111111 xorxor a5a5=0=0,则这个数后来就没有是11的位置了,最终得到结果为00,说明1111111111能被xorxor出。

 

 

个人谈一谈对线性基的理解:

  很多情况下,只有有关异或运算和求最值,就可以用到线性基。线性基有很多很好的性质,比如说如果有很多个数,我们可以构出这些数的线性基,那么这个线性基可以通过互相xorxor,能够构出原来的数可以相互xorxor构出的所有的数。所以可以大大减少判断的时间和次数。同时线性基的任何一个非空子集都不会使得其xorxor和为0,证明也很简单,反证法就可以说明。这个性质在很多题目中可以保证算法合法性,比如:BZOJ2460BZOJ2460。

 

  构造的方法有点像贪心,从大到小保证高位更大。也比较好理解。就是这几行代码:

  

1

2

3

4

5

6

7

8

9

10

11

12

13

for(int i=1;i<=n;i++) {    

 

        for(int j=62;j>=0;j--) {

 

             if(!(a[i]>>j)) continue;//对线性基的这一位没有贡献           

 

               if(!p[j]) { p[j]=a[i]; break; }//选入线性基中                   

 

               a[i]^=p[j];

 

             }

 

       }

 

线性基性质:

  1. 线性基能相互异或得到原集合的所有相互异或得到的值。

  2. 线性基是满足性质1的最小的集合

  3. 线性基没有异或和为0的非空子集。

下面是例题:

1.       Codeforces 1101G:https://codeforces.com/contest/1101/problem/G 

题意:将一个序列分成许多段,每段的值为元素异或和,求使得不存在非空子集异或和为0的最大段数。

题解:

将前缀异或和求出来,然后求这个前缀异或和的线性基的大小即可。

代码:

#include<bits/stdc++.h>
#define LL long long
#define MAXN 200010
using namespace std;
LL a[MAXN],s[MAXN],y[35];
LL XXJ(LL n)
{
	LL i,j,ans;
	memset(y,0,sizeof(y));
	for(i=1;i<=n;i++)
	{
		for(j=30LL;j>=0LL;j--)
		{
			if(((s[i]>>j)&1LL)==1LL)
			{
				if(y[j]==0LL)
				{
					y[j]=s[i];
					break;
				}
				else s[i]^=y[j];
			}
		}
	}
	ans=0LL;
	for(j=0LL;j<=30LL;j++)if(y[j]!=0LL)ans++;
	return ans;
}
int main()
{
	LL n,i;
	scanf("%lld",&n);
	for(i=1;i<=n;i++)scanf("%lld",&a[i]);
	s[0]=0LL;for(i=1;i<=n;i++)s[i]=s[i-1]^a[i];
	if(s[n]==0LL)printf("-1");
	else printf("%lld",XXJ(n));
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值