Codeforces 889C Maximum Element(DP + 计数)
题目链接 Maximum Element
题意 现在有这一段求序列中最大值的程度片段:
(假定序列是一个1-n的排列)
int fast_max(int n, int a[]) {
int ans = 0;
int offset = 0;
for (int i = 0; i < n; ++i)
if (ans < a[i]) {
ans = a[i];
offset = 0;
} else {
offset = offset + 1;
if (offset == k)
return ans;
}
return ans;
}
显然这段程序是错误的……有很多可以X掉这段程序的排列
求这样的排列有多少个。
题目是让我们求符合这样条件的排列个数:
1、存在某个数,他比前面的数都大并且小于$n$;
2、他比他后面$k$个数都要大。
假设“中间这个数”为$cnt$
假设$D(i)$为满足$p(i) = i$的这样的排列个数
我们可以把$D(i)$的求解分成两个过程。
1、计算$cnt$等于$i - 1$的排列个数
2、计算$cnt$不等于$i - 1$的排列个数
首先如果$i <= k + 1$,则$D(i) = 0$
当这个序列的$cnt$为$i - 1$时,只要满足$i - 1$和$i$之间的数大于等于$k$个即可。
于是对于$i - 1$这个数的位置的选择,我们有$i - k - 1$种。
然后呢,除了$i - 1$和$i$这两个数,其他数的位置随意(因为$i$排在最后,所以排在$i - 1$前的数字都比$i - 1$要小)
所以当前这种情况对答案的贡献为$(i - k - 1) * (i - 2)!$
考虑另外一种情况。
当$cnt$不等于$i - 1$的时候,一定有$cnt < i - 1$
设$i - 1$所在位置为$pos$,我们把$i - 1$之前的$pos - 1$个数离散化成一个$1$到$pos - 1$的排列
然后在这个排列的最后加上$pos$,就构成了一个$1$到$pos$并且以$pos$结尾的排列
那么如果这个排列是符合要求的,那么整个排列也是符合要求的。
于是我们枚举$i - 1$的位置$pos$,满足条件的位置为$i - k <= pos <= i - 1$
我们在剩下的$i - 2$个数中选出$pos - 1$个放到前$pos - 1$个位置,然后乘上$D(pos)$。
然后还要乘上$(i - pos - 1)!$,因为$i - 1$到$i$之间的数都是随意乱放的……
于是当前这种情况对答案的贡献为
于是我们终于推出了D(n)的公式
最后的答案怎么计算呢
我们假设$n$的位置为$pos$
那么当$p(pos) = n$的时候,前pos个数的方案数为$D(pos) * C(n - 1, pos - 1)$
后$n - pos$个数的方案数为$(n - pos)!$
所以当$p(pos) = n$的时候对答案的贡献为$D(pos) * C(n - 1, pos - 1)*(n - pos)!$
枚举$pos$,累加即可。
#include <bits/stdc++.h> using namespace std; #define rep(i, a, b) for (int i(a); i <= (b); ++i)
#define dec(i, a, b) for (int i(a); i >= (b); --i) const int N = 1e6 + 10;
const int mod = 1e9 + 7; int n, k;
int f[N], s[N];
int fac[N], inv[N];
int ans = 0; inline int Pow(int a, int b, int mod){
int ret(1);
for (; b; b >>= 1, a = (1ll * a * a) % mod) if (b & 1) (ret = 1ll * ret * a) % mod;
return ret;
} void init(){
fac[0] = 1;
rep(i, 1, 1e6 + 1) fac[i] = 1ll * fac[i - 1] * i % mod;
rep(i, 1, 1e6 + 1) inv[i] = Pow(fac[i], mod - 2, mod);
} inline void up(int &a, int b) { a = (0ll + a + b) % mod;}
inline void mulup(int &a, int b){ a = 1ll * a * b % mod;} int main(){ scanf("%d%d", &n, &k);
init(); rep(i, k + 2, n){
f[i] = i - k - 1;
up(f[i], s[i - 1] - s[i - k - 1]);
mulup(f[i], fac[i - 2]);
s[i] = (0ll + s[i - 1] + 1ll * f[i] * inv[i - 1] % mod) % mod;
} rep(i, 1, n) up(ans, (int)1ll * f[i] * fac[n - 1] % mod * inv[i - 1] % mod);
printf("%d\n", ans);
return 0;
}
Codeforces 889C Maximum Element(DP + 计数)的更多相关文章
- 【CF886E】Maximum Element DP
[CF886E]Maximum Element 题意:小P有一个1-n的序列,他想找到整个序列中最大值的出现位置,但是他觉得O(n)扫一遍太慢了,所以它采用了如下方法: 1.逐个遍历每个元素,如果这个 ...
- Codeforces Round #445 Div. 1 C Maximum Element (dp + 组合数学)
题目链接: http://codeforces.com/contest/889/problem/C 题意: 给你 \(n\)和 \(k\). 让你找一种全排列长度为\(n\)的 \(p\),满足存在下 ...
- Codeforces 886E Maximum Element 组合数学 + dp
我们定义dp[ i ]表示长度为 i 的序列, 最后没有一个==k的时候返回的方案数, 也就是最后强制返回 i 的方案数. 我们能得到dp方程 dp[ i ] = sum(dp[ i - j - ...
- Codeforces 702A Maximum Increase(dp)
题目链接:http://codeforces.com/problemset/problem/702/A 题意: 给你N个数,a[0], a[1], a[2], ....., a[n-1],让你找出最长 ...
- Codeforces - 102222A - Maximum Element In A Stack - 模拟
https://codeforc.es/gym/102222/problem/F 注意到其实用unsigned long long不会溢出. #include<bits/stdc++.h> ...
- Tetrahedron(Codeforces Round #113 (Div. 2) + 打表找规律 + dp计数)
题目链接: https://codeforces.com/contest/166/problem/E 题目: 题意: 给你一个三菱锥,初始时你在D点,然后你每次可以往相邻的顶点移动,问你第n步回到D点 ...
- 【CodeForces】889 C. Maximum Element 排列组合+动态规划
[题目]C. Maximum Element [题意]给定n和k,定义一个排列是好的当且仅当存在一个位置i,满足对于所有的j=[1,i-1]&&[i+1,i+k]有a[i]>a[ ...
- 动态规划(DP计数):HDU 5116 Everlasting L
Matt loves letter L.A point set P is (a, b)-L if and only if there exists x, y satisfying:P = {(x, y ...
- Codeforces 484B Maximum Value(高效+二分)
题目链接:Codeforces 484B Maximum Value 题目大意:给定一个序列,找到连个数ai和aj,ai%aj尽量大,而且ai≥aj 解题思路:类似于素数筛选法的方式,每次枚举aj,然 ...
随机推荐
- destoon去除编辑器替换图片删除原图功能,删除信息删除相关图片功能
去除这些功能会造成大量垃圾图片,但是客户存在大量复制内容,其中图片一样,为了防止客户替换其中一个图片或者删除信息 造成其他复制信息图片丢失 去除文章模型级联图片功能. 对应模块class.php se ...
- 《linux设备驱动开发详解》笔记——8阻塞与非阻塞IO
8.1 阻塞与非阻塞IO 8.1.0 概述 阻塞:访问设备时,若不能获取资源,则进程挂起,进入睡眠状态:也就是进入等待队列 非阻塞:不能获取资源时,不睡眠,要么退出.要么一直查询:直接退出且无资源时, ...
- opencv和numpy的安装
近日,学姐让我们切割图片,查了一下资料,发现我需要安装opencv和numpy.但是在安装过程中却出现了很多小问题,我在此结合自和自己的安装经验和网上查找的资料,做一个笔记. 1.opencv的安装 ...
- Latex:插入伪代码
*本文属于转载. *转载链接:https://blog.csdn.net/lwb102063/article/details/53046265 目录 clrscode algorithm algori ...
- WPF触控程序的开发(一)——有用的资源
迟来的一篇博文,每次都要撞到月末,这个月实在太忙了,除了在公司上班,还接了个单子,用wpf做一个触屏软件,类似iphone的相册功能.先说搭建开发环境吧,我是不可能去买个平板来的,再说基于win7的程 ...
- C#通过http post方式调用需要证书的webservice
前一段时间做花旗银行的项目,用到花旗的接口是websevice,由于很多原因直接在项目中引用webservice不成功,于是就用了http post方式请求,把请求信息(xml格式)组装之后发送到服务 ...
- Python学习-day7 类 部分socket
这周还是继续关于类的学习,在面向对象的学习过程中又学习了网络编程,并且提交了编写FTP的作业. 复习一下类的相关概念和定义 类 属性 实例变量:内存中 ...
- python学习-- Django传递数据给JS
var List = {{ List|safe }};//safe 必须存在
- 优秀的缓存请求库,快速请求接口和图片:WTRequestCenter
WTRequestCenter 方便缓存的请求库无需任何import和配置,目前实现了基础需求如果有其他需要请在issue 上提出,谢谢! 使用方法 Usage 注意:所有的请求都是缓存的 GET 请 ...
- 面向对象编程(二)封装--构造方法,this关键字,static关键字,方法重载
面向对象三大特点:封装.继承.多态 封装概念 ① 将东西包装在一起,然后以新的完整形式呈现出来: 将方法和字段一起包装到一个单元中,单元以类的形式实现; ② 信息隐藏,隐藏对象的实现细节,不让 ...