【GDOI2016模拟4.22】数字方阵

Description

Anica 做了一个很奇怪的梦：她梦见了一个无限大的平板，平板上填着无限行和无限列的整数。有趣的是，每个整数在那神奇的平板上只出现有限的次数。
机智的Anica很快便发现了这其中数字的规律，每一行第一列的数字表示当前的行号，其它非第一列的数字，为该位置左边一列的数字加上其数字的翻转数字之和。为了方便描述，我们定义A[i，j]表示平板上第i行第j列的数字：
1. A[i，1] ＝ i，
2. A[i，j] = A[i，j-1] + Rev( A[i，j-1] ), 其中Rev(x)表示x在10进制下的翻转，例如Rev(213) ＝ 312， Rev(406800) = 008604 ＝ 8604 。
然而，Anica对这堆数字不太感兴趣。她看了一下平板的背面，并发现了一盏神灯，同时神灯也发现了Anica，这时一个精灵冒了出来：
精灵：“Anica，如果你能答对我Q个问题，你将获得一份神秘的礼物，否则你将无法无法醒来，哈哈哈哈～～”
精灵：“每个问题将给你两个整数A和B，请问区间[A，B]的数字总共在平板上出现了多少次？”
Anica：“我不想要你的礼物，我只想要醒过来，我还有很多事情要做呢！”

Input

第一行给出一个整数Q (1 <= Q <=10^5)，表示问题的数量。
接下来Q行，每行两个整数A，B (1 <= A,B <= 10^10)，表示每个询问的区间[A,B]。

Output

输出Q行，第i行一个整数，表示第i个询问的答案。

Sample Input

样例输入1:
2
1 10
5 8
样例输入2:
3
17 144
121 121
89 98
样例输入3:
1
1 1000000000

Sample Output

样例输出1:
18
8
样例输出2:
265
25
10
样例输出3:
1863025563

Data Constraint

50%的分数保证1 <= A,B <= 10^6。

Solution

对于50%：发现数据比较小，暴力。直接暴力第一列的数，每次f[x+rev(x)]+=f[x]。对f求前缀和，每次询问就是f[r]-f[l-1]。比较简单。
对于100%：
同样是暴力，但是要特别一点。
可以发现，总共的数的数量不会超过$19^5$，等会证明。那么，如何暴力这些数呢？
如果
x长为L=a1 a2 a3.. aL
rev(x)=aL aL-1 .. a1
x+rev(x)=a1+aL a2+aL-1 a3+aL-2 aL+a1
得到一个对称的数。L最大为10，那么一半就是5。分别为0~18中一个共19个，所以大约有$19^5$个数。
爆搜这些数，算出它在第二列时有多少种可能从第一列转移过来，存到f里面。接着同50分一样，最后二分一下就行了。
注意：存数用hash，表开大一点。最好所有数组与变量开longlong。

code

#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <iostream>
#define fo(i,a,b) for(ll i=a;i<=b;i++)
#define ll long long
#define mo 6412345
using namespace std;
struct note
{
    ll x,y;
};
note a[mo];
ll b[6],hsh[mo+5],has[mo+5],len,k=0,c[19],tot=0;
ll rev(ll x)
{
    ll jy=0;
    for(;x;x/=10) jy=jy*10+x%10;
    return jy;
}
ll hash(ll x)
{
    ll y;
    y=x%mo+1;
    while (hsh[y]!=0 && hsh[y]!=x) y=(y+1)%mo;
    return y;
}
void dg(int j)
{
    if (j>len)
    {
        ll x=1,y=0;
        fo(i,1,len) 
        {
            ll jy=c[b[i]];  
            if (i==1 && b[i]<10) jy--;
            if (i==len && k%2==1) 
            {
                if(b[i]%2==1) jy=0;
                else jy=1;
            }
            x*=jy;
            if (i!=len) y=y*10+b[i];
        }
        if (x==0) return;
        if (k%2==0) y=y*10+b[len];
        for(int i=len;i>=1;i--) y=y*10+b[i];
        ll jy=hash(y);
        if (hsh[jy]==0) a[++tot].x=y,a[tot].y=x,hsh[jy]=y,has[jy]=tot;
        else a[has[jy]].y+=x;
        return; 
    }
    ll kk;kk=j==1?1:0;
    fo(i,kk,18)
    {
        b[j]=i;dg(j+1);
    }
}
bool cmp(note x,note y){return x.x<y.x;}
ll find(ll l,ll r,ll k)
{
    ll mid;
    for(;l+1<r;)
    {
        mid=(l+r)/2;
        if (a[mid].x<=k) l=mid;else r=mid;
    }   
    if (a[r].x<=k) l=r;return l;
}
int main()
{
    freopen("table.in","r",stdin);freopen("table.out","w",stdout);
    fo(i,0,9) 
    fo(j,0,9) c[i+j]++;
    for(len=1;len<=5;len++)
    {
        k++;dg(1);k++;dg(1);
    }
    sort(a+1,a+1+tot,cmp);
    fo(i,1,tot) 
    {
        ll x=a[i].x;x=x+rev(x);
        if (x<=a[tot].x) a[find(0,tot,x)].y+=a[i].y;
    }
    fo(i,1,tot) a[i].y+=a[i-1].y;
    int ak;scanf("%d",&ak);
    while (ak-->0)
    {
        ll l,r;scanf("%lld%lld",&l,&r);
        ll x=find(0,tot,r);ll y=find(0,tot,l);
        if (a[y].x==l) y--;
        printf("%lld\n",a[x].y-a[y].y+r-l+1);
    }
    fclose(stdin);fclose(stdout);
}