ACM-思维-D.Mr. Panda and Geometric Sequence

传送门 http://codeforces.com/gym/101775/problem/D
训练赛的时候被我误认为是数位dp了，结果没来得及模改。。。。
看了网上的题解才意识到可以进行预处理，根本不需要数位dp，55555
大致就是，预处理出1e15中所有的相关数值，然后排序去重，二分查找位置，r减去l-1就是对应的数量

#include<bits/stdc++.h>
using namespace std;
#define ll unsigned long long
#define mem(a,b) memset(a,b,sizeof a)
const int mod=1e9+7;
const int maxn=1e5;//3个确定 1e15
int p,k;
vector<ll>num;
ll qpow[20];
int sum(ll n){
    int cnt=0;
    while(n>0){
        n/=10;
        cnt++;
    }
    return cnt;
}
void init(){//预处理出所有15位以内符合条件的数值
    num.clear();
    qpow[0]=1;
    for(int i=1;i<=16;i++)qpow[i]=qpow[i-1]*10;
    //printf("debug\n");
    for(int i=1;i<=maxn;i++){//    r=j/i
        for(int j=i+1;j<=maxn/i;j++){
            if(__gcd(i,j)>1)continue;
            for(int k=1;k<=(maxn/i)/j;k++){//构造前三个数值
                ll x=1LL*k*i*i,y=1LL*k*i*j,z=1LL*k*j*j;
                ll sumx=sum(x),sumy=sum(y),sumz=sum(z);
                ll sumall=sumz+sumx+sumy;
                if(sumall>15)break;
                ll now=z;
                ll res=x*qpow[sumy+sumz]+y*qpow[sumz]+z;
                num.push_back(res);
                while(1){
                    if(now%i!=0)break;
                    now=now*j/i;//*j/i
                    int ge=sum(now);
                    if(sumall+ge>15)break;
                    sumall=sumall+ge;
                    res=res*qpow[ge]+now;
                    num.push_back(res);
                }
            }
        }
    }
    //排序并去重(新操作，学一波)
    sort(num.begin(),num.end());
    num.resize(unique(num.begin(),num.end())-num.begin());

}
ll erfen(ll n){//二分查找<=n的最后一个值,最后的答案就是r的位置-(l-1)的位置
    return upper_bound(num.begin(),num.end(),n)-num.begin();
}
int main(){
    int t;
    init();
    cin>>t;
    int tp=1;
    while(t--){
       ll l,r;
       cin>>l>>r;
       cout<<"Case #"<<tp++<<": "<<erfen(r)-erfen(l-1)<<"\n";
    }
    return 0;
}