原题戳这里

绝对是一道好题

需要注意到两个东西

**

1.符合条件的数可以拆成一堆\(11...11\)相加的形式,比如\(1145=1111+11+11+11+1\)

2.\(1,11,111,1111,...\)模\(p\)会出现循环,循环节长度不超过\(p\)

**

还有就是\(11...11\)最多为\(9\)个,然后就可以\(dp\)了

首先需要统计长度为\(1-n\)的全\(1\)串中有多少个模\(p\)为\(r\),记为\(cnt[r]\),这个可以\(O(p)\)的预处理出来

接着设\(f[k][i][j]\)表示在\(cnt[1]-cnt[i]\)中已经选了\(k\)个,且它们的和模\(p\)为\(j\)的方案数,那么转移如下:

\[f[k+l][i][(j+l\times i)\% p]+=\binom{cnt[i]+l-1}{l}f[k][i-1][j],0\leqslant k+l\leqslant 8
\]

为什么是\(8\)而不是\(9\)呢?因为我们至少要选一个长度为\(n\)的全\(1\)串,为了方便,我们在\(dp\)之前就把它选出来,相当于占用了\(9\)个中的\(1\)个

还有就是那个组合是可重组合

细节看代码吧:

#include <bits/stdc++.h>

using namespace std;

#define $SHOW(x) cout << #x" = " << x << endl

#define ll long long
#define MOD 999911659
#define MAXN 500 ll n, P, all, cnt[MAXN + 5], inv[10]; // all代表长度为n的全1串模p的值
int st, len, pos[MAXN + 5], f[10][MAXN + 5][MAXN + 5]; void add(int &x, int y) {
x = (x + y) % MOD;
if(x < 0) x += MOD;
} int Mul(int x, int y) {
return 1LL * x * y % MOD;
} int fpow(int x, int p) {
int ret = 1;
while (p) {
if (p & 1) ret = Mul(ret, x);
x = Mul(x, x);
p >>= 1;
}
return ret;
} int C(ll n, ll m) {
if (n < m) return 0;
int ret = 1;
for (ll x = n, y = 1; y <= m; --x, ++y)
ret = Mul(ret, Mul(x % MOD, inv[y]));
return ret;
} int main() {
cin >> n >> P;
if (n <= P) { // 预处理循环节
int sum = 0;
for (int i = 1; i <= n; ++i) {
sum = (sum * 10 + 1) % P;
cnt[sum]++;
}
all = sum;
}
else {
int sum = 0;
for (int i = 1; i <= P + 1; ++i) {
sum = (sum * 10 + 1) % P;
if (cnt[sum]) {
st = pos[sum], len = i - pos[sum];
break;
}
pos[sum] = i, cnt[sum]++;
}
for (int i = 0; i < P; ++i)
if (cnt[i] && pos[i] >= st) {
cnt[i] = (n - st + 1) / len;
if (pos[i] - st + 1 <= (n - st + 1) % len) cnt[i]++;
if ((pos[i] - st + 1) % len == (n - st + 1) % len) all = i;
}
}
for (int i = 1; i <= 8; ++i) inv[i] = fpow(i, MOD - 2);
f[0][0][all] = 1; // 初始化
for (int r = 0; r < P; ++r) {
for (int i = 0; i <= 8; ++i)
for (int k = 0; k < P; ++k)
for (int j = 0; j + i <= 8; ++j)
if(cnt[r]) add(f[i + j][r + 1][(k + j * r % P) % P], Mul(C(cnt[r] + j - 1, j), f[i][r][k]));
else f[i][r + 1][k] = f[i][r][k]; // 注意这里,要把值继承过来
}
int ans = 0;
for (int i = 0; i <= 8; ++i)
add(ans, f[i][P][0]);
printf("%d\n", ans);
return 0;
}

[SDOI2010]代码拍卖会——DP的更多相关文章

  1. Luogu2481 SDOI2010 代码拍卖会 DP、组合

    传送门 神仙DP 注意到\(N \leq 10^{18}\),不能够直接数位DP,于是考虑形成的\(N\)位数的性质. 因为低位一定不会比高位小,所以所有满足条件的\(N\)位数一定是不超过\(9\) ...

  2. SDOI2010代码拍卖会 (计数类DP)

    P2481 SDOI2010代码拍卖会 $ solution: $ 这道题调了好久好久,久到都要放弃了.洛谷的第五个点是真的强,简简单单一个1,调了快4个小时! 这道题第一眼怎么都是数位DP,奈何数据 ...

  3. BZOJ 1974: [Sdoi2010]auction 代码拍卖会( dp )

    在1, 11, 111……中选<=8个, + 11..(n个1)拼出所有可能...这些数mod p至多有p中可能, 找出循环的处理一下. 那么dp就很显然了...dp(i, j, k)表示前i种 ...

  4. 【BZOJ-1974】auction代码拍卖会 DP + 排列组合

    1974: [Sdoi2010]auction 代码拍卖会 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 305  Solved: 122[Submit ...

  5. [SDOI2010]代码拍卖会

    题目描述 随着iPig在P++语言上的造诣日益提升,他形成了自己一套完整的代码库.猪王国想参加POI的童鞋们都争先恐后问iPig索要代码库.iPig不想把代码库给所有想要的小猪,只想给其中的一部分既关 ...

  6. bzoj 1974: [Sdoi2010]代码拍卖会

    Description 随着iPig在P++语言上的造诣日益提升,他形成了自己一套完整的代 码库.猪王国想参加POI的童鞋们都争先恐后问iPig索要代码库.iPi g不想把代码库给所有想要的小猪,只想 ...

  7. 洛谷 P2481 [SDOI2010]代码拍卖会

    洛谷 这大概是我真正意义上的第一道黑题吧! 自己想出了一个大概,状态转移方程打错了一点点,最后还是得看题解. 一句话题意:求出有多少个\(n\)位的数,满足各个位置上的数字从左到右不下降,且被\(p\ ...

  8. [BZOJ1974][SDOI2010]代码拍卖会[插板法]

    题意 询问有多少个数位为 \(n\) 的形如 \(11223333444589\) 的数位值不下降的数字在\(\mod p\) 的意义下同余 \(0\). $n\leq 10^{18} ,p\leq ...

  9. luogu P2481 [SDOI2010]代码拍卖会

    luogu 题目中的那个大数一定是若干个1+若干个2+若干个3...+若干个9组成的,显然可以转化成9个\(\underbrace {111...1}_{a_i个1}(0\le a_1\le a_2\ ...

随机推荐

  1. Spring 视图层如何显示验证消息提示

    1.示例 <p th:if="${#fields.hasErrors('name')}" th:errors="*{name}" ></p&g ...

  2. [转帖]SSH远程登录配置文件sshd_config详解

    SSH远程登录配置文件sshd_config详解 2016年06月02日 17:42:25 Field_Yang 阅读数 61386   版权声明:本文为博主原创文章,遵循CC 4.0 by-sa版权 ...

  3. Hanlp-地名识别调试方法详解

    HanLP收词特别是实体比较多,因此特别容易造成误识别.下边举几个地名误识别的例子,需要指出的是,后边的机构名识别也以地名识别为基础,因此,如果地名识别不准确,也会导致机构名识别不准确. 类型1 数字 ...

  4. Vue.js + Element.ui 从搭建环境到打包部署

    一.搭建环境 由于新的node已经集成了npm,所以直接安装node,前往node官网下载最新版本的node,根据自己的操作系统选择相应的包,按照步骤一步步走就可以,这里不做过多介绍. 安装好后可以打 ...

  5. (五)mybatis开发dao层

    目录 SqlSession 是线程不安全的 原始 dao 开发方法 Mapper 代理方法 关于代理对象 SqlSession 是线程不安全的 SqlSession 是 线程不安全 的: 对于它,我们 ...

  6. logstash grok

    input { file { path => "/opt/service/test-service/logs/catalina-error*.log" type => ...

  7. 3037 插板法+lucas

    先说下lucas定理 1)Lucas定理:p为素数,则有: (2)证明: n=(ak...a2,a1,a0)p = (ak...a2,a1)p*p + a0 =  [n/p]*p+a0 (注意 这里( ...

  8. 轻松搭建CAS 5.x系列(3)-连接数据库,使用数据库表中的帐号做登录

    概要说明 前面的CAS SEVER中的登录帐号名是配置中写死的,实际情况中不太可能用这些方法.通常情况下,数据库的帐号名密码都是在数据库表中的,这样可以对登录帐号进行增删改的处理. 如果您对搭建固定帐 ...

  9. Faster RCNN学习笔记

    感谢知乎大神的分享 https://zhuanlan.zhihu.com/p/31426458 Ross B. Girshick在2016年提出了新的Faster RCNN,在结构上,Faster R ...

  10. Git忽略已追踪文件或文件夹

    今天拉取代码,用vs生成后发现obj文件夹下自动生成的文件被同事提交了,这个本应该加入到ignore的 我就需要把这个文件夹加入到gitignore, 不过已经追踪的文件和文件夹,直接添加到gitig ...