Faker_fan
Faker_fan

hdu6800

航电多校第三场

　　这道题的思维很是巧妙啊，虽然赛时大概有了想法，但是当初没学会树套树，还是过不了的啊，感觉题解有点太难懂了，再一篇博客，分享一下自己标注过的std；

　　首先是这题一共只有两只手，而且按顺序跑n个点，也就是说，假如dp[i][j]代表着你左手在i，右手在j的时候的最优解，其实dp[i][j]=dp[j][i];所以说当最大的点为i时，你只需要关心另外一只手放在哪个点了就行了，也就是说dp[i+1][j]只和dp[i][*]有关系，和dp[i-1][*]只有间接关系。

　　然后为了dp好弄一点，dp存的值不是当前答案时多少，而是我们通过交换手省下了多少的步长。

　　首先除了dp[i][i-1]以外，dp[i][j]都是等于dp[i-1][j]的，为什么呢，因为从dp[i-1][j]到dp[i][j]，我们不会换手，所以说我们不会省下任何步长，所以我们每跑一个点，只需要更新dp[i][i-1]即可。

　　然后dp[i][i-1]怎么更新呢，就变成了二维平面上，找一个j点与这个点的曼哈顿距离加上dp[i][j]最小，这个可以通过树套树loglog时间解决

#include 
#include <string.h>
#include 
#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;
typedef long long LL;
const int maxn = (int)1e5 + 1;

namespace Memory {
    const int maxm = maxn * 18 + 1, maxl = maxm * 8 + 1;
    int pool[maxm], *pool_tail;
    LL info[maxl], *info_tail;

    inline void reset() {
        pool_tail = pool;
        info_tail = info;
    }
    inline int *ask32(int len) {
        int *ret = pool_tail;
        pool_tail += len;
        return ret;
    }
    inline LL *ask64(int len) {
        LL *ret = info_tail;
        info_tail += len;
        return ret;
    }
}

int n, qtot;
LL dp[maxn];
pair<int, int> seq[maxn], que[maxn];

inline void upd_min(LL &x, LL y) {
    x > y && (x = y);
}

inline bool cmp_y(int const &u, int const &v) {
    return que[u].second < que[v].second;
}

struct Segment {
    int ytot, *yque;
    LL *info;
} seg[maxn << 1 | 1];

inline int seg_idx(int L, int R) {//线段树L-R的hash映射
    return (L + R) | (L < R);
}


void seg_init_inner(Segment &rt, int L, int R) {
    if(L == R)
        return;
    int M = (L + R) >> 1;
    seg_init_inner(rt, L, M);
    seg_init_inner(rt, M + 1, R);
    LL *cur = rt.info + (seg_idx(L, R) << 2);
    LL *lft = rt.info + (seg_idx(L, M) << 2);
    LL *rht = rt.info + (seg_idx(M + 1, R) << 2);
    cur[0] = min(lft[0], rht[0]);
    cur[1] = min(lft[1], rht[1]);
    cur[2] = min(lft[2], rht[2]);
    cur[3] = min(lft[3], rht[3]);
}

void seg_init_outer(int L, int R) {
    
    static int ord[maxn];
    
    if(L < R) {
        int M = (L + R) >> 1;
        seg_init_outer(L, M);
        seg_init_outer(M + 1, R);
        inplace_merge(ord + L, ord + M + 1, ord + R + 1, cmp_y);//按que[u].second 有序合并
    }
    else
    {
        ord[L] = L;
    }
    
    Segment &rt = seg[seg_idx(L, R)];
    
    rt.ytot = 1;
    for(int i = L + 1; i <= R; ++i)
        rt.ytot += cmp_y(ord[i - 1], ord[i]);//L到R中到底有多少个不同的y值
    
    rt.yque = Memory::ask32(rt.ytot);
    rt.info = Memory::ask64(((rt.ytot - 1) << 1 | 1) << 2); ((n-1)*2+1)*4;
    
    
    for(int i = L, j = 0; i <= R; ++j) {
        int u = ord[i++], y = que[u].second, xL = que[u].first, xR = xL;
        
        
        for(int v; i <= R && !cmp_y(u, v = ord[i]); ++i) { //如果u+1和u同高，继续u++
            int tmp = que[v].first;
            if(tmp < xL) {
                xL = tmp;
            } else if(tmp > xR) {
                xR = tmp;
            }
            //扩充当前的左右边界
        }
        //这里的j是按高度不同才增一次的
        rt.yque[j] = y;
        
        LL *val = rt.info + (seg_idx(j, j) << 2);//4*j
        val[0] = -xR - y;//各个象限对应的值
        val[1] = -xR + y;
        val[2] = xL - y;
        val[3] = xL + y;
    }
    //内层的线段树
    //内层的线段树
    seg_init_inner(rt, 0, rt.ytot - 1);
}

inline void seg_update_inner(Segment &rt, int y, LL val[]) {
    
    for(int L = 0, R = rt.ytot - 1; L <= R; ) {
        LL *cur = rt.info + (seg_idx(L, R) << 2);
        upd_min(cur[0], val[0]);
        upd_min(cur[1], val[1]);
        upd_min(cur[2], val[2]);
        upd_min(cur[3], val[3]);
        if(L == R)
            break;
        
        int M = (L + R) >> 1;
        
        if(y <= rt.yque[M]) {
            R = M;
        } else {
            L = M + 1;
        }
        
        
    }
    
}

inline void seg_update_outer(int pos) {
    int x = que[pos].first, y = que[pos].second;
    LL v = dp[pos], val[5] = {
        v - x - y,
        v - x + y,
        v + x - y,
        v + x + y
    };
    for(int L = 0, R = qtot - 1; L <= R; ) {
        
        seg_update_inner(seg[seg_idx(L, R)], y, val);
        
        if(L == R)
            break;
        
        int M = (L + R) >> 1;
        
        if(pos <= M) {
            
            R = M;
            
        } else {
            L = M + 1;
        }
    }
}

inline LL seg_query_left_inner(Segment &rt, int x, int y) {
    LL ret = LLONG_MAX;
    for(int L = 0, R = rt.ytot - 1; L <= R; ) {
        LL *cur = rt.info + (seg_idx(L, R) << 2);
        if(y <= rt.yque[L]) {
            upd_min(ret, cur[1] - y);
            break;
        } else if(y >= rt.yque[R]) {
            upd_min(ret, cur[0] + y);
            break;
        }
        
        int M = (L + R) >> 1;
        if(y <= rt.yque[M]) {
            LL *rht = rt.info + (seg_idx(M + 1, R) << 2);
            upd_min(ret, rht[1] - y);
            R = M;
        } else {
            LL *lft = rt.info + (seg_idx(L, M) << 2);
            upd_min(ret, lft[0] + y);
            L = M + 1;
        }
        
    }
    if(ret != LLONG_MAX)
        ret += x;
    return ret;
}

inline LL seg_query_right_inner(Segment &rt, int x, int y) {
    LL ret = LLONG_MAX;
    for(int L = 0, R = rt.ytot - 1; L <= R; ) {
        LL *cur = rt.info + (seg_idx(L, R) << 2);
        if(y <= rt.yque[L]) {
            upd_min(ret, cur[3] - y);
            break;
        } else if(y >= rt.yque[R]) {
            upd_min(ret, cur[2] + y);
            break;
        }//已经可以确定在哪个象限有答案了
        
        int M = (L + R) >> 1;
        
        if(y <= rt.yque[M]) {
            
            LL *rht = rt.info + (seg_idx(M + 1, R) << 2);
            //在右上方的点的贡献也可以算了
            upd_min(ret, rht[3] - y);
            R = M;//从新区间L-M再来
        } else {
            LL *lft = rt.info + (seg_idx(L, M) << 2);
            upd_min(ret, lft[2] + y);
            L = M + 1;
        }
    }
    if(ret != LLONG_MAX)
        ret -= x;
    return ret;
}

inline LL seg_query_outer(int pos) {
    int x = que[pos].first, y = que[pos].second;
    LL ret = dp[pos];
    for(int L = 0, R = qtot - 1; L < R; ) {
        int M = (L + R) >> 1;
        
        if(pos <= M) {
            upd_min(ret, seg_query_right_inner(seg[seg_idx(M + 1, R)], x, y));
            R = M;
        } else {
            upd_min(ret, seg_query_left_inner(seg[seg_idx(L, M)], x, y));
            L = M + 1;
        }
        
    }
    
    return ret;
}

void solve() {
    scanf("%d", &n);
    for(int i = 0; i < n; ++i) {
        scanf("%d%d", &seq[i].first, &seq[i].second);
        que[i] = seq[i];
        dp[i] = 0;
    }
    sort(que, que + n);
    qtot = unique(que, que + n) - que;

    Memory::reset();
    seg_init_outer(0, qtot - 1);//线段树初始化
    

    LL ans = 0, sum = 0;
        
    int las = lower_bound(que, que + qtot, seq[0]) - que;
    
    for(int i = 1; i < n; ++i) {
        LL adt = abs(seq[i].first - seq[i - 1].first) + abs(seq[i].second - seq[i - 1].second);
        int pos = lower_bound(que, que + qtot, seq[i]) - que;
        LL best = seg_query_outer(pos) - adt;//以换手的代价跑到pos，省了多少的代价
        
        if(dp[las] > best) {//比已经省下的多
            
            ans=min(ans , dp[las] = best );
            
            seg_update_outer(las);
        }
        
        sum += adt;
        las = pos;
    }
    printf("%lld\n", ans + sum);
}
int main() {
    int T;
    scanf("%d", &T);
    for(int Case = 1; Case <= T; ++Case) {
        // printf("Case #%d:\n", Case);
        solve();
    }
    return 0;
}

 

你可能感兴趣的:(hdu6800)

  • hdu6800 Faker_fan
    航电多校第三场这道题的思维很是巧妙啊，虽然赛时大概有了想法，但是当初没学会树套树，还是过不了的啊，感觉题解有点太难懂了，再一篇博客，分享一下自己标注过的std；首先是这题一共只有两只手，而且按顺序跑n个点，也就是说，假如dp[i][j]代表着你左手在i，右手在j的时候的最优解，其实dp[i][j]=dp[j][i];所以说当最大的点为i时，你只需要关心另外一只手放在哪个点了就行了，也就是说dp[i
  • LeetCode[Math] - #66 Plus One Cwind javaLeetCode题解AlgorithmMath
    原题链接：#66 Plus One   要求： 给定一个用数字数组表示的非负整数，如num1 = {1, 2, 3, 9}, num2 = {9, 9}等，给这个数加上1。 注意： 1. 数字的较高位存在数组的头上，即num1表示数字1239 2. 每一位（数组中的每个元素）的取值范围为0~9   难度：简单   分析： 题目比较简单，只须从数组
  • JQuery中$.ajax()方法参数详解 AILIKES JavaScriptjsonpjqueryAjaxjson
    url: 要求为String类型的参数，（默认为当前页地址）发送请求的地址。 type: 要求为String类型的参数，请求方式（post或get）默认为get。注意其他http请求方法，例如put和    delete也可以使用，但仅部分浏览器支持。 timeout: 要求为Number类型的参数，设置请求超时时间（毫秒）。此设置将覆盖$.ajaxSetup()方法的全局
  • JConsole & JVisualVM远程监视Webphere服务器JVM Kai_Ge JVisualVMJConsoleWebphere
        JConsole是JDK里自带的一个工具，可以监测Java程序运行时所有对象的申请、释放等动作，将内存管理的所有信息进行统计、分析、可视化。我们可以根据这些信息判断程序是否有内存泄漏问题。 　　使用JConsole工具来分析WAS的JVM问题，需要进行相关的配置。 　　首先我们看WAS服务器端的配置. 　　1、登录was控制台https://10.4.119.18
  • 自定义annotation 120153216 annotation
    Java annotation 自定义注释@interface的用法 一、什么是注释      说起注释，得先提一提什么是元数据(metadata)。所谓元数据就是数据的数据。也就是说，元数据是描述数据的。就象数据表中的字段一样，每个字段描述了这个字段下的数据的含义。而J2SE5.0中提供的注释就是java源代码的元数据，也就是说注释是描述java源
  • CentOS 5/6.X 使用 EPEL YUM源 2002wmj centos
    CentOS 6.X 安装使用EPEL YUM源1. 查看操作系统版本[root@node1 ~]# uname -a Linux node1.test.com 2.6.32-358.el6.x86_64 #1 SMP Fri Feb 22 00:31:26 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux [root@node1 ~]#
  • 在SQLSERVER中查找缺失和无用的索引SQL 357029540 SQL Server
    --缺失的索引 SELECT  avg_total_user_cost * avg_user_impact * ( user_scans + user_seeks ) AS PossibleImprovement ,          last_user_seek ,    
  • Spring3 MVC 笔记（二） —json+rest优化 7454103 Spring3 MVC
    接上次的 spring mvc 注解的一些详细信息！                          其实也是一些个人的学习笔记  呵呵！
  • 替换“\”的时候报错Unexpected internal error near index 1 \ ^ adminjun java“\替换”
    发现还是有些东西没有刻子脑子里,,过段时间就没什么概念了,所以贴出来...以免再忘...   在拆分字符串时遇到通过 \ 来拆分，可是用所以想通过转义 \\ 来拆分的时候会报异常   public class Main {          /*
  • POJ 1035 Spell checker(哈希表) aijuans 暴力求解--哈希表
    /* 题意：输入字典，然后输入单词，判断字典中是否出现过该单词，或者是否进行删除、添加、替换操作，如果是，则输出对应的字典中的单词 要求按照输入时候的排名输出 题解：建立两个哈希表。一个存储字典和输入字典中单词的排名，一个进行最后输出的判重 */ #include <iostream> //#define using namespace std; const int HASH =
  • 通过原型实现javascript Array的去重、最大值和最小值 ayaoxinchao JavaScriptarrayprototype
    用原型函数（prototype）可以定义一些很方便的自定义函数，实现各种自定义功能。本次主要是实现了Array的去重、获取最大值和最小值。 实现代码如下：   <script type="text/javascript"> Array.prototype.unique = function() { var a = {}; var le
  • UIWebView实现https双向认证请求 bewithme UIWebViewhttpsObjective-C
              什么是HTTPS双向认证我已在先前的博文 ASIHTTPRequest实现https双向认证请求 中有讲述，不理解的读者可以先复习一下。本文是用UIWebView来实现对需要客户端证书验证的服务请求，网上有些文章中有涉及到此内容，但都只言片语，没有讲完全，更没有完整的代码，让人困扰不已。但是此知
  • NoSQL数据库之Redis数据库管理(Redis高级应用之事务处理、持久化操作、pub_sub、虚拟内存) bijian1013 redis数据库NoSQL
    3.事务处理         Redis对事务的支持目前不比较简单。Redis只能保证一个client发起的事务中的命令可以连续的执行，而中间不会插入其他client的命令。当一个client在一个连接中发出multi命令时，这个连接会进入一个事务上下文，该连接后续的命令不会立即执行，而是先放到一个队列中，当执行exec命令时，redis会顺序的执行队列中
  • 各数据库分页sql备忘 bingyingao oraclesql分页
    ORACLE 下面这个效率很低 SELECT * FROM ( SELECT A.*, ROWNUM RN FROM (SELECT * FROM IPAY_RCD_FS_RETURN order by id desc) A ) WHERE RN <20; 下面这个效率很高 SELECT A.*, ROWNUM RN FROM (SELECT * FROM IPAY_RCD_
  • 【Scala七】Scala核心一：函数 bit1129 scala
    1. 如果函数体只有一行代码，则可以不用写{},比如 def print(x: Int) = println(x) 一行上的多条语句用分号隔开，则只有第一句属于方法体，例如   def printWithValue(x: Int) : String= println(x); "ABC"   上面的代码报错，因为，printWithValue的方法
  • 了解GHC的factorial编译过程 bookjovi haskell
    GHC相对其他主流语言的编译器或解释器还是比较复杂的，一部分原因是haskell本身的设计就不易于实现compiler，如lazy特性，static typed，类型推导等。 关于GHC的内部实现有篇文章说的挺好，这里，文中在RTS一节中详细说了haskell的concurrent实现，里面提到了green thread，如果熟悉Go语言的话就会发现，ghc的concurrent实现和Go有点类
  • Java-Collections Framework学习与总结-LinkedHashMap BrokenDreams LinkedHashMap
            前面总结了java.util.HashMap，了解了其内部由散列表实现，每个桶内是一个单向链表。那有没有双向链表的实现呢？双向链表的实现会具备什么特性呢？来看一下HashMap的一个子类——java.util.LinkedHashMap。       
  • 读《研磨设计模式》-代码笔记-抽象工厂模式-Abstract Factory bylijinnan abstract
    声明： 本文只为方便我个人查阅和理解，详细的分析以及源代码请移步 原作者的博客http://chjavach.iteye.com/ package design.pattern; /* * Abstract Factory Pattern * 抽象工厂模式的目的是： * 通过在抽象工厂里面定义一组产品接口，方便地切换“产品簇” * 这些接口是相关或者相依赖的
  • 压暗面部高光 cherishLC PS
    方法一、压暗高光&重新着色 当皮肤很油又使用闪光灯时，很容易在面部形成高光区域。 下面讲一下我今天处理高光区域的心得： 皮肤可以分为纹理和色彩两个属性。其中纹理主要由亮度通道（Lab模式的L通道）决定，色彩则由a、b通道确定。 处理思路为在保持高光区域纹理的情况下，对高光区域着色。具体步骤为：降低高光区域的整体的亮度，再进行着色。 如果想简化步骤，可以只进行着色（参看下面的步骤1
  • Java VisualVM监控远程JVM crabdave visualvm
    Java VisualVM监控远程JVM    JDK1.6开始自带的VisualVM就是不错的监控工具. 这个工具就在JAVA_HOME\bin\目录下的jvisualvm.exe, 双击这个文件就能看到界面   通过JMX连接远程机器, 需要经过下面的配置: 1. 修改远程机器JDK配置文件 (我这里远程机器是linux).    
  • Saiku去掉登录模块 daizj saiku登录olapBI
    1、修改applicationContext-saiku-webapp.xml <security:intercept-url pattern="/rest/**" access="IS_AUTHENTICATED_ANONYMOUSLY" />  <security:intercept-url pattern=&qu
  • 浅析 Flex中的Focus dsjt htmlFlexFlash
    关键字：focus、 setFocus、 IFocusManager、KeyboardEvent 焦点、设置焦点、获得焦点、键盘事件 一、无焦点的困扰——组件监听不到键盘事件 原因：只有获得焦点的组件（确切说是InteractiveObject）才能监听到键盘事件的目标阶段；键盘事件（flash.events.KeyboardEvent）参与冒泡阶段，所以焦点组件的父项（以及它爸
  • Yii全局函数使用 dcj3sjt126com yii
    由于YII致力于完美的整合第三方库，它并没有定义任何全局函数。yii中的每一个应用都需要全类别和对象范围。例如，Yii::app()->user;Yii::app()->params['name'];等等。我们可以自行设定全局函数，使得代码看起来更加简洁易用。(原文地址) 我们可以保存在globals.php在protected目录下。然后，在入口脚本index.php的，我们包括在
  • 设计模式之单例模式二（解决无序写入的问题） come_for_dream 单例模式volatile乱序执行双重检验锁
                    在上篇文章中我们使用了双重检验锁的方式避免懒汉式单例模式下由于多线程造成的实例被多次创建的问题，但是因为由于JVM为了使得处理器内部的运算单元能充分利用，处理器可能会对输入代码进行乱序执行（Out Of Order Execute）优化，处理器会在计算之后将乱序执行的结果进行重组，保证该
  • 程序员从初级到高级的蜕变 gcq511120594 框架工作PHPandroidhtml5
    软件开发是一个奇怪的行业，市场远远供不应求。这是一个已经存在多年的问题，而且随着时间的流逝，愈演愈烈。 我们严重缺乏能够满足需求的人才。这个行业相当年轻。大多数软件项目是失败的。几乎所有的项目都会超出预算。我们解决问题的最佳指导方针可以归结为——“用一些通用方法去解决问题，当然这些方法常常不管用，于是，唯一能做的就是不断地尝试，逐个看看是否奏效”。 现在我们把淫浸代码时间超过3年的开发人员称为
  • Reverse Linked List hcx2013 list
    Reverse a singly linked list.   /** * Definition for singly-linked list. * public class ListNode { * int val; * ListNode next; * ListNode(int x) { val = x; } * } */ p
  • Spring4.1新特性——数据库集成测试 jinnianshilongnian spring 4.1
    目录 Spring4.1新特性——综述 Spring4.1新特性——Spring核心部分及其他 Spring4.1新特性——Spring缓存框架增强 Spring4.1新特性——异步调用和事件机制的异常处理 Spring4.1新特性——数据库集成测试脚本初始化 Spring4.1新特性——Spring MVC增强 Spring4.1新特性——页面自动化测试框架Spring MVC T
  • C# Ajax上传图片同时生成微缩图(附Demo) liyonghui160com
        1.Ajax无刷新上传图片,详情请阅我的这篇文章。（jquery + c# ashx）         2.C#位图处理  System.Drawing。         3.最新demo支持IE7,IE8,Fir
  • Java list三种遍历方法性能比较 pda158 java
    从c/c++语言转向java开发，学习java语言list遍历的三种方法，顺便测试各种遍历方法的性能，测试方法为在ArrayList中插入1千万条记录，然后遍历ArrayList，发现了一个奇怪的现象，测试代码例如以下： package com.hisense.tiger.list; import java.util.ArrayList; import java.util.Iterator;
  • 300个涵盖IT各方面的免费资源（上）——商业与市场篇 shoothao seo商业与市场IT资源免费资源
    A.网站模板+logo+服务器主机+发票生成 HTML5 UP:响应式的HTML5和CSS3网站模板。 Bootswatch:免费的Bootstrap主题。 Templated:收集了845个免费的CSS和HTML5网站模板。 Wordpress.org|Wordpress.com:可免费创建你的新网站。 Strikingly:关注领域中免费无限的移动优
  • localStorage、sessionStorage uule localStorage
    W3School 例子   HTML5 提供了两种在客户端存储数据的新方法： localStorage - 没有时间限制的数据存储 sessionStorage - 针对一个 session 的数据存储   之前，这些都是由 cookie 完成的。但是 cookie 不适合大量数据的存储，因为它们由每个对服务器的请求来传递，这使得 cookie 速度很慢而且效率也不
按字母分类： ABCDEFGHIJKLMNOPQRSTUVWXYZ其他