BZOJ2844:albus就是要第一个出场——题解

https://www.lydsy.com/JudgeOnline/problem.php?id=2844

已知一个长度为n的正整数序列A(下标从1开始), 令 S = { x | 1 <= x <= n }, S 的幂集2^S定义为S 所有子集构成的集合。
定义映射 f : 2^S -> Z
f(空集) = 0
f(T) = XOR A[t] (对于一切t属于T)
现在albus把2^S中每个集合的f值计算出来, 从小到大排成一行, 记为序列B(下标从1开始)。
给定一个数, 那么这个数在序列B中第1次出现时的下标是多少呢?

这话真的不是人说出来的。

简单翻译一下题就是序列所有的子集(可为空)的异或和经过排序,给一个q,求q在排序后序列中的位置下标。

参考:https://blog.sengxian.com/algorithms/linear-basis

如果是无重的,我们按照类似HDU3949:XOR的做法二进制分解q即可做(看代码应该更好理解吧)。

但是有重集合就很尴尬。

不过我们有对于一个数x,它出现次数一定是2^(n-size)(size为线性基大小)。

详细证明可以看参考,这里给出我的想法。

显然对于原线性基有且仅有一种方法构造出来x,且有且仅有一种方法构造出来非线性基内但却在a序列中的数。

那么对于非线性基内但却在a序列中的数有(n-size)个,它们和线性基内的数异或为0的方法有(n-size),排列组合就有2^(n-size)个,再与仅有一种方法构造出来x,就有2^(n-size)个了。

于是我知道了它去重前的排名,它排名前面的数都乘上2^(n-size)后+1即为答案。

PS:自己曾经也纠结过第一次算ans是否需要减一,但后来发现集合可为空……emmm这不值就有0了吗,对于0也满足这个性质啊,但是ans算排名的时候忽略了0啊。所以相当于ans-1+1=ans了啊。

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<cmath>
#include<cctype>
#include<algorithm>
using namespace std;
const int N=100010;
const int p=10086;
const int BASE=30;
inline int read(){
    int X=0,w=0;char ch=0;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}
int qpow(int k,int n){
    int res=1;
    while(n){
    if(n&1)res=res*k%p;
    k=k*k%p;
    n>>=1;
    }
    return res;
}
int a[N],b[BASE+4],cnt;
vector<int>mp;
int main(){
    int n=read();
    for(int i=1;i<=n;i++){
    a[i]=read();
    for(int j=BASE;j>=0;j--){
        if(a[i]>>j&1){
        if(b[j])a[i]^=b[j];
        else{
            b[j]=a[i];cnt++;
            break;
        }
        }
    }
    }
    int q=read(),ans=0;
    for(int i=0;i<=BASE;i++)
    if(b[i])mp.push_back(i);
    for(int i=0;i<mp.size();i++){
    if(q>>mp[i]&1)ans+=1<<i;
    ans%=p;
    }
    ans=ans*qpow(2,n-cnt)%p+1;
    printf("%d\n",ans%p);
    return 0;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者:luyouqi233。               +

+欢迎访问我的博客:http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++

posted @ 2018-04-13 11:23  luyouqi233  阅读(251)  评论(0编辑  收藏  举报