BZOJ3834: [Poi2014]Solar Panels

一开始看错题了，后来重新看了题…emmmm还是不会啊啊啊啊，写了个貌似复杂度很对的暴力交上去过了？？？

（先讲一下我这个比较搞笑？的做法）
求a~b,c~d的数之间的最大的gcd，尝试枚举一下？
因为一个数x，他的$> \sqrt x$的约数可以通过枚举 <x√ <script id="MathJax-Element-62" type="math/tex">< \sqrt x</script>的约数求得，我们从1到$\sqrt b$枚举$d$，对于$xfrom\ a\ to\ b$，可以覆盖所有 <x√ <script id="MathJax-Element-66" type="math/tex"><\sqrt x</script>的约数，同时我们可以通过$b/d-(a-1)/d$算得一个区间，这个区间内的数都是a到b某个数的约数
于是我们算到了a~b的所有约数是什么
我们对a~b的算得$\sqrt n$个区间，再对c~d的算得$\sqrt n$个区间，再求这两类区间的交集最大的值
跑的好像挺快的qwq

然鹅看了题解发现我好菜
枚举gcd后，直接看$b/i-(a-1)/i$和$d/i-(c-1)/i$是否都为正就行了，因为这里都是带下取整的，所以i只有$O(\sqrt 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;
}