【BZOJ】2726: [SDOI2012]任务安排-斜率优化&CDQ二分

传送门：bzoj2726

---

题解

斜率优化。
按照“费用提前计算”的规则将每个任务的$S$花费提前处理，列出$dp$转移：
$d{p}_i=d{p}_j-(S+T_i)C_j+T_i{C}_i+S{C}_n$
其中$C_i=\sum _{j=1}^i{F}_j,T_i=\sum _{j=1}^i{T}_j$

转成斜率的形式：
$F_j=(S+T_i)C_j-T_i{C}_i-S{C}_n+F_i$

因为$C_j$是单调的，其实可以在凸包上二分，cdq分治也可以做。

---

代码

#include
using namespace std;
const int N=3e5+100;
typedef long long ll;
typedef double db;

int n,S,top,stk[N];
int c[N],t[N];ll f[N];

struct P{
	int x;ll y;int k;ll b;int id;
	bool operator<(const P&ky)const{
	    return k
inline void rd(yyy &x)
{
	cp=getchar();int pr=0;x=0;
	for(;!isdigit(cp);cp=getchar()) if(cp=='-') pr=1;
	for(;isdigit(cp);cp=getchar()) x=(x<<3)+(x<<1)+(cp^48);
	if(pr) x=-x;
}

inline db slp(int A,int B)
{
	if(p[A].x==p[B].x) return (p[B].y>=p[A].y)?(1e18):(-1e18);
	return (p[B].y-p[A].y)/(db)(p[B].x-p[A].x);
}

void sol(int l,int r)
{
	if(l==r) {p[l].y=f[l];return;}
	int mid=(l+r)>>1,k,i=l,j=mid+1;
	for(k=l;k<=r;++k)
	 if(p[k].id<=mid) rep[i++]=p[k];
	 else rep[j++]=p[k];
	for(k=l;k<=r;++k) p[k]=rep[k];
	sol(l,mid);
    top=0;
    for(k=l;k<=mid;++k){
    	for(;top>1 && slp(stk[top-1],stk[top])>=slp(stk[top],k);--top);
    	stk[++top]=k;
    }
    for(k=mid+1,i=1;k<=r;++k){
    	for(;ir) || (i<=mid && p[i].x