USACO2010FEB gold ice——STL的使用

要懂的知识

$pair$
$set$
$map$
迭代器
$upper\_bound$
$lower\_bound$

题目

题目描述
Bessie 在一个冰封的湖面上游泳，湖面可以表示为二维的平面，坐标范围是-1,000,000,000..1,000,000,000。湖面上的N(1 <= N <= 20,000)个位置有石块(编号分别为1到N)，其它位置是冰面。由于Bessie滑冰技术不够好，她通过推动自己旁边的石块，依靠反作用力向某一个方向前进，在碰到一个新的石块之前，Bessie是不会停下来的。(当然，最后会停留在某块石块的前一个格子里)由于Bessie无法计算复杂的角度，她只能够向东南西北四个方向前进。很显然，Bessie不能够穿越石块，因此，Bessie仅仅可以向三个方向滑。
滑冰不是没有风险，Bessie滑向某个方向后必须能碰到某个石块，因此她必须很小心。考虑下面的一个情况，Bessie希望到达在她东面的目标位置(x=5,y=1)，(. = 冰块,* = 石头, B = Bessie, G = 目的位置)如果她直接向东滑，那么她会滑过目标位置，因为她通过撞上某块石头来停下来，一个能够到达目标位置的方案是这样的：

在$(a)$中，Bessie
只能朝向北，东，或者南三个方向，但是只有在北面能撞上石头从而停下来，在$(b)$中，类似地，她只能向东走。
对于输入，石头$i$位于坐标为$X_i,Y_i$的位置，$(-1,000,000,000<= X_i <= 1,000,000,000; -1,000,000,000 <= Y_i <= 1,000,000,000)，$没有任何两块石头位于同一个位置，Bessie从Bx,By的位置出发(出发点一定与某个石头相邻)，Bessie的目标位置是$Gx,Gy(-1,000,000,000 <= Gx <= 1,000,000,000; -1,000,000,000 <= Gy <=1,000,000,000).$
Bessie 并不介意长时间滑冰，但是，不停地推石头依靠反作用力前进很累。FJ 非常关心Bessie的健康，因此他希望知道Bessie最少要推多少次石头才能到达终点。
输入格式
第一行: 五个用空格隔开的整数:$N, Bx, By, Gx, and Gy$
第二行到第$N+1$行: 第$i+1$行用两个空格隔开的整数来描述第i个石头
输出格式
第一行: 一个整数表示Bessie至少要推多少次石头才能够到达终点

不用STL的算法

这道题目是一个典型的BFS。
需要注意的是，我们需要把石子数做一个前缀和，进行二分查找石子位置就行了。
但是各种分情况，大概需要写12个不同的二分。虽然这些二分代码差不多，但是代码量太大了。查错也是一个难题。

用map和set优化

我们不用传统的数组来存放这些点的位置。

这样开两个映射，建立一个整数到一个集合的映射。$mapx[t]$代表横坐标为$t$的所有的点 的纵坐标的集合。
($ps$:这里$mapx[t]$可以直接当做一个$set$来用)
当我们要找在某一个位置找同行或者同列的石头的时候，可以使用$lower\_bound$函数.
($ps$:这函数的用处在于，在一个集合里找到，>=某个数的最小数，并返回它在集合中的地址)
如果要找坐标比当前位置大的石子，那就直接使用这个函数；如果要找坐标比当前位置小的石子，那就找到之后，位置-1就行了.

---

这里的pos结构体是用来存一个点的所有信息，x和y是坐标，w是最短距离。
POS是一个点的坐标。$dis$这个map把一个坐标映射到bool类型上，可以构造一个类似于以坐标为下标的数组，用来标记某个点是否被访问过。
q是广搜队列。

code

/*
by 柴犬首相
STL的使用
#290 ice 
2018.3.15 
*/ 
#include<bits/stdc++.h>
using namespace std;
inline int read()
{
    int num=0;
    char c=' ';
    bool flag=true;
    for(;c>'9'||c<'0';c=getchar())
    if(c=='-')
    flag=false;
    for(;c>='0'&&c<='9';num=num*10+c-48,c=getchar());
    return flag ? num : -num;
}//快读 
struct pos
{
    int x,y,w;
}B,G;
//这是一个点的所有信息，
//x和y代表是这个点的坐标
//w是代表源点B到这个点的最短距离 
typedef pair<int,int> POS;
//这个pair是存储坐标的 （用来当做map的下标） 
map<POS,bool>dis;
//这是一个以坐标为下标的数组
//用来存放某一个点是否访问过 
queue<pos>q;
//BFS队列 
map<int,set<int> >mapx;
map<int,set<int> >mapy;
//这是一个整数到集合的映射
//mapx[i]代表以横坐标为i的所有的点
//的纵坐标的集合 
int n;
void init()
{
    n=read();
    B.x=read();B.y=read();
    G.x=read();G.y=read();
    for (int i=1;i<=n;i++)
    {
        int a=read();
        int b=read();
        mapx[a].insert(b);
        mapy[b].insert(a);
    }
}//读入，自行理解！！ 
void bfs()
{
    q.push(B);//初始点 
    set<int>t;//临时集合，用来简化代码 
    set<int>::iterator it;
    //这是一个迭代器，存放的是元素在集合中的地址 
    //并且这个指针只会指向集合中的元素 
    while (q.size()>0)
    {
        pos now=q.front();//队首元素，当前位置 
        //思考1:q.pop()为什么不能放在这里？ 
        t=mapx[now.x];//用来代替，简化下面的代码 
        it=t.lower_bound(now.y);//找到目标点 
        if (it!=t.end()&&*it-1>now.y&&!dis[make_pair(now.x,*it-1)])
        //如果it==t.end()，那就代表，找不到比它大的 
        {
            if (now.x==G.x&&*it-1==G.y) return;//如果到目标点，那就弹出来 
            //思考2：这些判断为什么不能放在循环节开始的时候？ 
            dis[make_pair(now.x,*it-1)]=1;//已经访问过 
            q.push((pos){now.x,*it-1,now.w+1});//压入队列 
        }//向上 
        if (it!=t.begin()&&*(--it)+1<now.y&&!dis[make_pair(now.x,*it+1)])
        //这里有点不一样，因为是向下走，所以目标点不是it，
        //而是集合中it之前的那个元素 
        {
            if (now.x==G.x&&*it+1==G.y) return;//如果到目标点，那就弹出来 
            dis[make_pair(now.x,*it+1)]=1;
            q.push((pos){now.x,*it+1,now.w+1});
        }//向下 
        t=mapy[now.y];
        it=t.lower_bound(now.x);
        if (it!=t.end()&&*it-1>now.x&&!dis[make_pair(*it-1,now.y)])
        {
            if (now.y==G.y&&*it-1==G.x) return;//如果到目标点，那就弹出来 
            dis[make_pair(now.x,*it-1)]=1;
            q.push((pos){*it-1,now.y,now.w+1});
        }//向右 
        if (it!=t.begin()&&*(--it)+1<now.x&&!dis[make_pair(*it+1,now.y)])
        {
            if (now.y==G.y&&*it+1==G.x) return;//如果到目标点，那就弹出来 
            dis[make_pair(*it+1,now.y)]=1;
            q.push((pos){*it+1,now.y,now.w+1});
        }//向左 
        q.pop();//弹出 
    }
}
int main()
{
    init();
    if (G.x==B.x&&G.y==B.y) printf("%d\n",0);
    else
    {
        bfs();
        printf("%d\n",q.front().w+1); 
    }
    return 0;
}