应该可以看出来是做了 n n n次异或前缀和。那么 a i a_i ai对 b j b_j bj的贡献相当于从 ( 1 , i ) (1,i) (1,i)走到 ( j , n ) (j,n) (j,n)的路径方案数 ( j + n − i − 1 j − 1 ) \binom{j+n-i-1}{j-1} (j−1j+n−i−1),有贡献当且仅当 ( n − i ) ∩ ( j − 1 ) = 0 (n-i)\cap (j-1)=0 (n−i)∩(j−1)=0 。将 a a a序列反转就能得到 a i a_i ai对 b j b_j bj有贡献当且仅当 i ∩ j = 0 i\cap j=0 i∩j=0。
方便起见,我们让 n = 2 k n=2^k n=2k,可以通过在 a a a序列末尾补零实现。我们令 i = U − i ( U = 2 k − 1 ) i=U-i(U=2^k-1) i=U−i(U=2k−1),那么 b i b_i bi就等于 a i a_i ai的子集和。因为是异或所以做 F M T FMT FMT的时候不需要考虑 − 1 -1 −1的系数,因此 a i a_i ai也等于 b i b_i bi的子集和。
然而对于 b i b_i bi一段前缀的值我们是不知道的。这就比较麻烦了。
我们不妨将问题做一些泛化:
1.1 1.1 1.1 有两个自然数 n n n, t t t满足 n ≤ 2 t n\le 2^t n≤2t(不满足 t t t最小了)
1.2 1.2 1.2 已知 a n , a n + 1 , . . . , a 2 t − 1 a_n,a_{n+1,}...,a_{2^t-1} an,an+1,...,a2t−1(不限制为 0 0 0了)以及 b 2 t − n , b 2 t − n + 1 , . . . , b 2 t − 1 b_{2^t-n},b_{2^t-n+1},...,b_{2^{t}-1} b2t−n,b2t−n+1,...,b2t−1的值,并且满足 a i a_i ai, b i b_i bi互为子集异或和,问是否存在合法解。
考虑到 a a a, b b b两边剩下 n n n, 2 t − n 2^t-n 2t−n个数不知道,一定存在一边 ≤ 2 t − 1 \le 2^{t-1} ≤2t−1。
如果我们能把其中一边全部解出来,那么另一边做一个子集和就能得到答案。
假设 n ≤ 2 t − 1 n\le 2^{t-1} n≤2t−1。只需要消除 a a a序列 [ 2 t − 1 , 2 t − 1 ] [2^{t-1},2^{t}-1] [2t−1,2t−1]对 b b b序列 [ 2 t − 1 , 2 t − 1 ] [2^{t-1},2^{t}-1] [2t−1,2t−1]的贡献,然后解 a a a序列 [ 0 , 2 t − 1 − 1 ] [0,2^{t-1}-1] [0,2t−1−1]和 b b b序列 [ 2 t , 2 t − 1 ] [2^{t},2^t-1] [2t,2t−1]这两个部分即可。因为此时 b b b序列是 a a a序列的子集和,所以 a a a序列作逆变换还是可以得到 b b b。 n > 2 t − 1 n>2^{t-1} n>2t−1的情况类似。
复杂度 O ( n log n ) O(n\log n) O(nlogn)。
ps:求子集的异或和可以用下面这一段代码实现。
void conv(int *a,int *b,int m){
int n=1<<m;for(int i=0;i<n;i++)b[i]=a[i];
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
if(j>>i&1)b[j]^=b[j^(1<<i)];
}
}
}
如果有bug记得联系博主。
#include<bits/stdc++.h>
#define fi first
#define se second
#define ll long long
#define pb push_back
#define db double
#define inf 0x3f3f3f3f3f3f3f3f
using namespace std;
const int N=1e6+5;
int n,b[N],a[N],c[N],m;
void conv(int *a,int *b,int m){
int n=1<<m;for(int i=0;i<n;i++)b[i]=a[i];
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
if(j>>i&1)b[j]^=b[j^(1<<i)];
}
}
}
void solve(int *a,int *b,int m,int n){
if(n==(1<<m)){
conv(b,a,m);
return;
}
if(n==0){
conv(a,b,m);
return;
}
if(n<=(1<<m-1)){
conv(a+(1<<m-1),c,m-1);
for(int i=(1<<m)-n;i<(1<<m);i++){
b[i]^=c[i-(1<<m-1)];
}solve(a,b+(1<<m-1),m-1,n);
conv(a,b,m);
}
else {
conv(b+(1<<m-1),c,m-1);
for(int i=n;i<(1<<m);i++){
a[i]^=c[i-(1<<m-1)];
}solve(a+(1<<m-1),b,m-1,n-(1<<m-1));
conv(b,a,m);
}
}
int main(){
ios::sync_with_stdio(false);
cin.tie(0),cout.tie(0);
cin>>n;int m=0;while((1<<m)<n)m++;
for(int i=(1<<m)-1;i>=(1<<m)-n;i--)cin>>b[i];
solve(a,b,m,n);for(int i=n-1;i>=0;i--)cout<<a[i]<<" ";
}