题解 - 牛客H Playing games (想法+Fwt)
题目链接:https://ac.nowcoder.com/acm/contest/146/H
题意
给出N个数字,现在要求取出最少的数字,使其数组剩下的值异或和为0
数据范围: 1 ≤ N , V a l ≤ 5 ∗ 1 0 5 1\le N,Val \le 5*10^5 1≤N,Val≤5∗105
思路
根据线性基的思想,我们不难想到,最多取出19个数字,一定能够使得剩下数字异或和为零。
这个我们可以直接从1扫到19,或者二分取出 M M M个数字。
那么我们现在的问题就是,怎么判断取出了 M M M个数字后其可以组成异或的值。
令 F [ i ] F[i] F[i]代表,组成异或值为i的方案数,那么不难得到,我们需要求的是
F [ i ] = ∑ j ∧ k = i ( F [ j ] + F [ k ] ) F[i] = \sum_{j\land k= i} (F[j]+F[k]) F[i]=∑j∧k=i(F[j]+F[k])
这个式子,我们要求 M M M次,那么我们只需要将其FWT一下,再将其得到的 f p o w ( F ′ , M ) fpow(F',M) fpow(F′,M),即可。
要检验时,就在将其UFWT变成正常的样子,看其 F [ s u m ] F[sum] F[sum]是否有值
代码
#include <bits/stdc++.h>
using namespace std;
#define rep(i,j,k) for(ll i = (ll)j;i <= (ll)k;i ++)
#define debug(x) cerr<<#x<<":"<<x<<endl
#define pb push_back
typedef long long ll;
typedef pair<ll,ll> pi;
const ll MAXN = (ll)1e6+7;
const ll MOD = (ll)1e9+7,inv2 = (MOD+1)/2;
ll fpow(ll a,ll x) {
ll res = 1;
while (x) {
if (x&1) res = 1LL*res*a%MOD;
x >>= 1;
a = 1LL*a*a%MOD;
}
return res;
}
void fwtxor(ll f[],ll n,ll op){
for(ll p=2;p<=n;p<<=1){
ll len=p>>1;
for(ll k=0;k<n;k+=p)
for(ll i=k;i<k+len;++i){
ll t=f[i+len];
f[i+len]=(f[i]-f[i+len]+MOD)%MOD;
f[i]=(f[i]+t)%MOD;
}
if(op==-1) for(ll i=0;i<n;++i) f[i]=(ll)f[i]*inv2%MOD;
}
}
ll f[MAXN],a[MAXN];
ll N,sum,num=1<<19;
bool OK(ll m) {
rep(i,0,num-1) {
f[i] = fpow(a[i],m);
}
fwtxor(f,num,-1);
f[sum] = (f[sum]+MOD)%MOD;
return f[sum]>0;
}
int main()
{
scanf("%lld",&N);
rep(i,1,N) {
ll tmp;
scanf("%lld",&tmp);
a[tmp] ++;
sum ^= tmp;
}
a[0] ++;
fwtxor(a,num,1);
ll l = 0,r = 19;
while (l <= r) {
ll m = l+r>>1;
if (OK(m)) r = m-1;
else l = m+1;
}
ll ans = N-l;
printf("%lld\n",ans);
}