数位dp记录

最新推荐文章于 2024-05-20 20:08:58 发布
最新推荐文章于 2024-05-20 20:08:58 发布
阅读量234 收藏 5
点赞数 3
文章标签： 算法 c++ 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/kyrietheshy/article/details/137234944
版权

学习了oiwiki的习题和数位dp，记录打卡一下。

都是比较经典的数位dp，感觉收益很大。

Ahoi2009 self 同类分布

经典的数位dp。

dp[20][200][200]，开三维dp，第一维表示位数的遍历，第二维表示当前所有数字的位数之和，第三位表示当前数字的模数。通过暴力去枚举所有有可能的和（9*len）。然后一起进行数位dp，从而求解出最终得数。

// #pragma GCC optimize(3)  //O2优化开启
#include<bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
typedef pair<int,int> PII;
const int mod=1e9+7;
const int MX=0x3f3f3f3f3f3f3f3f; 
// static char buf[100000],*pa=buf,*pd=buf;
// #define gc pa==pd&&(pd=(pa=buf)+fread(buf,1,100000,stdin),pa==pd)?EOF:*pa++
// inline int read()
// {
//     register int x(0);register char c(gc);
//     while(c<'0'||c>'9')c=gc;
//     while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=gc;
//     return x;
// }
int l,r;
int dp[20][200][200];
int a[50];
int query(int x){
    if(x==0)return 0;
    int len=0;
    while(x)a[++len]=x%10,x/=10;
    auto dfs=[&](auto self,int id,int f,int sum,int p,int sux)->int{
        if(id==0){
            return sum==p&&(sux==0);
        }
        if(f==0&&dp[id][sum][sux]!=-1)return dp[id][sum][sux];
        if(sum>p)return 0;
        int res=9;
        if(f)res=a[id];
        int cnt=0;
        for(int i=0;i<=res;i++){
            cnt+=self(self,id-1,f&(i==res),sum+i,p,(sux*10ll+i)%p);
        }
        return dp[id][sum][sux]=cnt;
    };    
    int ans=0;
    for(int i=1;i<=9*len;i++){
        memset(dp,-1,sizeof dp);
        ans+=dfs(dfs,len,1,0,i,0);
    }
    // cout<<ans<<"+++\n";
    return ans;

}
void icealsoheat(){
    cin>>l>>r;
    l--;
    cout<<query(r)-query(l);

}
signed main(){
    ios::sync_with_stdio(false);          //int128不能用快读！！！！！！
    cin.tie();
    cout.tie();
    int _yq;
    _yq=1;
    // cin>>_yq;
    while(_yq--){
        icealsoheat();
    }
}

洛谷 P3413 SAC#1 - 萌数

dp[1005][10][10][2]

我们需要记录当前位数的位置，的上一位以及上上位置的数据值，以及当前是否已经满足了对称的条件从而分别作为数位dp的四维，来进行数位dp。

// #pragma GCC optimize(3)  //O2优化开启
#include<bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
typedef pair<int,int> PII;
const int mod=1e9+7;
const int MX=0x3f3f3f3f3f3f3f3f; 
// static char buf[100000],*pa=buf,*pd=buf;
// #define gc pa==pd&&(pd=(pa=buf)+fread(buf,1,100000,stdin),pa==pd)?EOF:*pa++
// inline int read()
// {
//     register int x(0);register char c(gc);
//     while(c<'0'||c>'9')c=gc;
//     while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=gc;
//     return x;
// }
string l,r;
int dp[1005][11][11][3];
int a[10005];
int len;
int ui;
int query(string &x){
    // int len=0;
    len=0;
    for(int i=0;i<x.size();i++){
        a[++len]=x[i]-'0';
    }
    auto dfs=[&](auto self,int id,int f,int f0,int le,int lle,int p)->int{
        if(id>len){
            return p;
        }
        if(le!=-1&&!f&&dp[id][le][lle][p]!=-1)return dp[id][le][lle][p];
        int res=9;
        if(f)res=a[id];
        int cnt=0;
        for(int i=0;i<=res;i++){
            int yu=-1;
            if(f0|i)yu=i;
            cnt=(cnt+self(self,id+1,f&&(i==res),f0||(i!=0),yu,le,p||((i==le)&&f0)||((i==lle)&&f0)))%mod;
        }
        
        if(lle!=-1&&f0){
            dp[id][le][lle][p]=cnt;
        }
        return cnt;
    };
    

    memset(dp,-1,sizeof dp);

    int ans=dfs(dfs,1,1,0,-1,-1,0);

    return ans;

}
void icealsoheat(){
    cin>>l>>r;
    int an=0;
    for(int i=0;i<l.size();i++){
        if(i>0&&l[i]==l[i-1]){
            an=1;
        }
        if(i>1&&l[i]==l[i-2])an=1;
    }
    cout<<((query(r)-query(l)+an)%mod+mod)%mod;


}
signed main(){
    ios::sync_with_stdio(false);          //int128不能用快读！！！！！！
    cin.tie();
    cout.tie();
    int _yq;
    _yq=1;
    // cin>>_yq;
    while(_yq--){
        icealsoheat();
    }
}

CF55D Beautiful numbers

我觉得很有代表性的一道题。

最开始想过用vector和map去维护1到9的数位dp，但最后以失败告终。看了题解才发现，我们只需要处理出1到9的最小公倍数lcm就行了。

三维dp代表了位数，当前数取模2520的值，以及lcm。

此时位dp[18][2520][2020]，依旧会爆内存，但是我们可以发现，lcm一共有的数字不超过50个，我们需要提前预处理一下，可以把第三维优化成50。

与此同时，我们也不需要每一次都清空。我之前的数位dp都是正序遍历，并且每一次把dp值手动清空。但是这道题告诉我数位dp还是逆序遍历的好，逆序遍历在这道题中是可以不用每次都清空dp的值的，能大大减少时间复杂度（不然会t）。

#pragma GCC optimize(3)  //O2优化开启
#include<bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
typedef pair<int,vector<int>> PII;
const int mod=2520;
// const int MX=0x3f3f3f3f3f3f3f3f; 
// static char buf[100000],*pa=buf,*pd=buf;
// #define gc pa==pd&&(pd=(pa=buf)+fread(buf,1,100000,stdin),pa==pd)?EOF:*pa++
// inline int read()
// {
//     register int x(0);register char c(gc);
//     while(c<'0'||c>'9')c=gc;
//     while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=gc;
//     return x;
// }
 
string l,r;
int dp[20][2525][50];
int hh[2525];
int a[20];
void yu(){
    int idx=0;
    for(int i=1;i<=2520;i++){
        if(2520%i==0){
            hh[i]=++idx;
        }
    }
}
 
int LCM(int a,int b){
    return a*b/__gcd(a,b);
}
 
int query(string &x){
    int len=0;
    for(int i=x.size()-1;i>=0;i--){
        a[++len]=x[i]-'0';
    }
    // memset(dp,-1,sizeof dp);
    auto dfs=[&](auto self,int id,int f,int sum,int lcm)->int{
        if(!id)return sum%lcm==0;
        if(dp[id][sum][hh[lcm]]!=-1&&f==0)return dp[id][sum][hh[lcm]];
        int res=9;
        if(f)res=a[id];
        int ans=0;
        for(int i=0;i<=res;i++){
            int u=i?LCM(lcm,i):lcm;
            ans+=self(self,id-1,f&&(i==res),(sum*10+i)%mod,u);
        }
 
        if(f==0)dp[id][sum][hh[lcm]]=ans;
        return ans;
 
    };
    int ans= dfs(dfs,len,1,0,1);
    // cout<<ans<<"+++\n";
    return ans;
}
void icealsoheat(){
    cin>>l>>r;
    int an=0;
    vector<int>p(10,0);
    vector<int>ff(10,0);
    for(int i=0;i<l.size();i++){
        int x=l[i]-'0';
        ff[x]=1;
        for(int j=1;j<=9;j++){
 
                p[j]=(p[j]*10+x)%j;
 
        }
    }
    an=1;
    for(int i=1;i<=9;i++){
        if(p[i]&&ff[i]){
            an=0;
            break;
        }
    }
 
    // cout<<(query(r)-query(l)+an)<<"\n";
    printf("%lld\n",(query(r)-query(l)+an));
    // cin>>l;
    // query(l);
 
 
}
signed main(){
    ios::sync_with_stdio(false);          //int128不能用快读！！！！！！
    cin.tie();
    cout.tie();
    int _yq;
    _yq=1;
    cin>>_yq;
    memset(dp,-1,sizeof dp);
    yu();
    while(_yq--){
        icealsoheat();
    }
}

CF628D Magic Numbers

和之前的题很像，直接上代码

#pragma GCC optimize(3)  //O2优化开启
#include<bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
typedef pair<int,int> PII;
const int mod=1e9+7;
const int MX=0x3f3f3f3f3f3f3f3f; 
// static char buf[100000],*pa=buf,*pd=buf;
// #define gc pa==pd&&(pd=(pa=buf)+fread(buf,1,100000,stdin),pa==pd)?EOF:*pa++
// inline int read()
// {
//     register int x(0);register char c(gc);
//     while(c<'0'||c>'9')c=gc;
//     while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=gc;
//     return x;
// }
int m,d;
string L,R;
int dp[2005][2005];
void icealsoheat(){
    cin>>m>>d;
    cin>>L>>R;
    L=' '+L;
    R=' '+R;
    int cnt=0;
    auto dfs=[&](auto self,int id,int f,string &s,int sum)->int{
        if(id>=s.size()){
            return sum==0;
        }
        if(f==0&&dp[id][sum]!=-1)return dp[id][sum];
        int res=9;
        if(f)res=s[id]-'0';
        int bns=0;
        for(int i=0;i<=res;i++){
            if(i==d&&(id&1))continue;
            if(i!=d&&!(id&1))continue;
            int sux=(sum*10+i)%m;
            bns=(bns+self(self,id+1,f&(i==res),s,sux))%mod;
        }
 
        return dp[id][sum]=bns;
    };
 
    int ans=0;
    memset(dp,-1,sizeof dp);
    ans=dfs(dfs,1,1,R,0);
    // cout<<ans<<"++++\n";
    memset(dp,-1,sizeof dp);
    ans=(ans-dfs(dfs,1,1,L,0))%mod+mod;
    ans%=mod;
    // cout<<ans;
    bool ff=0;
    int summ=0;
    for(int i=1;i<L.size();i++){
        summ*=10;
        summ+=L[i]-'0';
        summ%=m;
        if(L[i]-'0'!=d&&!(i&1)){
            ff=1;
        }
        else if((i&1)&&L[i]-'0'==d)ff=1;
    }
    if(!ff&&summ==0)ans++;
    ans%=mod;
    cout<<ans;
 
}
signed main(){
    ios::sync_with_stdio(false);          //int128不能用快读！！！！！！
    cin.tie();
    cout.tie();
    int _yq;
    _yq=1;
    // cin>>_yq;
    while(_yq--){
        icealsoheat();
    }
}

寒冰也是热
关注
  • 3
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
数位DP 详解+模版
glorious_dream的博客
02-21 786
首先清楚数位DP要解决什么样的问题： 求出在给定区间 [A,B] 内，符合条件 f(i) 的数 i 的个数。条件 f(i) 一般与数的大小无关，而与数的组成有关。由于数是按位dp，数的大小对复杂度的影响很小。 用记忆化搜索来实现。 先来看模板 dfs里面关键的几个变量（示题而定） 1.pos 位数 这个很好理解，就是位数，往下搜的时候把位数改变即可 2.limit 最高位标记 什么意思？ 比如要搜索1-555之间的数，如果最高位填5了，那么第二位的取值只有0-5，而最高位是1-4的时候，
leetcode530-Leetcode:Leetcode刷题记录
06-30
线性数结构 链表 2号 2号 串串 No.647回文 堆 32号 42号 94号 155号 滚动 二叉树 第113号： 第1028号： 第404号： 第530号： 第538号： 第968号： 第617号： dfs/bfs 17号 37号 39号 40号 46号 47号 126号 No.543 No...
数位DP模板
RBS的专栏
08-06 957
数位DP问题数位DP指一类问题：给定正整数区间[s, e]，问符合条件的数一共有多少个。例如hdu2089以及hdu3555等。 其中，hdu2089的条件为数字中不含4且不含62，hdu3555的条件为数字中包含49。一般而言，这类条件中至少有一个是与整个数的数值无关的，而是与每一位的数字有关。例如上面2道题的条件就与整个数的数值无关。 因此这类问题实际上是数位有关的问题，另一方面，这类问题通
数位dp
cccccccj
11-17 1万+
原文地址： https://blog.csdn.net/brazy/article/details/77427699 觉得写得很不错~~ ============================= 数位dp是一种计数用的dp，一般就是要统计一个区间[le,ri]内满足一些条件数的个数。所谓数位dp，字面意思就是在数位上进行dp咯。数位还算是比较好听的名字，数位的含义：一个数有个位、十位、百位、千位...
数位dp总结 之 从入门到模板
热门推荐
努力
08-03 7万+
for(int i=le;i<=ri;i++) if(right(i)) ans++; 基础篇 数位dp是一种计数用的dp，一般就是要统计一个区间[le,ri]内满足一些条件数的个数。所谓数位dp，字面意思就是在数位上进行dp咯。数位还算是比较好听的名字，数位的含义：一个数有个位、十位、百位、千位......数的每一位就是数位啦！ 之所以要引入数位的概念完全就是为了dp。
数位DP - 带49的数
伏城之外的博客
12-02 1528
带49的数，数位DP
数位DP总结
DIDCJS的博客
04-03 1469
数位DP基于动态的思想，记录状态以至于不用重复的计算，能够处理大规模的数 其时间复杂度为O（状态数*转移数） //状态数是dp数组的大小，转移数是for循环大小（这个时间我也不是很想得通，某大佬讲的。。） 基本模板如下： 基本的三种题型： 1）： 求1~n数字中包含某些数字特征的数量和（如数字子串中是否有13等） 2）：求1~n数字中各个位数相加能够整除某个数的数量和（如123各
数位dp入门详解
CillyB的博客
02-19 4653
基础篇 数位dp是一种计数用的dp，一般就是要统计一个区间[le,ri]内满足一些条件数的个数。所谓数位dp，字面意思就是在数位上进行dp咯。数位还算是比较好听的名字，数位的含义：一个数有个位、十位、百位、千位......数的每一位就是数位啦！ 之所以要引入数位的概念完全就是为了dp。数位dp的实质就是换一种暴力枚举的方式，使得新的枚举方式满足dp的性质，然后记忆化就可以了。
数位DP.
qq_45896050的博客
02-10 1924
数位DP
动态规划——数位dp
qq_57150526的博客
01-05 3919
数位是指把一个数字按照个、十、百、千等等一位一位地拆开，关注它每一位上的数字。如果拆的是十进制数，那么每一位数字都是 0~9，其他进制可类比十进制。
数位dp介绍
08-13 501
不了解dp的可以先看一下dp 数位dp含义： 数位：一个数有个位，十位，百位，千位等等，数的每一位都是数位。 数位dp归为计数dp，是在数位上进行操作的dp。 数位dp的实质是一种快速枚举的方式，它满足dp的性质，然后进行记忆化搜索。 用途： 有两个数，两个数范围很大(例如1e9，甚至更大)，求这两个数符合限定条件的个数。纯暴力不行，就要用数位dp。 例子：求从0到n,（n...
javalruleetcode-leetcode-diary:leetcode刷题记录
06-29
*寻找两个有序数组的中位数 #4 *最长回文子串 #5 Z 字形变换 #6 整数反转 #7 字符串转换整数 (atoi) #8 回文数 #9 *正则表达式匹配 #10 回溯 dp 盛最多水的容器 #11 整数转罗马数字 #12 罗马数字转整数 #13 最长公共...
上班时间刷leetcode-leetcode:记录leetcode的刷题过程和解题思路
06-30
两个有序数组中位数问题 利用递归减少问题规模 较难 花了较长的时间来解决该题 Date : 2018-09-07 Q5 最长回文子串问题 利用了DP的思想 Date : 2018-09-10 Q6 按一定规则转化字符串问题 无难点 按部就班即可 Date : ...
leetcode76-Hihocode:Hihocode经典算法问题
06-30
非法二进制数（数位dp） 1469 福字 计算几何 1040 矩形判断 树 1049 后续遍历 1062 最近公共祖先·一(LRC) 1325 平衡树·Treap AC 图 1066 无间道之并查集(并查集) 1307 穿越禁区（并查集） 1502 最大子矩阵（元素和...
leetcode中国-leetcode:力码
06-29
计数问题　(数位dp, 朴素算法) P1047 校门外的树　(线段树, 区间合并, 分块思想) P1161 开灯 (异或) P1321 单词覆盖还原　(想想这种方法,怎么这么简单!!) P2415 集合求和 (数学) P1518 两只塔姆沃斯牛 The Tamworth ...
【力扣一轮】数组-旋转矩阵
诚。
05-14 233
这时需要用到循环不变量，它是一个循环，每次处理的时候需要处理的方式是一样的，而不是有着条件的限制。2.对于进入再次循环的条件，没有考虑到，具体体现在for循环中用到了n，而n是作为退出条件改变了。对于这道题，循环不变量就是边元素的处理，只处理当前边的起始元素，不处理当前边的末尾元素。3.对于奇数矩阵和偶数矩阵的处理，没有考虑到中间值。对于奇数矩阵还是偶数矩阵，只是多加一个元素指定值。以上4点，导致看了4次卡尔老师的代码，才完成。1.对于循环退出条件，没有考虑到，需要设置。调试之后，发现自己欠缺考虑有几点。
day09
2202_76101487的博客
05-17 722
正值少年的日常打卡！！！
【算法】【二叉树，DFS，哈希集合，分类讨论】力扣1110. 删点成林
未来就在脚下。
05-20 491
本题通过DFS遍历二叉树，结合分类讨论的方法，逐步删除指定节点并生成新的森林。该算法有效地处理了节点删除后树结构的调整问题，并通过哈希集合优化删除判断的时间复杂度，最终实现了高效的解决方案。
力扣62 不同路径 Java版本
最新发布
m0_47066863的博客
05-20 185
机器人每次只能向下或者向右移动一步。机器人试图达到网格的右下角（在下图中标记为 “Finish” ）。一个机器人位于一个 m x n 网格的左上角 （起始点在下图中标记为 “Start” ）。从左上角开始，总共有 3 条路径可以到达右下角。题目数据保证答案小于等于 2 * 109。输入：m = 3, n = 7。输入：m = 3, n = 2。输入：m = 7, n = 3。输入：m = 3, n = 3。问总共有多少条不同的路径？
数位dp模板
05-27
数位 DP 是一类动态规划算法，用于计算某个数的某种性质。通常情况下，数位 DP 需要将数分解成某种形式，然后逐位进行计算。下面是一个简单的数位 DP 模板： ```python def digit_dp(pos, sum, limit): if pos == -1: return sum if not limit and dp[pos][sum] != -1: return dp[pos][sum] up = 9 if limit else digit[pos] res = 0 for i in range(up + 1): res += digit_dp(pos - 1, sum + i, limit and (i == up)) if not limit: dp[pos][sum] = res return res ``` 其中，`pos` 表示当前处理的位数，`sum` 表示当前已经处理的位数上的数的性质，`limit` 表示是否达到上界。在每一次递归中，我们枚举当前位数上可以取到的数，然后递归计算下一位的结果。在计算结果的过程中，需要根据 `limit` 的值来判断是否达到上界，以便得到正确的结果。另外，还需要记录已经计算过的状态，以便避免重复计算。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值