洛谷 P4071 [SDOI2016]排列计数 题解
P4071 [SDOI2016]排列计数
题目描述
求有多少种长度为 n 的序列 A,满足以下条件:
1 ~ n 这 n 个数在序列中各出现了一次
若第 i 个数 A[i] 的值为 i,则称 i 是稳定的。序列恰好有 m 个数是稳定的
满足条件的序列可能很多,序列数对 \(10^9+7\) 取模。
输入格式
第一行一个数 T,表示有 T 组数据。
接下来 T 行,每行两个整数 n、m。
输出格式
输出 T 行,每行一个数,表示求出的序列数
输入输出样例
输入 #1
5
1 0
1 1
5 2
100 50
10000 5000
输出 #1
0
1
20
578028887
60695423
说明/提示
测试点 1 ~ 3: $ T = 1000\(,\)n \leq 8\(,\)m \leq 8$;
测试点 4 ~ 6: $ T = 1000\(,\)n \leq 12\(,\)m \leq 12$;
测试点 7 ~ 9: $ T = 1000\(,\)n \leq 100\(,\)m \leq 100$;
测试点 10 ~ 12:$ T = 1000\(,\)n \leq 1000\(,\)m \leq 1000$;
测试点 13 ~ 14:$ T = 500000\(,\)n \leq 1000\(,\)m \leq 1000$;
测试点 15 ~ 20:$ T = 500000\(,\)n \leq 1000000\(,\)m \leq 1000000$。
【思路】
组合数学 / 错排
【题目大意】
a[i]在i的位置上面是稳定的
求有m个数是稳定的序列有多少个
【题目分析】
求符合条件序列的个数
符合条件序列可以先只看那稳定的数
一个序列中只有m个数是稳定的
其他的都是稳定的
那么稳定的数组合方式就是n个数里面取m个
因为如果稳定那一个数只对应一个位置
所以不存在顺序这一说,所以就是 \(C_n^m\)
知道了稳定数的组合方式
拿在看看除了这些稳定数之外数的组合方式
其他的数都是不在自己的位置上面的
也就是错排
直接用错排求出n-m(这里是减号下同下下同)个人错排方法的数量就好了
因为一种稳定数对应n-m个人的错排方式
所以数量数就是 \(C_n^m\) * (n-m)个人的错排方式
【存在的问题】
1.因为n和m的数据范围都是小于等于1e6
不是很小,而且T很大,
所以每次T 不能都单独求C和错排次数了
这样一定会超时
2.因为这道题中有取模运算
而递推求组合数只能过2000所以用这个阶乘求组合数就是必然的了
那么就要用到除法
而取模运算中不能用除法
【优化】
1.针对问题1
多次求解会超时
那就先预处理出来所有的阶乘和i个人的错排方式
到时候O(1)查询就好了
2.针对问题2
既然不能用除法那就用乘法
用逆元来代替除法就可以啦
求逆元因为模数是质数
所以可以用费马小定理求
费马小定理求逆元详见
这里
【完整代码】
#include<iostream>
#include<cstdio>
#define int long long
using namespace std;
const int k = 1e9 + 7;
const int Max = 1000005;
int d[Max];
int f[Max];
int inv[Max];
int p(int a,int b)
{
int ans = 1;
while(b)
{
if(b & 1)
ans = ans * a % k;
b >>= 1;
a = a * a % k;
}
return ans;
}
signed main()
{
d[0] = 1,d[2] = 1;
for(register int i = 3;i <= 1000000;++ i)
d[i] = (i - 1) * (d[i - 1] + d[i - 2]) % k;
f[0] = 1,f[1] = 1;
for(register int i = 2;i <= 1000000;++ i)
f[i] = f[i - 1] * i % k;
inv[0] = p(f[0],k - 2);
for(register int i = 1;i <= 1000000;++ i)
inv[i] = p(f[i],k - 2) % k;
int t;
cin >> t;
while(t --)
{
int n,m;
scanf("%lld%lld",&n,&m);
//c(n,m) * d(n-m)
//n!/m!(n - m)! * d[n-m]
printf("%lld\n",(f[n] * inv[m] % k * inv[n-m] % k * d[n-m]) % k);
}
return 0;
}
洛谷 P4071 [SDOI2016]排列计数 题解的更多相关文章
- 洛谷——P4071 [SDOI2016]排列计数(错排+组合数学)
P4071 [SDOI2016]排列计数 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列 ...
- 洛谷P4071 [SDOI2016] 排列计数 [组合数学]
题目传送门 排列计数 题目描述 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m ...
- 洛谷 P4071 [SDOI2016]排列计数
洛谷 这是一道组合数学题. 对于一个长为n的序列,首先我们要选m个使之稳定\(C^{m}_{n}\). 且要保证剩下的序列不稳定,即错排\(D_{n-m}\). 所以答案就是:\[ANS=C^{m}_ ...
- P4071 [SDOI2016]排列计数 题解
分析: 线性求逆元:https://blog.csdn.net/qq_34564984/article/details/52292502 代码: #include<cstdio> usin ...
- 洛谷 P2606 [ZJOI2010]排列计数 解题报告
P2606 [ZJOI2010]排列计数 题目描述 称一个\(1,2,...,N\)的排列\(P_1,P_2...,P_n\)是\(Magic\)的,当且仅当对所以的\(2<=i<=N\) ...
- BZOJ4517 & 洛谷4071:[SDOI2016]排列计数——题解
https://www.lydsy.com/JudgeOnline/problem.php?id=4517 https://www.luogu.org/problemnew/show/P4071 求有 ...
- ●洛谷P2606 [ZJOI2010]排列计数
题链: https://www.luogu.org/problemnew/show/P2606题解: 组合数(DP),Lucas定理 首先应该容易看出,这个排列其实是一个小顶堆. 然后我们可以考虑dp ...
- 洛谷P2606 [ZJOI2010]排列计数
题目描述 称一个1,2,...,N的排列P1,P2...,Pn是Magic的,当且仅当2<=i<=N时,Pi>Pi/2. 计算1,2,...N的排列中有多少是Magic的,答案可能很 ...
- 洛谷P2606 [ZJOI2010]排列计数(数位dp)
题目描述 称一个1,2,...,N的排列P1,P2...,Pn是Magic的,当且仅当2<=i<=N时,Pi>Pi/2. 计算1,2,...N的排列中有多少是Magic的,答案可能很 ...
随机推荐
- mongoDB对没有字段的记录新增字段
db.getCollection('ClassRecordOneDetail').updateMany({'synState':{'$exists':false}},{'$set':{'synStat ...
- C#实现RSA加密解密
RSA介绍 RSA公钥加密算法是1977年由Ron Rivest、Adi Shamirh和LenAdleman在(美国麻省理工学院)开发的。RSA取名来自开发他们三者的名字。 RSA的缺点: 产生密钥 ...
- vue自定义组件的递归
1.设置一个div显示数值(这个数值不要设置在组件内部) 2.设置组件的名字 <!DOCTYPE html> <html> <head> <title> ...
- 浏览器Notwork XHR被隐藏了
图片中红色区域内容被隐藏 解决方式,点击此处
- css 和常用快捷键
一.css概述: 1.规则:CSS 规则由选择器,以及一条或多条声明两个部分构成. 2.选择器:选择器通常是您需要改变样式的 HTML 元素. 3.声明:声明是您要设置的样式(每条声明由一个属性和一个 ...
- Spark实现二次排序
一.代码实现 package big.data.analyse.scala.secondsort import org.apache.log4j.{Level, Logger} import org. ...
- day 03 作业 预科
目录 作业 1.简述变量的组成 2.简述变量名的命名规范 3.简述注释的作用 4.使用turtle库构造一幅图,贴在markdown文档中 作业 1.简述变量的组成 变量由变量名.赋值符号.变量值所组 ...
- MySQL Replication--复制基本工作原理
复制工作原理(1) master将改变记录到二进制日志(binary log)中(这些记录叫做二进制日志事件,binary log events):(2) slave将master的binary lo ...
- pandas之数据处理
首先,数据加载 pandas提供了一些用于将表格型数据读取为DataFrame对象的函数,期中read_csv和read_table这两个使用最多. 1.删除重复元素 使用duplicated()函数 ...
- golang读写文件之Scan和Fprintf
1. 标准输入输出 os提供了标准输入输出: Stdin = NewFile(uintptr(syscall.Stdin), "/dev/stdin") Stdout = NewF ...