链接:



吉哥系列故事——完美队形I

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)

Total Submission(s): 1454    Accepted Submission(s): 420

Problem Description
  吉哥这几天对队形比较感兴趣。

  有一天,有n个人按顺序站在他的面前,他们的身高分别是h[1], h[2] ... h[n],吉哥希望从中挑出一些人,让这些人形成一个新的队形,新的队形若满足以下三点要求,则称之为完美队形:

  

  1、挑出的人保持他们在原队形的相对顺序不变;

  2、左右对称,假设有m个人形成新的队形,则第1个人和第m个人身高相同,第2个人和第m-1个人身高相同,依此类推,当然,如果m是奇数,中间那个人可以任意;

  3、从左到中间那个人,身高需保证递增,如果用H表示新队形的高度,则H[1] < H[2] < H[3] .... < H[mid]。



  现在吉哥想知道:最多能选出多少人组成完美队形?
 
Input
  第一行输入T,表示总共有T组数据(T <= 20);

  每组数据先输入原先队形的人数n(1<=n <= 200),接下来一行输入n个整数,表示按顺序从左到右原先队形位置站的人的身高(50 <= h <= 250,不排除特别矮小和高大的)。
 
Output
  请输出能组成完美队形的最多人数,每组数据输出占一行。
 
Sample Input
2
3
51 52 51
4
51 52 52 51
 
Sample Output
3
4
 
Source
 
Recommend
liuyiding



算法:LCIS 【最长公共上升子序列分析



正解思路:


主要路线:先输入数组 a[] , 然后另外用一个数组 b[] 反序存 a[], 然后求 a 和 b 的 LCIS

关键:如何解决 中间的人 在哪儿?如何避免 过多 重叠?


PS:由于是从 1 开始记录的,所以序列 a[i] 对应于序列 b[n-i+1]

用 dp[i][j] 表示 a 从 1到 i , b[] 从 1 到 j ,两段序列的 LCIS

解决方法:每次枚举只重叠一个人。。。重叠的人的编号 i 从 1 到 n-1【这样就一定保证了经过后面的处理,不会出现过多的重叠】
                  那么每次求 LCIS 时 a[] 对应从 1 到 i
                                                  b[] 对应从 1 到 n-i+1
                  那么每次更新 ans 最多只可能重叠了一人。【也就是题目中说的最终完美队列中有奇数个人】
                  然后再根据所求的最大 dp[i][j] 的编号 i 和 j 判断是否重叠了
                  如果对应于当前最大的 dp[i][j] :i < n-j+1 说明当前最高的人没有重叠, ans = max(ans, dp[i][j]*2)
                                                                   否则最高的人重叠了, ans = max(ans, dp[i][j]*2-1)

PS:下面贴上我的心酸 WA 史和思路历程


错误思想分析:

误区一:写着写着,发现没有测试出错误

昨晚比赛的时候,先是想着定义 ans 为全局变量,然后用 LCIS 返回在 b[] 中找到的最高的人的编号 index
然后再对应于它在 a[] 中的编号 index1 = n-index+1 ,看 a 中 index1 前面有没有一个这样高的人
如果有,那么说明最高的人没有重叠 ans = ans*2
如果 a[] 中 1 到 index1-1 中没有一个这样高的人那么 ans = ans*2-1.....不知道哪里错了,求数据求指导,附上WA了的代码




误区二:

今天早上看了下别人发的博客都是边求 LCIS 边更新 ans 的。
于是乎,就没有管重叠了多少人直接 LCIS(n,n)
void LCIS(int n)
{
ans = 1;
for(int i = 1; i <= n; i++)
{
int tmp = 0;
for(int j = 1; j <= n; j++) //想来有了下面的下标处理是不会出现过多重叠的,错的很严重。。。
{
if(a[i] > b[j])
tmp = max(tmp,dp[j]);
else if(a[i] == b[j])
dp[j] = max(dp[j],tmp+1);
if(i < (n-j+1)) //判断是否重叠
ans = max(ans, dp[j]*2);
else ans = max(ans, dp[j]*2-1);
}
}
}

这样的代码对于下面的数据样例:


1
5
1 3 2 1 3

正确的完美序列该是 1 3 1 或者 1 2 1 都可以,那么应该输出的是 3
但是求出最后的的 dp[] 的最大值确是3
序列: 1 3 2
1 3 
            3 1 23 1
然后到了下面判断重叠的情况的时候直接  dp[]*2-1 = 5 肯定是比原来的正确的 ans 要大的,然后就错了。。。

改正: 只要把第二重循环 j < n 改成 j < n-i+1 即可以过,保证最多只重叠一人

第一种思路开始想着错误的应该和这个差不多,但是实在是没有 找出错误样例,希望路过的大神指点下=_=!

LCIS模板:

PS: 不熟悉的,自己可以先做下 hdu 1424 看下:

const int maxn = 500+50;
int dp[maxn];
int a[maxn],b[maxn]; /*********************************************
求序列 a (长度为 n)和序列 b (长度为 m)的 LCIS
下标从 1 开始
返回 LCIS 长度
*********************************************/
int LCIS(int n,int m)
{
for(int i = 1; i <= n; i++)
{
int tmp = 0;
for(int j = 1; j <= m; j++)
{
if(a[i] > b[j])
{
tmp = max(tmp,dp[j]);
}
else if(a[i] == b[j])
dp[j] = max(dp[j],tmp+1);
}
} int ans = 0;
for(int i = 1; i <= m; i++)
ans = max(ans, dp[i]);
return ans;
}


code:

code1: 用一次 LCIS ,边求LCIS 边处理重叠的人,更新 ans
2013-08-04 13:58:02 Accepted 4512 15MS 232K 988 B C++ free斩

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std; const int maxn = 210;
int a[maxn], b[maxn];
int dp[maxn];
int n;
int ans; void LCIS(int n,int m)
{
ans = 1;
for(int i = 1; i <= n; i++)
{
int tmp = 0;
for(int j = 1; j <= (n-i+1); j++) //注意:最多中间只重复一人,方便后面的重叠处理
{
if(a[i] > b[j])
tmp = max(tmp,dp[j]);
else if(a[i] == b[j])
dp[j] = max(dp[j],tmp+1);
if(i < (n-j+1)) //判断是否重叠
ans = max(ans, dp[j]*2);
else ans = max(ans, dp[j]*2-1);
}
}
} int main()
{
int T;
scanf("%d", &T);
while(T--)
{
scanf("%d", &n);
for(int i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
b[n-i+1] = a[i];
}
memset(dp,0,sizeof(dp));
LCIS(n, n);
printf("%d\n", ans);
}
return 0;
}

code2: 枚举重叠的人进行 LCIS 2*N 次,思路一样效率稍微低一点,不过比较好理解

2013-08-04 10:40:36 Accepted 4512 125MS 232K 1064 B C++ free斩

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std; const int maxn = 210;
int a[maxn], b[maxn];
int dp[maxn];
int n; int LCIS(int n,int m)
{
memset(dp,0,sizeof(dp));
for(int i = 1; i <= n; i++)
{
int tmp = 0;
for(int j = 1; j <= m; j++)
{
if(a[i] > b[j])
tmp = max(tmp,dp[j]);
else if(a[i] == b[j])
dp[j] = max(dp[j],tmp+1);
}
}
int ans = 0;
for(int i = 1; i <= m; i++)
ans = max(ans, dp[i]);
return ans;
} int main()
{
int T;
scanf("%d", &T);
while(T--)
{
scanf("%d", &n);
for(int i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
b[n-i+1] = a[i];
}
int ans = 1;
for(int i = 1; i < n; i++)
{
ans = max(ans, 2*LCIS(i, n-i)); //不重叠
ans = max(ans, 2*LCIS(i+1, n-i)-1); //最坏恰好重叠一人
}
printf("%d\n", ans);
}
return 0;
}




hdu 4512 吉哥系列故事——完美队形I【LCIS经典应用】的更多相关文章

  1. HDU 4512 吉哥系列故事——完美队形(LCIS)

    Problem Description 吉哥这几天对队形比较感兴趣. 有一天,有n个人按顺序站在他的面前,他们的身高分别是h[1], h[2] ... h[n],吉哥希望从中挑出一些人,让这些人形成一 ...

  2. hdu 4512 吉哥系列故事——完美队形I LCIS

    题目链接 给一个数列, 在里面选出一些数组成一个对称的数列, 数的顺序不能打乱. 使得左半边是一个严格递增的数列, 右边递减, 并且a[i] = a[n-i+1]. 就是一个对称的LCIS.. #in ...

  3. HDU 4512——吉哥系列故事——完美队形I——————【LCIS应用】

    吉哥系列故事——完美队形I Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Tot ...

  4. hdu 4512 吉哥系列故事——完美队形I_LCIS

    题目链接 题意: 假设有n个人按顺序站在他的面前,他们的身高分别是h[1], h[2] ... h[n],吉哥希望从中挑出一些人,让这些人形成一个新的队形,新的队形若满足以下三点要    求,则就是新 ...

  5. hdu 4512 吉哥系列故事——完美队形I(最长公共上升自序加强版)

    首先要去学习最长公共上升子序列LCIS.然后是对n*n效率的算法进行优化,这里要注意的是能够求出来的序列中间能够有一个最高的.刚開始将输入的数组进行逆置,写下来发现这可能存在问题. 只是详细是什么也没 ...

  6. LCIS HDOJ 4512 吉哥系列故事——完美队形I

    题目传送门 题意:中文题面 分析:LCIS应用:设置b[]为a[]的反转,然后LCIS,若相等的是自己本身,则+1, 否则+2 代码: #include <cstdio> #include ...

  7. HDU 4513 吉哥系列故事――完美队形II(Manacher)

    题目链接:cid=70325#problem/V">[kuangbin带你飞]专题十六 KMP & 扩展KMP & Manacher V - 吉哥系列故事――完美队形I ...

  8. HDU 4513 吉哥系列故事――完美队形II

    http://acm.hdu.edu.cn/showproblem.php?pid=4513 吉哥系列故事——完美队形II Time Limit: 3000/1000 MS (Java/Others) ...

  9. HDU 4513 吉哥系列故事——完美队形II manacher

    吉哥系列故事——完美队形II Problem Description 吉哥又想出了一个新的完美队形游戏! 假设有n个人按顺序站在他的面前,他们的身高分别是h[1], h[2] ... h[n],吉哥希 ...

随机推荐

  1. jQuery国际化插件 jQuery.i18n.properties 【轻量级】

    jQuery.i18n.properties是一款轻量级的jQuery国际化插件,能实现Web前端的国际化. 国际化英文单词为:Internationalization,又称i18n,“i”为单词的第 ...

  2. python \uxxxx转中文,Python列表中的字典 \uxxxx转中文,

    import json a = [{u', u'roleFlag': 7}] print json.dumps(a).decode("unicode_escape") 输出结果是: ...

  3. Java并发编程(四):线程安全性

    一.定义 当多个线程访问某个类时,不管运行时环境采用何种调度方式或者这些进程将如何交替执行,并且在主调代码中不需要额外的同步或协同,这个类都能表现出正确的行为,那么就称这个类是线程安全的. 二.线程安 ...

  4. 两DD-WRT组建WDS设置

    今天把买了非常久的两台CISCO(已刷DD-WRT)路由器组成WDS,參考了三篇文章   DD-WRT v24-sp2的WDS中继设置 http://www.cnblogs.com/zjoch/arc ...

  5. iOS 图片的属性

    UIViewContentModeScaleToFill UIViewContentModeScaleAspectFit UIViewContentModeScaleAspectFill UIView ...

  6. asp.net写日志权限问题

    asp.net网站程序写不了日志,都是这个原因. 程序池的身份标识,设为内置的网络服务即可.

  7. PHP学习笔记(4)GD库画五角星

         <?php //加header头,不然浏览器乱码 header("content-type: image/png"); //创建画布资源 $img = imagec ...

  8. 小型web服务器thttpd的学习总结(上)

    1.软件的主要架构 软件的文件布局比较清晰,主要分为6个模块,主模块是thttpd.c文件,这个文件中包含了web server的主要逻辑,并调用了其他模块的函数.其他的5个模块都是单一的功能模块,之 ...

  9. Angular js ie 7,8 兼容性

    Angularjs  官网有云: 1)在html 里面 ,有ng-app 的标签里需要定义个id ,id='ng-app'; 2)ie 7及以下版本需要json2.js或json3.js,主要用来解析 ...

  10. javascript中的函数作用域和声明提前

    在一些类C的编程语言中,花括号内的每一段代码都具有各自作用域,并且变量在声明他们的代码段之外是不可见的,这个概念叫做块级作用域. javascript中没有块级作用域的概念,有的是函数作用域的概念:变 ...