[USACO18DEC]Sort It Out P
初看本题毫无思路,只能从特殊的 \(K = 1\) 出发。
但是直接考虑构造一组字典序最小的方案还是不好构造,可以考虑先手玩一下样例。通过自己手玩的样例可以发现,貌似没有被选出来的数在原排列中都是递增的子序列。比如说 \(1 \ 6 \ 5 \ 3 \ 4 \ 2\) 可以证明 \(2 \ 5 \ 6\) 是字典序最小的方案,那么剩下的数按照在原排列中的顺序依次写下是 \(1 \ 3 \ 4\) 是一个递增序列。可以思考一下为什么会有这样的结论出现。
不难发现不论我们怎么操作,原来序列中顺序对的相对位置必然是不会变动的。比如说上方的 \(1 \ 3 \ 4\) 不论怎么变换 \(1\) 都会在 \(3\) 之前,\(3\) 都会在 \(4\) 之前。并且你会发现将剩下的所有数按照字典序依次进行操作是能让这个排列重新排好的。于是我们可以得到一个结论,原排列中的任意一个递增序列保留下来剩下的作为一个操作集合都能时这个排列重新排好。
既然如此,为了让操作集合字典序最小,相反我们需要让保留集合字典序最大。因此保留集合就一定需要选择最长递增子序列中字典序最大的那个。可以考虑每次贪心地选取能拼成最长递增子序列中的最大的元素即可,具体实现可以先求出 \(dp_i\) 表示以 \(i\) 位置结尾的最长递增子序列长度。最后将每个位置按照 \(dp\) 值为第一关键字,按照 \(a_i\) 为第二关键字排序即可。
那么回来思考原问题,字典序第 \(K\) 小的方案怎么求。不难发现实际上我们还是要求保留集合第 \(K\) 大的方案。于是可以考虑从低到高位逐步确定选择的集合,那么我们就需要求出 \(dp_i\) 表示以 \(i\) 开头的最长递增子序列的数量。可以先考虑一个朴素的求法,令 \(f_i\) 为以 \(i\) 开头的最长递增子序列的长度,那么会有转移:
\]
实际上又因为 \(f_i = \max\limits_{a_i < a_j, i < j} f_j + 1\),实际上我们的 \(dp\) 就是求权值在一段后缀中 \(f\) 最大值的数量,这个可以直接在权值线段树上实现,于是我们就将求 \(dp_i\) 的复杂度降至 \(O(n \log n)\)。
那么最后确定第 \(K\) 大的保留集合时,将所有 \(dp\) 值相同的位置压入到一个 \(vector\) 当中并按照权值为第二关键字排序,依次确定每个位即可。
需要注意的是 \(dp_i\) 的数量可能会爆 \(long long\) 但因为我们只需要和 \(K\) 比大小所以需要随时将 \(dp\) 值与 \(1e18\) 取 \(\min\)。
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define ls (p << 1)
#define rs (p << 1 | 1)
#define mid (l + r >> 1)
#define rep(i, l, r) for (int i = l; i <= r; ++i)
#define dep(i, l, r) for (int i = r; i >= l; --i)
const int N = 100000 + 5;
const int inf = 1000000000000000001;
struct tree {
int mx, cnt;
}dp[N], t[N << 2];
int n, k, len, ans, a[N], book[N];
vector <int> G[N];
int read() {
char c; int x = 0, f = 1;
c = getchar();
while (c > '9' || c < '0') { if(c == '-') f = -1; c = getchar();}
while (c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
tree up(tree x, tree y) {
tree ans;
if(x.mx > y.mx) ans = x;
else if(x.mx < y.mx) ans = y;
else ans.mx = x.mx, ans.cnt = min(x.cnt + y.cnt, inf);
return ans;
}
void update(int p, int l, int r, int x, int y, tree k) {
if(l >= x && r <= y) { t[p] = k; return;}
if(mid >= x) update(ls, l, mid, x, y, k);
if(mid < y) update(rs, mid + 1, r, x, y, k);
t[p] = up(t[ls], t[rs]);
}
tree query(int p, int l, int r, int x, int y) {
if(l >= x && r <= y) return t[p];
tree ans; ans.mx = ans.cnt = 0;
if(mid >= x) ans = up(ans, query(ls, l, mid, x, y));
if(mid < y) ans = up(ans, query(rs, mid + 1, r, x, y));
return ans;
}
bool cmp(int x, int y) {
return a[x] > a[y];
}
signed main() {
n = read(), k = read();
rep(i, 1, n) a[i] = read();
dep(i, 1, n) {
dp[i] = query(1, 1, n, a[i], n);
++dp[i].mx, dp[i].cnt = (dp[i].mx == 1 ? dp[i].cnt + 1 : dp[i].cnt);
update(1, 1, n, a[i], a[i], dp[i]);
G[dp[i].mx].push_back(i), len = max(len, dp[i].mx);
}
rep(i, 1, len) sort(G[i].begin(), G[i].end(), cmp);
int P = 0;
dep(i, 1, len) {
for (int j = 0; j < G[i].size(); ++j) if(G[i][j] > P && a[G[i][j]] > a[P]) {
if(dp[G[i][j]].cnt < k) k -= dp[G[i][j]].cnt;
else { P = G[i][j]; break;}
}
book[a[P]] = 1; if(P != n + 1 && !ans) ans = i;
}
printf("%lld\n", n - ans);
rep(i, 1, n) if(!book[i]) printf("%lld\n", i);
return 0;
}
值得一提的是,在这种最优性问题毫无思路时,可以通过手玩样例发现一些策略和性质。
[USACO18DEC]Sort It Out P的更多相关文章
- 题解-USACO18DEC Sort It Out
Problem 洛谷5156 题意概要:给定一个长为\(n\)的排列,可以选择一个集合\(S\)使这个集合内部元素排到自己在整个序列中应该在的位置(即对于集合\(S\)内的每一个元素\(i\),使其排 ...
- 洛谷P5156 [USACO18DEC]Sort It Out
这题就是让你求字典序第k小的最短乱序子序列 转换一下,其实就是字典序第k大的最长上升子序列 就统计一下以i结尾的最长上升子序列\(f[i]\),长度为i的上升子序列的开头组成的集合\(v[i]\),转 ...
- p5156 [USACO18DEC]Sort It Out
传送门 分析 我们发现对于没有发现的点相对位置不会发生改变 于是我们可以吧问题转化为求一个lis 于是我们字典序第k小的答案就是字典序第k大的lis 代码 #include<iostream&g ...
- [USACO18DEC]Sort It Out(树状数组)
[Luogu5156] 题解 求字典序第 k 小的满足题意的集合,取反一下,就是求序列中字典序第 k 大的最长上升子序列 [51nod1376] 最长递增子序列的数量 置 \(f_{i}\)表示以权值 ...
- [算法]——归并排序(Merge Sort)
归并排序(Merge Sort)与快速排序思想类似:将待排序数据分成两部分,继续将两个子部分进行递归的归并排序:然后将已经有序的两个子部分进行合并,最终完成排序.其时间复杂度与快速排序均为O(nlog ...
- [算法]——快速排序(Quick Sort)
顾名思义,快速排序(quick sort)速度十分快,时间复杂度为O(nlogn).虽然从此角度讲,也有很多排序算法如归并排序.堆排序甚至希尔排序等,都能达到如此快速,但是快速排序使用更加广泛,以至于 ...
- shell之sort命令
1 sort的工作原理 sort将文件的每一行作为一个单位,相互比较,比较原则是从首字符向后,依次按ASCII码值进行比较,最后将他们按升序输出. [rocrocket@rocrocket progr ...
- 详细解说 STL 排序(Sort)
0 前言: STL,为什么你必须掌握 对于程序员来说,数据结构是必修的一门课.从查找到排序,从链表到二叉树,几乎所有的算法和原理都需要理解,理解不了也要死记硬背下来.幸运的是这些理论都已经比较成熟,算 ...
- SQL Tuning 基础概述06 - 表的关联方式:Nested Loops Join,Merge Sort Join & Hash Join
nested loops join(嵌套循环) 驱动表返回几条结果集,被驱动表访问多少次,有驱动顺序,无须排序,无任何限制. 驱动表限制条件有索引,被驱动表连接条件有索引. hints:use_n ...
随机推荐
- HDU 1052:Tian Ji -- The Horse Racing(贪心)
Tian Ji -- The Horse Racing Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Jav ...
- RotateRect(旋转矩形)的倾斜旋转变换矫正
在Opencv中的图像处理中,经常要用到minAreaRect()函数求最小外接矩形,该函数的返回值就是一个RotatedRect类对象. RotatedRect类定义如下: class CV_EXP ...
- 前端后端通信初步尝试(javascript - flask)
在某项目中,需要使用python flask做后端功能开发,web提供功能入口. 此时需要使用Ajax通信. 由于以前从未接触过网络传输,记录了一些基础知识. 资料参考<HTML5+CSS3+J ...
- SpringBootAdmin-使用踩坑
一.版本选择 目前项目中使用的springcloud版本为 Hoxton.SR3,使用的spring-cloud-alibaba版本为2.2.1.RELEASE, 使用的springboot版本为2. ...
- “老”的Flexbox和“新”的Flexbox
本文由大漠根据Chris Coyier的<"Old" Flexbox and "New" Flexbox>所译,整个译文带有我们自己的理解与思想,如 ...
- Python网络编程之网络基础
Python网络编程之网络基础 目录 Python网络编程之网络基础 1. 计算机网络发展 1.1. OSI七层模型 1.2. 七层模型传输数据过程 2. TCP/IP协议栈 2.1 TCP/IP和O ...
- 【Java】回形数
回形数 键盘读入一个整数n(1-20),以n为矩阵大小,把1,2,3,4,5-按顺时针螺旋的形式填入. import java.util.Scanner; public class HuiXingSh ...
- JavaScript 中BOM的常用操作
JavaScript BOM操作 1.获取浏览器窗口尺寸 var width=window,innerWidth //获取可视窗口宽度 var height=window.innerHeight // ...
- JQuery选择器的使用和分类
jQuery选择器 id选择器格式 $("#box") //获取标签里的id是box的标签 类选择器格式 $(".a") //获取标签里的类名是a的标签 标签选 ...
- ClassCastException: java.util.Date cannot be cast to java.sql.Date
解决办法 /** * 单个方法,作用,根据输入的day:yyyy-mm-dd格式的字符日期,将数据库中的该天所有数据更新为0 * 0表示假期 * @param day * @throws SQLExc ...