题目描述

随着iPig在P++语言上的造诣日益提升,他形成了自己一套完整的代码库。猪王国想参加POI的童鞋们都争先恐后问iPig索要代码库。iPig不想把代码库给所有想要的小猪,只想给其中的一部分既关系好又肯出钱的小猪,于是他决定举行了一个超大型拍卖会。

在拍卖会上,所有的N头小猪将会按照和iPig的好感度从低到高,从左到右地在iPig面前站成一排。每个小猪身上都有9猪币(与人民币汇率不明),从最左边开始,每个小猪依次举起一块牌子,上面写上想付出的买代码库的猪币数量(1到9之间的一个整数)。大家都知道,如果自己付的钱比左边的猪少,肯定得不到梦寐以求的代码库,因此从第二只起,每只猪出的钱都大于等于左边猪出的价钱。最终出的钱最多的小猪(们)会得到iPig的代码库真传,向着保送PKU(Pig Kingdom University)的梦想前进。

iPig对自己想到的这个点子感到十分满意,在去现场的路上,iPig就在想象拍卖会上会出现的场景,例如一共会出现多少种出价情况之类的问题,但这些问题都太简单了,iPig早已不敢兴趣了,他想要去研究更加困难的问题。iPig发现如果他从台上往下看,所有小猪举的牌子从左到右将会正好构成一个N位的整数,他现在想要挑战的问题是所有可能构成的整数中能正好被P整除的有多少个。由于答案过大,他只想要知道答案mod 999911659就行了。

输入输出格式

输入格式:

输入文件auction.in有且仅有一行:两个数N(1≤N≤10^18)、P(1≤P≤500),用一个空格分开。

输出格式:

输入文件auction.out有且仅有一行:一个数,表示答案除以999911659的余数。

输入输出样例

输入样例#1: 复制

2 3
输出样例#1: 复制

15

说明

样例解释

方案可以是:12 15 18 24 27 33 36 39 45 48 57 66 69 78 99,共15种。

数据规模

题解

  这题太神仙了……题解看都看不懂……最后基本只能硬生生的理解了……

  首先,我们考虑数列,原数列是一个不降的序列

  考虑如下数列,$1,1,2,3,4$

  这样我们是相当于竖着分割的

  那么怎么转换为横着分割呢?我们可以记录大于等于$1$的数的个数,为$5$,大于等于$2$的数的个数,为$3$……

  那么最后原数列可以转化为$11111,111,11,1$(每个数用相同个数的$1$表示),然后我们惊奇的发现他们的和和原来的$n$位数是一样的,也就是说,他们构成的整数取模之后也是一样的!

  这就相当于在每一个位置放相当于权值大小的石头,我们一开始是竖着分,而第二种方法是横着分割

  如果竖着分割,总共有$n$个数,如果横着割,把所有模$p$同余的看成一类,那么总共只有$p$个数(因为0,1,11,111这样下去模$p$的值肯定能构成一个循环节,所以只需要计算$p$次,剩下的可以直接计算)

  那么我们为什么不转化为横着分割呢?设$cnt[i]$表示模$p$为$i$的数的个数,那么题目就变成从$cnt[i]$中取$9$个使下标之和被$p$整除

  那么就可以转化为dp了,设$f[i][j][k]$表示考虑到第$i$个数,选了$k$个,他们的和模$p$为$j$,那么状态转移方程就是$f[i+1][(j+l*i)%mod][k+l]=(f[i][j][k]*C_{cnt_i}^l+f[i+1][(j+l*i)%mod][k+l])%mod$

  然后要先考虑加上一个$11111$($n$个$1$),因为我们dp的时候是允许有前导$0$的,所以得先强制至少为$1$才行

  细节有点多,都写在注解里了

 //minamoto
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int mod=;
ll ans,n,p,cnt[],beg,len,pos[],A[],c[][],f[][][],a;
int main(){
scanf("%lld%d",&n,&p);
ll sum=;
if(n<=p){
//直接把循环节跑出来
for(int i=;i<=n;++i) sum=(sum*+)%p,++cnt[sum];
a=sum;
}else{
//否则去找循环节,数出每一个余数的出现次数
for(int i=;i<=p+;++i){
sum=(sum*+)%p;
if(cnt[sum]){
beg=pos[sum],len=i-pos[sum];
break;
}
++cnt[sum],pos[sum]=i;
}
for(int i=;i<p;++i)
if(cnt[i]&&pos[i]>=beg){
cnt[i]=(n-beg+)/len;
if(pos[i]-beg+<=(n-beg+)%len) ++cnt[i];
if((pos[i]-beg+)%len==(n-beg+)%len) a=i;
}
}
A[]=,A[]=;
for(int i=;i<=;++i)
A[i]=(mod-mod/i)*A[mod%i]%mod;
for(int i=;i<p;++i){
c[i][]=;
if(cnt[i])
for(int j=;j<=;++j){
c[i][j]=cnt[i]*c[i][j-]%mod*A[j]%mod;
cnt[i]=(cnt[i]+)%mod;
//C(n,m-1)->C(n,m)
//C(n,m-1)/m*(n-m+1)=C(n,m)
//cnt[i]=n-m+1
}
}
f[][a][]=;
//默认加n个1,因为dp的时候可以有前缀0
for(int i=;i<p;++i)
for(int j=;j<p;++j)
for(int k=;k<;++k)
for(int l=;l<=k;++l)
(f[i+][j][k]+=f[i][(j-(l*i%p)+p)%p][k-l]*c[i][l]%mod)%mod;
for(int i=;i<=;++i) (ans+=f[p][][i])%=mod;
printf("%lld\n",ans);
return ;
}

BZOJ 1974 [Sdoi2010] auction 代码拍卖会(数位dp)的更多相关文章

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

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

  2. BZOJ 1974 [Sdoi2010]auction 代码拍卖会 ——动态规划

    把每一位上不递减的数拆成1+11+11111+11111+..... 然后就可以把巨大的N从复杂度中消掉,因为随着长度的增加1...111%p出现了循环节. 然后就是从n个数中选出几个使他们结果为0( ...

  3. Bzoj1974 [Sdoi2010]auction 代码拍卖会

    Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 375  Solved: 151 Description 随着iPig在P++语言上的造诣日益提升,他形成 ...

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

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

  5. bzoj 1026 [SCOI2009]windy数(数位DP)

    1026: [SCOI2009]windy数 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 4550  Solved: 2039[Submit][Sta ...

  6. bzoj 1026: [SCOI2009]windy数【数位dp】

    忘记limit不能记WA了一发-- 典型数位dp,变成work(r)-work(l-1),然后dfs的时候记录w当前位置,la上一个数选的什么,lm当前位是否有上限,ok当前位是否可以不考虑差大于等于 ...

  7. bzoj 3209 花神的数论题 —— 数位DP

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3209 算是挺简单的数位DP吧,但还是花了好久才弄明白... 又参考了博客:https://b ...

  8. bzoj 1902: Zju2116 Christopher lucas定理 && 数位DP

    1902: Zju2116 Christopher Time Limit: 1 Sec  Memory Limit: 64 MBSubmit: 172  Solved: 67[Submit][Stat ...

  9. BZOJ 3209 花神的数论题 数位DP+数论

    题目大意:令Sum(i)为i在二进制下1的个数 求∏(1<=i<=n)Sum(i) 一道非常easy的数位DP 首先我们打表打出组合数 然后利用数位DP统计出二进制下1的个数为x的数的数量 ...

随机推荐

  1. Thread之一:线程生命周期及五种状态

    <Thread之一:线程生命周期及五种状态> <juc线程池原理(四): 线程池状态介绍> 一.线程的生命周期及五种基本状态 关于Java中线程的生命周期,首先看一下下面这张较 ...

  2. 分治思想的应用:C++实现快速排序和随机化的快速排序

    分治思想的应用:C++实现快速排序和随机化的快速排序 原创 2014年09月08日 14:04:49 标签: 快速排序 / 随机化快速排序 / 排序算法 / 数据结构 947 1. 快速排序时冒泡排序 ...

  3. spring学习六

    1: @Valid 注解    @NotNull(message="名字不能为空") private String userName; @Max(value=120,message ...

  4. 在IIS下配置自定义的报错页面

    这里介绍在IIS中配置自定义出错页面的方法,主要以404为例,其他状态可类推 1.远程桌面连接IIS所在的服务器,进入控制面板>系统和安全>管理工具,双击打开IIS管理器,选择需要配置的网 ...

  5. JVM Class Loading过程

    转自:<Java Performance>第三章 VM Class Loading The Hotspot VM supports class loading as defined by ...

  6. java 截取替换掉括号 包括括号中的内容

    public static void main(String[] args) { String company = "华厦世纪(厦门)地产"; // System.out.prin ...

  7. python之Dict和set类型

    Dict就是一种key:value的表格: >>> d = { 'Adam':95, 'Lisa':85, 'Bart':59, 'Paul':75 } >>> p ...

  8. The Apache Tomcat installation at this directory is version 8.5.24 Tomcat 8.0 installation is expect

    在一台新电脑上搭建Java开发环境,JDK 是1.8,Tomcat下载了Tomcat 8.5.24,已经配置好了Java和Tomcat的环境变量,开发工具是Eclipse MARS,准备在Eclips ...

  9. 项目一:第九天 1、前台客户登录 2、Jquery citypicker省市区三级联动插件 4、业务受理(在线下单)

    1. 前台客户登录 2. Jquery citypicker省市区三级联动插件 3. 百度地图介绍 4. 业务受理(在线下单) 1 实现前台系统登录功能 1.1 Md5加密 admin(明文)---- ...

  10. C/C++读写csv文件

    博客转载自:http://blog.csdn.net/u012234115/article/details/64465398 C++ 读写CSV文件,注意一下格式即可 #include <ios ...