agc015E - Mr.Aoki Incubator(dp)
题意
平面上有$n$个点,每个点都有一个位置$x_i$,和向右的速度$v_i$
现在要求你对其中的一些点进行染色,当一个点被染色后,在无限距离内与它相遇的点也会被染色
问在可能的$2^n$种染色方案中,有多少种染色方案可以使得最后的点全部被染色
Sol
非常好的dp题。
首先考虑若点$i$被染色后,哪些点会被染色
最显然的两种情况
1. $x_j < x_i$且$v_j > v_i$
2. $x_j > x_i$且$v_j < v_i$
当然还有其他的情况,因为题目中的“染色”是一个连锁反应,也就是说$a$被染色,$a$又遇到了$b$,那么$b$也会被染色
实际上,当我们把所有点按速度从大到小排序。
对于$x_i$,向左找到第一个$x_l < x_i$的最小的$l$(此时同时保证了$v_l$最大)
向右找到第一个$x_r > x_i$的最大的$r$(同理保证了$v_r$最小)
那么区间$[l, r]$内所有的点都会被染色
现在问题转化成了,每个点能覆盖一段区间$[l, r]$,问把区间$[l, N]$覆盖的方案数
预处理出每个点能覆盖的区间$[l, r]$。
这些区间有很多好的性质
满足$l_1 \leqslant l_2 \leqslant \dots \leqslant l_n, r_1 \leqslant r_2 \leqslant \dots \leqslant r_n$
证明显然,拿l来说,对于按速度排序以后的排列中两个位置$x_k, x_{k + 1}$
若$x_{k+1} < x_k$,此时$l_{k + 1} > l_k$
若$x_{k + 1} > x_k$,此时$l_{k + 1} = l_k$
右端点同理
注意一个细节,由于我们是按从大到小排序,这样处理起来不是很方便,直接把所有区间reverse一下
设$f[i]$表示用前$i$个区间覆盖$[1, n]$的方案
$s[i]$为其前缀和
前缀和优化dp即可
统计方案那里不是特别懂,路过的大佬可以说一下思路qwq
到了考场上还是写树状数组吧,无脑一点
#include<bits/stdc++.h>
#define Pair pair<int, int>
#define fi first
#define se second
#define MP(x, y) make_pair(x, y)
#define LL long long
using namespace std;
const int MAXN = 1e6, mod = 1e9 + ;
inline int read() {
char c = getchar(); int x = , f = ;
while(c < '' || c > '') {if(c == '-') f = -; c = getchar();}
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x * f;
}
int N, date[MAXN], L[MAXN], R[MAXN], f[MAXN];
struct Node {
int x, v, id, l, r;
bool operator < (const Node &rhs) const {
return v > rhs.v;
}
}a[MAXN];
int comp(const Node &a, const Node &b) {
return a.r == b.r ? a.l < b.l : a.r < b.r;
}
void GetL(int *op) {
vector<Pair> q;
for(int i = ; i <= N; i++) q.push_back(MP(a[i].x, a[i].id));
sort(q.begin(), q.end());
for(int i = ; i <= N; i++) {
while(!q.empty()) {
if(a[i].x <= q.back().fi) op[q.back().se] = i, q.pop_back();
else break;
}
}
}
void GetR(int *op) {
vector<Pair> q;
for(int i = ; i <= N; i++) q.push_back(MP(a[i].x, a[i].id));
sort(q.begin(), q.end());
reverse(q.begin(), q.end());
for(int i = N; i >= ; i--) {
while(!q.empty()) {
if(a[i].x >= q.back().fi) op[q.back().se] = i, q.pop_back();
else break;
}
}
}
int dp() {
for(int i = ; i <= N; i++) a[i].r = N - L[i] + , a[i].l = N - R[i] + ;
sort(a + , a + N + , comp);
//for(int i = 1; i <= N; i++) printf("%d %d\n", a[i].l, a[i].r);
static int s[MAXN], f[MAXN];//f[i]:前i只猫覆盖所有区间的方案
int ans = ;
for(int i = , j = ; i <= N; i++) {
while(a[j].r < a[i].l - ) j++;
f[i] = (s[i - ] + mod - s[j - ] ) % mod;
if(a[i].l == ) f[i]++;
if(a[i].r == N) (ans += f[i]) %= mod;
s[i] = (s[i - ] + f[i]) % mod;
}
return ans % mod;
}
main() {
N = read();
for(int i = ; i <= N; i++) a[i].x = read(), a[i].v = read(), a[i].id = i, date[i] = a[i].v;
sort(date + , date + N + );
int num = unique(date + , date + N + ) - date - ;
for(int i = ; i <= N; i++) a[i].v = lower_bound(date + , date + num + , a[i].v) - date;
sort(a + , a + N + );
GetL(L); GetR(R);
//for(int i = 1; i <= N; i++) printf("%d %d\n", L[i], R[i]);
printf("%d", dp() % mod); }
agc015E - Mr.Aoki Incubator(dp)的更多相关文章
- 【AGC015E】Mr.Aoki Incubator DP
题目描述 数轴上有\(n\)个人,每个人的位置是\(x_i\),速度是\(v_i\). 最开始有一些人感染了传染病. 如果某一时刻一个正常人和一个被感染的人处于同一位置,那么这个正常人也会被感染. 问 ...
- AGC015E Mr.Aoki Incubator
atcoder luogu 首先可以考虑给一个人\(A\)染色.其他人被染色,要么被本来在后面的速度更快的人染色,要么被在前面的更慢的人染色.然后假设一个速度比最开始那个人慢的人\(B\)最后被染色了 ...
- AGC 015 E - Mr.Aoki Incubator
E - Mr.Aoki Incubator 链接 题意: 数轴上有N个黑点,每个点都有一个方向向右的正速度v.当两个点在同一个位置上重合时,若其中一个是红色,另一个也变成红色.保证没有相同速度或初始坐 ...
- AtCoder Grand Contest 015 E - Mr.Aoki Incubator
题目传送门:https://agc015.contest.atcoder.jp/tasks/agc015_e 题目大意: 数轴上有\(N\)个点,每个点初始时在位置\(X_i\),以\(V_i\)的速 ...
- 多校联训 DS 专题
CF1039D You Are Given a Tree 容易发现,当 \(k\) 不断增大时,答案不断减小,且 \(k\) 的答案不超过 \(\lfloor\frac {n}{k}\rfloor\) ...
- AtCoder Grand Contest 015
传送门 A - A+...+B Problem 题意:n个数最大值a,最小值b,求和的可能数量. #include<cstdio> #include<algorithm> us ...
- 【AtCoder】AGC015
AGC015 A - A+...+B Problem #include <bits/stdc++.h> #define fi first #define se second #define ...
- 做题记录 To 2019.2.13
2019-01-18 4543: [POI2014]Hotel加强版:长链剖分+树形dp. 3653: 谈笑风生:dfs序+主席树. POJ 3678 Katu Puzzle:2-sat问题,给n个变 ...
- AtCoder Grand Contest 015 题解
A - A+...+B Problem 常识 Problem Statement Snuke has N integers. Among them, the smallest is A, and th ...
随机推荐
- JNA 传参char[] 和结构体等
近日项目中需要用java调用c/c++编写的dll库,所有了解到jna这个东东,下面是使用的一些经验: 一.java使用Jna需要两个jar包,eg:jna-3.5.1.jar和platform-3. ...
- ArrayList反转
Collections.reverse(arrayList); 下面看例子: import java.util.ArrayList;import java.util.Collections; publ ...
- 7.27实习培训日志-Oracle SQL(三)
Oracle SQL(三) 视图 特性 简单视图 复杂视图 关联的表数量 1个 1个或多个 查询中包含函数 否 是 查询中包含分组数据 否 是 允许对视图进行DML操作 是 否 CREATE [OR ...
- doxygen+ graphviz 开源工具生成源码调用树的wiki
当拿到一含有大量代码的工程怎么看?!这时一个好的代码分析工具非常有用,网上有很多开源工具,但资料都参差不齐,偶然发现doxygen+ graphviz这两工具非常棒,使用工具直接生成函数调用链图,帮助 ...
- webconfig配置详解--转
花了点时间整理了一下ASP.NET Web.config配置文件的基本使用方法.很适合新手参看,由于Web.config在使用很灵活,可以自定义一些节点.所以这里只介绍一些比较常用的节点. <? ...
- sql server 中raiserror的使用
server数据库中raiserror的作用就和asp.net中的throw new Exception一样,用于抛出一个异常或错误.这个错误可以被程序捕捉到. raiserror('错误的描述',错 ...
- 洛谷P1092 虫食算
P1092 虫食算 题目描述 所谓虫食算,就是原先的算式中有一部分被虫子啃掉了,需要我们根据剩下的数字来判定被啃掉的字母.来看一个简单的例子: http://paste.ubuntu.com/2544 ...
- Mysql-2-数据库基础
(1)mysql是一个小型关系型数据库管理系统. (2)mysql是一个快速.多线程.多用户.健壮的SQL数据库服务器.与其他数据库管理系统比,mysql有以下的优势: mysql是一个关系数据库管理 ...
- C# Stack堆栈的使用方法
堆栈(Stack)代表了一个后进先出的对象集合.当您需要对各项进行后进先出的访问时,则使用堆栈.当您在列表中添加一项,称为推入元素,当您从列表中移除一项时,称为弹出元素. Stack 类的方法和属性 ...
- SMTP服务器配置
Windows Server 2012/2012 R2:安装和配置 SMTP 服务器 安装 SMTP 服务器 以下是安装 SMTP 服务器功能的步骤: 1. 打开“服务器管理器”:单击键盘上的 Win ...