嵌套矩形（nyoj 16)

矩形嵌套
时间限制：3000 ms  |  内存限制：65535 KB 
难度：4
描述 
有n个矩形，每个矩形可以用a,b来描述，表示长和宽。矩形X(a,b)可以嵌套在矩形Y(c,d)中当且仅当a<c,b<d或者b<c,a<d（相当于旋转X90度）。例如（1,5）可以嵌套在（6,2）内，但不能嵌套在（3,4）中。你的任务是选出尽可能多的矩形排成一行，使得除最后一个外，每一个矩形都可以嵌套在下一个矩形内。
输入
第一行是一个正正数N(0<N<10)，表示测试数据组数，
每组测试数据的第一行是一个正正数n，表示该组测试数据中含有矩形的个数(n<=1000)
随后的n行，每行有两个数a,b(0<a,b<100)，表示矩形的长和宽
输出
每组测试数据都输出一个数，表示最多符合条件的矩形数目，每组输出占一行
样例输入
1
10
1 2
2 4
5 8
6 10
7 9
3 1
5 8
12 10
9 7
2 2
样例输出
5

 

 

分析：看了刘汝佳的讲解，思路是这样的。如果我们将每个矩形看做是图上面的一个点，那么矩形A可以套住矩形B看做是图上的由A指向B的一条边，那么我们就可以将这个题目转化为求DAG（有向无环图）上的最长路了，这就是DAG上的动态规划，转移方程是d(i)=max{ d(j) + 1 | (i,j) 属于 E }。看下面代码：

 

代码1：

 

 
#include<stdio.h>
#include<string.h>
#include<stdlib.h> 

typedef struct
{
	int w;
	int l;
}M;

M a[1010];//用结构体数组来记录矩形的信息，方便对矩形进行排序，为建立map做准备; 

int map[1010][1010], max, d[1010],n;

void change(int *i, int *j)
{
	int t;
	t = *i;
	*i = *j;
	*j = t;
}

int compare(void const*p, void const*q)
{
	int x1 = ((M *)p)->w;
	int x2 = ((M *)q)->w;
	int y1 = ((M *)p)->l;
	int y2 = ((M *)q)->l;
	if(x1 != x2)
	{
		return x2 - x1;//如果长不相等，按长的大小排序 
	}
	else 
	{
		return y2 - y1;//长相等按宽的大小排序; 
	}
}

void init_map()
{
	int i, j;
	for(i= 0; i < n; i++)
	{
		for(j = 0; j < n; j++)
		{
			map[i][j] = 0;
		}
	}
}

int dp(int i)
{
	int j;
	if(d[i] > 0)//d[i] > 0说明先前已经从这个点搜索过已经记下了他的最长路径; 
	{
		return d[i];
	}
	d[i] = 1;
	for(j = 0; j < n;j++)
	{
		if(map[i][j] == 1 && d[i] < dp(j) + 1)
		{
			d[i] = d[j] + 1;
		}
	}
	return d[i];
}

int main()
{
	int N,i,j;
	scanf("%d", &N);
	while(N--)
	{
		max = 0;
		memset(a, 0, sizeof(a));
		memset(d, 0, sizeof(d));
		scanf("%d", &n);
		init_map();
		for(i = 0; i < n; i++)
		{
			scanf("%d%d", &a[i].w, &a[i].l);
			if(a[i].w < a[i].l)//把矩形按照统一的大的在前(长)，小的在后（宽）的顺序排列， 为后面的排列做准备; 
			{
				change(&a[i].w, &a[i].l);
			}
		}
		qsort(a,n, sizeof(M), compare);//把结构体数组按照从大到小的顺序排列； 
		for(i = 0; i < n; i++)
		{
			for(j = 0; j < i; j++)
			{
				if(a[i].w < a[j].w && a[i].l < a[j].l) 
				{
					//printf("%d %d %d %d %d %d\n", i, j, a[i].w, a[i].l, a[j].w,a[j].l);
					map[i][j] = 1;
				}
			}
		}
		for(i = n - 1; i >= 0; i--)//特别注意因为这个图可能不是连通图，所以不能只记忆搜索一次，而要把所有点都搜一遍; 
		{
		    dp(i); //从最小的出发，如果是连通图那么只要去递归一次就能求出答案，剩下的只要进入递归函数就会因为d[i] > 0 退出递归; 
		}
	/*	for(i = 0; i < n; i++)
		{
			printf("%d\n", d[i]);
		}*/
		for(i = 0; i < n; i++)
		{
			if(d[i]> max)
			{
				max = d[i];
			}
		}
		printf("%d\n", max);		
	}
	return 0;
}       

代码2:只是把邻接矩阵改成了vector;

 
 
#include<stdio.h>
#include<vector> 
#include<string.h>
#include<stdlib.h> 
using namespace std;

typedef struct
{
	int w;
	int l;
}M;

M a[1010];//用结构体数组来记录矩形的信息，方便对矩形进行排序，为建立map做准备; 


vector<int> map[1010];
int  maxn, d[1010],n;

void change(int *i, int *j)
{
	int t;
	t = *i;
	*i = *j;
	*j = t;
}

int compare(void const*p, void const*q)
{
	int x1 = ((M *)p)->w;
	int x2 = ((M *)q)->w;
	int y1 = ((M *)p)->l;
	int y2 = ((M *)q)->l;
	if(x1 != x2)
	{
		return x2 - x1;//如果长不相等，按长的大小排序 
	}
	else 
	{
		return y2 - y1;//长相等按宽的大小排序; 
	}
}

void init_map()
{
	int i;
	for(i= 0; i < n; i++)
	{
		map[i].clear();
	}
}

int dp(int i)
{
	int j;
	if(d[i] > 0)//d[i] > 0说明先前已经从这个点搜索过已经记下了他的最长路径; 
	{
		return d[i];
	}
	d[i] = 1;
	for(j = 0; j < map[i].size(); j++)
	{
		if(d[i] < dp(map[i][j]) + 1)//注意：这里是map[i][j]很容易写成j;
		{
			d[i] = d[map[i][j]] + 1;
		}
	}
	return d[i];
}

int main()
{
	int N,i,j;
	scanf("%d", &N);
	while(N--)
	{
		maxn = 0;
		memset(a, 0, sizeof(a));
		memset(d, 0, sizeof(d));
		scanf("%d", &n);
		init_map();
		for(i = 0; i < n; i++)
		{
			scanf("%d%d", &a[i].w, &a[i].l);
			if(a[i].w < a[i].l)//把矩形按照统一的大的在前(长)，小的在后（宽）的顺序排列， 为后面的排列做准备; 
			{
				change(&a[i].w, &a[i].l);
			}
		}
		qsort(a,n, sizeof(M), compare);//把结构体数组按照从大到小的顺序排列； 
		for(i = 0; i < n; i++)
		{
			for(j = 0; j < i; j++)
			{
				if(a[i].w < a[j].w && a[i].l < a[j].l) 
				{
					//printf("%d %d %d %d %d %d\n", i, j, a[i].w, a[i].l, a[j].w,a[j].l);
					map[i].push_back(j);
				}
			}
		}
		for(i = n - 1; i >= 0; i--)//特别注意因为这个图可能不是连通图，所以不能只记忆搜索一次，而要把所有点都搜一遍; 
		{
		    dp(i); //从最小的出发，如果是连通图那么只要去递归一次就能求出答案，剩下的只要进入递归函数就会因为d[i] > 0 退出递归; 
		}
		/*for(i = 0; i < n; i++)
		{
			printf("%d\n", d[i]);
		}*/
		for(i = 0; i < n; i++)
		{
			if(d[i]> maxn)
			{
				maxn = d[i];
			}
		}
		printf("%d\n", maxn);		
	}
	return 0;
}                        

 

 

代码3：

递推：

 
#include<stdio.h>
#include<vector> 
#include<string.h>
#include<stdlib.h> 
using namespace std;

typedef struct
{
	int w;
	int l;
}M;

M a[1010];//用结构体数组来记录矩形的信息，方便对矩形进行排序，为建立map做准备; 


vector<int> map[1010];
int  maxn, d[1010],n;

void change(int *i, int *j)
{
	int t;
	t = *i;
	*i = *j;
	*j = t;
}

int compare(void const*p, void const*q)
{
	int x1 = ((M *)p)->w;
	int x2 = ((M *)q)->w;
	int y1 = ((M *)p)->l;
	int y2 = ((M *)q)->l;
	if(x1 != x2)
	{
		return x2 - x1;//如果长不相等，按长的大小排序 
	}
	else 
	{
		return y2 - y1;//长相等按宽的大小排序; 
	}
}

void init_map()
{
	int i;
	for(i= 0; i < n; i++)
	{
		map[i].clear();
	}
}


int main()
{
	int N,i,j, k;
	scanf("%d", &N);
	while(N--)
	{
		maxn = 0;
		memset(a, 0, sizeof(a));
		memset(d, 0, sizeof(d));
		scanf("%d", &n);
		init_map();
		for(i = 0; i < n; i++)
		{
			scanf("%d%d", &a[i].w, &a[i].l);
			if(a[i].w < a[i].l)//把矩形按照统一的大的在前(长)，小的在后（宽）的顺序排列， 为后面的排列做准备; 
			{
				change(&a[i].w, &a[i].l);
			}
		}
		qsort(a,n, sizeof(M), compare);//把结构体数组按照从大到小的顺序排列； 
		for(i = 0; i < n; i++)
		{
			for(j = 0; j < i; j++)
			{
				if(a[i].w < a[j].w && a[i].l < a[j].l) 
				{
					//printf("%d %d %d %d %d %d\n", i, j, a[i].w, a[i].l, a[j].w,a[j].l);
					map[i].push_back(j);
				}
			}
		}
		for(i = 0; i < n; i++)
		{
			j = map[i].size();
			if(j == 0)
			{
				d[i] = 1;
			}
			else
			{
				for(k = 0; k < j; k++)
				{
					if(d[i] < d[map[i][k]] + 1)
					{
						d[i] = d[map[i][k]] + 1;
					}
				}
			}
		}
		for(i = 0; i < n; i++)
		{
			if(d[i]> maxn)
			{
				maxn = d[i];
			}
		}
		printf("%d\n", maxn);		
	}
	return 0;
}