[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;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章