给出N个数的集合A ,找出一个子集S满足S中的所有元素异或后为0(mod 1000000007),求∑|S|
算每个元素出现在满足要求的集合的次数,也就是每个元素对答案的贡献。所以先求出线性基d1
对于线性基以外的元素,必然可以由线性基的异或和表示,所以完全为线性基外元素的异或和也仍然可以由线性基的异或和表示。所以,线性基外元素的贡献为
2
n
−
r
−
1
2^{n-r-1}
2n−r−1,共计n-r个。
而对于线性基内的元素,需要判断该元素是否可代替,也就是说该元素是否能插入剩余n-1个元素的线性基,只有插入时才会有贡献。因为不是所有线性基的元素都是用来表示其他数组中其他元素,有的线性基元素只是因为别人表示不了他,才被划进线性基里的。此时贡献为 2 n − r − 1 2^{n-r-1} 2n−r−1
因为在判断线性基内元素时要求一个新的线性基,而每次重新插入必然是超时的,现在考虑如何快速的求解。
从超时的插入中我们可以看出,非线性基元素是一定要插入的,所以不妨先求出非线性基元素的线性基d2,再在这个基础上插入线性基元素,构成所求的线性基d3,再向d3中插入该元素判断。
#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
#define ll long long
#define Pair pair<int,int>
#define re return
#define getLen(name,index) name[index].size()
#define mem(a,b) memset(a,b,sizeof(a))
#define Make(a,b) make_pair(a,b)
#define Push(num) push_back(num)
#define rep(index,star,finish) for(register int index=star;index<finish;index++)
#define drep(index,finish,star) for(register int index=finish;index>=star;index--)
using namespace std;
const int maxn=1e5+7;
const int mod=1e9+7;
template<class T> void _deb(const char *name,T val){
cout<<name<<val<<endl;
}
class linerBase{
private:
int len;
ll d[maxn];
public:
void ini(){
len=0;
mem(d,0);
}
bool insert(ll x){
for(register int i=61;i>=0;i--){
if(x&(1LL<<i)){
if(!d[i]){
d[i]=x;
break;
}
x=x^d[i];
}
}
if(x){
len++;
re false;
}else{
re true;
}
}
inline int size(){
re len;
}
};
ll Pow[maxn];
ll qpow(ll a, ll n)
{
if(Pow[n])
re Pow[n];
int record=n;
ll res = 1;
while(n)
{
if(n & 1)
res = (res * a) % mod;
n >>= 1;
a = (a * a) % mod;
}
return Pow[record]=res % mod;
}
linerBase d1,d2,d3;
int lenR,lenOt;
ll store[maxn],r[maxn],ot[maxn];
set<ll> LB;
int main(){
Pow[0]=1;
int n;
while(scanf("%d",&n)!=EOF){
lenR=0;
d1.ini(),d2.ini();
rep(i,0,n){
scanf("%lld",&store[i]);
}
rep(i,0,n){
ll &num=store[i];
if(d1.insert(num)){
d2.insert(num);
}else{
r[lenR]=num;
lenR++;
}
}
ll amount=(n-d1.size());
for(register int i=0;i<lenR;i++){
ll &tar=r[i];
d3=d2;
for(register int j=0;j<lenR;j++){
ll &ele=r[j];
if(tar != ele){
d3.insert(ele);
}
}
if(d3.insert(tar))
amount++;
}
ll ans=amount*qpow(2,n-d1.size()-1)%mod;
printf("%lld\n",ans);
}
re 0;
}
除此之外,还有更加快速的方法判断。在判断线性基内元素时,枚举所有可替代元素。但是我还没看懂,看懂了再补。