dp学习笔记8 数位dp好难QAQ

题目
先贴个大佬的模板

typedef long long ll;
int a[20];
ll dp[20][state];//不同题目状态不同
ll dfs(int pos,/*state变量*/,bool lead/*前导零*/,bool limit/*数位上界变量*/)//不是每个题都要判断前导零
{
    //递归边界，既然是按位枚举，最低位是0，那么pos==-1说明这个数我枚举完了
    if(pos==-1) return 1;/*这里一般返回1，表示你枚举的这个数是合法的，那么这里就需要你在枚举时必须每一位都要满足题目条件，也就是说当前枚举到pos位，一定要保证前面已经枚举的数位是合法的。不过具体题目不同或者写法不同的话不一定要返回1 */
    //第二个就是记忆化(在此前可能不同题目还能有一些剪枝)
    if(!limit && !lead && dp[pos][state]!=-1) return dp[pos][state];
    /*常规写法都是在没有限制的条件记忆化，这里与下面记录状态是对应，具体为什么是有条件的记忆化后面会讲*/
    int up=limit?a[pos]:9;//根据limit判断枚举的上界up;这个的例子前面用213讲过了
    ll ans=0;
    //开始计数
    for(int i=0;i<=up;i++)//枚举，然后把不同情况的个数加到ans就可以了
    {
        if() ...
        else if()...
        ans+=dfs(pos-1,/*状态转移*/,lead && i==0,limit && i==a[pos]) //最后两个变量传参都是这样写的
        /*这里还算比较灵活，不过做几个题就觉得这里也是套路了
        大概就是说，我当前数位枚举的数是i，然后根据题目的约束条件分类讨论
        去计算不同情况下的个数，还有要根据state变量来保证i的合法性，比如题目
        要求数位上不能有62连续出现,那么就是state就是要保存前一位pre,然后分类，
        前一位如果是6那么这意味就不能是2，这里一定要保存枚举的这个数是合法*/
    }
    //计算完，记录状态
    if(!limit && !lead) dp[pos][state]=ans;
    /*这里对应上面的记忆化，在一定条件下时记录，保证一致性，当然如果约束条件不需要考虑lead，这里就是lead就完全不用考虑了*/
    return ans;
}
ll solve(ll x)
{
    int pos=0;
    while(x)//把数位都分解出来
    {
        a[pos++]=x%10;//个人老是喜欢编号为[0,pos),看不惯的就按自己习惯来，反正注意数位边界就行
        x/=10;
    }
    return dfs(pos-1/*从最高位开始枚举*/,/*一系列状态 */,true,true);//刚开始最高位都是有限制并且有前导零的，显然比最高位还要高的一位视为0嘛
}
int main()
{
    ll le,ri;
    while(~scanf("%lld%lld",&le,&ri))
    {
        //初始化dp数组为-1,这里还有更加优美的优化,后面讲
        printf("%lld\n",solve(ri)-solve(le-1));
    }
}

题意: f(x) = Sigma (ai * 2^(n-i))，n为数字a的位数，ai是数字a的第i位。给定l,r.求[0,r]中有多少个数满足f(i) <= f(l)
思路: 先求出来f(l)，会发现这个f(x)范围不大，可以开在数组的第二维。
那么我们这里要注意一个问题，这里的状态不能定义为枚举到第i位，和为sum时符合条件的数的数量。因为这样会因l而异，只有sum<=l的才行，l变了就寄了。多组测试样例又不能每次都memset。
我们可以将状态改为枚举到第i位，剩下的数凑出sum时符合条件的数的数量。这个就与l无关了，将初始条件设为f(l),当sum < 0时return，防止越界。
时间复杂度: O(能过)
代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int cas = 1;
int l,r;
int wh;
int a[20];
int f[20][20500]; //枚举了前i位当前和为j的方案数
int check(int x)
{
	string s = to_string(x);
    n = s.size();
    int res = 0;
    for(int i=0;i<s.size();++i) res += (1<<(n-i-1))*(s[i]-'0');
    return res;
}
int dfs(int cur,int sum,bool limit,bool lead)
{
	if(sum < 0) return 0;
	if(cur==-1) return 1;
	if(!lead && !limit && ~f[cur][sum]) return f[cur][sum];
	int up = limit?a[cur]:9;
	int ans = 0;
	for(int i=0;i<=up;++i)
	{
		ans += dfs(cur-1,sum-i*(1<<(cur)),limit&&i==up,lead&&i==0);
	}
	if(!lead&&!limit) f[cur][sum] = ans;
	return ans;
}
int fun(int x)
{
    int pos;
    for(pos=0;x;x/=10) a[pos++] = x%10;
    return dfs(pos-1,wh,1,1);	
}
void solve()
{
   cin>>l>>r;
   wh = check(l);
   printf("Case #%d: %d\n",cas++,fun(r));
}
signed main(void)
{  
   mem(f,-1);
   // T = 1;
   // OldTomato; cin>>T;
   read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}