一开始看错题了,后来重新看了题…emmmm还是不会啊啊啊啊,写了个貌似复杂度很对的暴力交上去过了???
(先讲一下我这个比较搞笑?的做法)
求a~b,c~d的数之间的最大的gcd,尝试枚举一下?
因为一个数x,他的
>x√
的约数可以通过枚举
<x√
<script id="MathJax-Element-62" type="math/tex">< \sqrt x</script>的约数求得,我们从1到
b√
枚举
d
,对于
于是我们算到了a~b的所有约数是什么
我们对a~b的算得
n√
个区间,再对c~d的算得
n√
个区间,再求这两类区间的交集最大的值
跑的好像挺快的qwq
然鹅看了题解发现我好菜
枚举gcd后,直接看
b/i−(a−1)/i
和
d/i−(c−1)/i
是否都为正就行了,因为这里都是带下取整的,所以i只有
O(n√)
种取值,直接分块就行了qwq
code:
#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
inline void up(int &x,const int &y){if(x<y)x=y;}
inline void down(int &x,const int &y){if(x>y)x=y;}
const int maxn = 110000;
int n;
struct node{int l,r;}A[maxn],B[maxn],t1[maxn],t2[maxn]; int an,bn;
int div(node s[],int &sn,const int a,const int b)
{
int n1=0,n2=0;
int u=sqrt(b); if((u+1)*(u+1)==b) u++;
for(int i=1;i<=u;i++) if(b/i-(a-1)/i)
{
t1[++n1]=(node){i,i};
int l=a/i,r=b/i;
if(i*l<a) l++;
t2[++n2]=(node){l,r};
}
sn=0; int now1=1,now2=n2;
while(now1<=n1&&now2)
{
if(t1[now1].r<=t2[now2].r) s[++sn]=t1[now1++];
else s[++sn]=t2[now2--];
while(sn>1&&s[sn].l<=s[sn-1].r) down(s[sn-1].l,s[sn].l),up(s[sn-1].r,s[sn].r),sn--;
}
while(now1<=n1)
{
s[++sn]=t1[now1++];
while(sn>1&&s[sn].l<=s[sn-1].r+1) down(s[sn-1].l,s[sn].l),up(s[sn-1].r,s[sn].r),sn--;
}
while(now2)
{
s[++sn]=t2[now2--];
while(sn>1&&s[sn].l<=s[sn-1].r+1) down(s[sn-1].l,s[sn].l),up(s[sn-1].r,s[sn].r),sn--;
}
}
int main()
{
scanf("%d",&n);
while(n--)
{
int a,b,c,d; scanf("%d%d%d%d",&a,&b,&c,&d);
div(A,an,a,b);
div(B,bn,c,d);
int ans=1;
for(int i=an;i>=1&&bn;i--)
{
while(bn&&B[bn].l>A[i].r) bn--;
if(B[bn].r>=A[i].l) { ans=min(B[bn].r,A[i].r); break; }
}
printf("%d\n",ans);
}
return 0;
}