前几天被数位DP搞崩了,开始重新学(......)
① HDU 3555 Bomb
[http://acm.hdu.edu.cn/showproblem.php?pid=3555]
Description:计算1~n之间的数中含有‘49’的数的个数
一种解法:记录之前‘49’是否出现和前一位数的状态,板子题
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn = 34;
int T;
LL n,dp[maxn][maxn][2],num[maxn];
LL dfs(int pos,bool limit,int pre,bool ok){
if(pos==-1) return ok;
if(!limit&&dp[pos][pre][ok]!=-1) return dp[pos][pre][ok];
int up = limit?num[pos]:9;
LL res = 0;
for(int i=0;i<=up;i++) res+=dfs(pos-1,limit&&(i==num[pos]),i,ok||(pre==4&&i==9));
if(!limit) dp[pos][pre][ok] = res;
return res;
}
LL solve(LL x){
int pos = 0;
while(x){
num[pos++] = x%10;
x/=10;
}
return dfs(pos-1,true,0,false);
}
int main(){
scanf("%d",&T);
memset(dp,255,sizeof(dp));
while(T--){
scanf("%lld",&n);
printf("%lld\n",solve(n));
}
return 0;
}
另一种解法:只记录‘49’出现的状态,板子题
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn = 34;
int T;
LL n,dp[maxn][3],num[maxn];
LL dfs(int pos,bool limit,int sta){
if(pos==-1) return sta==2;
if(!limit&&dp[pos][sta]!=-1) return dp[pos][sta];
int up = limit?num[pos]:9;
LL res = 0;
for(int i=0;i<=up;i++){
int tsta = 0;
if(sta==2) tsta = 2;
else if(sta==1&&i==9) tsta = 2;
else if(i==4) tsta = 1;
res+=dfs(pos-1,limit&&(i==num[pos]),tsta);
}
if(!limit) dp[pos][sta] = res;
return res;
}
LL solve(LL x){
int pos = 0;
while(x){
num[pos++] = x%10;
x/=10;
}
return dfs(pos-1,true,0);
}
int main(){
scanf("%d",&T);
memset(dp,255,sizeof(dp));
while(T--){
scanf("%lld",&n);
printf("%lld\n",solve(n));
}
return 0;
}
② HDU 3652 B-number
Description:计算1~n之间的数中,能被13整除并且数中含有‘13’的数的个数
[http://acm.hdu.edu.cn/showproblem.php?pid=3652]
记录两个状态,当前在该位之前的和对13的余数和含有‘13’的状态直接套板子就行了
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn = 32;
int n;
int num[maxn],dp[maxn][13][3];
int dfs(int pos,bool limit,int mod,int sta_B){
if(pos==-1) return sta_B==2&&mod==0;
if(!limit&&dp[pos][mod][sta_B]!=-1) return dp[pos][mod][sta_B];
int up = limit?num[pos]:9;
int res = 0;
for(int i=0;i<=up;i++){
int sta = 0;
if(sta_B==2) sta = 2;
else if(sta_B==1&&i==3) sta = 2;
else if(i==1) sta = 1;
res+=dfs(pos-1,limit&&i==num[pos],(mod+i*(int)(pow(10,pos)))%13,sta);
}
if(!limit) dp[pos][mod][sta_B] = res;
return res;
}
int solve(int x){
int pos = 0;
while(x){
num[pos++] = x%10;
x/=10;
}
return dfs(pos-1,true,0,0);
}
int main(){
memset(dp,255,sizeof(dp));
while(scanf("%d",&n)!=EOF){
printf("%d\n",solve(n));
}
return 0;
}
③ HDU 2089 不要62
Description:计算l~r之间不含4且不含62的数的个数
记录是否含有‘4’和含有‘62’的状态
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 10;
int n,m;
int num[maxn],dp[maxn][2][3];
int dfs(int pos,int limit,bool have_4,int sta_62){
if(pos==-1) return !have_4 && sta_62!=2;
if(!limit&&dp[pos][have_4][sta_62]!=-1) return dp[pos][have_4][sta_62];
int up = limit?num[pos]:9;
int res = 0;
for(int i=0;i<=up;i++){
int sta = 0;
if(sta_62==2||sta_62==1&&i==2) sta = 2;
else if(i==6) sta = 1;
res+=dfs(pos-1,limit&&i==num[pos],have_4||i==4,sta);
}
if(!limit) dp[pos][have_4][sta_62] = res;
return res;
}
int solve(int x){
int pos = 0;
while(x){
num[pos++] = x%10;
x/=10;
}
return dfs(pos-1,true,false,0);
}
int main(){
memset(dp,255,sizeof(dp));
while(scanf("%d %d",&n,&m)!=EOF && n+m) printf("%d\n",solve(m)-solve(n-1));
return 0;
}
④HDU 3709 Balanced Number
【http://acm.hdu.edu.cn/showproblem.php?pid=3709 】
Description:计算L~R之间平很数有多少个,平衡数指的是一个数存在一个中点数位,其左边的各数位值乘上离中点数位的距离的和等于右边各数位值乘上离中点数位的距离的和
记录平衡点的位置和当前的权值和,对每一位都进行计算然后求和,注意要减掉全为0的情况
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn = 20;
int T;
LL l,r;
int digit[maxn];
LL dp[maxn][1700][maxn];
LL dfs(int pos,bool limit,int mid,int sum){
if(pos==-1) return sum==0;
if(sum<0) return 0;
if(!limit && dp[pos][sum][mid]!=-1) return dp[pos][sum][mid];
int up = limit?digit[pos]:9;
LL res = 0;
for(int i=0;i<=up;i++) res+=dfs(pos-1,limit&&i==up,mid,sum+i*(pos-mid));
if(!limit) dp[pos][sum][mid] = res;
return res;
}
LL solve(LL x){
int pos = 0;
while(x){
digit[pos++] = x%10;
x/=10;
}
LL res = 0;
for(int i=0;i<pos;i++) res+=dfs(pos-1,true,i,0);
return res-(pos-1);
}
int main(){
scanf("%d",&T);
memset(dp,255,sizeof(dp));
while(T--){
scanf("%I64d %I64d",&l,&r);
printf("%I64d\n",solve(r)-solve(l-1));
}
return 0;
}
⑤HDU4737 F(x)
【http://acm.hdu.edu.cn/showproblem.php?pid=4734 】
Description:若数x表示为:
A
n
A
n
−
1
A
n
−
2
.
.
.
.
.
.
A
2
A
1
A_n A_{n-1}A_{n-2}......A_2A_1
AnAn−1An−2......A2A1那么
F
(
x
)
=
A
n
⋅
2
n
−
1
+
A
n
−
1
⋅
2
n
−
2
+
.
.
.
.
.
.
+
A
2
⋅
2
+
A
1
⋅
1
F(x)=A_n·2^{n-1}+A_{n-1}·2^{n-2}+......+A_2·2+A_1·1
F(x)=An⋅2n−1+An−1⋅2n−2+......+A2⋅2+A1⋅1
现在给出A和B,要求计算出在0~B之间有多少个数x使得F(x)<=F(A)
由于F(A)已知在变化,所以不能记录当前权值和然后和F(A)进行比较,如果这样做的话dp数组每次都要初始化,会超时.
所以记录的是剩余权值和,这样就和F(A)没有关系了,只要初始化dp数组一遍就好了
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 15;
int T,a,b;
int digit[maxn],dp[maxn][5000];
int dfs(int pos,bool limit,int maxx){
if(maxx<0) return 0;
if(pos==-1) return maxx>=0;
if(!limit && dp[pos][maxx]!=-1) return dp[pos][maxx];
int up = limit?digit[pos]:9;
int res = 0;
for(int i=0;i<=up;i++) res+=dfs(pos-1,limit&&i==digit[pos],maxx-(i<<pos));
if(!limit) dp[pos][maxx] = res;
return res;
}
int solve(int x,int y){
int pos = 0;
while(x){
digit[pos++] = x%10;
x/=10;
}
int maxx = 0;
for(int i=0;y;i++,y/=10) maxx+=(y%10)<<i;
return dfs(pos-1,true,maxx);
}
int main(){
scanf("%d",&T);
memset(dp,255,sizeof(dp));
for(int kase=1;kase<=T;kase++){
scanf("%d %d",&a,&b);
printf("Case #%d: %d\n",kase,solve(b,a));
}
return 0;
}
⑥HDU5179 beautiful number
【http://acm.hdu.edu.cn/showproblem.php?pid=5179 】
Description:
给出一个数字a1 a2 a3 a4 … an (n<=9)问区间[L,R]之间有多少数满足数位中后一位(小的那一位)不大于前一位并且当前位能被之后所有位整除。
DP记录当前位置和前一位的数是多少就好了,然后在dfs的时候判断一下能不能继续往下搜
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 11;
int digit[maxn];
int T,dp[maxn][maxn];
int dfs(int pos,bool limit,int pre){
if(pos==-1) return 1;
if(!limit&&dp[pos][pre]!=-1) return dp[pos][pre];
int up = limit?digit[pos]:9;
int res = 0;
for(int i=0;i<=up;i++) if(pre==0||(i!=0&&i<=pre&&pre%i==0)) res+=dfs(pos-1,limit&&i==digit[pos],i);
if(!limit) dp[pos][pre] = res;
return res;
}
int solve(int x){
int pos = 0;
while(x){
digit[pos++] = x%10;
x/=10;
}
return dfs(pos-1,true,0);
}
int main(){
memset(dp,255,sizeof(dp));
scanf("%d",&T);
while(T--){
int l,r;
scanf("%d %d",&l,&r);
printf("%d\n",solve(r)-solve(l-1));
}
return 0;
}
⑦HDU4507 吉哥系列故事——恨7不成妻
Description:
如果一个数满足和7有关,则:
1、整数中某一位是7;
2、整数的每一位加起来的和是7的整数倍;
3、这个整数是7的整数倍;
要求计算[L,R]之间和7无关的数的平方和
DP要记录三个数据:符合条件的个数,符合条件的数的和,符合条件的数的平方和
根据公式ans为平方和,sum为和,num为个数,当前枚举到数位pos,当前数为i,当前为now,dfs下一位为pre
则满足:
a
n
s
[
n
o
w
]
=
a
n
s
[
p
r
e
]
+
(
2
⋅
i
⋅
1
0
p
o
s
⋅
s
u
m
[
p
r
e
]
)
2
+
n
u
m
⋅
(
i
⋅
1
0
p
o
s
)
2
ans[now]=ans[pre]+(2·i·10^{pos}·sum[pre])^2+num·(i·10^{pos})^2
ans[now]=ans[pre]+(2⋅i⋅10pos⋅sum[pre])2+num⋅(i⋅10pos)2
然后就是鬼畜的取余操作
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long int LL;
const int maxn = 20;
const LL mod = 1e9+7;
LL pow_10[maxn];
int digit[maxn];
struct RES{
LL ans,sum,num;
}dp[maxn][7][7];
RES dfs(int pos,bool limit,int sum_7,int mul_7){
if(pos==-1) return (RES){0,0,sum_7&&mul_7};
if(!limit&&dp[pos][sum_7][mul_7].num!=-1) return dp[pos][sum_7][mul_7];
int up = limit?digit[pos]:9;
RES res = (RES){0,0,0};
for(int i=0;i<=up;i++){
if(i==7) continue;
RES tmp = dfs(pos-1,limit&&i==digit[pos],(sum_7+i)%7,(mul_7*10+i)%7); //这里不能是mul_7+i*pow_10[pos] 会有取余精度损失
if(tmp.num!=0){
res.num=(res.num+tmp.num)%mod;
res.sum=(res.sum+(tmp.sum+i*pow_10[pos]%mod*tmp.num)%mod)%mod;
res.ans=(res.ans+(tmp.ans+(2*i*(pow_10[pos]*tmp.sum%mod))%mod)%mod+(pow_10[pos]*pow_10[pos]%mod*tmp.num%mod*i*i%mod))%mod;
}
}
if(!limit) dp[pos][sum_7][mul_7] = res;
return res;
}
LL solve(LL x){
int pos = 0;
while(x){
digit[pos++] = x%10;
x/=10;
}
return dfs(pos-1,true,0,0).ans;
}
int T;
int main(){
scanf("%d",&T);
for(int i=0;i<maxn;i++) for(int j=0;j<7;j++) for(int k=0;k<7;k++) dp[i][j][k].ans=dp[i][j][k].sum=dp[i][j][k].num=-1;
pow_10[0] = 1;
for(int i=1;i<maxn;i++) pow_10[i] = pow_10[i-1]*10%mod;
while(T--){
LL l,r;
scanf("%I64d %I64d",&l,&r);
printf("%I64d\n",((solve(r)-solve(l-1))%mod+mod)%mod);
}
return 0;
}