Slope Trick

问题引入

CF713C
给出$n$个整数$a_1,a_2,...,a_n$，一次操作可以选择一个整数$a_i$使其变为$a_i=a_{i-1}+1$或$a_i=a_{i-1}-1$，问使得$a_1<a_2<...,a_n$所需要的最小操作次数。
$$\begin{gathered} 1\le n\le 3000 \\ 1\le a_i\le 10^9 \end{gathered}$$

解法

$a_i<a_{i+1}\leftrightarrow a_i\le a_{i+1}-1$。
首先令$a_i=a_i-i$，问题等价于使$a_1\le a_2\le ...\le a_n$需要的最小操作次数。
设$f_{i,x}$为前$i$个数满足$a_1\le a_2\le ...\le a_i\le x$所需要的最小操作次数。
有$$\begin{gathered} f_{1,x}=mi{n}_{y\le x}|y-a_1| \\ {f}_{i,x}=mi{n}_{y\le x}(f_{i-1,y}+|y-a_i|) \end{gathered}$$
$f_{i,x}$中$x$的取值可以只考虑集合 { a 1 , a 2 , . . . , a n } \{a_1,a_2,...,a_n\} {a1​,a2​,...,an​}内的，因此这个$dp$可以在$O(n^2)$时间内解决，但下面可以把它优化到$O(nlogn)$。

slope trick

对$f_{i,x}$以$x$为横坐标，$f_{i,x}$为纵坐标，在$(x,f_{i,x})$和$(x+1,f_{i,x+1})$之间连一条直线，我们可以得到一个分段函数，这个分段函数在$x$为整数处的斜率都为整数。

首先$f_{i,x}$是一个单调递减的函数。因为若$x_1<x_2$，则$a_1\le a_2\le ...\le a_i\le x_1\le x_2$。

对$f_{i,x}$维护一个集合 S i = { s i , i − 1 , s i , i − 2 , . . . , s i , 0 } S_i=\{s_{i,i-1},s_{i,i-2},...,s_{i,0}\} Si​={si,i−1​,si,i−2​,...,si,0​}和$f_{i,s_0}$的答案。
其含义为对于$x\in [s_{i,j},s_{i,j-1})(i\le j\le 0)$，有$f_{i,x+1}-f_{i,x}=-j$（即函数$f_{i,x}$在点$x$处的斜率为$-j$，特别的$s_{i,i}=-\infty ,s_{i,-1}=\infty$）。
由$f_{i,x}$的单调性有$s_{i,i-1}\le s_{i,i-2}\le ...\le s_{i,0}$

现在假设求出了$S_{i-1}$和$f_{i-1,s_0}$，现在求$S_i$，考虑三种情况：
$1.x<a_i$，$f_{i-1,x}$在点$x$处的斜率为$-j(j\ge 0)$（即$f_{i-1,x+1}-f_{i-1,x}=-j$）
$$\begin{gathered} f_{i,x}=f_{i-1,x}+(a_i-x) \\ {f}_{i,x+1}=f_{i-1,x+1}+(a_i-(x+1))=f_{i-1,x}-j+a_i-x-1 \\ {f}_{i,x+1}-f_{i,x}=-j-1 \end{gathered}$$
故所有$x<a_i$的点的斜率减$1$，即有$s_{i,j}=s_{i-1,j-1}(s_{i-1,j-1}<a_i)$

$2.x\ge a_i$，$f_{i-1,x}$在点$x$处的斜率为$-j(j>0)$。
$$\begin{gathered} f_{i,x}=f_{i-1,x}+(x-a_i) \\ {f}_{i,x+1}=f_{i-1,x+1}+(x+1-a_i)=f_{i-1,x}-j+x+1-a_i \\ {f}_{i,x+1}-f_{i,x}=-j+1=-(j-1) \end{gathered}$$
故当$x\ge a_i$且$x$的斜率不为$0$时，它的斜率加$1$，有$s_{i,j}=s_{i-1,j+1}(s_{i-1,j+1}\ge a_i\&\&j+1>0)$

$3.x\ge a_i$，$f_{i-1,x}$在点$x$处的斜率为$0$。
$$\begin{gathered} f_{i,x}=f_{i-1,x}+(x-a_i) \\ {f}_{i,x+1}=f_{i-1,x+1}+(x+1-a_i)=f_{i-1,x}+x+1-a_i \\ {f}_{i,x+1}-f_{i,x}=1 \end{gathered}$$
此时由于$a_1\le a_2\le ...\le a_i\le x\le x+1$，因此$f_{i,x+1}$取$f_{i,x}$答案反而更优，即$f_{i,x+1}=min(f_{i,x+1},f_{i,x)}=f_{i,x}$。
因此如果$s_{i-1,0}\ge a_i$，则其不需要加入到$S_i$中。

基于上述结论，再分情况讨论：
$1.s_{i-1,0}<a_i$，保持集合$S_{i-1}$的元素不动。
S i − 1 = { s i − 1 , i − 2 , s i − 1 , i − 3 , . . . , s i − 1 , 0 } S_{i-1}=\{s_{i-1,i-2},s_{i-1,i-3},...,s_{i-1,0}\} Si−1​={si−1,i−2​,si−1,i−3​,...,si−1,0​}转为 S i = { s i , i − 1 , s i , i − 2 , . . . , s i , 1 } S_i=\{s_{i,i-1},s_{i,i-2},...,s_{i,1}\} Si​={si,i−1​,si,i−2​,...,si,1​}，由于$a_i-1$处的斜率变为$-1$，但$a_i$处的斜率仍然为$0$，故再加入$s_{i,0}=a_i$。
$f_{i,s_{i,0}}=f_{i,a_i}=f_{i-1,a_i}+|a_i-a_i|=f_{i-1,a_i}$。

$2.s_{i-1,0}=a_i$，保持集合$S_{i-1}$的元素不动。
$s_{i-1,0}$的斜率不变， S i − 1 = { s i − 1 , i − 2 , s i − 1 , i − 3 , . . . , s i − 1 , 0 } S_{i-1}=\{s_{i-1,i-2},s_{i-1,i-3},...,s_{i-1,0}\} Si−1​={si−1,i−2​,si−1,i−3​,...,si−1,0​}转为 S i = { s i , i − 1 , s i , i − 2 , . . . , s i , 2 , s i , 0 } S_i=\{s_{i,i-1},s_{i,i-2},...,s_{i,2},s_{i,0}\} Si​={si,i−1​,si,i−2​,...,si,2​,si,0​}，缺少$s_{i,1}$，加入$s_{i,1}=s_{i,0}=a_i$即可，因为已经不存在斜率为$-1$的整点，用$[s_{i,1},s_{i,0})=[a,a)$表示这部分为空集。
$f_{i,s_{i,0}}=f_{i,a_i}=f_{i-1,a_i}+|a_i-a_i|=f_{i-1,a_i}$。

$3.s_{i-1,0}>a_i$。
此时需要弹出$s_{i-1,0}$，设$s_{i,k}<a_i,s_{i,k-1}\ge a_i$
S i − 1 = { s i − 1 , i − 2 , s i − 1 , i − 3 , . . . , s i , k , s i − 1 , k − 1 , . . . , s i − 1 , 1 } S_{i-1}=\{s_{i-1,i-2},s_{i-1,i-3},...,s_{i,k},s_{i-1,k-1},...,s_{i-1,1}\} Si−1​={si−1,i−2​,si−1,i−3​,...,si,k​,si−1,k−1​,...,si−1,1​}转为 S i = { s i , i − 1 , s i , i − 2 , . . . , s i , k + 1 , s i , k − 2 , . . . , s i , 2 , s i , 0 } S_i=\{s_{i,i-1},s_{i,i-2},...,s_{i,k+1},s_{i,k-2},...,s_{i,2},s_{i,0}\} Si​={si,i−1​,si,i−2​,...,si,k+1​,si,k−2​,...,si,2​,si,0​}，加入$s_{i,k}=s_{i,k-1}=a_i$。
$f_{i,s_{i,0}}=f_{i-1,s_{i-1,0}}+s_{i-1,0}-a_i$。

由于$s_{i,i-1}\le s_{i,i-2}\le ...\le s_{i,0}$，整个过程只需要维护一个优先队列，从而时间复杂度可以优化到$O(nlogn)$。

#include <bits/stdc++.h>
typedef unsigned long long ull;
typedef long long ll;
#define rep(i, l, r) for (int i = l; i <= r; i++)
#define nep(i, r, l) for (int i = r; i >= l; i--)
void sc(int &x) { scanf("%d", &x); }
void sc(int &x, int &y) { scanf("%d%d", &x, &y); }
void sc(int &x, int &y, int &z) { scanf("%d%d%d", &x, &y, &z); }
void sc(ll &x) { scanf("%lld", &x); }
void sc(ll &x, ll &y) { scanf("%lld%lld", &x, &y); }
void sc(ll &x, ll &y, ll &z) { scanf("%lld%lld%lld", &x, &y, &z); }
void sc(char *x) { scanf("%s", x); }
void sc(char *x, char *y) { scanf("%s%s", x, y); }
void sc(char *x, char *y, char *z) { scanf("%s%s%s", x, y, z); }
void out(int x) { printf("%d\n", x); }
void out(ll x) { printf("%lld\n", x); }
void out(int x, int y) { printf("%d %d\n", x, y); }
void out(ll x, ll y) { printf("%lld %lld\n", x, y); }
void out(int x, int y, int z) { printf("%d %d %d\n", x, y, z); }
void out(ll x, ll y, ll z) { printf("%lld %lld %lld\n", x, y, z); }
void ast(ll x,ll l,ll r){assert(x>=l&&x<=r);}
using namespace std;
const int N=5e3+5,mod=998244353;
const ll inf=1e18;
int n;
priority_queue<int>q;
void sol(int cas)
{
    sc(n);
	ll ans=0;
	rep(i,1,n)
	{
		int x;sc(x);
		x-=i;
		q.push(x);
		if(q.top()>x)
		{
			ans+=q.top()-x;
			q.pop();
			q.push(x);
		}
	}
	out(ans);
}
int main()
{
  // freopen("1.in", "r",stdin);
  // freopen("1.out", "w", stdout);
  srand(time(0));
  int t=1,cas=0;
//   scanf("%d",&t);
  while(t--)
  {
    sol(++cas);
  }
  return 0;
}

abc217H - Snuketoon

最初你在点$0$，每一秒你可以从点$x$走到$x+1,x,x-1$。
有$N$次攻击。
第$i$次攻击发生在时刻$T_i$，设此时在点$x$，如果$D_i=0$，则你会受到$max(0,X_i-x)$点伤害，$D_i=1$则你会受到$max(0,x-X_i)$点伤害。
问$N$次攻击后你受到的总伤害最小是多少。
$$\begin{gathered} 1\le N\le 200000 \\ 1\le T_1\le T_2\le ...\le T_N\le 10^9 \\ 0\le D_i\le 1 \\ -10^9\le X_i\le 10^9 \end{gathered}$$

设$f_{i,x}$表示在$T_i$时刻你在点$x$受到的最小伤害。
有$f_{i,x}=cost(D_i,X_i,x)+\underset{x-(T_i-T_{i-1})\le y\le x+(T_i-T_{i-1})}{\min }{f}_{i-1,y}$

与上个问题不同的是，$f_{i,x}$是一个下凸函数，形如下图：

因此这里需要维护两个优先队列。
一个为$L_i=l_{i,i-1},l_{i,i-2},...,l_{i,0}$，另一个为$R_i=r_{i,0},r_{i,1},...,r_{i,i-1}$。
表示$x\in [l_{i,j},l_{i,j-1})$时$x$处的斜率为$-j$。
$x\in [l_{i,0},r_{i,0})$时$x$处斜率为$0$。
$x\in [r_{i,j},r_{i,j+1})$时$x$处斜率为$j$。
对于$x-(T_i-T_{i-1})\le y\le x+(T_i-T_{i-1})$这个限制。
如图，我们对红色的图像左平移$T_i-T_{i-1}$得到黄色图像，对红色的图像右平移$T_i-T_{i-1}$得到粉色图像。最后我们要取的图像是棕色图像。
即$g_{i,x}=\underset{x-(T_i-T_{i-1})\le y\le x+(T_i-T_{i-1})}{\min }{f}_{i-1,y}$只需要将$L_{i-1}$中的点坐标都减去$T_i-T_{i-1}$，将$R_{i-1}$中的点都加上$T_i-T_{i-1}$就可以得到。
得到$g_{i,x}$后，再将其加上分段函数$Cost(D_i,X_i)$，对$L_i$和$R_i$进行形如上一题中的讨论即可，另外这里只需要维护$ans=f_{i,l_{i,0}}=f_{i,l_{i,0}+1}=...=f_{i,r_{i,0}}$的值即可。
注意这里斜率为$0$的点的斜率也要加$1$或减$1$。
另外一个地方是初始时$<0$处的斜率为$-\infty$，$>0$处的斜率为$\infty$。
所以初始时还要额外插入较多的$0$。

#include <bits/stdc++.h>
typedef unsigned long long ull;
typedef long long ll;
#define rep(i, l, r) for (int i = l; i <= r; i++)
#define nep(i, r, l) for (int i = r; i >= l; i--)
void sc(int &x) { scanf("%d", &x); }
void sc(int &x, int &y) { scanf("%d%d", &x, &y); }
void sc(int &x, int &y, int &z) { scanf("%d%d%d", &x, &y, &z); }
void sc(ll &x) { scanf("%lld", &x); }
void sc(ll &x, ll &y) { scanf("%lld%lld", &x, &y); }
void sc(ll &x, ll &y, ll &z) { scanf("%lld%lld%lld", &x, &y, &z); }
void sc(char *x) { scanf("%s", x); }
void sc(char *x, char *y) { scanf("%s%s", x, y); }
void sc(char *x, char *y, char *z) { scanf("%s%s%s", x, y, z); }
void out(int x) { printf("%d\n", x); }
void out(ll x) { printf("%lld\n", x); }
void out(int x, int y) { printf("%d %d\n", x, y); }
void out(ll x, ll y) { printf("%lld %lld\n", x, y); }
void out(int x, int y, int z) { printf("%d %d %d\n", x, y, z); }
void out(ll x, ll y, ll z) { printf("%lld %lld %lld\n", x, y, z); }
void ast(ll x,ll l,ll r){assert(x>=l&&x<=r);}
using namespace std;
const int N=2e5+5,mod=998244353;
const ll inf=1e18;
int n;
ll ans,sl,sr;
priority_queue<ll>L;
priority_queue<ll,vector<ll>,greater<ll>>R;
void sol(int cas)
{
	sc(n);
	int pre=0;
	L.push(0);
	R.push(0);
	for(int i=0;i<=n;i++) L.push(0),R.push(0);
	rep(i,1,n)
	{
		int t,d,x;sc(t,d,x);
		sl-=t-pre;
		sr+=t-pre;
		pre=t;
		if(d==0)//Xi-x
		{
			if(x<=R.top()+sr)
				L.push(x-sl);
			else
			{
				ans+=x-(R.top()+sr);
				ll y=R.top()+sr;
				R.pop();
				R.push(x-sr);
				L.push(y-sl);
			}
		}
		else//x-Xi
		{
			if(x>=L.top()+sl)
				R.push(x-sr);
			else
			{
				ans+=(L.top()+sl)-x;
				ll y=L.top()+sl;
				L.pop();
				L.push(x-sl);
				R.push(y-sr);
			}
		}
	}
	out(ans);
}
int main()
{
  // freopen("1.in", "r",stdin);
  // freopen("1.out", "w", stdout);
  srand(time(0));
  int t=1,cas=0;
//   scanf("%d",&t);
  while(t--)
  {
    sol(++cas);
  }
  return 0;
}
/*
befor submit code check:
freopen
size of N
mod
debug output
*/