参照给出的题解算法:
ki | gcd(xi , v) 该条件成立,即xi %ki ==0 && v %ki ==0 , 若 xi % ki != 0 则无解,当有解时,只需寻找a中的 v,使得 v % ki = 0 ,且 v + xi <= si ,同时最大化 v xor xi;
首先我们考虑 v % ki = 0 ,且 v + xi <= si ,因为ki,v的范围1至1e5,可以建立1e5个独立的数据结构,每个数据集collection[k] , 代表当前数组a内能够被k整除的所有数。
那么这个数据结构如何选取,首先必须是动态申请内存,因为每个数据集都可能达到1e5的存储量,再者,考虑到第三个条件,最大化 v xor xi,我们用trie树存储,可以做到每次log(1e5)的查询,为了使得查询方便trie树每个节点要存储一个当前所在子树最小值。
注意c++的初始化清零,不是递归的深层清零。所以注意,结构体内的变量,特别是指针手动清零,不然指针可能指向随机方向。
使用指针版本的trie要注意根节点使用指针比较好,这样在节点往下走时,能将第一步与下面做到一致
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn = 1e5 + 5;
struct node{
int min_;
node *bit[2];
node():min_(maxn){
bit[0] = (NULL),bit[1] = (NULL);
}
}*a[maxn];
void add(int x ,int u){
if(a[x] ==NULL) a[x] = new node();
node* p = a[x];
p->min_ = min(p->min_ , u);
for(int i = 18; i>=0 ; i--){
int ibit = ((u>>i) & 1);
if(p->bit[ibit] == NULL){
p->bit[ibit] = new node();
}
p = p->bit[ibit];
p->min_ = min(p->min_ , u);
}
}
// xi + vi <= si , max -> xi^v
int find_vi(int xi ,int si ,int ki){
node* p = a[ki];
int lim = si - xi , flag ;
if( p == NULL || lim < p->min_ ) return -1;
int ans = 0;
for(int i=18;i>=0;i--){
int now = 0 , xi_i = (((xi>>i)&1)<<i);
if(p->bit[0] == NULL){
p = p->bit[1];
continue;
}
if(p->bit[1] == NULL){
p = p->bit[0];
continue;
}
if(xi_i > 0){
ans += xi_i , p = p->bit[0];
continue;
}
if(xi_i == 0 && p->bit[1]->min_ <= lim){
ans += xi_i, p = p->bit[1];
} else p=p->bit[0];
}
return p->min_;
}
bool vis[maxn];
vector<int> st[maxn];
void init(){
for(int i=1;i<maxn ; i++)
for(int j = i ; j<maxn ; j+=i){
st[j].push_back(i);
}
}
int main()
{
init();
int n;
cin>>n;
for(int i=1;i<=n;i++){
int Q;
scanf("%d",&Q);
if(Q == 1){
int val;
scanf("%d",&val);
if(vis[val]) continue;
vis[val] = true;
for(int j = 0; j< st[val].size() ; j++){
add(st[val][j] , val);
}
} else {
int xi , si , ki ;
scanf("%d %d %d",&xi,&ki,&si);
if(xi % ki == 0)
printf("%d\n",find_vi(xi , si , ki));
else printf("-1\n");
}
}
return 0;
}