百练1088 滑雪 dp

   动态规划，转化成图来做，即转化成有向图上的最长路问题。

   假设建图已好，存在邻接矩阵G中。d[ i ]表示从结点 i 出发的最长路长度，则

      d[i]=max { d[ j ] +1| ( i , j )属于E}

     E为边集合。 

    记忆化搜索程序（调用前d数组初始化为0）

    

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



 在主程序中，调用：
for(i=1;i<=n;i++)
    dp(i);
则最长路长度是d数组中最大的。

通常会是稀疏图，采用邻接表建有向图:

int n,m;
int first[max];
int u[max],v[max],w[max],next[max];
void read_graph()
{
   scanf("%d%d",&n,&m);//节点数，边数 
   for(int i=0;i<n;i++)
     first[i]=-1;
   for(int e=0;e<m;e++)//输入边 
   {
     scanf("%d%d%d",&u[e],&v[e],&w[e]);
     next[e]=first[u[e]];
     first[u[e]]=e;
   }
}   

这样建表后，若想判断结点 x和 y之间是否存在边，则可以：
int yes(int x,int y)
{
  int i,p;
  p=first[x];
  while(p!=-1)
  {
    if(v[p]==y) return 1;
    p=next[p];
  }
  return 0;
} 
这些问题解决后就差不多了


#include<stdio.h>
#include<stdlib.h>
#define N 20000
#define max 20002
int map[105][105],d[max];
int u[max],v[max],first[max],next[40005];
int ti[4]={-1,0,1,0},tj[4]={0,1,0,-1};
int h1,h2,m,n;
int yes(int x,int y)
{
  int i,p;
  p=first[x];
  while(p!=-1)
  {
    if(v[p]==y) return 1;
    p=next[p];
  }
  return 0;
} 
int dp(int i)
{
  int j;  
  if(d[i]>0)  return d[i];
  d[i]=1;
  for(j=1;j<=n*m;j++)
    if(yes(i,j))
    {
      if(d[i]<dp(j)+1) d[i]=dp(j)+1;
    }
  return d[i];
}  
int main()
{
  int i,j,k,e=0;
  scanf("%d%d",&n,&m);
  for(i=1;i<=n;i++)
     for(j=1;j<=m;j++)
       scanf("%d",&map[i][j]);
  for(i=0;i<=m+1;i++)
    map[0][i]=map[n+1][0]=N;
  for(i=0;i<=n+1;i++)
    map[i][0]=map[i][m+1]=N;
  
  memset(d,0,sizeof(d));  
  memset(first,-1,sizeof(first));
  
  for(i=1;i<=n;i++)
    for(j=1;j<=m;j++)
    { 
      h1=(i-1)*m+j;      
      for(k=0;k<4;k++)
        if(map[i][j]>map[i+ti[k]][j+tj[k]])
        {
          h2=(i+ti[k]-1)*m+j+tj[k];  
          u[e]=h1;
          v[e]=h2;
          next[e]=first[u[e]];
          first[u[e]]=e;
          e++;
        }
    }
    
  for(i=1;i<=n*m;i++)
    dp(i);
  k=d[1];
  for(i=2;i<=n*m;i++)
    if(k<d[i]) k=d[i];
  printf("%d",k); 
  system("pause"); 
  return 0;
}