数位dp

最新推荐文章于 2024-04-04 08:54:13 发布
最新推荐文章于 2024-04-04 08:54:13 发布
阅读量903 收藏
点赞数
分类专栏： ACM-ICPC 2016暑期特训1 动态规划——经典题例
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/XerxesSimon/article/details/68060430
版权

讲解：

http://blog.csdn.net/niuox/article/details/9864199

http://m.blog.csdn.net/article/details?id=52100392

第一题 http://acm.timus.ru/problem.aspx?space=1&num=1057

#include<bits/stdc++.h>
using namespace std;
int f[50][50],g[50];
int cal(int x,int k,int b)
{
    int len=0,res=0;
    while(x) g[++len]=x%b,x/=b;
    for(int i=len;i;i--)
    {
        if(g[i]>1) res+=f[i-1][k]+f[i-1][k-1];
        else if(g[i]==1) res+=f[i-1][k],k--;
        if(k<0) return res;
    }
    return res;
}
int main()
{
    int x,y,k,b;
    for(int i=0;i<33;i++)
    {
        f[i][0]=1;f[i][1]=i;
        for(int j=2;j<=i;j++) f[i][j]=f[i-1][j]+f[i-1][j-1];
    }
    scanf("%d%d%d%d",&x,&y,&k,&b);
    printf("%d\n",cal(y+1,k,b)-cal(x,k,b));
    return 0;
}

第二题 http://www.lydsy.com/JudgeOnline/problem.php?id=1026 

#include <iostream>  
#include <stdio.h>  
#include <stdlib.h>  
#include <string.h>  
#include <math.h>  
#include <map>  
#include <queue>  
#include <algorithm>  
using namespace std;  
  
int A[12];  
  
int f[12][10];  
  
//f[i][j]代表长度为i,最高位为j的windy数个数  
void init()  
{  
   memset(f,0,sizeof(f));  
   for(int i=0;i<10;i++) f[1][i] = 1;  
   for(int i=2;i<=10;i++)  
   {  
      for(int j=0;j<10;j++)  
      {  
         for(int k=0;k<10;k++)  
         {  
            if(abs(j-k)>1) f[i][j] += f[i-1][k];  
         }  
      }  
   }  
}  
//(0,a)范围内的windy数个数  
int calc(int a)  
{  
   int m = 0;  
   while(a)  
   {  
      A[m++] = a%10;  
      a/=10;  
   }  
   int ans = 0;  
   //先处理长度小于m的windy数的个数  
   for(int i=1;i<m;i++)  
   {  
      //题目要求不含前导0  
      for(int j=1;j<10;j++)  
      {  
         ans += f[i][j];  
      }  
   }  
   //长度等于m且最高位和原数不同且小于原数的windy数  
   for(int j=1;j<A[m-1];j++) ans += f[m][j];  
   //依次循环将最高位 变为和原数相同  
   for(int i=m-1;i>=1;i--)  
   {  
      for(int j=0;j<A[i-1];j++)  
      {  
         if(abs(j-A[i]) > 1) ans += f[i][j];  
      }  
      if(abs(A[i] - A[i-1])<=1) break;  
   }  
   return ans;  
}  
  
  
int main()  
{  
   #ifndef ONLINE_JUDGE  
      freopen("in.txt","r",stdin);  
   #endif  
   int a,b;  
   init();  
   while(scanf(" %d %d",&a,&b)!=EOF)  
   {  
      int ans = calc(b+1) - calc(a);  
      printf("%d\n",ans );  
   }  
   return 0;  
}


第三题 http://acm.hdu.edu.cn/showproblem.php?pid=2089

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+10;
typedef long long ll;
ll f[20][2],g[20];
ll dfs(int pos,int pre,int sta,bool limit)
{
    if(pos==-1) return 1;
    if(!limit&&f[pos][sta]!=-1) return f[pos][sta];
    int up=limit?g[pos]:9;
    ll res=0;
    for(int i=0;i<=up;i++)
    {
        if(pre==6&&i==2) continue;
        if(i==4) continue;
        res+=dfs(pos-1,i,i==6,limit&&i==g[pos]);
    }
    if(!limit) f[pos][sta]=res;
    return res;
}
ll cal(ll x)
{
    ll pos=0;
    while(x) g[pos++]=x%10,x/=10;
    return dfs(pos-1,-1,0,1);
}
int main()
{
    ll l,r;
    memset(f,-1,sizeof f);
    while(~scanf("%lld%lld",&l,&r),l|r)
    {
        printf("%lld\n",cal(r)-cal(l-1));
    }
    return 0;
}

第四题 http://acm.hdu.edu.cn/showproblem.php?pid=3555 

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+10;
typedef long long ll;
ll f[20][2],g[20];
ll dfs(int pos,int pre,int sta,bool limit)
{
    if(pos==-1) return 1;
    if(!limit&&f[pos][sta]!=-1) return f[pos][sta];
    int up=limit?g[pos]:9;
    ll res=0;
    for(int i=0;i<=up;i++)
    {
        if(pre==4&&i==9) continue;
        res+=dfs(pos-1,i,i==4,limit&&i==g[pos]);
    }
    if(!limit) f[pos][sta]=res;
    return res;
}
ll solve(ll x)
{
    ll pos=0;
    while(x) g[pos++]=x%10,x/=10;
    return dfs(pos-1,-1,0,1);
}
int main()
{
    ll x,t;
    scanf("%lld",&t);
    memset(f,-1,sizeof f);
    while(t--)
    {
        scanf("%lld",&x);
        printf("%lld\n",x+1-solve(x));
    }
    return 0;
}


第五题 http://acm.hdu.edu.cn/showproblem.php?pid=3709 

#include<bits/stdc++.h>// http://acm.hdu.edu.cn/showproblem.php?pid=3709
using namespace std;
const int N=1e3+10;
typedef long long ll;
ll f[20][20][2005],g[20];
ll dfs(int pos,int pre,int sum,bool limit)
{
    if(pos<0) return sum==0;
    if(sum<0) return 0;
    if(!limit&&f[pos][pre][sum]!=-1) return f[pos][pre][sum];
    int up=limit?g[pos]:9;
    ll res=0;
    for(int i=0;i<=up;i++)
    {
        res+=dfs(pos-1,pre,sum+i*(pos-pre),limit&&i==g[pos]);
    }
    if(!limit) f[pos][pre][sum]=res;
    return res;
}
ll cal(ll x)
{
    if(x<0) return 0;
    ll len=0,res=0;
    while(x) g[len++]=x%10,x/=10;
    for(int i=0;i<len;i++)
        res+=dfs(len-1,i,0,1);
    return res-len+1;//排除掉0,00,000....这些情况
}
int main()
{
    ll x,y,t;
    scanf("%lld",&t);
    memset(f,-1,sizeof f);
    while(t--)
    {
        scanf("%lld%lld",&x,&y);
        printf("%lld\n",cal(y)-cal(x-1));
    }
    return 0;
}


第六题 http://acm.hit.edu.cn/hoj/problem/view?id=1983 
#include<bits/stdc++.h>// http://acm.hit.edu.cn/hoj/problem/view?id=1983
using namespace std;
const int N=1e3+10;
typedef long long ll;
ll f[20][2525][55],g[20],ha[2525];
int gcd(int x,int y)
{
    return y?gcd(y,x%y):x;
}
ll dfs(int pos,int mod,int lcm,bool limit)
{
    if(pos<0) return mod%lcm==0;
    if(!limit&&f[pos][mod][ha[lcm]]!=-1) return f[pos][mod][ha[lcm]];
    int up=limit?g[pos]:9;
    ll res=0;
    for(int i=0;i<=up;i++)
    {
        res+=dfs(pos-1,(mod*10+i)%2520,i?(i*lcm/gcd(i,lcm)):lcm,limit&&i==g[pos]);
    }
    if(!limit) f[pos][mod][ha[lcm]]=res;
    return res;
}
ll cal(ll x)
{
    if(x<0) return 0;
    ll len=0;
    while(x) g[len++]=x%10,x/=10;
    return dfs(len-1,0,1,1);
}
int main()
{
    ll x,y,t;
    memset(f,-1,sizeof f);
    int id=0;
    for(int i=1;i*i<2520;i++)
    if(2500%i==0)
    {
        ha[i]=id++;
        if(i*i!=2520) ha[2520/i]=id++;
    }
    while(~scanf("%lld%lld",&x,&y))
    {
        printf("%lld\n",cal(y)-cal(x-1));
    }
    return 0;
}

第七题 http://acm.hdu.edu.cn/showproblem.php?pid=4507 

#include<bits/stdc++.h>// http://acm.hdu.edu.cn/showproblem.php?pid=4507
using namespace std;
const int N=1e3+10;
const int p=1e9+7;
typedef long long ll;
struct node
{
    ll n,s,sq;
}f[20][10][10];
ll g[20],pw[2525];
node dfs(int pos,int mod,int sum,bool limit)
{
    if(pos<0)
    {
        node tmp;
        tmp.n=(mod!=0&&sum!=0);
        tmp.s=tmp.sq=0;
        return tmp;
    }
    if(!limit&&f[pos][mod][sum].n!=-1)  return f[pos][mod][sum];
    int up=limit?g[pos]:9;
    node res,tmp;
    res.n = res.s = res.sq = 0;
    for(int i=0;i<=up;i++)
    {
        if(i==7) continue;
        tmp=dfs(pos-1,(mod*10+i)%7,(sum+i)%7,limit&&i==g[pos]);
        res.n=(res.n+tmp.n)%p;
        res.s=(res.s+tmp.s+((i * pw[pos])%p * tmp.n) % p) % p ;
        res.sq=(res.sq+tmp.sq+((2*i*pw[pos])%p*tmp.s)%p +
                (((i*i*pw[pos])%p*pw[pos])%p*tmp.n)%p)%p;
    }
    if(!limit) f[pos][mod][sum]=res;
    return res;
}
ll cal(ll x)
{
    if(x<0) return 0;
    ll len=0;
    while(x) g[len++]=x%10,x/=10;
    return dfs(len-1,0,0,1).sq;
}
int main()
{
    ll x,y,t;
    memset(f,-1,sizeof f);
    pw[0]=1;
    for(int i=1;i<=19;i++)
    {
        pw[i]=(pw[i-1]*10)%p;
    }
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&x,&y);
        printf("%lld\n",(cal(y)-cal(x-1)+p)%p);
    }
    return 0;
}

第八题 https://www.bnuoj.com/v3/problem_show.php?pid=52325   （他的oj现在好像过不去了）

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
//位数  数字 关系
ll f[20][10][5],sum[20],g[20];
ll dfs(ll pos,ll num,ll lis,ll flag)
{
	if(pos<=0) return 1;
	if(f[pos][num][lis]!=-1&&!flag) return f[pos][num][lis];
	ll x=flag?g[pos]:9;
	ll res=0;
	for(ll i=0;i<=x;i++)
	{
		if(i<=num&&lis==0)  res+=dfs(pos-1,i,0,flag&&(i==x));
		if(i==num&&lis==1)  res+=dfs(pos-1,i,1,flag&&(i==x));
		if(i>=num&&lis==2)  res+=dfs(pos-1,i,2,flag&&(i==x));
	}
	if(!flag) f[pos][num][lis]=res;
	return res;
}
ll solve(ll x)
{
	ll res=0,pos=0;
	memset(g,0,sizeof g);
	while(x)
	{
		g[++pos]=x%10;
		x/=10;
	}
	pos--;  x=g[pos+1];
	for(ll i=1;i<x;i++)
		res+=dfs(pos,i,0,0)+dfs(pos,i,2,0)-dfs(pos,i,1,0);
	res+=dfs(pos,x,0,1)+dfs(pos,x,2,1)-dfs(pos,x,1,1);
	return res+sum[pos];
}
int main()
{
	ll t,tt=0,x=0,y,z;
	memset(f,-1,sizeof f);
	for(ll i=1;i<20;i++)
	{
		x=x*10+9;
		sum[i]=solve(x);
	}
	scanf("%lld",&t);
	while(t--)
	{
		scanf("%lld%lld",&x,&y);
		printf("%lld\n",solve(y)-solve(x-1));
	}
}


第九题 http://gdutcode.sinaapp.com/problem.php?cid=1057&pid=6

#include<bits/stdc++.h>
using namespace std;
const int N=20000+10;
const int inf=0x3f3f3f3f;
typedef long long ll;
ll f[20][20][10][2][2][2],g[20],num[20];
//代表当前第 位上一位是什么，增过，减过，            长度为 len   回文?
ll rec(int pos,int pre,int up,int down,int flag,int q,int len,int ispa)
{
    if(pos<0) return up&&down&&ispa;
    if(~f[pos][len][pre][up][down][ispa]&&!flag&&!q)
    return f[pos][len][pre][up][down][ispa];
    ll x=flag?g[pos]:9;
    ll res=0;
    for(int i=0;i<=x;i++)
    {
        num[pos]=i;
        if(q) res+=rec(pos-1,i,0,0,i<x?0:flag,q&&i==0,len-(q&&i==0),ispa);
        else if(i==pre)
        {
            if(ispa&&pos<len/2)
            res+=rec(pos-1,i,up,down,i<x?0:flag,q&&i==0,len,i==num[len-pos-1]);
            else res+=rec(pos-1,i,up,down,i<x?0:flag,q&&i==0,len,ispa);
        }
        else if(i>pre)
        {
            if(!down) continue;
            if(ispa&&pos<len/2)
            res+=rec(pos-1,i,1,down,i<x?0:flag,q&&i==0,len,i==num[len-pos-1]);
            else res+=rec(pos-1,i,1,down,i<x?0:flag,q&&i==0,len,ispa);
        }
        else if(i<pre)
        {
            if(up) continue;
            if(ispa&&pos<len/2)
            res+=rec(pos-1,i,up,1,i<x?0:flag,q&&i==0,len,i==num[len-pos-1]);
            else res+=rec(pos-1,i,up,1,i<x?0:flag,q&&i==0,len,ispa);
        }
    }
    if(!flag&&!q) f[pos][len][pre][up][down][ispa]=res;
    return res;
}
ll cal(ll x)
{
    int len=0;
    while(x) g[len++]=x%10,x/=10;
    return rec(len-1,0,0,0,1,1,len,1);
}
int main()
{
    ll t,l,r;
    memset(f,-1,sizeof f);
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&l,&r);
        printf("%lld\n",cal(r)-cal(l-1));
    }
    return 0;
}


XerxesSimon
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
最短路四种算法模板
qq_43416494的博客
08-27 153
1.floyed算法求多源最短路 #include<bits/stdc++.h> using namespace std; const int maxn = 1e3 + 10; #define INF 0x3f3f3f3f int mp[maxn][maxn];//记录图之间的路径 int main() { int n,m;//n表示节点数，m表示边 while(s...
ccf201604-1 折点计数
青鱼的博客
09-13 226
#include<bits/stdc++.h> using namespace std; const int maxn = 1e3 + 10; int a[maxn]; int main(){ int n,ans = 0; scanf("%d",&n); for(int i = 0;i < n;++i) scanf("%d",&a[i]); fo...
算法---树状数组
gaoqiandr的博客
11-04 402
对于算法竞赛中遇到的区间问题的求解，如区间最值问题，区间修改单点求值问题等，一般都会超时，本文采用树状数组来求解，简单易懂还不超时更可以AC，属于模板类型。
c语言中1e-10啥意思,【1e+10是多大的数？那个e是什么含义？】-1e-
热门推荐
weixin_32869687的博客
05-21 3万+
概述：本道作业题是柯汾汲同学的课后练习，分享的知识点是1e，指导老师为万老师，涉及到的知识点涵盖：【1e+10是多大的数？那个e是什么含义？】-1e-物理，下面是柯汾汲作业题的详细。题目：【1e+10是多大的数？那个e是什么含义？】-1e-物理楼上的说得都不是太正确,这里的e应该不是通常说的自然底数我觉得楼主遇到的只是科学计数法中的数1e+10实际上是指1乘以10的10次方例如科学计数法中：234...
第15届蓝桥杯省赛模拟第二期题解及代码实现
最新发布
Li的博客
04-04 1080
问题描述输入一个仅包含小写英文字母的字符串，请问这个字符串中的最后一个元音是什么。在英文中，a,e,i, o,u共5个字母是元音字母，其它字母不是元音字母。输入格式输入一行包含一个字符串，仅由小写英文字符组成，字符串中至少包含一个元音字母。输出格式输出一行包含一个字符，表示答案。最后一个首先想到倒着遍历.本题从字符串的尾部开始遍历，找到第一个元音字母即停止。cin>>s;string t;t=c;cout<<t;return 0;
数位dp算法详细图文解答
02-13
数位动态规划（Digital Dynamic Programming, 简称数位DP）是一种用于解决与数字的组成有关的问题的算法。它通常涉及到对数字的每一位进行操作，通过递推和记忆化搜索来优化计算效率。数位DP的核心思想是将问题分解...
数位DP.pdf
07-12
数位DP整理的模板和一些题目的题解。QAQ没钱下载资源了只好上传点东西。
数位dp——某次队内分享任务
07-21
数位动态规划(Digit Dynamic Programming, 简称数位DP)是一种在处理与数字相关的动态规划问题时，通过分析每个数字位上的状态来优化算法的技术。它通常用于解决涉及数字序列特性的数学问题，例如计数特定类型的数...
数位DP入门
08-22
数位DP入门思想，包括详细的思考过程以及一些习题的讲解
4数位DP1
08-08
然后通过三重循环更新`dp`数组，分别对序列长度、和以及当前数位进行迭代。 总结来说，这个题目涉及到的IT知识点包括： 1. 动态规划（DP）及其记忆化搜索 2. 深度优先搜索（DFS）用于位枚举 3. 状态转移和状态设计...
信息学奥赛一本通 1286：怪盗基德的滑翔翼（evd）
everwide的博客
11-17 955
【题目描述】 怪盗基德是一个充满传奇色彩的怪盗，专门以珠宝为目标的超级盗窃犯。而他最为突出的地方，就是他每次都能逃脱中村警部的重重围堵，而这也很大程度上是多亏了他随身携带的便于操作的滑翔翼。 有一天，怪盗基德像往常一样偷走了一颗珍贵的钻石，不料却被柯南小朋友识破了伪装，而他的滑翔翼的动力装置也被柯南踢出的足球破坏了。不得已，怪盗基德只能操作受损的滑翔翼逃脱。 假设城市中一共有N幢建筑排成一条线，每幢建筑的高度各不相同。初始时，怪盗基德可以在任何一幢建筑的顶端。他可以选择一个方向逃跑，但是不能中途改变方向（因
22.4.10
aidaqiudeaichao的博客
04-10 379
1，滑雪（记忆化搜索） 1，滑雪（记忆化搜索） 题意： 直接dfs会爆tle，所以用dp来做； 闫式dp分析法： 第一步： 状态表示： f[ i, j ] 表示从当前位置开始滑能滑到的最远路径； 属性：max 第二步： 状态计算： 先划分集合： 很明显，根据集合定义可以划分为四个部分： 然后计算子集： 拿向上滑举例： i , j → （滑到） i-1 , j → .. .. . . .只有第一步滑到 i-1 ,j 是确定的，剩下的都确定，又因为求的是最大值，把每个集合都减1对答案没有
算法竞赛入门（C++）
zhr_the_great的博客
07-06 2010
C++版算法竞赛入门笔记
FJNU - 第二十一届低编赛题解
过路人丶的博客
11-26 1158
题目链接:点击打开链接 A: 题意:题意很清楚了，就不解释了。 解题思路:很容易发现这个题目其实跟01排列顺序没有关系长度才是关键，就比如要求步数相同的情况下1001和0000,1111答案数是一样的，这个想一想就很好理解了。 接下来代码有注释说明。 代码: #include #include #include #include #define M 100005 using
C++中数值的表示方法
峰升的专栏
12-08 2878
1. 1e3  表示10的3次方   1e-3  则表示10的3次方的倒数 2.整型常量的表示方法： (1）十进制：数字为0-9，如123; (2)  八进制：以0开头，数字为0-7，如0101表示八进制101Q; (3)  十六进制：以0x或0X开头，数字为0-9，A-F如0x41  41H  0xFA  0FAH;
数据结构课程设计【C++实现】
nefu_ljw的博客
12-12 5900
POJ 3984 迷宫问题 #include <iostream> #include <cstdio> #include <queue> using namespace std; const int N=1e3+10; char a[N][N]; bool vis[N][N]; int dir[4][2]={{1,0},{-1,0},{0,1},{0,-1}};...
c语言1e3和1e3,数据格式1e3是多少
weixin_39988888的博客
05-20 9906
第一题：浮点型数据有两种1.使用小数点表示的数。正确的小数形式有3种：整数位+小数点；小数点+小数位；整数位+小数点+小数位例如： .0 0. 0.0 三个数都是正确的。2.指数形式表示的数。正确的指数形式是 xen 其中x可以是整数也可以是小数(带小数点的)，e或者E代表底数10，n代表幂值。例如： 3e5 等价于 3乘以10的5次幂12.34e-3 等价于 12.34...
最佳淘汰算法
梦醒了
03-29 5567
4272: 最佳淘汰算法 Time Limit: 12 Sec  Memory Limit: 128 MB Submit: 412  Solved: 64 Description Lyf在做操作系统作业的时候，遇到了一题要他实现下Opt算法（全称页面置换最佳淘汰算法）。可是Lyf发现最近好忙啊，于是他想请他的学弟,学妹们帮他来完成。 作业里的内容是这样的： 提出最佳页面淘汰算法
网络流（最大流、最小费用最大流、有上下界的网络流）
梦醒了
08-09 5101
最大流的算法： Ford-Fulkerson算法  求最大流的过程，就是不断找到一条源到汇 的路径，然后构建残余网络，再在残余网络 上寻找新的路径，使总流量增加，然后形成 新的残余网络，再寻找新路径…..直到某个 残余网络上找不到从源到汇的路径为止，最 大流就算出来了。   每次寻找新流量并构造新残余网络的过程， 就叫做寻找流量的“增广路径”，也叫“增 广”  现在假设每条边的容量
数位DP详解：区间统计与C语言基础
文中通过不同的例题解释了数位DP的基本思想、方法，并提供了几个具体实例，如HDU2089等，帮助读者理解和应用这一技术。" 在信息学竞赛和编程问题中，数位DP是一种非常实用的策略，它主要用于处理与数字的各个位相关...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值