[bzoj 4725--POI2017]Reprezentacje ró?nicowe

给定一个数列a：
当n<=2时,a[n]=n
当n>2，且n是奇数时，a[n]=2a[n-1]
当n>2，且n是偶数时，a[n]=a[n-1]+r[n-1]
其中r[n-1]=mex(|a[i]-a[j]|)(1<=i<=j<=n-1)，mex{S}表示最小的不在S集合里面的非负整数。
数列a的前若干项依次为：1,2,4,8,16,21,42,51,102,112,224,235,470,486,972,990,1980。
可以证明，对于任意正整数x，只存在唯一一对整数(p,q)满足x=a[p]-a[q]，定义为repr(x)。比如repr(17)=(6,3)，repr(18)=(16,15)。
现有n个询问，每次给定一个正整数x，请求出repr(x)。

这道题首先可以发现a数组增长地非常的快，所以可以分成两种情况，一种为p-q可能为1，另一种为p必为q+1。
打个表后就可以发现第56之后的数只有与前一个的数相减才可能小于10^9，所以先判断一下是否为第一种情况，如果不是，那么就只可能是第二种情况。
那第二种情况怎么处理，可以先处理出前56个数中两两之间的所有差，然后因为p必为q+1，r数组是递增的，所以二分一下，再瞎jb乱搞就可以了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define inf 1000000000
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=x*10+ch-'0',ch=getchar();
    return x*f;
}
inline void write(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
int top,sta[2100];
int a[56]={1,2,4,8,16,21,42,51,102,112,224,235,470,486,972,990,1980,2002,4004,4027,8054,8078,16156,16181,32362,32389,64778,64806,129612,129641,259282,259313,518626,518658,1037316,1037349,2074698,2074734,4149468,4149505,8299010,8299049,16598098,16598140,33196280,33196324,66392648,66392693,132785386,132785432,265570864,265570912,531141824,531141876,1062283752,1062283805};
inline int erfen(int limit,int l,int r,bool id)
{
    int ans=top;
    while(l<=r)
    {
        int mid=(l+r)>>1,wy=a[mid];if(id==0)wy=sta[mid];
        if(limit>=wy)ans=mid,l=mid+1;
        else r=mid-1;
    }
    return ans;
}
int main()
{
    //freopen("4725.in","r",stdin);
    //freopen("4725.out","w",stdout);
    for(int i=1;i<56;i++)
    {
        for(int j=0;j<i;j++)
        {
            int wy=a[i]-a[j];
            if(wy>inf)continue;
            sta[++top]=wy;
        }
    }sort(sta+1,sta+top+1);
    int T=read();
    while(T--)
    {
        int n=read();bool bk=false;
        for(int i=1;i<56;i++)
        {
            int ul=erfen(a[i]-n,0,i-1,1);
            if(a[ul]+n==a[i] && ul<56){printf("%d %d\n",i+1,ul+1),bk=true;break;}
        }
        if(bk==false)
        {
            int ul=erfen(n,53,top,0)-53,p=55+(n-53-ul)*2;
            printf("%d %d\n",p+1,p);
        }
    }
    return 0;
}