bzoj 2844: albus就是要第一个出场（线性基）

2844: albus就是要第一个出场

Time Limit: 6 Sec   Memory Limit: 128 MB
Submit: 1635   Solved: 680
[ Submit][ Status][ Discuss]

Description

已知一个长度为n的正整数序列A（下标从1开始）， 令 S = { x | 1 <= x <= n }, S 的幂集2^S定义为S 所有子

集构成的集合。定义映射 f : 2^S -> Zf(空集) = 0f(T) = XOR A[t] , 对于一切t属于T现在albus把2^S中每个集

合的f值计算出来， 从小到大排成一行， 记为序列B（下标从1开始）。 给定一个数， 那么这个数在序列B中第1

次出现时的下标是多少呢？

Input

第一行一个数n, 为序列A的长度。接下来一行n个数， 为序列A， 用空格隔开。最后一个数Q， 为给定的数.

Output

共一行， 一个整数， 为Q在序列B中第一次出现时的下标模10086的值.

Sample Input

3

1 2 3

1

Sample Output

3

题意：

给定n个数以及一个数q，将这n个数的所有子集（可以为空）的异或值从小到大排序得到序列B（这样B内就有2^n 个 元素），求出q在序列B中第一次出现的下标，答案对10086取膜（默认q一定存在）

思路：

这题要解决两个问题：

①B序列去重之后，q在B中是第几个

②B序列中每个数字出现多少次

对于问题①，我们已经知道了如何求出去重后B数组的第k小，而这题其实是反过来的给你一个数q，求出它是第几小

其实也很简单，把操作也反过来就行了，你甚至不用求出行最简式矩阵

步骤：1. 求出线性基后线性基内所有元素从小到大排序，并且只保留最高位

2. 之后将q二进制拆分，初始化rank=0，如果q的第i位为1，就看线性基中是否存在一个数x满足x = (1<<i)，如果存在，那么rank = rank^(x在线性基中排名)，最后的rank就是q的排名-1，为什么是排名-1呢因为还有个0没有算

对于问题②，定理：B数列中每个数字出现次数都是2^(n-cnt)！其中cnt是线性基的大小

有了这个定理问题②其实就解决了，这样原问题答案就是rank*2^(n-cnt)+1，一个快速幂就可以解决

证明如下：

还记得这道题么：Codeforces Round #448 (Div. 2): C. Square Subsets

这道题是求n个数有多少个集合满足异或为0，如果包含空集答案就是2^(n-cnt)

而线性基有个性质：线性基内任意集合异或结果唯一，而所有数异或0还是本身

所以每个数字出现的个数就等于(1*异或值为0的集合个数) = 2^(n-cnt)

证毕

#include<stdio.h>
#include<algorithm>
using namespace std;
#define mod 10086
int a[100005], p[35], jz[35];
int Pow(int x, int y)
{
	int ans = 1;
	while(y)
	{
		if(y%2)
			ans = ans*x%mod;
		x = x*x%mod;
		y /= 2;
	}
	return ans;
}
int main(void)
{
	int n, i, j, cnt, K, q;
	scanf("%d", &n);
	for(i=1;i<=n;i++)
		scanf("%d", &a[i]);
	for(i=1;i<=n;i++)
	{
		for(j=30;j>=0;j--)
		{
			if(a[i]&(1<<j))
			{
				if(p[j]==0)
				{
					p[j] = a[i];
					break;
				}
				else
					a[i] ^= p[j];
			}
		}
	}
	scanf("%d", &q);
	cnt = 0, K = 0;
	for(i=0;i<=30;i++)
	{
		if(p[i])
			jz[cnt++] = i;
	}
	for(i=0;i<=cnt-1;i++)
	{
		if(q&(1<<jz[i]))
			K ^= (1<<i);
	}
	printf("%d\n", (K%mod*Pow(2, n-cnt)+1)%mod);
	return 0;
}