(DAG)NestedRectangle

//这个算法是书上的提示，比我的版本好多了。简洁

#include<iostream>
#include<fstream>
#include<cstdlib>
using namespace std;


ifstream fin("C:\\data12_1.in");


struct rect
{
rect(int a=0,int b=0):length(a),width(b){}
int length;
int width;
};


rect arr[100];
int n;
int G[100][100];
int d[100];


void Init()
{
int a,b;
int temp;
n=0;
while(fin>>a>>b)
{
if(a<b)
{
temp=a;
a=b;
b=temp;
}
rect rt(a,b);
arr[n++]=rt;
}
memset(G,0,sizeof(G));
memset(d,0,sizeof(d));
for(int i=0;i<n;++i)
{
for(int j=i+1;j<n;++j)
{
if(arr[i].length>arr[j].length&&arr[i].width>arr[j].width)
G[i][j]=1;
else if(arr[i].length<arr[j].length&&arr[i].width<arr[j].width)
G[j][i]=1;
}
}
}

//利用了全局数组和有向无环图的性质，用递归简化代码。。。
int dp(int i)
{
int& ans=d[i];
if(ans>0)
return ans;
ans=1;
for(int j=0;j<n;++j)
   if(G[i][j]==1)
   ans>?=dp(j)+1;
return ans;
}


void print(int i)
{
cout<<"Rect("<<arr[i].length<<","<<arr[i].width<<")\t";
for(int j=0;j<n;++j)
{
if(G[i][j]==1&&d[i]==d[j]+1)
{
print(j);
break;
}
}
}


int main()
{
    Init();
    dp(0);
int max=0;
int pos;
for(int i=0;i<n;++i)
{
if(max < dp(i))
{
max=dp(i);
pos=i;
}
}
cout<<"最长路径为："<<max<<endl;
print(pos);
system("pause");
return 0;
}