NYOJ 201 作业题 动态规划

作业题

时间限制： 3000 ms  |  内存限制： 65535 KB

难度： 3

描述

小白同学这学期有一门课程叫做《数值计算方法》，这是一门有效使用数字计算机求数学问题近似解的方法与过程，以及由相关理论构成的学科……

今天他们的Teacher S，给他们出了一道作业题。Teacher S给了他们很多的点，让他们利用拉格朗日插值公式，计算出某严格单调函数的曲线。现在小白抄下了这些点，但是问题出现了，由于我们的小白同学上课时走了一下神，他多抄下来很多点，也就是说这些点整体连线不一定还是严格递增或递减的了。这可怎么处理呢。为此我们的小白同学制定了以下的取点规则：

1、取出尽可能多的满足构成严格单调曲线的点，作为曲线上的点。

2、通过拉格朗日插值公式，计算出曲线的方程

但是，他又遇到了一个问题，他发现他写下了上百个点。[- -!佩服吧],这就很难处理了(O_O).。由于拉格朗日插值公式的计算量与处理的点数有关，因此他请大家来帮忙，帮他统计一下，曲线上最多有多少点，以此来估计计算量。

已知：没有任何两个点的横坐标是相同的。

输入

本题包含多组数据：
首先，是一个整数T，代表数据的组数。
然后，下面是T组测试数据。对于每组数据包含两行：
第一行：一个数字N(1<=N<=999),代表输入的点的个数。
第二行：包含N个数对X(1<=x<=10000),Y(1<=Y<=10000),代表所取的点的横纵坐标。

输出

每组输出各占一行，输出公一个整数，表示曲线上最多的点数

样例输入

2
2
1 2 3 4
3
2 2 1 3 3 4

样例输出

2
2

     
     
      
      思路：因为横坐标没有重复的，所以不用考虑垂直上升或下降
     
     
     
     
      
      因此，只需将横坐标升序排列后，求出纵坐标的最长递增序列长度与最长递减序列长度
     
     
     
     
      
      二者中最大的就是答案
     
     
     
     
      
      

       
       01.
       
       #include<stdio.h>
      
      
      
      

       
       02.
       
       #include<string.h>
      
      
      
      

       
       03.
       
       #include<algorithm>
      
      
      
      

       
       04.
       
       using namespace std;
      
      
      
      

       
       05.
       
       int dp1[10005],dp2[10005];
      
      
      
      

       
       06.
       
       struct point  //定义结构体，存储坐标
      
      
      
      

       
       07.
       
       {
      
      
      
      

       
       08.
       
       int x,y;
      
      
      
      

       
       09.
       
       }a[10005];  
      
      
      
      

       
       10.
       
       int max(int x,int y)  //  求最大值
      
      
      
      

       
       11.
       
       {
      
      
      
      

       
       12.
       
       return x>y?x:y;
      
      
      
      

       
       13.
       
       }
      
      
      
      

       
       14.
       
       bool cmp(point b,point c)  //对x从小到大排序
      
      
      
      

       
       15.
       
       {
      
      
      
      

       
       16.
       
       return b.x<c.x;
      
      
      
      

       
       17.
       
       }
      
      
      
      

       
       18.
       
       int main()
      
      
      
      

       
       19.
       
       {
      
      
      
      

       
       20.
       
       int n,t,i,j,p,q;
      
      
      
      

       
       21.
       
       scanf("%d",&t);
      
      
      
      

       
       22.
       
       while(t--)
      
      
      
      

       
       23.
       
       {
      
      
      
      

       
       24.
       
       memset(dp1,0,sizeof(dp1));  //储存最长递增序列长度
      
      
      
      

       
       25.
       
       memset(dp2,0,sizeof(dp2));  //储存最长递减序列长度
      
      
      
      

       
       26.
       
       getchar();
      
      
      
      

       
       27.
       
       scanf("%d",&n);
      
      
      
      

       
       28.
       
       for(i=0;i<n;i++)
      
      
      
      

       
       29.
       
       scanf("%d%d",&a[i].x,&a[i].y);
      
      
      
      

       
       30.
       
       sort(a,a+n,cmp);  //排序
      
      
      
      

       
       31.
       
       for(i=1;i<n;i++)
      
      
      
      

       
       32.
       
       for(j=0;j<i;j++)
      
      
      
      

       
       33.
       
       {
      
      
      
      

       
       34.
       
       if(a[i].y>a[j].y)  //递增序列
      
      
      
      

       
       35.
       
       dp1[i]=max(dp1[i],dp1[j]+1);
      
      
      
      

       
       36.
       
       if(a[i].y<a[j].y)   //递减序列
      
      
      
      

       
       37.
       
       dp2[i]=max(dp2[i],dp2[j]+1);
      
      
      
      

       
       38.
       
       }
      
      
      
      

       
       39.
       
       int max1=dp1[0],max2=dp2[0];
      
      
      
      

       
       40.
       
       for(i=1;i<n;i++)  //求最大值
      
      
      
      

       
       41.
       
       {
      
      
      
      

       
       42.
       
       if(dp1[i]>max1)
      
      
      
      

       
       43.
       
       max1=dp1[i];
      
      
      
      

       
       44.
       
       if(dp2[i]>max2)
      
      
      
      

       
       45.
       
       max2=dp2[i];
      
      
      
      

       
       46.
       
       }
      
      
      
      

       
       47.
       
       printf("%d\n",max(max1,max2)+1);  
      
      
      
      

       
       48.
       
       }
      
      
      
      

       
       49.
       
       }