1049: [HAOI2006]数字序列

Description

现在我们有一个长度为n的整数序列A。但是它太不好看了，于是我们希望把它变成一个单调严格上升的序列。但是不希望改变过多的数，也不希望改变的幅度太大。
Input

第一行包含一个数n，接下来n个整数按顺序描述每一项的键值。
Output

第一行一个整数表示最少需要改变多少个数。第二行一个整数，表示在改变的数最少的情况下，每个数改变的绝对值之和的最小值。
Sample Input
4
5 2 3 5
Sample Output
1
4

【数据范围】
90%的数据n<=6000。
100%的数据n<=35000。
保证所有数列是随机的。

这个算法的复杂度真的是O（n^2）吗？怎么我觉得不对呢

和机房的小伙伴讨论之后，觉得这个复杂度应该是O（n^3）的（虽然跑得很快）

先转换一下，b[i]=a[i]-i，然后就是求b的最长不下降了，这样好算代价一些

cost[i]=cost[j]+w(j,i)(f[j]=f[i]-1,j<i)

计算w的时候有一个定理，肯定有一种最优方案是把左半边变成b[j]右半边变成b[i]

证明：http://pan.baidu.com/share/link?uk=2651016602&shareid=1490516411

 const

     maxn=;

     inf=;

 var

     a,f,b,first,next,last:array[..maxn]of longint;

     g:array[..maxn]of int64;

     n,tot,max:longint;

 procedure insert(x,y:longint);

 begin

     inc(tot);

     last[tot]:=y;

     next[tot]:=first[x];

     first[x]:=tot;

 end;

 procedure find(x:longint);

 var

     l,r,mid:longint;

 begin

     l:=;

     r:=n;

     while l<>r do

       begin

         mid:=(l+r)>>;

         if b[mid]>a[x] then r:=mid

         else l:=mid+;

       end;

     if max<l then max:=l;

     f[x]:=l;

     b[l]:=a[x];

     insert(l,x);

 end;

 procedure init;

 var

     i:longint;

 begin

     read(n);

     for i:= to n do

       begin

         read(a[i]);

         dec(a[i],i);

       end;

     inc(n);

     a[n]:=inf-;

     a[]:=-inf;

     for i:= to n do

       b[i]:=inf;

     insert(,);

 end;

 function calc(l,r:longint):int64;

 var

     i:longint;

     s:int64;

 begin

     s:=;

     for i:=l to r do

       inc(s,abs(a[i]-a[r]));

     calc:=s;

     for i:=l to r do

       begin

         s:=s-abs(a[i]-a[r])+abs(a[i]-a[l]);

         if calc>s then calc:=s;

       end;

 end;

 procedure work;

 var

     i,j:longint;

     num:int64;

 begin

     for i:= to n do

       begin

         find(i);

         g[i]:=inf*inf;

         j:=first[f[i]-];

         while j<> do

           begin

             if a[last[j]]<=a[i] then

               begin

                 num:=calc(last[j],i);

                 if g[i]>g[last[j]]+num then g[i]:=g[last[j]]+num;

               end;

             j:=next[j];

           end;

       end;

     writeln(n-max);

     write(g[n]);

 end;

 begin

     init;

     work;

 end.

pascal第三（前两名是怎么回事......）

1049: [HAOI2006]数字序列 - BZOJ的更多相关文章

【BZOJ 1049】 1049: [HAOI2006]数字序列（LIS+动态规划）
1049: [HAOI2006]数字序列 Description 现在我们有一个长度为n的整数序列A.但是它太不好看了,于是我们希望把它变成一个单调严格上升的序列.但是不希望改变过多的数,也不希望改变 ...
bzoj 1049 [HAOI2006]数字序列
[bzoj1049][HAOI2006]数字序列 Description 现在我们有一个长度为n的整数序列A.但是它太不好看了,于是我们希望把它变成一个单调严格上升的序列.但是不希望改变过多的数,也不 ...
【BZOJ】1049: [HAOI2006]数字序列（lis+特殊的技巧）
http://www.lydsy.com/JudgeOnline/problem.php?id=1049 题意:给一个长度为n的整数序列.把它变成一个单调严格上升的序列.但是不希望改变过多的数,也不希 ...
bzoj 1049: [HAOI2006]数字序列【dp+二分+瞎搞】
第一问明显就是用b[i]=a[i]-i来做最长不下降子序列然后第二问,对于一对f[i]=f[j]+1的(i,j),中间的数一定要改的,并且是等于b[i]或者b[j],我不会证,然后因为是随机数据,所 ...
【BZOJ1049】 [HAOI2006]数字序列
BZOJ1049 [HAOI2006]数字序列 dp好题? 第一问第一问我会做!令\(b_i=a_i-i\),求一个最长不下降子序列. \(n-ans\)就是最终的答案. 第二问好难啊.不会.挖坑 ...
洛谷 P2501 [HAOI2006]数字序列解题报告
P2501 [HAOI2006]数字序列题目描述现在我们有一个长度为n的整数序列A.但是它太不好看了,于是我们希望把它变成一个单调严格上升的序列.但是不希望改变过多的数,也不希望改变的幅度太大. ...
洛谷P2501 bzoj1049 [HAOI2006]数字序列
题目链接 bzoj 洛谷题解第一问: 假如 \(i < j\) 如果 \(j\)能从\(i\)转移过来显然中间空隙必须足够例如:\(50\) \(53\) \(53\) \(52\) 就 ...
BZOJ1049 [HAOI2006]数字序列0
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000作者博客:http://www.cnblogs.com/ljh2000-jump/转 ...
BZOJ1049: [HAOI2006]数字序列
题目:http://www.lydsy.com/JudgeOnline/problem.php?id=1049 题解: ydc的题解:http://pan.baidu.com/share/link?u ...

随机推荐

apk瘦身（包大小优化）
最近太忙好久没有写随笔,放假才有空写写随笔. 最近在项目中一直在做包大小优化随着业务需求的增加,包大小是忍不住的往上涨为了提高应用市场下载转化率,一直在优化包大小首先分析你的apk中占比的主要 ...
LeetCode16 3Sum Closest
题意: Given an array S of n integers, find three integers in S such that the sum is closest to a given ...
hdu 4424 并查集
思路:将边从大到小排序,判断向哪边连,能使总和最大. #include<map> #include<set> #include<cmath> #include< ...
js 定义函数的几种方法以及如何调用
/*1.方法调用模式: 先定义一个对象,然后在对象的属性中定义方法,通过myobject.property来执行方法,this即指当前的myobject 对象.*/ var car = { carId ...
oracle 数据库分割字符串返回结果集函数
CREATE OR REPLACE FUNCTION "UFN_SPLIT" ( p_list varchar2, p_sep varchar2 := ',' ...
java如何获得数据库表中各字段的字段名
public class TestDemo { public static Connection getConnection() { Connection conn = null; try { Cla ...
springboot2.0添加logback
程序代码 import org.slf4j.Logger; import org.slf4j.LoggerFactory; public class PreFilter extends ZuulFil ...
Win7系统中如何查看当前文件被哪一个程序占用了
https://superuser.com/questions/117902/find-out-which-process-is-locking-a-file-or-folder-in-windows ...
CentOS 文件及目录等
1.在linux中一切皆是文件,只是类型不同,通过ls -l看到的一个字母表示文件的类型 -:普通文件. d:目录文件. l:链接文件. b:块设备文件. c:字符设备文件. p:管道文件. 2.文件 ...
[sj系统] phabricator系统升级
sj是公司产品bug管理 wiki一体化系统 ,很刁. 安装: http://www.cnblogs.com/chorulex/p/5381558.html 升级: https://phabricat ...

1049: [HAOI2006]数字序列 - BZOJ

1049: [HAOI2006]数字序列 - BZOJ的更多相关文章

随机推荐

热门专题