BZOJ3834: [Poi2014]Solar Panels

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=3834

题意：求max（gcd（i，j））  其中 a<=i<=b,c<=j<=d

题解：

先考虑怎么判断是否存在两个给定范围内的数使得gcd（i，j）>=k

这很好办，我们只要判断 b/k>(a-1)/k &&d/k>(c-1)/k

那我们枚举是个好方法，但是超时，我们发现上面的四个值每个都只有根号级别的数量，那我们分块即可。

注意用块内的最后一个数来更新答案。

代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 200000+5

#define maxm 200000+5

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)

#define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)

#define mod 1000000007

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    int T=read();
    while(T--)
    {
        int a=read(),b=read(),c=read(),d=read(),ans;
        if(b>d)swap(a,c),swap(b,d);
        if(b>=c){printf("%d
",b);continue;}
        a--;c--;
        for(int i=1,j;i<=b;i=j+1)
        {
            j=min(b/(b/i),d/(d/i));
            if(i<=a)j=min(j,a/(a/i));
            if(i<=c)j=min(j,c/(c/i));
            if(b/j>a/j&&d/j>c/j)ans=j;
        }
        printf("%d
",ans);
    }

    return 0;

}