蒟蒻的数位DP专题总结

最新推荐文章于 2024-09-13 18:48:07 发布
最新推荐文章于 2024-09-13 18:48:07 发布
阅读量79 收藏
点赞数
文章标签： php
原文链接：http://www.cnblogs.com/zxhl/p/4954444.html
版权

BZOJ  1026: [SCOI2009]windy数：

题目链接： http://www.lydsy.com/JudgeOnline/problem.php?id=1026

          dp[11][11][2]：dep,pre,f

           要求的性质就是相邻数字差至少是2。
          递归函数的状态设计如下dep，pre，f，分别表示已经枚举到第dep位，他的前一位（更高的位）是pre，f表示大小关系是否已经确定.

///1085422276

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define pb push_back
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }return x*f;
}
//****************************************
const int  N=500000+50;
#define mod 1000000007
#define inf 10000007

int dp[11][11][2],vis[11][11][2],d[11];
int dfs(int dep,int pre,int f) {
    if(dep<0) return 1;
    else if(pre>=0&&vis[dep][pre][f]) return dp[dep][pre][f];
    else {
        int re=0;
        if(pre == -1) {
            re+=dfs(dep-1,-1,f||d[dep]>0);
            if(f) for(int i=1;i<=9;i++) re+=dfs(dep-1,i,1);
            else {
                for(int i=1;i<=d[dep];i++) {
                    re+=dfs(dep-1,i,i<d[dep]);
                }
            }
        }
        else  {
            if(f) {
                for(int i=0;i<10;i++) {
                    if(abs(i-pre)>1) {
                        re+=dfs(dep-1,i,1);
                    }
                }
            }
            else {
                for(int i=0;i<10;i++) {
                    if(i<=d[dep]&&abs(i-pre)>1) {
                        re+=dfs(dep-1,i,i<d[dep]);
                    }
                }
            }
        }
        if(pre>=0) vis[dep][pre][f]=1,dp[dep][pre][f]=re;
        return re;
    }
}
int cal(int n) {
   mem(dp);mem(vis);
   int len=0;
   while(n) {
    d[len++]=n%10;n/=10;
   }
   return dfs(len-1,-1,0);
}
int main() {
    int a,b;
    while(scanf("%d%d",&a,&b)!=EOF) {
        cout<<cal(b)-cal(a-1)<<endl;
    }
   return 0;
}
BZOJ1026

 BZOJ 3209: 花神的数论题：

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=3209

             求的是1到N内进制种类数的个数的乘积，显然是数学题嘛。。。。

             仔细想想：我们如果只考虑二进制上01的个数，也就是如果枚举1的个数，将二进制上含有i个1的个数用数位dp求解出来，再算快速幂，就可以得到答案了

///1085422276

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define pb push_back
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }return x*f;
}
//****************************************
const int  N=60+50;
#define mod 10000007
#define inf 10000007
#define maxn 10000

ll dp[65][65][3],vis[65][65][3],d[N],c[N][N],ans[N];;
ll dfs(int dep,int left,int f) {
   if(vis[dep][left][f]) return dp[dep][left][f];
   if(dep+1<left) return 0;if(left==0) return 1;
   ll& re=dp[dep][left][f];
   vis[dep][left][f]=1;
   if(!f) {
       if(d[dep]==0) {
        re+=dfs(dep-1,left,f);
       }
       else {
        re+=dfs(dep-1,left,1);
        re+=dfs(dep-1,left-1,0);
       }
   }
   else {
     re+=c[dep+1][left];
   }//vis[dep][left][f]=1,dp[dep][left][f]=re;
   return re;
}
ll quick_pow(ll n,ll m) {
   ll ret=1;
   while(m) {
     if(m&1) ret=(ret*n)%mod;
     n=(n*n)%mod;
     m>>=1;
   }
   return ret;
}
ll cal(ll n) {
    int len=0;mem(vis);mem(dp);
    while(n) {
        d[len++]=n%2;
        n/=2;
    }
    for (int i=1;i<=len;i++) {
        ans[i]=dfs(len-1,i,0);
    }
    ll A=1;
    for(int i=2;i<=len;i++)A=(A*quick_pow(i,ans[i]))%mod;
    return A;
}
int main(){
       ll n;
       for(int i=0;i<=60;i++) {
        c[i][0]=c[i][i]=1;
        for(int j=1;j<i;j++) {
            c[i][j]=c[i-1][j]+c[i-1][j-1];
        }
       }
       while(scanf("%lld",&n)!=EOF) {
        printf("%lld\n",cal(n));
       }
  return 0;
}
BZOJ3209

BZOJ 1799: [Ahoi2009]self 同类分布：

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1799

              因为最多18位，最大也就是18*9，我们枚举位数和，去数位dp就可以了

///1085422276

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define pb push_back
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }return x*f;
}
//****************************************
const int  N=600+50;
#define mod 10000007
#define inf 10000007
#define maxn 10000

ll dp[20][200][200][2],vis[20][200][200][2];
int d[N];
ll dfs(int dep,int T,int f,int sum,int M) {
   if(dep<0) return T == 0 && sum == M;
   if(sum>M)return 0;
   if(vis[dep][T][sum][f]) return dp[dep][T][sum][f];
   ll& re=dp[dep][T][sum][f];
   vis[dep][T][sum][f]=1;
   if(f) {
     for(int i=0;i<=9;i++) {
        re+=dfs(dep-1,(T*10+i)%M,f,sum+i,M);
     }
   }
   else {
    for(int i=0;i<=d[dep];i++) {
        re+=dfs(dep-1,(T*10+i)%M,i<d[dep],sum+i,M);
    }
   }
   return re;
}

ll cal(ll x) {
   int len=0;
   ll A=0;
   while(x) d[len++]=x%10,x/=10;
   for(int i=1;i<=9*len;i++) {
        mem(dp),mem(vis);
        A+=dfs(len-1,0,0,0,i);
   }
    return A;
}
int main() {
    ll l,r;
    while(scanf("%I64d%I64d",&l,&r)!=EOF) {
        printf("%I64d\n",cal(r)-cal(l-1));
        //cout<<cal(r)-cal(l-1)<<endl;
    }
    return 0;
}
BZOJ1799

 HDU 4734 F(x):

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4734

              我们先求出F(A),题意只要求满足F(x)<=F(A) 的就行，所以我们从构造数位dp  将F(A)设为上限，

              F(A)可以发现最大就是9*2^9多一点，我们再枚举一个数位大小的时候，将剩余上限减去使用当前数位的F()贡献

              最后也就是枚举完，剩余>=0就可行

///1085422276

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define pb push_back
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }return x*f;
}
//****************************************
const int  N=6000+50;
#define mod 10000007
#define inf 10000007
#define maxn 10000

int dp[10][N],vis[10][N],d[20];

int dfs(int dep,int sum,int f) {
    if(sum<0) return 0;
    if(dep<0) return sum>=0;
    if(f&&vis[dep][sum]) return dp[dep][sum];
    if(f) {
        int& re=dp[dep][sum];
        vis[dep][sum]=1;
        for(int i=0;i<=9;i++) {
            re+=dfs(dep-1,sum-i*(1<<dep),f);
        }
        return re;
    }
    else {
        int re=0;
        for(int i=0;i<=d[dep];i++) {
            re+=dfs(dep-1,sum-i*(1<<dep),i!=d[dep]);
        }
        return re;
    }
}
int cal(int A,int B) {
    if(B==0) return 1;
    int sum=0;int len=0;//mem(dp),mem(vis);
    while(A) sum+=(A%10)*(1<<len),len++,A/=10;
    len=0;
    while(B) d[len++]=B%10,B/=10;
    return dfs(len-1,sum,0);
}

int main() {
    int T=read();
    int A,B;
    for(int i=1;i<=T;i++) {
       A=read(),B=read();
        printf("Case #%d: ",i);
        printf("%d\n",cal(A,B));
    }
    return 0;
}
HDU4734

 HDU 3555 Bomb：

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3555

             求1到N内有49相连的数字个数，经典数位；

///1085422276

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define pb push_back
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }return x*f;
}
//****************************************
const int  N=6000+50;
#define mod 10000007
#define inf 10000007
#define maxn 10000

ll dp[60][60],vis[50][60],d[60];
ll dfs(int dep,int pre,int f) {
    if(dep<0) return 1;
    if(f&&vis[dep][pre]) return dp[dep][pre];
    if(f) {
            ll& re=dp[dep][pre];
           vis[dep][pre]=1;
        for(int i=0;i<=9;i++) {
            if(pre!=4||(pre==4&&i!=9)) {
                re+=dfs(dep-1,i,f);
            }
        }return re;
    }
    else {
        ll re=0;
        for(int i=0;i<=d[dep];i++) {
            if(pre!=4||(pre==4&&i!=9)) 
            re+=dfs(dep-1,i,i<d[dep]);
        } return re;
    }
}
ll cal(ll x) {
    int len=0;mem(dp),mem(vis);
    while(x) d[len++]=x%10,x/=10;
    return dfs(len-1,0,0);
}
int main() {

    int T=read();
    while(T--) {
        ll n;
        scanf("%I64d",&n);
        cout<<n-cal(n)+1<<endl;
    }
    return 0;
}
HDU3555

 HDU 3709 Balanced Number：

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3709

             我们枚举平衡点就OK，注意0 的次数重复计算了，减去len-1就是答案

///1085422276

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define pb push_back
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }return x*f;
}
//****************************************
const int  N=6000+50;
#define mod 10000007
#define inf 10000007
#define maxn 10000

ll dp[20][6000][20],vis[20][6000][20],d[N];
ll dfs(int dep,int D,int mid,int f) {
   if(dep<0) return D==3000;
   if(f&&vis[dep][D][mid]) return dp[dep][D][mid];
   if(f) {
     vis[dep][D][mid]=1;
     ll& re=dp[dep][D][mid];
     for(int i=0;i<=9;i++) {
        re+=dfs(dep-1,D-(dep-mid)*i,mid,f);
     }
     return re;
   }
   else {
    ll re=0;
    for(int i=0;i<=d[dep];i++) {
        re+=dfs(dep-1,D-(dep-mid)*i,mid,i<d[dep]);
    }
    return re;
   }

}
ll cal(ll x) {
   if(x==0) return 1;
   if(x<0) return 0;
   int len=0;ll re=0;mem(dp),mem(vis);
   while(x) d[len++]=x%10,x/=10;
   for(int i=0;i<len;i++) {
     re+=dfs(len-1,3000,i,0);
   }
   return re-len+1;
}
int main() {

    int T=read();
    while(T--) {
            ll l,r;
      scanf("%I64d%I64d",&l,&r);
      cout<<cal(r)-cal(l-1)<<endl;
    }
    return 0;
}
HDU3709

转载于:https://www.cnblogs.com/zxhl/p/4954444.html

weixin_30390075
关注
dfs求解数位dp基础题
avenLJY的博客
05-13 490
#数位dp基础题(dfs求解） 本人蒟蒻，第一篇博客，来数位dp几道水题 【例1】洛谷2602 [ZJOI2010]数字计数 大意：求l到r所有数中数字0的个数，数字1的个数……数字9的个数 代码： #include<bits/stdc++.h> using namespace std; #define ll long long ll dp[20][20][20],bit[20],a...
数位DP题目
qq_61130688的博客
10-28 261
f(x)定义为x二进制下1的个数，定义一个m位的01串a，看 [0,L] 有多少数字满足对于任意 i 属于 [0, m-1]，f(x+i) mod2 = ai。
二进制1的个数
[L, R)
03-17 896
给定一个整数x(x >= 1 && x 比如x= 5 1的二进制表示为1 2的二进制表示为10 3的二进制表示为11 4的二进制表示为100 5的二进制表示为101 所以1的总数为7 对于这么大的数据，将所有结果打表保存，内存吃不消，打表的时间也要很久。 我们换个思路，做一张特殊的表sum[33], sum[x]代表:1到2^x - 1内所有1的总数。 求sum的时候我们
蒟蒻的笔记本一、数位dp
diker 的博客
10-12 201
本博客搬自本人的洛谷博客 https://www.luogu.org/blog/diker2018/post-di-bi-ji-ben-yi-shuo-wei-dp 蒟蒻的我学了一天终于稍微明白一点原理了orz。 1.对数位dp的一些理解 真如其名，数位dp就是在数位上做动态规划。 数位dp一般应用于： 求出在给定区间[A,B]内，符合条件P的数的个数。 条件P一般与数的大小无关，而与数的组成有关...
数位 dp 总结
anchengao1503的博客
08-06 177
T1：不要 62 题干： 　　杭州人称那些傻乎乎粘嗒嗒的人为 $62$（音：$laoer$）。 　　杭州交通管理局经常会扩充一些的士车牌照，新近出来一个好消息，以后上牌照，不再含有不吉利的数字了，这样一来，就可以消除个别的士司机和乘客的心理障碍，更安全地服务大众。 　　不吉利的数字为所有含有 $62$ 或 $4$ 的号码。例如： $62315$、$73418$、$88914$ 都...
LuoguP4317(数位dp)
stay_accept的专栏
10-27 177
链接：点击打开链接 题意：话说花神这天又来讲课了。课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了。 花神的题目是这样的：设表示的二进制表示中的个数。给出一个正整数，花神要问你，也就是的乘积。 代码： #include <bits/stdc++.h> using namespace std; typedef long long ll; const ll mo...
windy数（数位DP）
ChrysanDog的博客
10-25 295
题目描述 windywindywindy定义了一种windywindywindy数。不含前导零且相邻两个数字之差至少为222的正整数被称为windywindywindy数。windywindywindy想知道，在AAA和BBB之间，包括AAA和BBB，总共有多少个windywindywindy数？ 输入格式 包含两个整数，AAA，BBB。 输出格式 一个整数 SampleSample...
数位DP入门
BroDrinkWater的博客
08-25 1104
考试考dp的时候时常会碰见有关数位dp的问题，每次考到就是一脸懵逼加吃惊，所以今天抽空看了一下有关数位dp的知识，网上有很多大神都说的很好，推荐看几篇blog。 入门经典 慢慢看，很不错数位DP的套路数位dp其实看了那么多篇blog感觉就是一个套路，一个记忆化搜索，方法无非是用两个端点的答案相减得到答案。dp主要考得是状态的设定和转移及其优化，对于数位dp来说，转移和优化其实是固定的，变的只是状
BZOJ 3209 花神的数论题 数位DP 数论
Chlience的博客
05-30 397
BZOJ 3209 花神的数论题 题目描述 Description 背景 众所周知，花神多年来凭借无边的神力狂虐各大 OJ、OI、CF、TC …… 当然也包括 CH 啦。 描述 话说花神这天又来讲课了。课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了。 花神的题目是这样的 设 sum(i) 表示 i 的二进制表示中 1 的个数。给出一个正整数 N ，花神要问你 ∏(Sum(i))∏(...
算法-数位dp
KonnyWen的博客
02-24 1235
算法-数位dp 先学学。未完待续。
P2602 [ZJOI2010]数字计数（非数位dp做法）
w2319_Willy的博客
12-12 331
题目描述 给定两个正整数a和b，求在[a,b]中的所有整数中，每个数码(digit)各出现了多少次。 输入格式： 输入文件中仅包含一行两个整数a、b，含义如上所述。 输出格式： 输出文件中包含一行10个整数，分别表示0-9在[a,b]中出现了多少次。 输入输出样例 输入 #1 1 99 输出 #1 9 20 20 20 20 20 20 20 20 20 说明/提示 数据规模与约定 对于 30% 的数据，保证 a ≤ b ≤ 106 对于 100% 的数据，保证 1≤ a ≤ b ≤ 1012 思路
[2019年国庆专题训练] dp专题训练
cosini的博客
10-06 640
文章目录T1.生日礼物题目描述样例输入输出题解参考代码T2.大理石题目描述样例输入输出题解参考代码T3.猜数字题目描述样例输入输出题解参考代码T4.加分二叉树题目描述样例输入输出题解参考代码T5.数字计数题目描述样例输入输出题解参考代码T6.逛公园题目描述样例输入输出题解参考代码T7.征途题目描述样例输入输出题解参考代码T8.股票交易题目描述样例输入输出题解参考代码T9.麻将题目描述样例输入输出题...
[SGU390]Tickets && 数位DP
shiyukun1998的专栏
12-23 563
超级难想的一道题 用一棵十叉数来统计  因为是各个数位之和 所以要把之前各个数位的和sum记录下来  因为每个数位（子树）并不一定完整的分成cnt组 最后可能有剩余 所以每一个数位（子树）前面肯能会有多余的空间 用rem记录 因为需要每组末的空间是在统计时不断改变的 所以dp时要同时记录组数和最后的空余 实现方面 沿用上一个记忆化搜索的思路 传入两个参数unlimit1
一码空传临时网盘PHP源码，支持提取码功能
爱酷码的博客
09-11 411
一码空传临时网盘源码V2.0免费授权，该源码提供了一个简单易用的无数据库版临时网盘解决方案。前端采用了layui开发框架，后端使用原生PHP编写，没有引入任何开发框架，保持了代码的简洁和高效。这个程序使用了一个无数据库配置读写类，并借鉴了网络上的config文件读写代码。用户可以通过提取码来提取文件，无需上传文件到服务器。该程序还配备了一个后台管理系统，而且不需要数据库支持。用户上传的文件会保存在upload文件夹中，并按照md5算法进行加密和重命名，确保每次加密结果都不相同。
Redis 多线程模型详解
lssffy的博客
09-11 955
在 Redis 的早期版本中，Redis 采用了典型的单线程模型，即所有的客户端请求都是由 Redis 的主线程依次处理的。具体来说，Redis 使用单线程的事件循环模型来处理客户端的请求和响应，使用了操作系统的selectpoll或epoll来管理多个客户端连接。为什么 Redis 选择单线程？简单性：单线程模型不需要考虑线程间数据竞争问题，避免了锁机制带来的复杂性和潜在的性能开销。避免上下文切换：多线程环境下，线程的上下文切换会引入额外的 CPU 负担，而单线程架构没有这个问题。性能足够。
[极客大挑战 2019]Http
2201_75556700的博客
09-10 583
3、访问Secret.php发现页面有提示，它说它不是来自这个网页。6、添加X-Forwarded-For: 127.0.0.1。5、修改User-Agent为：Syclover。2、查看页面源码发现一个。
SEAFARING靶场渗透
最新发布
kknifek的博客
09-13 446
??echo '<??
Netty笔记03-组件Channel
weixin_46425661的博客
09-12 772
本文主要讲解Netty中的组件Channel、ChannelFuture和CloseFuture
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值