因为两个人方案的对称性,可以将$k$除以$2$,转化为在$n-1$个间隔中设置若干断点,求第$k$小的增量。

对于选中的相邻的断点$(a,a+1)$和$(b,b+1)$,增量为$|x_a-x_{b+1}|$。

将绝对值拆开,用可持久化权值线段树优化建图,然后求$k$短路即可。

时间复杂度$O(n\log^2n+k\log k)$。

#include<cstdio>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
typedef long long ll;
typedef pair<ll,int>P;
const int MAXN=10010,N=400000,M=2000000;
const ll inf=1LL<<60;
int n,K,i,a[MAXN],b[MAXN],T0,T1,l[N],r[N],tot;ll base,ans;
namespace G{
int n,m,i,S,T,g[N],v[M],u[M],w[M],nxt[M],f[N],h[N],tot;ll d[N];bool is[M],vis[N];
struct Node{
int l,r,d;P v;
Node(){}
Node(int _l,int _r,int _d,P _v){l=_l,r=_r,d=_d,v=_v;}
}pool[M*4];
inline int build(P v){
pool[++tot]=Node(0,0,0,v);
return tot;
}
int merge(int a,int b){
if(!a||!b)return a+b;
if(pool[a].v>pool[b].v)swap(a,b);
int x=++tot;
pool[x]=pool[a];
pool[x].r=merge(pool[a].r,b);
if(pool[pool[x].l].d<pool[pool[x].r].d)swap(pool[x].l,pool[x].r);
pool[x].d=pool[x].r?pool[pool[x].r].d+1:0;
return x;
}
void cal(int x){
if(vis[x])return;
vis[x]=1;
d[x]=inf,f[x]=0;
if(x==T){d[x]=0;return;}
for(int i=g[x];i;i=nxt[i]){
cal(u[i]);
if(d[u[i]]+w[i]<d[x])d[x]=d[u[i]]+w[i],f[x]=i;
}
}
void dfs(int x){
if(!f[x]||vis[x])return;
vis[x]=1;
dfs(u[f[x]]);
h[x]=merge(h[x],h[u[f[x]]]);
}
inline void add(int x,int y,int z){
if(!x||!y)return;
v[++m]=x;u[m]=y;w[m]=z;nxt[m]=g[x];g[x]=m;
}
void solve(){
for(i=1;i<=n;i++)cal(i);
ans=d[S];
if(K==1)return;
K--;
for(i=1;i<=n;i++)vis[i]=0,is[f[i]]=1;
for(i=1;i<=m;i++)if(!is[i]&&d[u[i]]<inf)h[v[i]]=merge(h[v[i]],build(P(w[i]-d[v[i]]+d[u[i]],u[i])));
for(i=1;i<=n;i++)dfs(i);
priority_queue<P,vector<P>,greater<P> >q;
ll x,y;
y=h[S];
if(y)q.push(P(d[S]+pool[y].v.first,y));
while(!q.empty()&&K){
K--;
P t=q.top();q.pop();
ans=t.first;
x=t.second,y=pool[x].l;
if(y)q.push(P(ans-pool[x].v.first+pool[y].v.first,y));
y=pool[x].r;
if(y)q.push(P(ans-pool[x].v.first+pool[y].v.first,y));
y=h[pool[x].v.second];
if(y)q.push(P(ans+pool[y].v.first,y));
}
}
}
int ins(int x,int a,int b,int c,int X,int W){
int y=++tot;
if(a==b){
G::add(X,y,W);
return y;
}
int mid=(a+b)>>1;
if(c<=mid)l[y]=ins(l[x],a,mid,c,X,W),r[y]=r[x];
else l[y]=l[x],r[y]=ins(r[x],mid+1,b,c,X,W);
G::add(l[y],y,0),G::add(r[y],y,0);
return y;
}
void ask(int x,int a,int b,int c,int d,int X,int W){
if(!x)return;
if(c<=a&&b<=d){
G::add(x,X,W);
return;
}
int mid=(a+b)>>1;
if(c<=mid)ask(l[x],a,mid,c,d,X,W);
if(d>mid)ask(r[x],mid+1,b,c,d,X,W);
}
int main(){
scanf("%d%d",&n,&K);
K=(K+1)/2;
for(i=1;i<=n;i++){
scanf("%d",&a[i]);
b[i]=a[i];
base+=abs(a[i]-a[i-1]);
}
sort(b,b+n+1);
for(i=0;i<=n;i++)a[i]=lower_bound(b,b+n+1,a[i])-b;
G::S=1;
G::T=2;
tot=2;
G::add(1,2,0);
for(i=1;i<n;i++){
int x=++tot,w=-abs(b[a[i]]-b[a[i+1]]);
G::add(1,x,w+abs(b[a[i+1]]));
G::add(x,2,0);
ask(T0,0,n,0,a[i+1],x,w+b[a[i+1]]);
if(a[i+1]<n)ask(T1,0,n,a[i+1]+1,n,x,w-b[a[i+1]]);
T0=ins(T0,0,n,a[i],x,-b[a[i]]);
T1=ins(T1,0,n,a[i],x,b[a[i]]);
}
G::n=tot;
G::solve();
return printf("%lld",ans+base),0;
}

  

BZOJ3945 : 无聊的邮递员的更多相关文章

  1. 无聊的 邮递员 插头dp

    邮递员想知道,如果他每天都用不同路线走过10×20个点阵邮筒,他必须活过多少个世纪才能走遍所有方案? 7:00 改完T1,开始肝插头dp 7:10 放弃,颓博客 7:20 学习插头dp 7:21 放弃 ...

  2. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  3. 无聊的人用JS实现了一个简单的打地鼠游戏

    直入正题,用JS实现一个简单的打地鼠游戏 因为功能比较简单就直接裸奔JS了,先看看效果图,或者 在线玩玩 吧 如果点击颜色比较深的那个(俗称坏老鼠),将扣分50:如果点击颜色比较浅的那个(俗称好老鼠) ...

  4. [无聊的事-连载之连开100个VS] 折腾你的骚机吧,骚年们!!!

    [无聊的事-连载之连开100个VS] 折腾你的骚机吧,骚年们!!! 只有8G内存的,用这个:打开你的CMD命令提示符如果只装了Visual Studio 2010,用这条命令,打开50个VSfor / ...

  5. 一个无聊的python + opencv 示例

    opencv不用多说,先推荐一个给力的教程: https://opencv-python-tutroals.readthedocs.io/en/latest/py_tutorials/py_tutor ...

  6. Python无聊的总结

    在公司无聊的时候看了前辈写的python代码,突然发现一个比较好玩的python表达式: lambda x,y:x+y 咋一看,这个应该类似方法之类的,上网查了查,所以特此总结下 lambda:上代码 ...

  7. Python的高级特性10:无聊的@property

    @property装饰器其实有点无聊,单独拿出来作为一个知识点其实没必要,尽管它可以将方法变成属性,让get和set方法更好用,但是,它破坏了python的简洁(不是代码的简洁而是指语法上). 下面来 ...

  8. 洛谷P1629 邮递员送信

    题目描述 有一个邮递员要送东西,邮局在节点1.他总共要送N-1样东西,其目的地分别是2~N.由于这个城市的交通比较繁忙,因此所有的道路都是单行的,共有M条道路,通过每条道路需要一定的时间.这个邮递员每 ...

  9. 暑假热身 E. 无聊的LSY

    LSY大牛没事就爱玩游戏,包括很多很无聊的游戏.某日,LSY大牛又找到了一个无聊的游戏:每一局游戏的开始,LSY大牛将代表自己的棋子放在一个线性棋盘的最左端(第0个格子,可以认为向右端无限延伸),接着 ...

随机推荐

  1. Could not create SSL/TLS secure channel.

    解决办法: ServicePointManager.Expect100Continue = true;ServicePointManager.SecurityProtocol = SecurityPr ...

  2. Java+selenium之WebDriver模拟鼠标键盘操作(六)

    org.openqa.selenium.interactions.Actions类,主要定义了一些模拟用户的鼠标mouse,键盘keyboard操作.对于这些操作,使用 perform()方法进行执行 ...

  3. 美国主机BlueHost vs HostEase

    网站备案对于大部分个人站长而言,花费成本高.程序复杂,因此在挑选主机时,常选择免备案IDC服务商,如美国.香港主机,且前几日国内某免备案机房因个别网站涉及非法言论而配合相关部门采取停网整顿等,为此站长 ...

  4. What's news in Visual Studio 2017

    文字总结: 1.高级智能提示  在属性列表中输入 M C即可查询属性中包含字母m\c的属性 2.更快的导航查询,在Go To All中输入任意查询的字符,可快速查到任何包含关键字的文件 3.代码智能分 ...

  5. 用jQuery实现Ajax

    前置知识:ajax原理,json字符串进行信息传递. Ajax主要的功能是实现了浏览器端 异步 访问服务器:通过浏览器的XMLHttpRequest对象发出小部分数据,与服务端进行交互, 服务端返回小 ...

  6. [BZOJ1977][BeiJing2010组队]次小生成树

    题解: 首先要证明一个东西 没有重边的图上 次小生成树由任何一颗最小生成树替换一条边 但是我不会证啊啊啊啊啊啊啊 然后就很简单了 枚举每一条边看看能不能变 但有一个特殊情况就是,他和环上的最大值相等, ...

  7. 【AtCoder】ARC074

    ARC 074 C - Chocolate Bar 直接枚举第一刀横切竖切,然后另一块要求如果横切分成\(H / 2\)竖切分成\(W/2\)即可 #include <bits/stdc++.h ...

  8. Python_tuple部分功能介绍

    x.count():元素在元组内的个数 x.index():元素在元组内的位置

  9. Centos7X部署Zabbix监控

    一:yum安装LAMP环境 zabbix-server端防火墙配置(可以选择iptables -F清空) iptables -A INPUT -m state --state NEW -m tcp - ...

  10. Vim的合并行操作

    日常常用到多行合并的功能,记录如下: 第一种, 多行合并成一行,即: AAAAABBBBBCCCCC 合并为:AAAAA BBBBB CCCCC 方法1: normal状态下 3J 其中的3是范围,可 ...