[SCOI2013]数数

最新推荐文章于 2021-09-05 21:23:03 发布

quan_tum

最新推荐文章于 2021-09-05 21:23:03 发布

阅读量291

点赞数

分类专栏：动态规划——数位dp 文章标签：数位dp

本文链接：https://blog.csdn.net/quan_tum/article/details/82183857

版权

动态规划——数位dp 专栏收录该内容

9 篇文章 0 订阅

订阅专栏

这题不是很懂，要看详细的到这里
这题很明显是数位dp，数组的第二维都是表示之前几位是否达到最大值。
注意前导零…

#include<bits/stdc++.h>
#define ll long long
#define il inline
using namespace std;
const int N=100005,mo=20130427;
#define getchar()(p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
il int read(){
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=(x+(x<<2)<<1)+c-48;
    return x*f;
} 
int n,m,b,L[N],R[N];
ll c[N],S[N],a[N][2],s[N][2],ss[N][2],sl[N][2];
int solve(int *p,int l){
    memset(a,0,sizeof a);
    memset(s,0,sizeof s);
    memset(ss,0,sizeof ss);
    memset(sl,0,sizeof sl);
    a[l][0]=1;
    for(int i=l-1;~i;--i){
        int c=(i==l-1?0:b);
        a[i][0]=a[i+1][0];
        a[i][1]=(c-1+a[i+1][1]*b+a[i+1][0]*p[i])%mo;
        sl[i][0]=sl[i+1][0]+a[i+1][0];
        sl[i][1]=(c-1+sl[i][0]*p[i]+(sl[i+1][1]+a[i+1][1])*b)%mo;
        ss[i][0]=(ss[i+1][0]*b+p[i]*sl[i][0])%mo;
        ss[i][1]=(S[c]+ss[i+1][0]*b*p[i]+S[p[i]]*sl[i][0]+ss[i+1][1]*b%mo*b+S[b]*(sl[i+1][1]+a[i+1][1]))%mo;
        s[i][0]=(s[i+1][0]+ss[i][0])%mo;
        s[i][1]=(s[i+1][0]*p[i]+s[i+1][1]*b+ss[i][1])%mo;
    }
    return(s[0][0]+s[0][1])%mo;
}
int main(){
    b=read();
    c[0]=1;
    for(int i=0;i<N-1;++i) c[i+1]=(c[i]*b+1)%mo;
    S[0]=0;
    for(int i=0;i<b;++i) S[i+1]=(S[i]+i)%mo;
    n=read();
    for(int i=0;i<n;++i) L[n-i-1]=read();
    for(int i=0;i<n;++i){
        if(L[i]>0){--L[i];break;}
        L[i]=b-1;
    }
    if(!L[n-1]) --n;
    m=read();
    for(int i=0;i<m;++i) R[m-i-1]=read();
    return !printf("%d",(solve(R,m)-solve(L,n)+mo)%mo);
}