HDU 5714 拍照 [杂题] [离散化]

拍照
Time Limit: 3000MS Memory Limit: 65536KB 64bit IO Format: %I64d & %I64u

Description
小明在旅游的路上看到了一条美丽的河，河上有许多船只，有的船只向左航行，有的船只向右航行。小明希望拍下这一美丽的风景，并且把尽可能多的船只都完整地拍到一张照片中。

小明位于河的边上，并且可以在河边的任意位置进行拍照，照相机的视野恰好为90度角，只能以垂直于河边的方向进行拍照。河上的船只全都可看作是平行于河边的一条线段，跟河边的距离各不相同，有的正在向左移动，有的正在向右移动，但移动速度恰好都是一样的。小明可以等待恰当的时间让尽量多的船只都走进照相机的视野里，你不需要考虑船只之间会互相遮挡视野的情况。

Input
第一行为T，表示输入数据组数。

下面T组数据，对于每组数据：

第一行是一个数n(1\leq n\leq 10^{4})，表示船只的数量。

接下来n行，每行四个整数
x,y,z,d(-10^{6}\leq x < y \leq 10^{6},1\leq z\leq 10^{4})，表示船只的左端点位置、右端点位置、距离河边的距离，以及航行的方向。d为-1表示向左航行，1表示向右航行。
Output
对第i组数据，输出

Case #i:

然后输出一行，仅包含一个整数，表示最多可以拍到多少完整的船只。
Sample Input
3
2
1 3 1 1
2 4 1 -1
2
1 3 1 -1
2 4 1 1
1
1 4 1 1
Sample Output
Case #1:
2
Case #2:
1
Case #3:
0
Source
2016”百度之星” - 复赛（Astar Round3）

---

这道题考思维呐。。。
首先对于向左走的船只，其左边且向右走的船只必定会与之相遇，那么算法大致就是枚举向左走的船只。
对于任意一条船，可以转换成能完整观察到的区间[y-z,x+z]，其中y-z>x+z时这条船无论如何也不能被完整地观察到。
不过数据10^6，用离散化吧。。
最后注意一只船有两个坐标点，数组开两倍！

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 10005;
struct Boat
{
    int x,y,dis,d;
    inline void read() { scanf("%d%d%d%d",&x,&y,&dis,&d); }
}boat[maxn];
int node[maxn<<1],maxnode;
int suml[maxn<<1],sumr[maxn<<1];
int n;
void init()
{
    memset(node,0,sizeof(node));
    maxnode=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        boat[i].read();
        Boat t=boat[i];
        int L=t.y-t.dis,R=t.x+t.dis;
        if(L>R) continue;
        node[++maxnode]=L,node[++maxnode]=R;
    }
    sort(node+1,node+maxnode+1);
    maxnode=unique(node+1,node+maxnode+1)-node-1;
}
int work()
{
    memset(suml,0,sizeof(suml));
    memset(sumr,0,sizeof(sumr));
    int ret=0;
    for(int i=1;i<=n;i++)
    {
        Boat t=boat[i];
        int L=t.y-t.dis,R=t.x+t.dis;
        if(L>R) continue;
        L=lower_bound(node+1,node+maxnode+1,L)-node;
        R=upper_bound(node+1,node+maxnode+1,R)-node;
        if(~t.d) sumr[L]++,sumr[R]--;
        else suml[L]++,suml[R]--;
    }
    for(int i=1;i<=maxnode;i++)
        suml[i]+=suml[i-1],sumr[i]+=sumr[i-1];
    int maxcnt=0;
    for(int i=1;i<=maxnode;i++)
    {
        smax(maxcnt,sumr[i]);
        smax(ret,maxcnt+suml[i]);
    }
    return ret;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("photo.in","r",stdin);
    freopen("photo.out","w",stdout);
#endif
    int T,cas=0;
    scanf("%d",&T);
    while(T--)
    {
        init();
        printf("Case #%d:\n%d\n",++cas,work());
    }
    return 0;
}