洛谷 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 型机器人而言,将它 ...
随机推荐
- 关于 我的博客和Git-hub
欢迎大家到我的GitHub 热烈讨论 https://github.com/ljj-19951010 由于另一个博客忘了怎么登陆了,换用此博客(仅供个人学习使用,请勿传播) 如果想看 特别详细的教程请 ...
- SLAM名词介绍
gauge freedom:测量自由度 degrees-of-freedom(DoF) 自由度 wide-baseline matches:宽基线匹配 宽基线匹配:从描绘同一场景的两个或多个图像中建立 ...
- 痞子衡嵌入式:超级下载算法RT-UFL v1.0在Keil MDK下的使用
痞子衡主导的"学术"项目 <RT-UFL - 一个适用全平台i.MXRT的超级下载算法设计> v1.0 版发布近 4 个月了,部分客户已经在实际项目开发调试中用上了这个 ...
- CanalAdmin搭建Canal Server集群
CanalAdmin搭建Canal Server集群 一.背景 二.机器情况 三.实现步骤 1.下载canal admin 2.配置canalAdmin 3.初始化canal admin数据库 4.启 ...
- 一文读懂Android进程及TCP动态心跳保活
一直以来,APP进程保活都是 各软件提供商 和 个人开发者 头疼的问题.毕竟一切的商业模式都建立在用户对APP的使用上,因此保证APP进程的唤醒,提升用户的使用时间,便是软件提供商和个人开发者的永恒追 ...
- camera HSYNC:VSYNC
HSYNC:行锁存,换行信号VSYNC:祯锁存,换页信号 320×240的屏,每一行需要输入320个脉冲来依次移位.锁存进一行的数据,然后来个HSYNC 脉冲换一行:这样依次输入240行之后换行同时来 ...
- 整数转化 牛客网 程序员面试金典 C++ Python
整数转化 牛客网 程序员面试金典 C++ Python 题目描述 编写一个函数,确定需要改变几个位,才能将整数A转变成整数B. 给定两个整数int A,int B.请返回需要改变的数位个数. 测试样例 ...
- Tarjan算法离线 求 LCA(最近公共祖先)
本文是网络资料整理或部分转载或部分原创,参考文章如下: https://www.cnblogs.com/JVxie/p/4854719.html http://blog.csdn.net/ywcpig ...
- 2020美亚团体—Daniel篇
Daniel的桌上计算机的哈希值(SHA-256)是甚么? 通过取证大师计算 SHA-256值为 07DD40CF28603F421F3A09CD38F1C8AA40A2AC4BFB46ECF8299 ...
- 1.在项目中使用D3.js
在项目中使用D3.js D3.js(全称:Data-Driven Documents)是一个基于数据操作文档的JavaScript库.D3帮助您使用HTML.SVG和CSS使数据生动起来.D3对web ...