题目背景
SOURCE:NOIP2015-SHY4
题目描述
Alice 和 Bob 正在玩一个游戏,两个人从 1 轮流开始报数,如果遇到 7 的倍数或者遇到的这个数的十进制表示中含 7 ,则遇到的那个人需要喊“过”。
例如:
1 2 3 4 5 6 过 8 9 10 11 12 13 过 15 16 过 18 ……
游戏过后,Bob 提出了一个问题:在区间 [L,R] 里有多少数要喊“过”?
输入格式
第一行一个整数 N ,表示共有 N 组数据。
接下来 N 行,每行两个整数 L 和 R ,表示区间 [L,R] 。
输出格式
共 N 行,每行一个整数。分别表示每一组数据的答案。
样例数据 1
输入
3
5 10
7 30
2 100
输出
1
6
30
备注
【数据范围】
对 40% 的输入数据 :N≤30, L,R≤106
对 70% 的输入数据 :N≤300, L,R≤109
对 100% 的输入数据 :N≤3000, L,R≤1018
分析
很明显的数位dp,不用多说
代码
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#define ll long long
using namespace std;
int n;
ll l,r,f[20][520][2][7];
int a[25],t;
ll dp(int pos,bool flag,bool limit,int mod){
if(pos<=0){
if(flag) return 1;
if(mod%7==0) return 1;
return 0;
}
if(f[pos][flag][limit][mod]!=-1) return f[pos][flag][limit][mod];
int up=limit?a[pos]:9;
ll ans=0;
for(int i=0;i<=up;++i){
if(i==7) ans+=dp(pos-1,true,limit&&i==a[pos],(mod*10+i) % 7);
else ans+=dp(pos-1,flag,limit&&i==a[pos],(mod*10+i) % 7);
}
f[pos][flag][limit][mod]=ans;
return ans;
}
ll solve(ll x){
t=0;
while(x){
a[++t]=x%10;
x/=10;
}
memset(f,255,sizeof(f));
return dp(t,0,1,0);
}
int main(){
scanf("%d",&n);
int i,j,k;
for(i=1;i<=n;++i){
scanf("%I64d%I64d",&l,&r);
cout<<solve(r)-solve(l-1)<<'\n';
}
return 0;
}
为什么前面写那么短,因为重点是我的自我分析,好了,我要开始瞎bibi了,你们就别看了
考场上满心欢喜以为自己的数位dp写对了,只是在10^18这种大数据的时候跑不过去,觉得差不多小数据应该没毛病,就往下做了,完全没有想过自己的“正解”可能是错的,没有人能肯定自己写的一定正确,请记住这句话了!那就分段处理啊,至少不能让暴力的分白白流走,这个数据范围40%的分是随便暴力的,所以到最后你连 40分 基础分都搞不到手!
下面来看一看为什么我自以为的正解 WA 了
放一个考场上写的
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#define ll long long
using namespace std;
int n;
ll l,r,f[20][520][2][15];///
int a[25],t;
ll dp(int pos,int status,bool limit,int mod){
if(pos<=0){
if(status==0) return 0;
if((1<<6)&status) return 1;
if(mod%7==0) return 1;
return 0;
}
if(f[pos][status][limit][mod]!=-1) return f[pos][status][limit][mod];
mod%=7;///
int up=limit?a[pos]:9;
ll ans=0;
for(int i=0;i<=up;++i){
if(i==0) ans+=dp(pos-1,status,limit&&i==a[pos],mod*10+i);
else ans+=dp(pos-1,(1<<i-1)|status,limit&&i==a[pos],mod*10+i);
}
f[pos][status][limit][mod]=ans;
return ans;
}
ll solve(ll x){
t=0;
while(x){
a[++t]=x%10;
x/=10;
}
return dp(t,0,1,0);
}
int main(){
scanf("%d",&n);
int i,j,k;
for(i=1;i<=n;++i){
memset(f,255,sizeof(f));///
scanf("%I64d%I64d",&l,&r);
printf("%I64d\n",solve(r)-solve(l-1));
}
return 0;
}
错误点:
1. f 数组要在 solve 里面初始化,因为每次重新 solve 的时候 f 值都是不一样的(memset 可以直接赋0,-1,除此之外不要赋其他)
2. mod 要在前面取,不然 f 数组会溢出(之前做的那道题是因为最后一次的 mod 对检验有贡献,但这里又不是,灵活处理!)
3.检验这个数是否含有 7 根本不能像我这样搞,因为会TLE!,每次memset的时候你以为不耗时间啊,你开那么大的数组,哪怕没用到,但赋初值的时候也会考虑进去
4.还有那个前导0你根本不用管啊,两个相减后不就没有了吗
唉唉唉唉唉,错错错,莫莫莫,数位dp没学到家啊