51nod 1607 卷积和(枚举搜索)

1607 卷积和
基准时间限制：1 秒 空间限制：131072 KB 分值: 40 难度：4级算法题 收藏 关注
杰西最近在研究卷积，例如一个数字1234，那么经过杰西的变换后，就会变成1*4+2*3+3*2+4*1，例如一个数字是234，那么就会变成2*4+3*3+4*2。
形象化的，我们可以设f(x)表示x经过杰西变换后的值。用公式表示 x=dn−1dn−2…d2d1d0 那么 f(x)=∑n−1i=0di∗dn−1−i
例如f(1234)=20，f(234)=25。
现在他想到了一个问题。
如果他想知道对于所有i=L~R时的f(i)的值的和，该怎么做呢。
现在这个问题交给了你。
但是这个答案可能会非常大，因此杰西只想知道最终答案对1000000007取模后的答案。

Input
单组测试数据
一行两个整数L,R（1<=L<=R<=10^18）。
Output
一个数表示答案。
Input示例
3 233
Output示例
8730
System Message (题目提供者)

---

这题看着应该是类似数位dp的那种样子，但是又不是数位dp，因为没办法记录下前面十几位分别是啥，但是我们发现，一旦枚举的那一位小于边界了，你就能直接算出所有的答案，所以我们又可以直接进行搜索。
搜索的过程很麻烦，首先你要预处理出每个长度的所有答案的和，对长度i的时候，需要只要i-2长度的所有和，并且有9*10=90次，对面两端的和，是45 * 45 * 2并且要乘$10^{i-2}$次，因为中间i-2位每个位都有10次选择，开头漏了这里，一直错。

后面计算也是，一旦确定了某一位，对于无关紧要的位，一定要乘以10的次方，不然就漏了。
然后就是类似数位dp一样计算了，先把所有小于这一位的取值都算出来，后面不确定的位不能漏乘倍数，然后考虑这一位相同，搜下一位。就这样

---

代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker,"/STACK:102400000,102400000")

using namespace std;
#define   MAX           1000005
#define   MAXN          1000005
#define   maxnode       205
#define   sigma_size    26
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
const double inf   = 1e18;
const double eps   = 1e-4;
const LL    mod    = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

int a[25];
LL all[25];
LL sum[25];
LL cnt[25];
LL fact[25];
LL dfs(int len,int pos){
    if(len==0){
        LL tmp=0;
        for(int i=pos;i>0;i--){
            tmp+=a[i]*a[pos-i+1];
        }
        return tmp%mod;
    }
    LL tmp=0;
    int cnt=(pos==len)?1:0;
    for(int i=cnt;i<a[len];i++){
        int num=len-1;
        for(int j=pos;j>0;j--){
            int k=pos-j+1;
            if(j>len){
                if(k>len) tmp=(tmp+a[j]*a[k]*fact[num])%mod;
                else if(k==len) tmp=(tmp+a[j]*i*fact[num])%mod;
                else tmp=(tmp+a[j]*45*fact[num-1])%mod;
            }
            else if(j==len){
                if(k>len) tmp=(tmp+i*a[k]*fact[num])%mod;
                else if(k==len) tmp=(tmp+i*i*fact[num])%mod;
                else tmp=(tmp+i*45*fact[num-1])%mod;
            }
            else{
                if(k>len) tmp=(tmp+a[k]*45*fact[num-1])%mod;
                else if(k==len) tmp=(tmp+i*45*fact[num-1])%mod;
                else{
                    if(k==j) tmp=(tmp+285*fact[num-1])%mod;
                    else tmp=(tmp+45*45*fact[num-2])%mod;
                }
            }
        }
    }
    return (tmp+dfs(len-1,pos))%mod;
}

LL solve(LL x){
    if(x==0) return 0;
    int len=0;
    while(x){
        a[++len]=x%10;
        x/=10;
    }
    //cout<<sum[len-1]<<endl;
    return (dfs(len,len)+sum[len-1])%mod;
}

void init(){
    all[0]=0;cnt[0]=0;
    all[1]=285;cnt[1]=285;
    fact[0]=1,fact[1]=10;
    for(int i=2;i<20;i++){
        all[i]=(cnt[i-2]*90%mod+90*45*fact[i-2]%mod)%mod;
        cnt[i]=(cnt[i-2]*100%mod+90*45*fact[i-2]%mod)%mod;
        fact[i]=fact[i-1]*10%mod;
    }
    sum[0]=0;
    for(int i=1;i<20;i++) sum[i]=(sum[i-1]+all[i])%mod;
        //cout<<sum[2]<<endl;
}
LL ans[100005];
int main(){
    //freopen("in.txt","r",stdin);
    init();
    LL L,R;
    cin>>L>>R;
   /* for(int i=1;i<=100000;i++){
        ans[i]=solve(i);
    }
    LL tmp=0;
    for(int i=1;i<=100000;i++){
        int k=i;
        vector<int> v;
        while(k){
            v.push_back(k%10);
            k/=10;
        }
        for(int j=0;j<v.size();j++){
            tmp+=v[j]*v[v.size()-1-j];
        }
        cout<<tmp<<" "<<ans[i]<<" "<<i<<endl;
    }*/
    cout<<((solve(R)-solve(L-1))%mod+mod)%mod<<endl;
    return 0;
}