嵌套模型（DAG上的动态规划）—动态规划入门（算法经典入门）

矩形嵌套

时间限制： 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

[cpp]  view plain copy

1. #include <iostream>  
2. #include <string.h>  
3. #include <cstdio>  
4. using namespace std;  
5. int G[1010][1010];  
6. int d[1010];  
7. int n;  
8.   
9. int dp(int i) {  
10.     int &ans=d[i];  
11.     if(ans>0) return ans; //记忆化  
12.     ans=1;  
13.     for(int j=1; j<=n; j++) {  
14.         if(G[i][j]) {  
15.             if(ans < dp(j) + 1) {  
16.                 ans = dp(j)+1;  
17.             }  
18.         }  
19.     }  
20.     return ans;  
21. }  
22.   
23. void print_ans(int i){  
24.     printf("%d ",i);  
25.     for(int j=1;j<=n;j++) {  
26.         if(G[i][j]&&d[i]==d[j]+1){  
27.             print_ans(j);  
28.             break;  
29.         }  
30.     }  
31. }  
32.   
33. void print_map(){  
34.     for(int i = 1; i <= n; i++) {  
35.         for(int j = 1; j <= n; j++) {  
36.             cout << G[i][j] << ' ';  
37.         }  
38.         cout << endl;  
39.     }  
40. }  
41.   
42. int main()  
43. {  
44.     int N;  
45.     scanf("%d", &N);  
46.     while(N--){  
47.         int a[1010], b[1010];  
48.         memset(G, 0, sizeof(G));  
49.         memset(d, 0, sizeof(d));  
50.         scanf("%d", &n);  
51.         for(int i=1;i<=n;i++) {  
52.             scanf("%d%d", &a[i], &b[i]);  
53.         }  
54.         for(int i=1; i<=n; i++) {  
55.             for(int j=1;j<=n; j++) {  
56.                 if((a[i]>a[j]&&b[i]>b[j])||(a[i]>b[j]&&b[i]>a[j])) {  
57.                     G[i][j]=1;                                      ///建图,有向无环图。  
58.                 }  
59.             }  
60.         }  
61.         //print_map();  
62.         int max1 = 0;  
63.         //int max_i = 0;  
64.         for(int i=1; i<=n; i++) {  
65.             if(dp(i) > max1) {  
66.                 max1=dp(i);  
67.                // max_i = i;  
68.             }  
69.         }  
70.        // print_ans(max_i);  
71.        // cout << endl;  
72.         cout<<max1<<endl;  
73.     }  
74.     return 0;  
75. }  
76. /******************************* 
77. http://acm.nyist.edu.cn/JudgeOnline/status.php?pid=16 
78.  
79. *******************************/