[SCOI2013]数数

版权声明:欢迎转载(标记出处),写得不好还请多指教 https://blog.csdn.net/quan_tum/article/details/82183857

这题不是很懂,要看详细的到这里
这题很明显是数位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);
}

猜你喜欢

转载自blog.csdn.net/quan_tum/article/details/82183857