简单的数位dp

2024-08-28 10:18
文章标签 简单 dp 数位

本文主要是介绍简单的数位dp,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

数位dp的模板


hdu 2089


设dp[len][flag],flag = 1表示前一位是6,否则前一位不是6.

#include <stdio.h>
#include <iostream>
#include <map>
#include <set>
#include <list>
#include <stack>
#include <vector>
#include <math.h>
#include <string.h>
#include <queue>
#include <string>
#include <stdlib.h>
#include <algorithm>
#define LL __int64
#define eps 1e-12
#define PI acos(-1.0)
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 50010;int dig[10];
int dp[10][2];int dfs(int len, int flag, int up)
{if(len == 0)return 1;if(!up && dp[len][flag] != -1)return dp[len][flag];int n = up == 1?dig[len]:9;int res = 0;for(int i = 0; i <= n; i++){if(i == 4 || (flag == 1 && i == 2))continue;res += dfs(len-1,i==6?1:0,up&&(i==n) );}if(!up)dp[len][flag] = res;return res;
}int cal(int num)
{int cnt = 0;while(num){dig[++cnt] = num%10;num /= 10;}memset(dp,-1,sizeof(dp));return dfs(cnt,0,1);
}
int main()
{int a,b;while(~scanf("%d %d",&a,&b)){if(a == 0 && b == 0)break;printf("%d\n",cal(b) - cal(a-1) );}return 0;
}


hdu 3555


不能出现49,与上题类似

#include <stdio.h>
#include <iostream>
#include <map>
#include <set>
#include <list>
#include <stack>
#include <vector>
#include <math.h>
#include <string.h>
#include <queue>
#include <string>
#include <stdlib.h>
#include <algorithm>
#define LL __int64
#define eps 1e-12
#define PI acos(-1.0)
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 4010;int dig[100];
LL dp[100][4];LL dfs(int len, int flag, int up)
{if(len == 0)return flag == 2;if(!up && dp[len][flag] != -1)return dp[len][flag];int n = up ? dig[len] : 9;LL res = 0;for(int i = 0; i <= n; i++){int tflag = flag;if(flag == 0 && i == 4)tflag = 1;else if(flag == 1 && i == 9)tflag = 2;else if(flag == 1 && i != 4)tflag = 0;res += dfs(len-1,tflag,up&&(i==n));}if(!up)dp[len][flag] = res;return res;
}LL cal(LL num)
{int len = 0;while(num){dig[++len] = num % 10;num /= 10;}memset(dp,-1,sizeof(dp));return dfs(len,0,1);
}int main()
{int test;LL n;scanf("%d",&test);while(test--){scanf("%I64d",&n);printf("%I64d\n",cal(n));}return 0;
}


SCOI2009


要求两个相邻的数之差大于等于2,对于第一位是特殊的。所以dfs的时候加一个变量first表示前面是否全是前导0。

#include <stdio.h>
#include <iostream>
#include <map>
#include <set>
#include <list>
#include <stack>
#include <vector>
#include <math.h>
#include <string.h>
#include <queue>
#include <string>
#include <stdlib.h>
#include <algorithm>
#define LL __int64
#define eps 1e-12
#define PI acos(-1.0)
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 4010;int dig[15];
int dp[15][15]; //dp[i][j]表示第i位它的上一位即i+1位为jint dfs(int len, int pre, int up, int first)
{if(len == 0)return first == 0;if(!up && dp[len][pre] != -1 && !first)return dp[len][pre];int n = up ? dig[len] : 9;int res = 0;for(int i = 0; i <= n; i++){if(first) //当前这一位是第一位,当i==0时,第i-1为也是第一位,否则不是res += dfs(len-1,i,up&&i==n,first&&i==0);else if(abs(pre-i) >= 2)res += dfs(len-1,i,up&&i==n,first);}if(!up && !first)dp[len][pre] = res;return res;
}int cal(int num)
{int len = 0;while(num){dig[++len] = num%10;num /= 10;}memset(dp,-1,sizeof(dp));return dfs(len,0,1,1);
}int main()
{int a,b;while(~scanf("%d %d",&a,&b)){printf("%d\n",cal(b) - cal(a-1));}return 0;}






这篇关于简单的数位dp的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1114485

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

hdu4865(概率DP)

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。 思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] )  代码如下: #include <stdio.h>#include

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

uva 10387 Billiard(简单几何)

题意是一个球从矩形的中点出发,告诉你小球与矩形两条边的碰撞次数与小球回到原点的时间,求小球出发时的角度和小球的速度。 简单的几何问题,小球每与竖边碰撞一次,向右扩展一个相同的矩形;每与横边碰撞一次,向上扩展一个相同的矩形。 可以发现,扩展矩形的路径和在当前矩形中的每一段路径相同,当小球回到出发点时,一条直线的路径刚好经过最后一个扩展矩形的中心点。 最后扩展的路径和横边竖边恰好组成一个直

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 10154 DP 叠乌龟

题意: 给你几只乌龟,每只乌龟有自身的重量和力量。 每只乌龟的力量可以承受自身体重和在其上的几只乌龟的体重和内。 问最多能叠放几只乌龟。 解析: 先将乌龟按力量从小到大排列。 然后dp的时候从前往后叠,状态转移方程: dp[i][j] = dp[i - 1][j];if (dp[i - 1][j - 1] != inf && dp[i - 1][j - 1] <= t[i]