紫书p263 嵌套矩形如何输出所有路径（DAG上的动态规划）

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <stdlib.h>
#include <algorithm>

using namespace std;

const int maxn=1000+5;
int G[maxn][maxn],d[maxn],n,b[maxn],sum=0,path[maxn];//b数组保存起点，sum表示起点有多少个

int dp(int i)
{
int& ans = d[i];
if(ans > 0)
return ans;
ans =1;
for(int j = 1; j <= n; j++)
if(G[i][j])
ans = max(ans, dp(j)+1);
return ans;
}

void print_ans(int i,int k)
{

path[k]=i;
if(d[i]==1)
{
for(int i=0;i<=k;i++)
printf("%d ",path[i]);
printf("\n");
return ;
}
for(int j = 1; j <=n; j++)
if(G[i][j]&&d[i]==d[j]+1)
{
print_ans(j,k+1);

}
}

int main()
{
// freopen("in.txt","r",stdin);
memset(G,0,sizeof(G));
memset(d,0,sizeof(d));
int x[maxn],y[maxn],t,max_sum=0,sub;
scanf("%d",&n);
for(int i=1; i<=n; i++)
{
scanf("%d %d",&x[i],&y[i]);
if(x[i]<y[i])
{
t=x[i];
x[i]=y[i];
y[i]=t;
}
}
for(int i=1; i<=n; i++)
for(int j=i+1; j<=n; j++)
{
if(x[i]>x[j]&&y[i]>y[j])
G[i][j]=1;
else if(x[i]<x[j]&&y[i]<y[j])
G[j][i]=1;
}
for(int i=1; i<=n; i++)
dp(i);
for(int i=1; i<=n; i++)
{
if(d[i]>max_sum)
{
max_sum=d[i];
sub=i;
}
}
for(int i=1; i<=n; i++)
{
if(d[i]==max_sum)
b[++sum]=i;
}
for(int i=1;i<=sum;i++)
print_ans(b[i],0);
return 0;
}