首页 技术 正文
技术 2022年11月15日
0 收藏 408 点赞 3,780 浏览 3785 个字

Kuro is currently playing an educational game about numbers. The game focuses on the greatest common divisor (GCD), the XOR value, and the sum of two numbers. Kuro loves the game so much that he solves levels by levels day by day.

Sadly, he’s going on a vacation for a day, and he isn’t able to continue his solving streak on his own. As Katie is a reliable person, Kuro kindly asked her to come to his house on this day to play the game for him.

Initally, there is an empty array aa. The game consists of qq tasks of two types. The first type asks Katie to add a number uiui to aa. The second type asks Katie to find a number vv existing in aa such that ki∣GCD(xi,v)ki∣GCD(xi,v), xi+v≤sixi+v≤si, and xi⊕vxi⊕v is maximized, where ⊕⊕ denotes the bitwise XOR operation, GCD(c,d)GCD(c,d) denotes the greatest common divisor of integers cc and dd, and y∣xy∣x means xx is divisible by yy, or report -1 if no such numbers are found.

Since you are a programmer, Katie needs you to automatically and accurately perform the tasks in the game to satisfy her dear friend Kuro. Let’s help her!

Input

The first line contains one integer qq (2≤q≤1052≤q≤105) — the number of tasks the game wants you to perform.

qq lines follow, each line begins with an integer titi — the type of the task:

  • If ti=1ti=1, an integer uiui follow (1≤ui≤1051≤ui≤105) — you have to add uiui to the array aa.
  • If ti=2ti=2, three integers xixi, kiki, and sisi follow (1≤xi,ki,si≤1051≤xi,ki,si≤105) — you must find a number vv existing in the array aa such that ki∣GCD(xi,v)ki∣GCD(xi,v), xi+v≤sixi+v≤si, and xi⊕vxi⊕v is maximized, where ⊕⊕ denotes the XOR operation, or report -1 if no such numbers are found.

It is guaranteed that the type of the first task is type 11, and there exists at least one task of type 22.

Output

For each task of type 22, output on one line the desired number vv, or -1 if no such numbers are found.

Examples

Input

5
1 1
1 2
2 1 1 3
2 1 1 2
2 1 1 1

Output

2
1
-1

Input

10
1 9
2 9 9 22
2 3 3 18
1 25
2 9 9 20
2 25 25 14
1 20
2 26 26 3
1 14
2 20 20 9

Output

9
9
9
-1
-1
-1

Note

In the first example, there are 5 tasks:

  • The first task requires you to add 11 into aa. aa is now {1}{1}.
  • The second task requires you to add 22 into aa. aa is now {1,2}{1,2}.
  • The third task asks you a question with x=1x=1, k=1k=1 and s=3s=3. Taking both 11and 22 as vv satisfies 1∣GCD(1,v)1∣GCD(1,v) and 1+v≤31+v≤3. Because 2⊕1=3>1⊕1=02⊕1=3>1⊕1=0, 22is the answer to this task.
  • The fourth task asks you a question with x=1x=1, k=1k=1 and s=2s=2. Only v=1v=1satisfies 1∣GCD(1,v)1∣GCD(1,v) and 1+v≤21+v≤2, so 11 is the answer to this task.
  • The fifth task asks you a question with x=1x=1, k=1k=1 and s=1s=1. There are no elements in aa that satisfy the conditions, so we report -1 as the answer to this task.

题意:开始有个空数组,现在有两种操作:

(1,x):给数组加一个数a[]=x;

(2,x,k,s):在k|x的情况下,在数组中找一个a[i],满足a[i]<=s-x,而且k|a[i];现在需要找满足条件的a[],它异或x的值最大。

思路:最大异或,想到Trie树,树上可以贪心地走和x不相同的方向,使得异或最大,而且树上不难满足a[i]<=s-x的条件。

主要问题在于如何满足a[i]是k的倍数。如果操作1,假如a[i],把a[i]的每个的因子都走一遍Trie树,YY一下感觉复杂度很高,所以不敢写,而且不会用指针,空间是肯定要爆炸的。 然而标解就是如此:

时间上:O=1e5*lg1e5*18=2*1e7;ok!

空间上:和时间复杂度差不多大。

(不会写指针,所以比赛的时候只有暴力Trie写了发回溯来满足a[i]是k的倍数。。。第一次写指针,感觉可以入门了!奶思!

#include<bits/stdc++.h>
using namespace std;
const int maxn=;
vector<int>G[maxn];
int vis[maxn];
void read(int &x){
x=; char c=getchar();
while(c>''||c<'') c=getchar();
while(c>=''&&c<='') x=(x<<)+(x<<)+c-'',c=getchar();
}
struct Trie
{
struct node
{
int Min,val;
node *ch[];
node(){
Min=maxn;
ch[]=ch[]=NULL;
}
}*rt[maxn];
void init()
{
for(int i=;i<maxn;i++)
for(int j=i;j<maxn;j+=i)
G[j].push_back(i);
for(int i=;i<maxn;i++) rt[i]=new node;
}
void insert(int x)
{
int Len=G[x].size();
for(int i=;i<Len;i++){
node *cur=rt[G[x][i]];
cur->Min=min(cur->Min,x);
for(int j=;j>=;j--){
if(cur->ch[x>>j&]==NULL) cur->ch[x>>j&]=new node;
cur=cur->ch[x>>j&];
cur->Min=min(cur->Min,x);
} cur->val=x;
}
}
int query(int x,int k,int s)
{
if(x%k!=) return -;
node *cur=rt[k];
if(cur->Min>s-x) return -;
for(int i=;i>=;i--){
int tb=x>>i&;
if(cur->ch[tb^]!=NULL&&cur->ch[tb^]->Min<=s-x) cur=cur->ch[tb^];
else cur=cur->ch[tb];
}
return cur->val;
}
}T;
int main()
{
//cout<<100000*18*log(100000);=2e7
int N,i,j,opt,x,k,s;
T.init();
read(N); while(N--){
read(opt);
if(opt==){
read(x);
if(!vis[x]) vis[x]=,T.insert(x);
}
else{
read(x); read(k); read(s);
printf("%d\n",T.query(x,k,s));
}
}
return ;
}
相关推荐
python开发_常用的python模块及安装方法
adodb:我们领导推荐的数据库连接组件bsddb3:BerkeleyDB的连接组件Cheetah-1.0:我比较喜欢这个版本的cheeta…
日期:2022-11-24 点赞:878 阅读:9,083
Educational Codeforces Round 11 C. Hard Process 二分
C. Hard Process题目连接:http://www.codeforces.com/contest/660/problem/CDes…
日期:2022-11-24 点赞:807 阅读:5,558
下载Ubuntn 17.04 内核源代码
zengkefu@server1:/usr/src$ uname -aLinux server1 4.10.0-19-generic #21…
日期:2022-11-24 点赞:569 阅读:6,407
可用Active Desktop Calendar V7.86 注册码序列号
可用Active Desktop Calendar V7.86 注册码序列号Name: www.greendown.cn Code: &nb…
日期:2022-11-24 点赞:733 阅读:6,180
Android调用系统相机、自定义相机、处理大图片
Android调用系统相机和自定义相机实例本博文主要是介绍了android上使用相机进行拍照并显示的两种方式,并且由于涉及到要把拍到的照片显…
日期:2022-11-24 点赞:512 阅读:7,817
Struts的使用
一、Struts2的获取  Struts的官方网站为:http://struts.apache.org/  下载完Struts2的jar包,…
日期:2022-11-24 点赞:671 阅读:4,900