题目描述
Fish 是一条生活在海里的鱼，有一天他很无聊，就开始数数玩。他数数玩的具体规则是：
确定数数的进制B
确定一个数数的区间[L, R]
对于[L, R] 间的每一个数，把该数视为一个字符串，列出该字符串的每一个（连续的）子串对应的B进制数的值。
对所有列出的数求和。现在Fish 数了一遍数，但是不确定自己的结果是否正确了。由于[L, R] 较大，他没有多余精力去验证是否正确，你能写一个程序来帮他验证吗？
输入输出格式
输入格式：
输入包含三行。
第一行仅有一个数B，表示数数的进制。
第二行有N +1 个数，第一个数为N，表示数L 在B 进制下的长度为N，接下里的N个数从高位到低位的表示数L 的具体每一位。
第三行有M+ 1 个数，第一个数为M，表示数R 在B 进制下的长度为M，接下里的M个数从高位到低位的表示数R 的具体每一位。
输出格式：
输出仅一行，即按照Fish 数数规则的结果，结果用10 进制表示，由于该数可能很大，输出该数模上20130427的模数。
输入输出样例
输入样例#1：
10
3 1 0 3
3 1 0 3
输出样例#1：
120
说明
【样例解释】
[103, 103] 之间仅有数103，该数的所有子串包括1, 10, 103, 0, 03, 3，其和为120。
【数据范围与约定】
20% 数据，0 <= R <= L <= 10^5。
50% 数据，2 <= B <= 1000，1 <= N,M <= 1000。
100% 数据，2 <= B <= 10^5，1 <= N,M <= 10^5。
题解
我不会
题解的思路明白，但代码根本看不懂
米娜自己去看好了……我实在无能为力……
//minamoto
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#define ll long long
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[<<],*p1=buf,*p2=buf;
inline int read(){
#define num ch-'0'
char ch;bool flag=;int res;
while(!isdigit(ch=getc()))
(ch=='-')&&(flag=true);
for(res=num;isdigit(ch=getc());res=res*+num);
(flag)&&(res=-res);
#undef num
return res;
}
const int N=,mod=;
int n,m,B,L[N],R[N];
ll SB[N],S[N],a[N][],s[N][],ss[N][],sl[N][];
int solve(int *p,int l){
memset(a,,sizeof(a));
memset(s,,sizeof(s));
memset(ss,,sizeof(ss));
memset(sl,,sizeof(sl));
a[l][]=;
for(int i=l-;~i;--i){
int c=(i==l-?:B);
a[i][]=a[i+][];
a[i][]=(c-+a[i+][]*B+a[i+][]*p[i])%mod;
sl[i][]=sl[i+][]+a[i+][];
sl[i][]=(c-+sl[i][]*p[i]+
(sl[i+][]+a[i+][])*B)%mod;
ss[i][]=(ss[i+][]*B+p[i]*sl[i][])%mod;
ss[i][]=(S[c]+ss[i+][]*B*p[i]+S[p[i]]*sl[i][]+
ss[i+][]*B%mod*B+
S[B]*(sl[i+][]+a[i+][]))%mod;
s[i][]=(s[i+][]+ss[i][])%mod;
s[i][]=(s[i+][]*p[i]+s[i+][]*B+ss[i][])%mod;
}
return (s[][]+s[][])%mod;
}
int main(){
//freopen("testdata.in","r",stdin);
B=read();SB[]=;
for(int i=;i<N-;++i) SB[i+]=(SB[i]*B+)%mod;
S[]=;
for(int i=;i<B;++i) S[i+]=(S[i]+i)%mod;
n=read();
for(int i=;i<n;++i) L[n-i-]=read();
for(int i=;i<n;++i){
if(L[i]>){--L[i];break;}
L[i]=B-;
}
if(!L[n-]) --n;
m=read();
for(int i=;i<m;++i) R[m-i-]=read();
printf("%d\n",(solve(R,m)-solve(L,n)+mod)%mod);
return ;
}