BZOJ4725: [POI2017]Reprezentacje ró?nicowe

$n leq 1e5$，$x leq 1e9$。

1e9呵呵，暴力处理$a_n$的前几项直到1e9。然后处理出差的数列，每次在这里面找，找得到就回答，找不到，那有贡献的只有$a_i-a_{i-1},a_i>1e9$，并且是越来越大的，那在原有的里面二分一下看要补几个即可。

注意$Max=1e9+多一点$，否则会少算$a_{56}-a_{54}$等！

#include<stdio.h>
#include<string.h>
//#include<iostream>
#include<algorithm>
using namespace std;

#define LL long long
int qread()
{
    char c; int s=0,t=1; while ((c=getchar())<'0' || c>'9') (c=='-') && (t=-1);
    do s=s*10+c-'0'; while ((c=getchar())>='0' && c<='9'); return s*t;
}

//

int tot;
#define maxn 233333
int a[maxn];
struct Node
{
    int x,y,v;
    bool operator < (const Node &b) const {return v<b.v;}
}b[maxn]; int lb;

int main()
{
    int Max=1100000000,n=2;
    a[1]=1; a[2]=2;
    b[++lb]=(Node){2,1,1};
    for (n=3;;n++)
    {
        if (n&1) a[n]=a[n-1]<<1; else
        {
            for (int k=2;k<=lb;k++) if (b[k].v!=b[k-1].v+1) {a[n]=b[k-1].v+1; break;}
            if (!a[n]) a[n]=b[lb].v+1;
            a[n]+=a[n-1];
        }
        for (int i=1;i<n;i++) b[++lb]=(Node){n,i,a[n]-a[i]};
        sort(b+1,b+1+lb);
        if (a[n]>Max) break;
    }
    int m=qread();
    while (m--)
    {
        int x=qread(),y;
        if (b[y=lower_bound(b+1,b+1+lb,(Node){0,0,x})-b].v==x) printf("%d %d
",b[y].x,b[y].y);
        else
        {
            int p=x-(y-1);
            printf("%d %d
",n+p*2-1,n+p*2-2);
        }
    }
    return 0;
}