zoukankan      html  css  js  c++  java
  • POJ1631 Bridging signals LIS

      题目链接:http://poj.org/problem?id=1631

    先确定一方有序后,就是LIS问题了。

          

    仔细观察上图:红色航线是合法的,那么他们满足什么规律呢?
            C1   C2   C3    C4
    北岸红线的端点: 4     9    15    17
    南岸红线的端点: 2     8    12    17
            D1   D2   D3   D4
    不难看出无论线的斜率如何,都有这样的规律:
      C1<C2<C3<C4 且 D1<D2<D3<D4
    如果我们把输入数据排升序,问题变抽象为:
    在一个序列(D)里找到最长的序列满足DI<DJ<Dk……且i<j<k
    这样的话便是典型的最长非降子序列问题了。。。。
    法二:边界条件法(我前面总结提到的第4个方法)
    边界法其实就是把数据往小了缩,显然N=1是答案是1。N=2时呢?
    考虑这样一组数据:
      N=2
      C1 D1
      C2 D2
    当 C1<C2 时,如果D1>D2 那么一定会相交,反之则不会相交。
    当C1 >C2时,如果D1<D2 那么一定会相交,反之则不会相交。
      N=3
      C1 D1
      C2 D2
      C3 D3
      ……
    其实不用在推导N=3了,有兴趣的可以推导去。看N=2时就能得出:
    对于任意两条航线如果满足Ci<Cj 且Di<Dj 则两条航线不相交。这样的话要想在一个序列里让所有的航线都不相交那比然满足,C1<C2<C3…Cans且D1<D2<D3…<Dans ,也就是将C排序后求出最长的满足这个条件的序列的长度就是解。
    这样分析后显然是一个最长非降子序列问题。<摘自动态规划经典教程> ( 很不错的一本书,强烈推荐!)

     1 //STATUS:C++_AC_125MS_448KB
     2 #include <functional>
     3 #include <algorithm>
     4 #include <iostream>
     5 //#include <ext/rope>
     6 #include <fstream>
     7 #include <sstream>
     8 #include <iomanip>
     9 #include <numeric>
    10 #include <cstring>
    11 #include <cassert>
    12 #include <cstdio>
    13 #include <string>
    14 #include <vector>
    15 #include <bitset>
    16 #include <queue>
    17 #include <stack>
    18 #include <cmath>
    19 #include <ctime>
    20 #include <list>
    21 #include <set>
    22 #include <map>
    23 using namespace std;
    24 //define
    25 #define pii pair<int,int>
    26 #define mem(a,b) memset(a,b,sizeof(a))
    27 #define lson l,mid,rt<<1
    28 #define rson mid+1,r,rt<<1|1
    29 #define PI acos(-1.0)
    30 //typedef
    31 typedef __int64 LL;
    32 typedef unsigned __int64 ULL;
    33 //const
    34 const int N=40010;
    35 const int INF=0x3f3f3f3f;
    36 const int MOD=100000,STA=8000010;
    37 const LL LNF=1LL<<60;
    38 const double EPS=1e-8;
    39 const double OO=1e15;
    40 const int dx[4]={-1,0,1,0};
    41 const int dy[4]={0,1,0,-1};
    42 //Daily Use ...
    43 inline int sign(double x){return (x>EPS)-(x<-EPS);}
    44 template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
    45 template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
    46 template<class T> inline T Min(T a,T b){return a<b?a:b;}
    47 template<class T> inline T Max(T a,T b){return a>b?a:b;}
    48 template<class T> inline T Min(T a,T b,T c){return min(min(a, b),c);}
    49 template<class T> inline T Max(T a,T b,T c){return max(max(a, b),c);}
    50 template<class T> inline T Min(T a,T b,T c,T d){return min(min(a, b),min(c,d));}
    51 template<class T> inline T Max(T a,T b,T c,T d){return max(max(a, b),max(c,d));}
    52 //End
    53 
    54 int num[N],f[N];
    55 int T,n;
    56 
    57 int binary(int l,int r,int tar)
    58 {
    59     int mid;
    60     while(l<r){
    61         mid=(l+r)>>1;
    62         if(f[mid]<=tar)l=mid+1;
    63         else r=mid;
    64     }
    65     return l;
    66 }
    67 
    68 int main()
    69 {
    70  //   freopen("in.txt","r",stdin);
    71     int i,j,l,r,k;
    72     scanf("%d",&T);
    73     while(T--)
    74     {
    75         scanf("%d",&n);
    76         for(i=0;i<n;i++){
    77             scanf("%d",&num[i]);
    78         }
    79 
    80         l=1;r=2;
    81         f[1]=INF;
    82         for(i=0;i<n;i++){
    83             k=binary(l,r,num[i]);
    84             f[k]=num[i];
    85             r=Max(r,k+1);
    86         }
    87 
    88         printf("%d\n",r-1);
    89     }
    90     return 0;
    91 }
  • 相关阅读:
    新手小白如何购买阿里云服务器?
    nginx 泛域名解析 实现 多域名动态网站内容
    CentOS 安装MySQL 详细教程
    CentOS添加并挂载新硬盘
    docker下nginx+tomcat 文件的上传与访问
    关于人人开源renren-fast-vue 中npm install各种报错的解决方案
    《素面朝天,青春盛宴》
    手把手教你springboot中导出数据到excel中
    使Word文档和PDF文件不可被复制的方法
    博客园中使用MarkDown编辑器向博文中添加视频
  • 原文地址:https://www.cnblogs.com/zhsl/p/3089983.html
Copyright © 2011-2022 走看看