当前位置: 代码迷 >> 综合 >> 数据备份[APIO/CTSC 2007][WQS二分]
  详细解决方案

数据备份[APIO/CTSC 2007][WQS二分]

热度:59   发布时间:2023-11-19 10:03:21.0

文章目录

  • 题目
  • 思路
  • 代码

题目

Luogu
nnn 个点在 xxx 轴上,位置 pip_ipi?,求配对 kkk 对距离之和最小值
n,k≤105n,k\le 10^5nk105

思路

主要是复习 WQSWQSWQS 二分原理
f(x)f(x)f(x) 为配对 xxx 的最小值
显然答案递增,理解为凸函数
在这里插入图片描述
然后记每次划分会有附加权值 ccc
g(x,c)g(x,c)g(x,c) 为有划分权值前提下划分 xxx 个的最小值
那么有 g(x,c)=f(x)+cxg(x,c)=f(x)+cxg(x,c)=f(x)+cx
因为后面的 cxcxcx 并不会影响决策
(凸函数叠加一次函数还是凸函数,和决策无关)
移项可得:
f(x)=?cx+g(x,c)f(x)=-cx+g(x,c)f(x)=?cx+g(x,c)
画图可得:
在这里插入图片描述
二分即可
现在有个问题:
在这里插入图片描述
来自CXH
也就是共线情况
可以记录相同 g(x,c)g(x,c)g(x,c)xxx 的最小值即可
二分时候区间为 (L,R](L,R](L,R] 保证 RRR 始终合法即可

代码

#include<set>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<cstdio>
#include<climits>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define LL long long
#define ULL unsigned long long
using namespace std;
LL read(){
    LL f=1,x=0;char c=getchar();while(c<'0'||'9'<c){
    if(c=='-')f=-1;c=getchar();}while('0'<=c&&c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();return f*x;
}
#define MAXN 500000
#define INF 0x3f3f3f3f
int n,k;
LL p[MAXN+5],f[MAXN+5],g[MAXN+5];
int check(LL c){
    //f[i]:前i个最小值,g:f前提下最小划分段数memset(f,0x3f,sizeof(f));f[0]=0,f[1]=0;for(int i=2;i<=n;i++){
    f[i]=f[i-1],g[i]=g[i-1];if(f[i]>f[i-2]+p[i]-p[i-1]+c)f[i]=f[i-2]+p[i]-p[i-1]+c,g[i]=g[i-2]+1;else if(f[i]==f[i-2]+p[i]-p[i-1]+c)g[i]=min(g[i],g[i-2]+1);}return g[n];
}
int main(){
    n=read(),k=read();for(int i=1;i<=n;i++)p[i]=read();LL L=-p[n]-1,R=0;while(L+1<R){
    LL Mid=(L+R)>>1;if(check(Mid)<=k)R=Mid;else L=Mid;}check(R);printf("%lld\n",f[n]-k*R);return 0;
}
  相关解决方案