4845. 【NOIP2016提高A组集训第5场11.2】寻找

Description

“我有个愿望，我希望穿越一切找到你。”
这是个二维平面世界，平面上有n个特殊的果实，我从(0,0)点出发，希望得到尽量多的果实，但是出于某种特殊的原因，我的运动方式只有三种（假设当前我在(x,y)）：
1、我可以走到(x+1,y)
2、我可以走到(x,y+1)
3、我可以走到(x+1,y+1)
现在我需要你的帮助，帮我找出我最多能够得到多少个果实。

Data Constraint

对于70%的数据1<=n<=1000
对于100%的数据1<=n<=100000，-10^9<=x,y<=10^9

Solution

这道题我们可以直接将横坐标为第一关键字，纵坐标为第二关键字排序，将纵坐标离散化。然后对于一个点，我们用树状数组查询比当前纵坐标小的点的最大值是多少，并把这个点的答案加入树状数组。

代码

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=100005;
struct code{
    int a,b;
}a[maxn];
int n,i,t,j,k,l,num,x,y,f[maxn];
bool cmp(code x,code y){
    return x.a<y.a || x.a==y.a && x.b<y.b;
}
int lowbit(int x){
    return x&(-x);
}
int find(int x){
    if (!x) return 0;
    return max(f[x],find(x-lowbit(x)));
}
void insert(int x){
    if (x>n) return;
    f[x]=max(f[x],t);
    insert(x+lowbit(x));
}
int main(){
    //freopen("find.in","r",stdin);freopen("find.out","w",stdout);
    scanf("%d",&n);
    for (i=1;i<=n;i++){
        scanf("%d%d",&x,&y);
        if (x<0 || y<0) continue;
        a[++num].b=x;a[num].a=y;
    }
    n=num;num=0;
    sort(a+1,a+n+1,cmp);
    x=-1;
    for (i=1;i<=n;i++){
        if (a[i].a!=x) num++;
        x=a[i].a;a[i].a=num;
    }
    for (i=1;i<=n;i++)
        swap(a[i].a,a[i].b);
    sort(a+1,a+n+1,cmp);
    for (i=1;i<=n;i++){
        t=find(a[i].b)+1;
        insert(a[i].b);
    }
    t=find(n);
    printf("%d\n",t);
}