Brief Description:
题目要求统计出区间[a, b]中能够被自身的数位之和整除的数的个数。
Analysis:
摘自题解:
分别计算出[1, b]中符合条件的个数和[1, a-1]中符合条件的个数。
d[l][i][j][k]表示前l位和为i模j的结果为k的数的个数,那么就有方程
d[l+1][i+x][j][(k*10+x)%j] += d[l][i][j][k]
预处理出d[l][i][j][k],然后再逐位统计即可。
一开始我也是这么想的,然后用了自己一直用的卡上界的方法,复杂度是10^7级别的,但是题目有50组数据且每组数据都大得惊人,因此用遍各种常数优化都TLE到死。
赛后看别人的做法都是打表,不用DP过的。
至于为什么和题解方法一样,但是会TLE呢?这是因为标程卡上界的方法比较巧,它其实预处理了一部分内容,用了“半打表”,卡上界的方法还是见代码吧,注意打印成模板,以后卡上界会超时就按标程那么写!这样就不用担心超时了。注意程序里有很多细节问题。
标程如下:
#include <iostream>
#include <cstdio>
using namespace std;
int dp[10][82][82][82];
void Pre(){
int i,j,k,l,x;
for(i=1; i<=81; i++) dp[0][0][i][0] = 1;
for(l=0; l<9; l++)
for(i=0; i<=l*9; i++)
for(j=1; j<=81; j++)
for(k=0; k<j; k++)
for(x=0; x<=9; x++){
dp[l+1][i+x][j][(k*10+x)%j] += dp[l][i][j][k];
}
}
int Solve(int x){
if(x == 0) return 0;
int ret = 0;
if(x == 1000000000) { ret++; x--; }
int i,j,k;
int len,sum,t,tt,power=1,bit[11];
len = sum = 0;
t = tt = x;
while(t){
bit[++len] = t % 10;
t /= 10;
sum += bit[len];
}
if(x % sum == 0) ret++; // remember
for(i=1; i<=len; i++){
// basic coding
tt /= 10; power *= 10;
t = tt * power;
sum -= bit[i];
for(j=0; j<bit[i]; j++){ // basic coding
for(k=sum+j; k<=sum+j+9*(i-1); k++){
if(k == 0) continue;
int dd = t % k;
if(dd > 0) dd = k - dd;
ret += dp[i-1][k-sum-j][k][dd];
}
t += power / 10;
}
}
return ret;
}
int main()
{
int ta=1,cas;
Pre();
scanf("%d",&cas);
while(cas--){
int a,b;
scanf("%d%d",&a,&b);
printf("Case %d: %d\n",ta++,Solve(b)-Solve(a-1));
}
return 0;
}