Palindrome Function

题目飞机直达
题目:Palindrome Function
题意:求sigma(f(i,j)),L<=i<=R, l<=j<+r, 2<=l<=36, 1<=L<=R<=1e9; 数据组数可达1e5;
f(i, j)=: j, j进制下i是回文数字,
1, j进制下i不是回文数字.
思路:由于数据组数太大了, 所以考虑得出结果应该不费劲.这说明要么是结论题,要么需要大量的预处理,或者是要用到记忆化什么的.
任意进制下1e9以内的回文数会有多大呢?不可能超过1e5个吧, 因为回文数的半边一定是不相同的数, 那就好办了。可以把所有进制
下的回文数全部构造出来（就枚举一边就行了，想想就知道不会有重复）， 然后排个序， 对每个询问每个进制进行两次二分查找
就可以轻松得到答案。

#include <cstdio>#include <vector>#include <algorithm>using namespace std;typedef long long LL;int const maxa=1e9+0.2;int num[100];vector<int> G[40];void init(){    for (int i=1; i<=200000; i++){        //一开始i到100000 Wa了, 改大一些莫名就A了...至今想不明白到100000怎么会漏掉的呢, 如果有同学知道麻烦告知一声        for (int b=2; b<=36; b++){            int p=0, t=i;            while (t>0){                num[p++]=t%b;                t/=b;            }            LL a1=i, a2=i*b+num[0];//a2偶数位,  a1奇数位(最后一个数字只取一次)            for (int j=1; j<p; j++){                a1=a1*b+num[j];                a2=a2*b+num[j];            }            //a1, a2在b进制下一定是回文数            if (a1<=maxa) G[b].push_back(a1);            if (a2<=maxa) G[b].push_back(a2);        }    }    for (int b=2; b<=36; b++) {        sort(G[b].begin(), G[b].end());    }}int main(void){    init();    int T, ca=1, L, R, l, r;    for (scanf("%d", &T); T; T--){        scanf("%d%d%d%d", &L, &R, &l, &r);        LL ans=0;        for (int b=l; b<=r; b++){            int len=upper_bound(G[b].begin(), G[b].end(), R)-lower_bound(G[b].begin(), G[b].end(), L);            //得到b进制下,[L, R]的回文数            ans+=b*len+(R+1-L-len);        }        printf("Case #%d: %I64d\n", ca++, ans);    }    return 0;}