• 数位DP 上 day44


    今天算是摆了半天 又睡了四分之一天

    最后无法忍受没有训练的自己

    还是来写了几道数位dp

    这个的思考是很有迹可循的

    值得注意的就是很多细节

    又是五点了 该觉觉了zzz

    1081. 度的数量

    我也是第一次写数位dp 不会很正常把 呃呃

    大概的重点就是把他看成一个b进制数 我们找这个数有多少个符合条件的

    符合条件就相当于<=x 并且其中有k个1

    我们像树一样check 每一位 要是此位为>=2 那么我们 0 1 就可以算出所有的符合条件的数 直接break 就可以了

    last 是 前面有多少个1

    要是为0的话就直接下一位了

    要是为1的话我们又要算0的又要算1的符合的数 但是1作为开头 那我们必须去下一层

    还要last++

    最后我们还要算右边的值 当且仅当k==last&&i==0 要是k不等于last 

    要吗last>k 早就break了 我们不要那么多1 意思符合的数早就算完了 再多也是不合法的(1太多了

    要吗k>last 那么我们自然还有下一个分支 我们可以想出来肯定是0的

    1. #include
    2. using namespace std;
    3. const int N = 35;
    4. const int M = 1<<12;
    5. const int mod = 1e9+7;
    6. #define int long long
    7. #define endl '\n'
    8. #define Endl '\n'
    9. #define _ 0
    10. #define inf 0x3f3f3f3f3f3f3f3f
    11. #define fast ios::sync_with_stdio(false);cin.tie(nullptr);
    12. int f[N][N],x,y,k,b;
    13. void init(){
    14. for(int i=1;i<=N;i++){
    15. for(int j=0;j<=i;j++){
    16. if(j==0||i==j)f[i][j]=1;
    17. else f[i][j]=f[i-1][j-1]+f[i-1][j];
    18. }
    19. }
    20. }
    21. int dp(int x){
    22. if(!x)return 0;
    23. int last=0,res=0;
    24. vector<int>nums;
    25. while(x)nums.push_back(x%b),x/=b;
    26. for(int i=nums.size()-1;i>=0;i--){
    27. int r=nums[i];
    28. if(r){
    29. res+=f[i][k-last];
    30. if(r>1){
    31. if(k-last-1>=0)res+=f[i][k-last-1];
    32. break;
    33. }else{
    34. last++;
    35. if(last>k)break;
    36. }
    37. }
    38. if(i==0&&last==k)res++;
    39. }
    40. return res;
    41. }
    42. signed main(){
    43. fast
    44. init();
    45. cin>>x>>y>>k>>b;
    46. cout<<dp(y)-dp(x-1)<
    47. return ~~(0^_^0);
    48. }

    1082. 数字游戏

    和上题一样的树

    不过我们预处理需要用dp来求

    状态表示 我们需要一个首位的值 还有一个位数 

    说实话我本来以为要三维 可是发现 我的三维表示还是错的

    其实两位就够了

    f[i][j]表示i位数字首位为j的cnt

    状态计算 f[i][j]+=f[i-1][j-9];

    注意的是我们0为首的也要算 不然会有答案会有遗漏

    last注意初始化为较小的数 也可以是负无穷 加一行特判就是了

    然后我们每次都要更新last啊 要是我们后面的数比last还要小了 那后面就再也没有合法ans了

    要是比last大或者等于last 那就要继续做

    当我们枚举到最后一层时 我们必须加上右边的值 就是它本身 (因为其本身就是一个不降数就会有

    1. #include
    2. using namespace std;
    3. const int N = 15;
    4. const int M = 1<<12;
    5. const int mod = 1e9+7;
    6. #define int long long
    7. #define endl '\n'
    8. #define Endl '\n'
    9. #define _ 0
    10. #define inf 0x3f3f3f3f3f3f3f3f
    11. #define fast ios::sync_with_stdio(false);cin.tie(nullptr);
    12. int f[N][N],l,r;//f[i][j]表示i位数字首位为j的cnt
    13. void init(){
    14. for(int i=0;i<=9;i++)f[1][i]=1;
    15. for(int i=2;i<=N;i++){
    16. for(int j=0;j<=9;j++){
    17. for(int k=j;k<=9;k++){
    18. f[i][j]+=f[i-1][k];
    19. }
    20. }
    21. }
    22. }
    23. int dp(int n){
    24. if(!n)return 1;
    25. int last=0,res=0;
    26. vector<int>v;
    27. while(n)v.push_back(n%10),n/=10;
    28. for(int i=v.size()-1;i>=0;i--){
    29. int x=v[i];
    30. for(int j=last;j
    31. res+=f[i+1][j];
    32. }
    33. if(last>x)break;
    34. else last=x;
    35. if(!i)res++;
    36. }
    37. return res;
    38. }
    39. signed main(){
    40. fast
    41. init();
    42. while(cin>>l>>r)cout<<dp(r)-dp(l-1)<
    43. return ~~(0^_^0);
    44. }

    1083. Windy数

    预处理出 //f[i][j]表示i位首位为j的合法方案cnt  最开始就抽风想成了 windy cnt 这不是直接就算出来了吗 额额

    状态转移方程就很好写了

    f[i][j]+=f[i-1][k] [abs(k-j)>=2]

    last注意初始化成一个很大的或者很小的都可以

    让后重要的来了 我们不能让第一个进来就有前导0

    为什么呢 因为我们算的时候会把前导0算进去 比如13 我们要是算进013那就寄了

    当然你可以特判着做当然会更麻烦 我们可以就让第一个不允许有前导0

    那么后面自然都可以有0了

    我们最后再处理有前导0的情况那么就很简单了 因为我们f[i][0]就是存的这种情况分别加进去就是

    然后值得主义的一点是 要是我们这边把0当作了一个windy数的话

    我们这边加的时候还要特判最后一位 也可以是0

    1. #include
    2. using namespace std;
    3. const int N = 15;
    4. const int M = 1<<12;
    5. const int mod = 1e9+7;
    6. #define int long long
    7. #define endl '\n'
    8. #define Endl '\n'
    9. #define _ 0
    10. #define inf 0x3f3f3f3f3f3f3f3f
    11. #define fast ios::sync_with_stdio(false);cin.tie(nullptr);
    12. int f[N][N],l,r;//f[i][j]表示i位首位为j的合法方案cnt 不是windy数数量
    13. void init(){
    14. for(int i=0;i<=9;i++)f[1][i]=1;
    15. for(int i=2;i<=N;i++){
    16. for(int j=0;j<=9;j++){
    17. for(int k=0;k<=9;k++){
    18. if(abs(j-k)>=2)f[i][j]+=f[i-1][k];
    19. }
    20. }
    21. }
    22. }
    23. int dp(int n){
    24. if(!n)return 1;
    25. vector<int>v;
    26. while(n)v.push_back(n%10),n/=10;
    27. int res=0,last=-inf;
    28. for(int i=v.size()-1;i>=0;i--){
    29. int x=v[i];
    30. for(int j=(i==v.size()-1);j
    31. if(abs(j-last)>=2)res+=f[i+1][j];
    32. }
    33. if(abs(x-last)>=2)last=x;
    34. else break;
    35. if(!i)res++;
    36. }
    37. for(int i=1;isize();i++){
    38. for(int j=(i==1?0:1);j<=9;j++){
    39. res+=f[i][j];
    40. }
    41. }
    42. return res;
    43. }
    44. signed main(){
    45. fast
    46. init();
    47. cin>>l>>r;
    48. cout<<dp(r)-dp(l-1)<
    49. return ~~(0^_^0);
    50. }

  • 相关阅读:
    LeetCode·583.两个字符串的删除操作·动态规划
    基于SSM的在线家教管理系统
    MyBatisPlus笔记
    深度学习(十四)——优化器
    多目标遗传算法NSGA-II改进策略
    Springboot毕设项目办公物资管理系统6304w(java+VUE+Mybatis+Maven+Mysql)
    c# wpf template ItemsPanel 简单试验
    jdk8更新到333了,你确定不更新你的Java吗
    入门小白拥有服务器的建议
    YoLo V3 SPP u模型的讲解与总结
  • 原文地址:https://blog.csdn.net/qq_23852555/article/details/126066710