CF1163E Magical Permutation（线性基，构造）

虽然做起来有一点裸……但是就是想不到啊……

首先令 $d_i=p_i\oplus p_{i-1}$，那么 $d_i$ 都是 $S$ 中的数，$a_i=d_i\oplus d_{i-1}\oplus \cdots\oplus d_2$。也就是每个数都能被表示成 $S$ 的某个子集的异或和。

要用 $S$ 表示出 $1$ 到 $2^x-1$ 的所有数（不用考虑 $0$，因为每个数是可以重复用的，可以 $S_i\oplus S_i=0$）。怎么求出最大的 $x$？

其实就是建出线性基，然后最小的没有数的位就是 $x$ 了。为什么？当 $0$ 到 $x-1$ 都有数时是可以表示出所有 $0$ 到 $2^x-1$ 的，而当 $x$ 没有数时无法填第 $x$ 位。

（想用严谨一点的语言的……然而实力不允许……）

然后，求出 $x$ 后如何构造排列？

首先有最原始的想法：DFS，每次 $O(n)$ 枚举下一个数是这一个数异或上啥。肯定不可能过。

然后发现只需要保留线性无关的最大子集（可以在建线性基的过程中就完成），此时肯定还是能表示出全部的数（线性无关的定义）。同时数的个数降到了 $O(\log)$。

看起来还是不能过，但是……它的复杂度是对的！为什么？

DFS 过程中要判断 $vis[x\oplus S[i]]$，这就是在原始想法中耗费大量无用时间的原因。

而出现这样的情况，当且仅当出现了区间异或和为 $0$ 的情况。

但是对于这个线性无关的子集，会出现这样的情况，当且仅当当前枚举的 $S_i$ 等于目前的某一个后缀异或和。因为既然线性无关，所以所有不可重子集的异或和都不为 $0$。

由于 $S_i$ 只有 $O(\log)$ 种，所以也只会有 $O(\log)$ 次无法递归。

那么复杂度就对了。$O(n\log)$。下面的代码实现比较丑，所以复杂度似乎是两个 $\log$。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=266666;
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define ROF(i,a,b) for(int i=(a);i>=(b);i--)
#define MEM(x,v) memset(x,v,sizeof(x))
inline int read(){
    int x=0,f=0;char ch=getchar();
    while(ch<'0' || ch>'9') f|=ch=='-',ch=getchar();
    while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
    return f?-x:x;
}
int n,s[maxn],lb[20],vec[20],vl,ans,seq[maxn],sl;
bool vis[maxn],ok;
void insert(int x){
    int tmp=x;
    ROF(i,19,0) if((x>>i)&1){
        if(!lb[i]){
            lb[i]=x;
            vec[++vl]=tmp;
            return;
        }
        else x^=lb[i];
    }
}
bool check(int x){
    FOR(i,0,x-1) if(!lb[i]) return false;
    return true;
}
void dfs(int x){
    vis[x]=true;
    seq[++sl]=x;
    if(sl==(1<<ans)){ok=true;return;}
    FOR(i,1,vl) if(!vis[x^vec[i]]){
        dfs(x^vec[i]);
        if(ok) return;
    }
    vis[x]=false;
    sl--;
}
int main(){
    n=read();
    FOR(i,1,n) s[i]=read();
    ROF(_,19,0){
        MEM(lb,0);MEM(vec,0);vl=0;
        FOR(i,1,n) if(s[i]<=(1<<_)-1) insert(s[i]);
        if(check(_)){ans=_;break;}
    }
    printf("%d\n",ans);
    dfs(0);
    FOR(i,1,sl) printf("%d ",seq[i]);
}

View Code

 

转载于:https://www.cnblogs.com/1000Suns/p/11126028.html