题面传送门

题意:

给出 \(n,p\),求有多少 \(n\) 位数 \(X=a_1a_2a_3\dots a_n\) 满足:

  1. 该 \(n\) 位数不含前导零
  2. \(a_i \leq a_{i+1}\)
  3. \(X\) 为 \(p\) 的倍数。

答案对 \(998244353\) 取模。

\(1 \leq n \leq 10^{18}\),\(1 \leq p \leq 500\)。

CSP 之前做的了,隔了好久才把题解补了。。。

本题的突破口在于怎样处理 \(a_i \leq a_{i+1}\) 这个条件。

我们不妨进行一个转化:每次加一个全是 \(1\) 的后缀,最多可以加 \(9\) 次,这样就能保证得到的数一定满足条件 2。

又由于 \(X\) 不能含前导零,所以我们加的后缀中必须有一个是长度 \(n\) 的后缀。

于是我们有了优秀的 \(n^8\) 的做法,可以拿到 10 分的好成绩。

但仔细观察就可以发现,每个后缀 \(111...11\) 对 \(p\) 取模的余数呈周期分部。

例如当 \(p=12\) 的时候,各后缀模 \(p\) 的余数分别为:\(1,11,9,7,11,9,7,11,...\)

不难看出除了 \(1\) 以外每三个一循环。

有了这个发现,我们就可以将这 \(n\) 个后缀分为三类:

  1. 还没进循环节
  2. 在完整的循环节中
  3. 在最后多出的部分中

证明可用扩展欧拉定理。

记 \(cnt_m\) 为模 \(p\) 余 \(m\) 的后缀数。分类讨论可以在 \(\mathcal O(p)\) 的时间内求出 \(cnt\)。

接下来就可以 \(dp\) 了。\(dp[i][j][k]\) 表示选择了余数为 \(0\) 到 \(i\) 的后缀共 \(j\),它们的和模 \(p\) 余 \(k\) 的方案数。

转移的时候枚举选择多少个余数为 \(i\) 的后缀,记为 \(s\)。

由于可以重复选择,这一部分的方案数可以用隔板法求,\(\dbinom{cnt_i+s-1}{s}\)。

最后别忘了我们必须选择长度为 \(n\) 的后缀,所以答案为 \(\sum\limits_{i=0}^8dp[p-1][i][(p-pn)\mod m]\),其中 \(pn\) 为长度为 \(n\) 的后缀对 \(p\) 取模的结果。

别忘了特判 \(p=1\),卡了我很久。。。。。。

综上,这是一道非常不错的综合题。

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define fz(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define ffe(it,v) for(__typeof(v.begin()) it=v.begin();it!=v.end();it++)
#define fill0(a) memset(a,0,sizeof(a))
#define fill1(a) memset(a,-1,sizeof(a))
#define fillbig(a) memset(a,63,sizeof(a))
#define pb push_back
#define mp make_pair
#define int long long
typedef pair<int,int> pii;
typedef long long ll;
const ll MOD=999911659;
ll n;
int p,a[505];
int vis[505];
ll cnt[505];
ll dp[505][10][505];
inline ll qpow(ll x,ll e){
ll ans=1;
while(e){
if(e&1) ans=ans*x%MOD;
x=x*x%MOD;e>>=1;
}
return ans;
}
ll inv[505];
inline ll getc(ll x,ll y){
ll ans=1;
for(int i=x;i>=x-y+1;i--) ans=ans*(i%MOD)%MOD;
for(int i=1;i<=y;i++) ans=ans*inv[i]%MOD;
return ans;
}
inline void inc(ll &x,ll y){
x+=y;if(x>=MOD) x-=MOD;
}
signed main(){
scanf("%lld%lld",&n,&p);
for(int i=1;i<=10;i++) inv[i]=qpow(i,MOD-2);
int cur=1;vis[cur]=1%p;
int len1,len2;a[1]=1%p;
for(int i=2;i<=1000;i++){
cur=(cur*10+1)%p;a[i]=cur;
if(vis[cur]){
len1=vis[cur]-1;len2=i-vis[cur];
break;
}
vis[cur]=i;
}
if(n<=len1) for(int i=1;i<=n;i++) cnt[a[i]]++;
else{
ll cyc=(n-len1)/len2;
ll rem=(n-len1)%len2;
for(int i=1;i<=len1;i++) cnt[a[i]]++;
for(int i=len1+1;i<=len1+len2;i++) cnt[a[i]]+=cyc;
for(int i=len1+1;i<=len1+rem;i++) cnt[a[i]]++;
}
for(int i=0;i<9;i++){
dp[0][i][0]=getc(cnt[0]+i-1,i)%MOD;
}
for(int i=0;i<p-1;i++) for(int j=0;j<9;j++) for(int k=0;k<p;k++){
for(int l=0;l+j<9;l++){
inc(dp[i+1][l+j][(k+l*(i+1))%p],dp[i][j][k]*getc(cnt[i+1]+l-1,l)%MOD);
}
}
int gn;
if(n<=len1) gn=a[n];
else if((n-len1)%len2==0) gn=a[len1+len2];
else gn=a[len1+(n-len1)%len2];
ll ans=0;
for(int i=0;i<9;i++) inc(ans,dp[p-1][i][(p-gn)%p]);
printf("%lld\n",ans);
return 0;
}
/*
1000000000 499
1000000000000000000 1
2 1
*/

洛谷 P2481 [SDOI2010]代码拍卖会(背包+隔板法)的更多相关文章

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

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

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

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

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

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

  4. 洛谷 P2014 选课(树形背包)

    洛谷 P2014 选课(树形背包) 思路 题面:洛谷 P2014 如题这种有依赖性的任务可以用一棵树表示,因为一个儿子要访问到就必须先访问到父亲.然后,本来本题所有树是森林(没有共同祖先),但是题中的 ...

  5. 洛谷P4525 【模板】自适应辛普森法1与2

    洛谷P4525 [模板]自适应辛普森法1 与P4526[模板]自适应辛普森法2 P4525洛谷传送门 P4525题目描述 计算积分 结果保留至小数点后6位. 数据保证计算过程中分母不为0且积分能够收敛 ...

  6. BZOJ1925或洛谷2467 [SDOI2010]地精部落

    BZOJ原题链接 洛谷原题链接 先讲下关于波动数列的\(3\)个性质. 性质\(1\):对于数列中的每一对\(i\)和\(i + 1\),若它们不相邻,那么交换这两个数形成的依旧是一个波动数列. 性质 ...

  7. 洛谷 P4389: 付公主的背包

    题目传送门:洛谷 P4389. 题意简述: 有 \(n\) 个物品,每个物品都有无限多,第 \(i\) 个物品的体积为 \(v_i\)(\(v_i\le m\)). 问用这些物品恰好装满容量为 \(i ...

  8. 洛谷 P2467 [SDOI2010]地精部落

    洛谷 我讲的应该没有这个[https://www.luogu.org/blog/user55639/solution-p2467]清楚. 贴个代码算了: #include <bits/stdc+ ...

  9. 洛谷P2468 [SDOI2010]粟粟的书架

    来了来了,随便拽一道题写题解[大雾] 最近发现自己基础奇差于是开始复习之前学过的东西,正好主席树我几乎完全没学会,然后打开洛谷试炼场… 发现了这么一道二合一的题. 这道题其实分成两个部分,前50%是一 ...

随机推荐

  1. Hadoop面试题(四)——YARN

    1.简述hadoop1与hadoop2 的架构异同 1)加入了yarn解决了资源调度的问题. 2)加入了对zookeeper的支持实现比较可靠的高可用. 2.为什么会产生 yarn,它解决了什么问题, ...

  2. 如果你还不知道Apache Zookeeper?你凭什么拿大厂Offer!!

    很多同学或多或少都用到了Zookeeper,并知道它能实现两个功能 配置中心,实现表分片规则的统一配置管理 注册中心,实现sharding-proxy节点的服务地址注册 那么Zookeeper到底是什 ...

  3. 写了10000条Airtest截图脚本总结出来的截图经验,赶紧收藏!

    前言 今天想先给大家分享1个小白用户的Airtest从入门到放弃的故事: 小A是一个自动化的小白,在逛测试论坛的时候,偶然间发现了Airtest这个基于图像识别的UI自动化框架. 出于好奇,小A试用了 ...

  4. oo第一单元学习总结

    写在开头: 第一次接触面向对象思想和java语言,在学习以及完成作业的过程经历了一个比较痛苦的过程, 虽然在每次写作业时总是会有一些小小的抱怨,虽然写出的代码还是很差, 但是看到自己有所进步,还是感觉 ...

  5. js--数组的 fill() 填充方法详解

    前言 我们知道了很多了初始化数组的方法,但是初始化数组之后,数组中的每一项元素默认为 empty 空位占位,如何对数组这些空位添加默认的元素,ES6提供了 fill() 方法实现这一操作.本文总结数组 ...

  6. 鸿蒙轻内核M核的故障管家:Fault异常处理

    摘要:本文先简单介绍下Fault异常类型,向量表及其代码,异常处理C语言程序,然后详细分析下异常处理汇编函数实现代码. 本文分享自华为云社区<鸿蒙轻内核M核源码分析系列十八 Fault异常处理& ...

  7. 算法:N-皇后问题

    一.八皇后问题 八皇后问题是一个以国际象棋为背景的问题:如何能够在8 × 8 的国际象棋棋盘上放置八个皇后(Queen),使得任何一个皇后都无法直接吃掉其他的皇后.为了达到此目的,任两个皇后都不能处于 ...

  8. TensorFlow从入门到入坑(1)

    初识TensorFlow 一.术语潜知 深度学习:深度学习(deep learning)是机器学习的分支,是一种试图使用包含复杂结构或由多重非线性变换构成的多个处理层对数据进行高层抽象的算法. 深度学 ...

  9. Lambda-让人又爱又恨的“->"

    写在前边 聊到Java8新特性,我们第一反应想到的肯定是Lambda表达式和函数式接口的出现.要说ta到底有没有在一定程度上"优化"了代码的简洁性呢?抑或是ta在一定程度上给程序员 ...

  10. Java基础语法5-运算符

    运算符 基本运算符(算术.赋值.关系.逻辑等)不再赘述 位运算符 位运算符&.|.^.~.<<.>> &与 |或 ~非 ^异或 <<左移 >& ...