思路:挺难的一道题,是队友教了才会的。首先要知道一个很重要的结论:一个长度为n从小到大排列的子序列中,两两最小的异或值一定在n-1个相邻对中出现(证明略)。
然后就不难想到O(N^2)的暴力dp,但是过不了,需要用01trie加速。
思路如下图:
代码:
#include<bits/stdc++.h>
using namespace std;
#define debug cout<<"!!!"<<endl;
#define pii pair<int,int>
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
const int N = 3e5+5, inf=0x3f3f3f3f;
int n, k, t[N*32][2], dp[N*32], pre[N*32], tot=0, ans[N*32];
pii a[N]; //{val,id}
inline void ins(int x){
int p=0;
ROF(i,30,0){
int c=(x>>i)&1; //当前位数字
if(!t[p][c]) t[p][c]=++tot; //创建新结点
p=t[p][c];
}
}
inline void ask(int x,int id){
int p = 0, res = 0;
ROF(i,30,0){
int c=(x>>i)&1; //(a[i]^k)的当前位
if (!((k>>i)&1) && t[p][c^1]){
if (dp[t[p][c ^ 1]] > res){
res = dp[t[p][c ^ 1]];
pre[id] = ans[t[p][c ^ 1]];
}
}
if (!t[p][c]) break; //走不下去了,结束
p = t[p][c];
}
if (!t[p][0] && !t[p][1] && res < dp[p]){ //如果走到终点,还有一个答案要更新
res = dp[p];
pre[id] = ans[p];
}
res++, p=0; //res+1(因为包括自己),然后从头开始把a[i]所有结果更新一遍
if(dp[0]<=res) dp[0]=res, ans[0]=id; //更新原点的dp值
ROF(i,30,0){
int c=((x^k)>>i)&1; //a[i]的当前位
p = t[p][c];
if(dp[p] <= res) dp[p]=res, ans[p]=id; //更新dp值,记录坐标
}
}
signed main(){
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cin>>n>>k;
FOR(i,1,n) cin>>a[i].first, a[i].second=i, ins(a[i].first);
if(k==0) {cout<<n<<'\n'; FOR(i,1,n) cout<<i<<' '; return 0;} //特判
sort(a+1,a+n+1); //按val排序
ROF(i,n,1) ask(a[i].first^k, a[i].second);
if(dp[0]==1) {cout<<-1<<'\n'; return 0;}
cout<<dp[0]<<'\n';
int id=ans[0];
while(id){
cout<<id<<' ';
id=pre[id];
}
return 0;
}