题目描述
有一条横贯东西的大河,河有笔直的南北两岸,岸上各有位置各不相同的N个城市。北岸的每个城市有且仅有一个友好城市在南岸,而且不同城市的友好城市不相同。每对友好城市都向政府申请在河上开辟一条直线航道连接两个城市,但是由于河上雾太大,政府决定避免任意两条航道交叉,以避免事故。编程帮助政府做出一些批准和拒绝申请的决定,使得在保证任意两条航道不相交的情况下,被批准的申请尽量多。
输入输出格式
输入格式:
第1行,一个整数N,表示城市数。
第2行到第n+1行,每行两个整数,中间用一个空格隔开,分别表示南岸和北岸的一对友好城市的坐标。
输出格式:
仅一行,输出一个整数,表示政府所能批准的最多申请数。
输入输出样例
输入样例#1: 复制
7
22 4
2 6
10 3
15 12
9 8
17 17
4 2
输出样例#1: 复制
4
说明
50% 1<=N<=5000,0<=xi<=10000
100% 1<=N<=2e5,0<=xi<=1e6
最长上升子序列(LIS)的典型变形,熟悉的n^2的动归会超时。LIS问题可以优化为nlogn的算法。
定义d[k]:长度为k的上升子序列的最末元素,若有多个长度为k的上升子序列,则记录最小的那个最末元素。
注意d中元素是单调递增的,下面要用到这个性质。
首先len = 1,d[1] = a[1],然后对a[i]:若a[i]>d[len],那么len++,d[len] = a[i];
否则,我们要从d[1]到d[len-1]中找到一个j,满足d[j-1]<a[i]<d[j],则根据D的定义,我们需要更新长度为j的上升子序列的最末元素(使之为最小的)即 d[j] = a[i];
最终答案就是len
利用d的单调性,在查找j的时候可以二分查找,从而时间复杂度为nlogn。
#include <algorithm>
#include <iostream>
#include <cstdio>
using namespace std;
struct cyq
{
int x,y;
}a[200005];
int n,d[200005],ans;
bool cmp(cyq aa,cyq bb){
if(aa.x==bb.x) return aa.y<bb.y;
return aa.x<bb.x;
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y);
sort(a+1,a+1+n,cmp);
d[++ans]=a[1].y;
for(int i=2;i<=n;i++)
{
int point=upper_bound(d+1,d+ans+1,a[i].y)-d;
d[point]=a[i].y;
if(point>ans)
ans++;
}
printf("%d\n",ans);
return 0;
}