热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

【题解】AHOI2009同类分布

好开心呀~果然只有不看题解做出来的题目才会真正的有一种骄傲与满足吧ヾ(๑╹◡╹)ノ实际上这题只要顺藤摸瓜就可以了。首先按照数位dp的套路,有两维想必是省不掉:1.当前dp到到的位数;2.

  好开心呀~果然只有不看题解做出来的题目才会真正的有一种骄傲与满足吧ヾ(๑╹◡╹)ノ"

  实际上这题只要顺藤摸瓜就可以了。首先按照数位dp的套路,有两维想必是省不掉:1.当前dp到到的位数;2.0/1状态表示是否受限制(这一条是因为有数字上限)。然后根据这两个维度来接着往下想。第二个维度先撇开不看,我们只考虑如何从第 \(i - 1\) 位dp到第 \(i\) 位。在这里其实卡了有点久,因为如果除数与被除数都在改变,那么两维的转移是非常凉凉的。

  这个时候联想题目的特殊性质 ----- 当感觉无法优化转移 / 转移方式的时候,考虑状态的重新设计 & 题目的特别要求。然后很开心的发现:\(1e18\) 实际上各位数字的和最大都只有 \(162\)。那么岂不是乱搞也可以?所以我们固定除数 \(Q\) 为 \(\left ( 1, 162 \right )\) 当中的任意一个数,分别进行dp即可。此时的转移就简单了,因为除数固定,自然地追加一维表示余数。状态固定为 \(f[i][j][k][L]\),表示dp到第 \(i\) 位,要求第 \(\left ( 1, i \right )\) 位的数字之和加起来为 \(j\),且原数除以 \(Q\) 的余数为 \(k\),限制为\(L\left ( 0, 1 \right )\)的总个数。

  感觉这份代码写的还行,跑得也还行……能看。

#include 
using namespace std;
#define int long long
int a[20], Res, mul[20];
int f[20][165][165][2];

int read()
{
    int x = 0, k = 1;
    char c;
    c = getchar();
    while(c <'0' || c > '9') { if(c == '-') k = -1; c = getchar(); }
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * k;
}

#define Pre f[now][tot][rm][lim]
int DP(int now, int tot, int rm, bool lim)
{
    if(~Pre) return Pre;
    else Pre = 0;
    if(tot > now * 9) return 0;
    if(now == 1) 
    {
        if(tot > 9 || (tot > a[now] && lim)) return Pre = 0;
        return Pre = ((Pre = tot % Res == rm) ? 1 : 0);
    }
    for(int i = 0; i <= 9; i ++)
    {
        if(i > a[now] && lim) break;
        if(tot break;
        int q = (mul[now] * i) % Res; q = q % Res; 
        int L = (i == a[now] && lim);
        f[now][tot][rm][lim] += DP(now - 1, tot - i, (rm - q + Res) % Res, L);
    }
    return f[now][tot][rm][lim];
}
#undef Pre

int Solve(int x)
{
    int k = x, ans = 0, num = 0; 
    while(k) num ++, a[num] = k % 10, k /= 10;
    for(int i = 1; i <= 163; i ++)
    {
        Res = i; if(i > num * 9) continue;
        memset(f, -1, sizeof(f));
        ans += DP(num, i, 0, 1);
    }
    return ans;
}

signed main()
{
    int a = read(), b = read(); mul[1] = 1;
    for(int i = 2; i <= 20; i ++) mul[i] = mul[i - 1] * 10;
    printf("%lld\n", Solve(b) - Solve(a - 1));
    return 0;
}

 


推荐阅读
  • 套路题?感觉讲不清,先写建图把每个点拆成两个,A和B,S-Ai流量1费用0,Bi-T流量1费用0ÿ ... [详细]
  • 做了题还是忍不住要写一发题解,感觉楼下的不易懂啊。本题解使用latex纯手写精心打造。题意:求\(\frac{1}{x}\frac{1}{y}\frac ... [详细]
  • 题意:多模式串匹配,输出模式串的ID思路:典型AC自动机.用向量存储答案ID#include<cstdio>#include<cstring> ... [详细]
  • 对象与对象之间的成员变量是相互独立的.要想共用数据,则需要使用静态成员或静态方法#只要在类中声明静态成员变量,即使不定义对象,也可以为静态成员变量分配空间,进而可以使用静态成员变 ... [详细]
  • 2017年二级计算机c真题语言,2017全国计算机二级C考试真题
    2017全国计算机二级C考试真题【2017全国计算机二级C考试真题】一、改错题函数fun的功能是:将一副扑克牌编号为1,2,353,54,以某种特 ... [详细]
  • rk3399的gpu测试节点在:catsysdevicesplatformff9a0000.gpudevfreqff9a0000.gpuload如果没有使用gpu的话 ... [详细]
  • 实现方式1、首先设置一个Qt下的一个窗口基类;2、窗口基类继承自一个重写的QGLWidget类和一个osgViewer::Viewer类3、重新QGLWidget类& ... [详细]
  • 实验六提交版
    1.21.3part2共用体与结构体类型的区别?答:共用体与结构体的区别在于它们的表示方法不同。结构体内,结构体的各成员顺序排列存储,每个成员都有自己独立的存储位置,而共用体的情况 ... [详细]
  • 1.中点画圆算法(1)P为当前点亮象素,那么,下一个点亮的象素可能是P1(Xp+1,Yp)或P2(Xp+1,Yp+1)。(2)构造函数:F(X,Y) ... [详细]
  • ACM Trick点常用操作记录(持续更新)(敏感时间)
    敏感时间注意多维vector速度很慢如下图string类的连接#includeusingnamespacestd;intmain(){stringt ... [详细]
  • 0-0#includeintmain(intargc,charconst*argv[]){std::cout ... [详细]
  • Bzoj1185最小矩阵覆盖[旋转卡壳+凸包+处理[0]情况]
    题目链接题目大意:就是给你若干个点用一个最小的矩形把这些点覆盖掉就是给你若干个点用一个最小的矩形把这些点覆盖掉就是给你若干个点用一个最小的矩形把这些点覆盖掉解题思路& ... [详细]
  • 哈密顿圈|回溯-6原文:https://www.geesfo ... [详细]
  • [转]Makefile 使用总结
    2019独角兽企业重金招聘Python工程师标准1.Makefile简介Makefile是和make命令一起配合使用的.很多大型项目的编译都是通过Makefile来组织的,如 ... [详细]
  • 数据结构学习记录(六)树
    前言树:有且只有一个称为根的节点,有若干个互不相交的子树,这些子树本身也是一颗树。通俗理解:树由节点与边组成, ... [详细]
author-avatar
局外人
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有