BZOJ2757 : [SCOI2012]Blinker的仰慕者

BZOJ AC900题纪念~~

 

若K>0，则

设f[i][j]表示i位数字，积为j的数字的个数

g[i][j]表示i位数字，积为j的数字的和

DP+Hash预处理

查询时枚举LCP然后统计贡献

若K=0，则

设f[i][j][k][l]表示已知前i位，乘积是否不为0，是否等于x，是否有数字的数字的个数

g[i][j][k][l]表示已知前i位，乘积是否不为0，是否等于x，是否有数字的数字的和

数位DP即可

 

#include<cstdio>

typedef long long ll;

const int P=20120427,H=1000037,M=66062;

int a[19],len,i,j,k,t,T,g[H],nxt[M],ed;ll n,l,r,K,tmp,pow[19],fin,v[M];

inline int get(ll u){

  int x=u%H,i=g[x];

  for(;i;i=nxt[i])if(v[i]==u)return i;

  return v[++ed]=u,nxt[ed]=g[x],g[x]=ed;

}

inline int ask(ll u){

  for(int i=g[u%H];i;i=nxt[i])if(v[i]==u)return i;

  return 0;

}

namespace Subtask1{

int f[19][M],g[19][M];

void init(){

  f[0][get(1)]=1;

  for(i=0;i<18;i++)for(j=1;j<=ed;j++)for(k=1;k<=9;k++){

    if(v[j]*k>1000000000000000000LL)break;

    (f[i+1][t=get(v[j]*k)]+=f[i][j])%=P;

    (g[i+1][t]+=g[i][j]*10+f[i][j]*k)%=P;

  }

}

inline ll cal(ll x,ll K){

  ll ans=0,pre=0,mul=1;

  for(len=0,tmp=x;tmp;a[++len]=tmp%10,tmp/=10);

  for(i=len;i;i--){

    if(!a[i])break;

    for(j=1;j<a[i];j++)if((K%(mul*j))==0){

      tmp=K/mul/j;

      t=ask(tmp);

      (ans+=(pre*10+j)*pow[i-1]%P*f[i-1][t]%P+g[i-1][t])%=P;

    }

    (pre*=10)+=a[i],mul*=a[i];

  }

  for(K=ask(K),i=1;i<len;i++)(ans+=g[i][K])%=P;

  return ans;

}

}

namespace Subtask2{

ll f[19][2][2][2],g[19][2][2][2];

inline ll cal(ll x){

  for(len=0,tmp=x;tmp;a[++len]=tmp%10,tmp/=10);

  for(i=1,j=len;i<=len&&i<j;i++,j--)t=a[i],a[i]=a[j],a[j]=t;

  for(i=1;i<=len;i++)for(j=0;j<2;j++)for(k=0;k<2;k++)for(t=0;t<2;t++)f[i][j][k][t]=g[i][j][k][t]=0;

  for(i=0;i<=a[1];i++)f[1][i>0][i==a[1]][i>0]++,g[1][i>0][i==a[1]][i>0]+=i;

  for(i=1;i<len;i++)for(j=0;j<2;j++){

    for(k=0;k<=9;k++)(f[i+1][k>0][0][k>0]+=f[i][j][0][0])%=P,(g[i+1][k>0][0][k>0]+=g[i][j][0][0]*10%P+f[i][j][0][0]*k%P)%=P;

    for(k=0;k<=9;k++)(f[i+1][j&&k>0][0][1]+=f[i][j][0][1])%=P,(g[i+1][j&&k>0][0][1]+=g[i][j][0][1]*10%P+f[i][j][0][1]*k%P)%=P;

    for(k=0;k<=a[i+1];k++)(f[i+1][j&&k>0][k==a[i+1]][1]+=f[i][j][1][1])%=P,(g[i+1][j&&k>0][k==a[i+1]][1]+=g[i][j][1][1]*10%P+f[i][j][1][1]*k%P)%=P;

  }

  return g[len][0][0][1];

}

}

int main(){

  for(pow[0]=i=1;i<19;i++)pow[i]=pow[i-1]*10;

  Subtask1::init();

  scanf("%d",&T);

  while(T--){

    scanf("%lld%lld%lld",&l,&r,&K);

    if(K)fin=(Subtask1::cal(r+1,K)-Subtask1::cal(l,K))%P;else fin=(Subtask2::cal(r+1)-Subtask2::cal(l))%P;

    while(fin<0)fin+=P;

    printf("%lld\n",fin);

  }

  return 0;

}