【BZOJ4167】永远的竹笋采摘 分块+树状数组
【BZOJ4167】永远的竹笋采摘
题解:我们考虑有多少点对(a,b)满足a与b的差值是[a,b]中最小的。以为是随机数据,这样的点对数目可能很少,实测是O(n)级别的,那么我们已知了有这么多可能对答案造成贡献的点对,如何将它们求出来呢?
考虑分块,因为所有数大小在[1,n]中,我们可以对于每个块,预处理出整个块到所有数的最小差值。然后从右往左枚举每一个点,再枚举右面所有的块,如果这个块到当前数的差值比之前的要小,那就暴力进入块中扫一遍。与此同时,我们需要知道是否已经存在这样的点对,被当前的点对完全包含且差值更小。这个可以用树状数组搞定。
最后,我们得到所有的点对,问题就变成了在数轴上选取k个互不相交的线段,使得线段权值和最小。跑个DP就行了。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=60010;
int n,m,B,cnt,minn;
int s[maxn],v[maxn],p[maxn];
int cls[250][maxn],f[2][maxn],to[500000],next[500000],head[maxn],val[500000];
int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-')f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar();
return ret*f;
}
int z(int x)
{
return x>0?x:-x;
}
void updata(int x,int val)
{
for(int i=x;i<=n;i+=i&-i) s[i]=min(s[i],val);
}
int query(int x)
{
int i,ret=1<<30;
for(i=x;i;i-=i&-i) ret=min(ret,s[i]);
return ret;
}
void test(int a,int b)
{
int c=z(v[b]-v[a]);
a++,b++,minn=min(minn,c);
if(query(b)<=c) return ;
updata(b,c);
to[cnt]=a,val[cnt]=c,next[cnt]=head[b],head[b]=cnt++;
}
int main()
{
//freopen("bz4168.in","r",stdin);
n=rd(),m=rd(),B=ceil(sqrt(n));
int i,j,k,last;
for(i=0;i<n;i++) v[i]=rd();
memset(cls,0x3f,sizeof(cls));
memset(s,0x3f,sizeof(s));
for(i=0;i<n;i+=B)
{
for(j=i;j<i+B&&j<n;j++) p[v[j]]=1;
for(last=-1<<30,j=1;j<=n;j++) cls[i/B][j]=min(cls[i/B][j],j-last),last=p[j]?j:last;
for(last=1<<30,j=n;j>=1;j--) cls[i/B][j]=min(cls[i/B][j],last-j),last=p[j]?j:last;
for(j=i;j<i+B&&j<n;j++) p[v[j]]=0;
}
memset(head,-1,sizeof(head));
for(i=n-1;i>=0;i--)
{
minn=1<<30;
for(j=i+1;j<i/B*B+B&&j<n;j++) if(v[j]!=v[i]&&z(v[j]-v[i])<minn) test(i,j);
for(j=i/B+1;j*B<n;j++) if(cls[j][v[i]]<minn) for(k=j*B;k<j*B+B&&k<n;k++) if(v[k]!=v[i]&&z(v[k]-v[i])<minn) test(i,k);
}
for(k=1;k<=m;k++)
{
for(i=0;i<=n;i++) f[k&1][i]=1<<30;
for(i=1;i<=n;i++)
{
f[k&1][i]=f[k&1][i-1];
for(j=head[i];j!=-1;j=next[j]) f[k&1][i]=min(f[(k&1)^1][to[j]-1]+val[j],f[k&1][i]);
}
}
printf("%d\n",f[m&1][n]);
return 0;
}
【BZOJ4167】永远的竹笋采摘 分块+树状数组的更多相关文章
- 【BZOJ 3295】动态逆序对 - 分块+树状数组
题目描述 给定一个1~n的序列,然后m次删除元素,每次删除之前询问逆序对的个数. 分析:分块+树状数组 (PS:本题的CDQ分治解法见下一篇) 首先将序列分成T块,每一块开一个树状数组,并且先把最初的 ...
- 【bzoj2141】排队 分块+树状数组
题目描述 排排坐,吃果果,生果甜嗦嗦,大家笑呵呵.你一个,我一个,大的分给你,小的留给我,吃完果果唱支歌,大家乐和和.红星幼儿园的小朋友们排起了长长地队伍,准备吃果果.不过因为小朋友们的身高有所区别, ...
- 【bzoj3744】Gty的妹子序列 分块+树状数组+主席树
题目描述 我早已习惯你不在身边, 人间四月天 寂寞断了弦. 回望身后蓝天, 跟再见说再见…… 某天,蒟蒻Autumn发现了从 Gty的妹子树(bzoj3720) 上掉落下来了许多妹子,他发现 她们排成 ...
- 【分块+树状数组】codechef November Challenge 2014 .Chef and Churu
https://www.codechef.com/problems/FNCS [题意] [思路] 把n个函数分成√n块,预处理出每块中各个点(n个)被块中函数(√n个)覆盖的次数 查询时求前缀和,对于 ...
- Bzoj 3295: [Cqoi2011]动态逆序对 分块,树状数组,逆序对
3295: [Cqoi2011]动态逆序对 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2886 Solved: 924[Submit][Stat ...
- 【XSY2111】Chef and Churus 分块 树状数组
题目描述 有一个长度为\(n\)的数组\(A\)和\(n\)个区间\([l_i,r_i]\),有\(q\)次操作: \(1~x~y\):把\(a_x\)改成\(y\) \(2~x~y\):求第\(l\ ...
- BZOJ3787:Gty的文艺妹子序列(分块,树状数组)
Description Autumn终于会求区间逆序对了!Bakser神犇决定再考验一下他,他说道: “在Gty的妹子序列里,某个妹子的美丽度可也是会变化的呢.你还能求出某个区间中妹子们美丽度的逆序对 ...
- 2018.06.30 BZOJ4765: 普通计算姬(dfs序+分块+树状数组)
4765: 普通计算姬 Time Limit: 30 Sec Memory Limit: 256 MB Description "奋战三星期,造台计算机".小G响应号召,花了三小时 ...
- 【xsy2111】 【CODECHEF】Chef and Churus 分块+树状数组
题目大意:给你一个长度为$n$的数列$a_i$,定义$f_i=\sum_{j=l_i}^{r_i} num_j$. 有$m$个操作: 操作1:询问一个区间$l,r$请你求出$\sum_{i=l}^{r ...
随机推荐
- LeetCode OJ-- Reverse Words in a String
https://oj.leetcode.com/problems/reverse-words-in-a-string/ 给一个字符串 abc dd m,返回 m dd abc. 注意:输入中可能有 ...
- 转载——Java与WCF交互(二):WCF客户端调用Java Web Service
在上篇< Java与WCF交互(一):Java客户端调用WCF服务>中,我介绍了自己如何使用axis2生成java客户端的悲惨经历.有同学问起使用什么协议,经初步验证,发现只有wsHttp ...
- 更新到xcode10以后出现几个无奈的问题,谨已此篇告诫广大ioser升级请慎重
1.第一次用xcode 10 archive的时候遇到的电脑卡死不动的问题,期间鼠标键盘通通都动不了,只能强制关机来解决,于是又进行了一次可还是遇到相同的问题,无奈之下只能等待,大约20分钟左右(20 ...
- 关于超大binlog事件的问题
我手里维护了一个项目,其功能是用Java模拟一个MariaDB的slave库连接到主库,对从主库传输过来的binlog事件进行监听与分析 碰到一个问题是: 如果主库做了一个很大的修改操作(比方说直接d ...
- (21)Oracle表查询进阶
转到基本查询 一.多表查询 笛卡尔积:每张表的列数相加,行数相乘. 连接条件:得出笛卡尔积后需要用where条件筛选出正确的数据.连接条件至少需要n张表减1个 1.等值连接 连接条件为等号 selec ...
- 洛谷 P1014 Cantor表【蛇皮矩阵/找规律/模拟】
题目描述 现代数学的著名证明之一是Georg Cantor证明了有理数是可枚举的.他是用下面这一张表来证明这一命题的: 1/1 1/2 1/3 1/4 1/5 … 2/1 2/2 2/3 2/4 … ...
- Codeforces 482B Interesting Array(线段树区间更新)
题目链接 Interesting Array 区间更新.然后对于每一个约数重新求一遍区间的&值,不符合就跳出. #include <bits/stdc++.h> using nam ...
- smartsvn学习(-)
3. Day-to-day use 3.1 Everyday commands Most of the common commands you will need are in the SmartSV ...
- c# winform 最小化到托盘
STEP1.添加托盘图标控件NotifyIcon(直接从工具箱中拖动添加即可) STEP2.添加(重写)窗口尺寸变动函数Form1_Resize private void Main_SizeChang ...
- 导出excel(利用工具类导出excel)
/** * 添加导出功能 * @param creditPageResult * @param request * @param response */ @RequestMapping(value = ...