洛谷 P6276 - [USACO20OPEN]Exercise P(组合数学+DP)
废了,又不会做/ll
orz czx 写的什么神仙题解,根本看不懂(%%%%%%%%%
首先显然一个排列的贡献为其所有置换环的乘积。考虑如何算之。
碰到很多数的 LCM 之积只有两种可能,一是 Min-Max 容斥将 LCM 转化为 GCD,而是枚举质因子及其次数算贡献。但对于此题而言前者不是太可做(可能有复杂度不错(大概 \(n^2d(n)\)?)的解法,不过我没有细想所以也不太清楚),因此考虑后者。
考虑用类似于差分的思想,对于每个质因子 \(p\) 的每个次数 \(k\),我们考虑计算一下存在一个置换环大小为 \(p^k\) 的排列个数 \(c\),那么我们令答案乘以 \(p^c\),不难发现对于一个排列,如果其置换环大小中 \(p\) 次数最多的一项的次数为 \(k'\),那么它的贡献会在 \(k=1,2,3,\cdots,k'\) 处各被算一次,因此这样计算恰好可以算到所有置换环的贡献。于是限制问题转化为,如何求有多少个排列,满足其至少有一个置换环大小为 \(x=p^k\) 的倍数。
考虑从反面入手,即计算有多少个长度为 \(n\) 的排列不含任何长度为 \(x\) 的置换环,设其为 \(f_n\)。那么这个子问题的答案就是 \(n!-f_n\)。那么如何求 \(f_i\) 呢?我们再从反面入手,计算一下含有长度为 \(x\) 的置换环的排列个数(为什么要一来一回搞两次反面呢?因为直接做递推式中需要用到 \(f_i\)),那么我们枚举长度为 \(x\) 的置换环的长度之和 \(j\),设 \(g_j\) 表示有多少个长度为 \(j\) 的排列满足其每个置换环大小都是 \(x\) 的倍数,那么有 \(f_i=i!·\sum\limits_{j\in[1,i]}\dbinom{i}{j}g_jf_{i-j}\),组合数表示分配 \(j\) 个数给大小是 \(x\) 的倍数的置换环的方案数,一目了然。接下来考虑如何求 \(g_i\),我们枚举 \(1\) 所在的置换环的大小 \(j\),类似地有 \(\sum\limits_{j\in[1,i]}g_{i-j}\dbinom{i-1}{j-1}(j-1)!\)。注意到只有 \(x\mid i\) 时 \(g_i\ne 0\),因此有用的 \(f,g\) 各只有 \(\dfrac{n}{x}\) 个,也就保证了复杂度。这部分可能可以容斥,不过没有细想(
总复杂度 \(\sum\limits_{i=1}^n\dfrac{n^2}{i^2}=\mathcal O(n^2)\),因为 \(\sum\limits_{n>0}\dfrac{1}{n^2}<2\)。
const int MAXN=7500;
int n,mod,phi,c[MAXN+5][MAXN+5],fac[MAXN+5];
int qpow(int x,int e){
int ret=1;
for(;e;e>>=1,x=1ll*x*x%mod) if(e&1) ret=1ll*ret*x%mod;
return ret;
}
int pr[MAXN/5+5],prcnt=0,vis[MAXN+5];
void sieve(int n){
for(int i=2;i<=n;i++){
if(!vis[i]) pr[++prcnt]=i;
for(int j=1;j<=prcnt&&pr[j]*i<=n;j++){
vis[pr[j]*i]=1;if(i%pr[j]==0) break;
}
}
}
int f[MAXN+5],g[MAXN+5];
int calc(int x){
memset(f,0,sizeof(f));memset(g,0,sizeof(g));g[0]=1;
for(int i=x;i<=n;i+=x) for(int j=x;j<=i;j+=x) g[i]=(g[i]+1ll*c[i-1][j-1]*g[i-j]%phi*fac[j-1])%phi;
// for(int i=x;i<=n;i+=x) printf("%d %d\n",i,g[i]);
for(int i=n%x;i<=n;i+=x){
f[i]=fac[i];
for(int j=x;j<=i;j+=x) f[i]=(f[i]-1ll*f[i-j]*g[j]%phi*c[i][j]%phi+phi)%phi;
} return (fac[n]-f[n]+phi)%phi;
}
int main(){
scanf("%d%d",&n,&mod);phi=mod-1;sieve(n);
for(int i=(fac[0]=1);i<=n;i++) fac[i]=1ll*fac[i-1]*i%phi;
for(int i=0;i<=n;i++){
c[i][0]=1;
for(int j=1;j<=i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%phi;
} int res=1;
// for(int i=1;i<=n;i++) printf("%d\n",calc(i));
for(int i=1;i<=prcnt;i++) for(int j=pr[i];j<=n;j*=pr[i])
res=1ll*res*qpow(pr[i],calc(j))%mod;
printf("%d\n",res);
return 0;
}
洛谷 P6276 - [USACO20OPEN]Exercise P(组合数学+DP)的更多相关文章
- 洛谷P3158 [CQOI2011]放棋子 组合数学+DP
题意:在一个m行n列的棋盘里放一些彩色的棋子,使得每个格子最多放一个棋子,且不同颜色的棋子不能在同一行或者同一列.有多少祌方法? 解法:这道题不会做,太菜了qwq.题解是看洛谷大佬的. 设C是组合数, ...
- 洛谷P2606 [ZJOI2010]排列计数 组合数学+DP
题意:称一个1,2,...,N的排列P1,P2...,Pn是Magic的,当且仅当2<=i<=N时,Pi>Pi/2. 计算1,2,...N的排列中有多少是Magic的,答案可能很大, ...
- 洛谷P2507 [SCOI2008]配对 题解(dp+贪心)
洛谷P2507 [SCOI2008]配对 题解(dp+贪心) 标签:题解 阅读体验:https://zybuluo.com/Junlier/note/1299251 链接题目地址:洛谷P2507 [S ...
- 洛谷 P3177 [HAOI2015]树上染色 树形DP
洛谷 P3177 [HAOI2015]树上染色 树形DP 题目描述 有一棵点数为 \(n\) 的树,树边有边权.给你一个在 \(0 \sim n\)之内的正整数 \(k\) ,你要在这棵树中选择 \( ...
- 洛谷 P4072 [SDOI2016]征途 斜率优化DP
洛谷 P4072 [SDOI2016]征途 斜率优化DP 题目描述 \(Pine\) 开始了从 \(S\) 地到 \(T\) 地的征途. 从\(S\)地到\(T\)地的路可以划分成 \(n\) 段,相 ...
- 洛谷P1880 石子合并(区间DP)(环形DP)
To 洛谷.1880 石子合并 题目描述 在一个园形操场的四周摆放N堆石子,现要将石子有次序地合并成一堆.规定每次只能选相邻的2堆合并成新的一堆,并将新的一堆的石子数,记为该次合并的得分. 试设计出1 ...
- 洛谷P1063 能量项链(区间DP)(环形DP)
To 洛谷.1063 能量项链 题目描述 在Mars星球上,每个Mars人都随身佩带着一串能量项链.在项链上有N颗能量珠.能量珠是一颗有头标记与尾标记的珠子,这些标记对应着某个正整数.并且,对于相邻的 ...
- 洛谷P1282 多米诺骨牌 (DP)
洛谷P1282 多米诺骨牌 题目描述 多米诺骨牌有上下2个方块组成,每个方块中有1~6个点.现有排成行的 上方块中点数之和记为S1,下方块中点数之和记为S2,它们的差为|S1-S2|.例如在图8-1中 ...
- 洛谷 P5469 - [NOI2019] 机器人(区间 dp+拉格朗日插值)
洛谷题面传送门 神仙题,放在 D1T2 可能略难了一点( 首先显然对于 P 型机器人而言,将它放在 \(i\) 之后它会走到左边第一个严格 \(>a_i\) 的位置,对于 Q 型机器人而言,将它 ...
随机推荐
- 返回值优化 RVO
<深度探索C++对象模型>-- 2.3 返回值的初始化 & 在编译器层面做优化
- Excel一对多查找
很多人在Excel中用函数公式做查询的时候,都必然会遇到的一个大问题,那就是一对多的查找/查询公式应该怎么写?大多数人都是从VLOOKUP.INDEX+MATCH中入门的,纵然你把全部的多条件查找方法 ...
- Spark RDD编程(博客索引,日常更新)
本篇主要是记录自己在中解决RDD编程性能问题中查阅的论文博客,为我认为写的不错的建立索引方便查阅,我的总结会另立他篇 1)通过分区(Partitioning)提高spark性能https://blog ...
- Java版人脸检测详解下篇:编码
欢迎访问我的GitHub https://github.com/zq2599/blog_demos 内容:所有原创文章分类汇总及配套源码,涉及Java.Docker.Kubernetes.DevOPS ...
- Vue el 使用el-checkbox-group复选框进行单选框操作
el-checkbox-group这个组件与其他复选框不一样,我当初也是半天不知道怎么操作 页面使用v-model绑定 size就是等比例缩小放大,v-ror循环应该看的懂.重要的是@chage到我们 ...
- docker multi-stage 多阶段构建
多阶段构建 一.需求 二.普通构建 1.编写Dockerfile 2.构建镜像 三.多阶段(multi-stage)构建 1.编写Dockerfile 2.构建镜像 四.比较2个镜像的体积大小 我们在 ...
- 2021.10.26考试总结[冲刺NOIP模拟16]
T1 树上的数 \(DFS\)一遍.结构体存边好像更快? \(code:\) T1 #include<bits/stdc++.h> using namespace std; namespa ...
- linux下使用shell命令通过wpa_cli控制wpa_supplicant连接wifi
最近在调试wifi,已经把wpa_supplicant 工具编译打包好了,为了测试wif驱动及wifi模块是否ok,需要用shell命令临时启动wifi服务连接wifi热点测试. 首先板子启动用ifc ...
- jQuery淡入淡出效果
如果是通过鼠标点击事件来触发动画效果可以使用 $("#button").click(function(){ $("#div").stop().fadeToggl ...
- 【Python+postman接口自动化测试】(3)什么是接口测试?
什么是接口测试? 接口测试是测试系统组件间接口的一种测试.接口测试主要用于检测外部系统与系统之间以及内部各个子系统之间的交互点.测试的重点是要检查数据的交换.传递和控制管理过程,以及系统间的相互逻辑依 ...