UR #13 Ernd
考试的时候没有注意到可以将(a,b)放在二维平面上之后旋转坐标系,使得转移变成树状数组二维偏序
这样就算我想出来了第二个转移的斜率优化也没有什么卵用啊(摔西瓜
设g(i)表示当前站在第i个水果下面且第i个水果此时并没有记分的最大得分
设f(i)表示当前站在第i个水果下面且第i个水果此时已经记分的最大得分
g(i)=max(f(j)) (满足j能到达i)
f(i)=max(g(j)+(i-j+1)^2)(满足从i出发可以一路接水果走到j)
首先考虑g(i)的转移,我们注意到j能到达i当且仅当
abs(ai-aj)<=bi-bj
又注意到bi>=bj
我们对绝对值进行分类讨论后可以得到j可以到达的i都在j的10点半钟方向到13点半钟方向
之后我们将坐标系旋转,即将每个点的坐标变为(bi+ai,bi-ai)
那么j能到达的i都在j的右上方啦,那么对于这个转移就变成了一个二维偏序问题
对于一维排序,另一维用树状数组维护即可,时间复杂度O(nlogn)
之后我们考虑f(i)的转移,我们很容易发现这是可以分段的
对于每个i不能到达i+1分一段,那么这个转移每一段都是相互独立的
我们又惊奇的发现在第一个转移排序后,每一段内的相对顺序是不变的
注意到这是一个很经典的可以用斜率优化的式子
化简完后f(i)=max(g(j)-2*j+j*j-2*i*j)+(i+1)^2
对于斜率2*i,我们用很容易发现这个斜率是递增的
我们可以对于每一段用单调栈维护一个上凸包,这样这个转移就是O(n)的了
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std; typedef long long LL;
const int maxn=500010;
int n,cnt;
int a[maxn],b[maxn];
int pos[maxn];
LL dp[maxn];
LL t[maxn];
LL m1,m2,ans;
struct Point{
LL x,y,id;
Point(LL x=0,LL y=0,LL id=0):x(x),y(y),id(id){}
}p[maxn];
vector<Point>V[maxn];
int top[maxn];
typedef Point Vector;
Vector operator -(const Point &A,const Point &B){return Vector(A.x-B.x,A.y-B.y);}
LL Cross(const Vector &A,const Vector &B){return A.x*B.y-A.y*B.x;}
bool cmpy(const Point &A,const Point &B){return A.y<B.y;}
bool cmpx(const Point &A,const Point &B){
if(A.x==B.x&&A.y==B.y)return A.id<B.id;
if(A.x==B.x)return A.y<B.y;
return A.x<B.x;
}
void read(int &num){
num=0;char ch=getchar();
while(ch<'!')ch=getchar();
while(ch>='0'&&ch<='9')num=num*10+ch-'0',ch=getchar();
}
int lowbit(int x){return x&(-x);}
void modify(int x,LL v){for(int i=x;i<=cnt;i+=lowbit(i))t[i]=max(t[i],v);}
LL ask(int x){
LL mx=0;
for(int i=x;i>=1;i-=lowbit(i))mx=max(mx,t[i]);
return mx;
}
LL Get_ans(int id,LL k){
while(V[id].size()>1&&V[id][top[id]].y-2*k*V[id][top[id]].x<=V[id][top[id]-1].y-2*k*V[id][top[id]-1].x)top[id]--,V[id].pop_back();
if(!V[id].empty())return V[id][top[id]].y-2*k*V[id][top[id]].x+(k+1)*(k+1);
else return 0;
}
void push(int id,Point P){
while(V[id].size()>1&&Cross(P-V[id][top[id]],V[id][top[id]]-V[id][top[id]-1])<=0)top[id]--,V[id].pop_back();
top[id]++;V[id].push_back(P);
} int main(){
read(n);
for(int i=1;i<=n;++i){
read(a[i]);read(b[i]);
p[i]=Point(b[i]+a[i],b[i]-a[i],i);
}
for(int i=1;i<=n;++i){
if(i==1||abs(a[i]-a[i-1])>b[i]-b[i-1])pos[i]=pos[i-1]+1;
else pos[i]=pos[i-1];
}
for(int i=1;i<=n;++i)a[i]=b[i]-a[i];
sort(a+1,a+n+1);sort(p+1,p+n+1,cmpy);
for(int i=1;i<=n;++i){
if(i==1||a[i]!=a[i-1])p[i].y=++cnt;
else p[i].y=cnt;
}
sort(p+1,p+n+1,cmpx);
memset(top,-1,sizeof(top));
for(int i=1;i<=n;++i){
m1=ask(p[i].y);m2=Get_ans(pos[p[i].id],p[i].id);
dp[p[i].id]=max(m1+1,m2);
ans=max(dp[p[i].id],ans);
modify(p[i].y,dp[p[i].id]);
push(pos[p[i].id],Point(p[i].id,m1-2*p[i].id+p[i].id*p[i].id));
}printf("%lld\n",ans);
return 0;
}
UR #13 Ernd的更多相关文章
- uoj#187. 【UR #13】Ernd
http://uoj.ac/problem/187 每个点只能从时间,b+a,b-a三维都不大于它的点转移过来,将点按时间分成尽量少的一些段,每段内三维同时非严格单调,每段内的点可能因为连续选一段而产 ...
- UR #13 Yist
第一次打UR,打了一个半小时就弃疗了QAQ 这是我唯一一道考试的时候做出来的题目,其他两道连暴力都懒得写了 很容易发现对于每个要删除的点 我们找到左边第一个比他小的不用删除的点,右边第一个比他小的不用 ...
- UR#13 SRAND
总感觉这位大仙讲的很清楚:bztminamoto 题意 题目讲的是求 l~r 内所有数的次大质因子,这里设 f(x) 为 x 的次大质因子 我们差分一下就变成求两个前缀和信息了 按照套路,我们考虑 S ...
- uoj#188. 【UR #13】Sanrd(Min_25筛)
题面 传送门 题解 这是一道语文题 不难看出,题目所求即为\(l\)到\(r\)中每个数的次大质因子 我们考虑\(Min\_25\)筛的过程,设 \[S(n,j)=\sum_{i=1}^nsec_p( ...
- UOJ 188 【UR #13】Sanrd——min_25筛
题目:http://uoj.ac/problem/188 令 \( s(n,j)=\sum\limits_{i=1}^{n}[min_i>=p_j]f(j) \) ,其中 \( min_i \) ...
- UOJ #188. 【UR #13】Sanrd
Description 给定 \(\sum_{i=l}^r f[i]\) \(f[i]=\) 把 \(i\) 的每一个质因子都从小到大排列成一个序列(\(p_i^{c_i}\)要出现 \(c_i\) ...
- UOJ188. 【UR #13】Sanrd
传送门 Sol 设 \(f_i\) 表示 \(i\) 的次大质因子 题目就是要求 \[\sum_{i=l}^{r}f_i\] 考虑求 \(\sum_{i=1}^{n}f_i\) 所求的东西和质因子有关 ...
- 「uoj#188. 【UR #13】Sanrd」
题目 不是很能看懂题意,其实就是求\([l,r]\)区间内所有数的次大质因子的和 这可真是看起来有点鬼畜啊 这显然不是一个积性函数啊,不要考虑什么特殊的函数了 我们考虑Min_25筛的过程 设\(S( ...
- UOJ188. 【UR #13】Sanrd [min_25筛]
传送门 思路 也可以算是一个板题了吧qwq 考虑min_25筛最后递归(也就是DP)的过程,要枚举当前最小的质因子是多少. 那么可以分类讨论,考虑现在这个质因子是否就是次大质因子. 如果不是,那么就是 ...
随机推荐
- AutoPoco的使用
官方开发指导https://autopoco.codeplex.com/documentation 初步使用: SimpleUser是自己要批量创建的类 1)创建管理工厂 IGenerationSes ...
- 膜拜 2014-2 (献给L之三)
文/安然 深深的夜静静的想你细细的品味满满的甜蜜爱,心灵的对话让我流泪战栗谢谢你给我这么弥足珍贵的体会不能一生相随又有什么关系你一直都在我的爱就永不停息此生足以——献给心底的爱
- AE实现矢量图层标注属性
添加引用ESRI.ArcGIS.Carto 1.获取图层 IGeoFeatureLayer pFtrLayer = m_pLayer as IGeoFeatureLayer; 2.初始化标注属性集合 ...
- GDI+
1, 编译error的话一般是却 #include <comdef.h>#include <Windows.h> Windows.h内会包含Windows.h,但是因为在std ...
- iOS 非ARC基本内存管理系列 4-autorelease方法和@autoreleasepool
1.autorelease 基本用法 对象执行autorelease方法时会将对象添加到自动释放池中 当自动释放池销毁时自动释放池中所有对象作release操作 对象执行autorelease方法后自 ...
- 在HTML5中怎样实现Canvas阴影效果
该文章是由e良师益友技术部小陈原创作品,转载是请注明来源,谢谢! 今天我给大家介绍一下在HTML5中怎样实现Canvas阴影效果,我们知道现在HTML5的Canvas阴影也经常使用的,这个就是HTML ...
- 从 Java 代码逆向工程生成 UML 类图和序列图
from:http://blog.itpub.net/14780914/viewspace-588975/ 本文面向于那些软件架构师,设计师和开发人员,他们想使用 IBM® Rational® Sof ...
- css和jQuery ,字符串中重新定义单个字符样式背景
<table class="tb2"> <tr><td class="td1">融资登记企业<span>985& ...
- Keep two divs sync scroll and example
srcDiv has visible horizontal scrollbar.(style="overflow:auto;") targetDiv has no scrollba ...
- linux 输入子系统(2)----简单实例分析系统结构(input_dev层)
实例代码如下: #include <linux/input.h> #include <linux/module.h> #include <linux/init.h> ...