基于动态规划对友好城市问题O(nlogn)做法的深入研究

Problem（附链接）

题目描述

有一条横贯东西的大河，河有笔直的南北两岸，岸上各有位置各不相同的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

Solution

这道题应该不难想到用dp吧！
但我们最可能想到的应该是O(n^2)的朴素dp，这儿就不多讲解（我们今天的重点内容是O(n*logn)的dp），但还是贴一下代码。

Code–O(n^2)

#include<bits/stdc++.h>
using namespace std;

const int MAX=2e5+10; 

struct node
{
	int X1;
	int X2;
}line[MAX];

int n,f[MAX],ans;

bool cmp(const node &a,const node &b)
{
	return a.X1<b.X1;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d %d",&line[i].X1,&line[i].X2);
	sort(line+1,line+n+1,cmp);
	for(int i=1;i<=n;i++) f[i]=1;
	for(int i=2;i<=n;i++)
	{
		for(int j=1;j<=i-1;j++)
		{
			if(line[j].X2<line[i].X2) f[i]=max(f[i],f[j]+1);
		}
	}
	for(int i=1;i<=n;i++) ans=max(ans,f[i]);
	cout<<ans;
	return 0;
}

好了，现在我们该想想应该如何改进这个做法。

我们还是回到问题的本质，当你把这些城市的南岸城市（北岸也可以）坐标从小到大排好序后，我们分析一下选择哪些北岸城市的坐标可以构成“合法航线”？

来看张图
当A1A2(A1<B1，即B1在A1的右边)连好后，要使B1B2与A1A2成为合法航线组，那就要让B2也在A2的右边（即A2<B2)。

等等！这不就是让北岸的坐标构成一个最长上升子序列吗？
没错，这道题的答案就是这个最长上升子序列(LIS)的长度！
这个长度还不好求吗？就是用O(nlogn)做法的模板呀！
（不会用O(nlogn)做法求最长上升子序列长度的朋友请参看我的另一篇文章） 【链接】
好了，看看代码吧！

Code–O(nlogn)

#include<bits/stdc++.h>
using namespace std;

const int MAX=2e5+10; 

struct node
{
	int X1;
	int X2;
}line[MAX];

int n,f[MAX],t[MAX],cnt;

bool cmp(const node &a,const node &b)
{
	return a.X1<b.X1;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d %d",&line[i].X1,&line[i].X2);
	sort(line+1,line+n+1,cmp);
	t[++cnt]=line[1].X2;
	for(int i=2;i<=n;i++)
	{
		int tmp=lower_bound(t+1,t+cnt+1,line[i].X2)-t;
		if(tmp<=cnt) t[tmp]=line[i].X2;
		else t[++cnt]=line[i].X2;
	}
	cout<<cnt;
	return 0;
}

代码中使用了一个二分查找函数----lower_bound( )，不会的朋友可以参看我的另一篇文章。【链接】

总结与思考

现在回过头再来看一下这道题，其实并不难对吧？

那当初为什么没想到O(nlogn)的做法呢？其实就是没有将问题的本质化为一个求最长上升子序列长度的模板。

所以一个OI的应该配备一种基本素养：
那就是应针对竞赛题目的要求构建数学模型，将问题抽象，进而构造出计算机可以接受的算法，最后才是编写出计算机能够执行的程序！

学习厌倦了？点我有更多精彩哦！