hdu1176 免费馅饼动态规划二维数组实现

免费馅饼

Time Limit: 1000MS Memory Limit: 32768KB

Submit Statistic Discuss

Problem Description

都说天上不会掉馅饼，但有一天gameboy正走在回家的小径上，忽然天上掉下大把大把的馅饼。说来gameboy的人品实在是太好了，这馅饼别处都不掉，就掉落在他身旁的10米范围内。馅饼如果掉在了地上当然就不能吃了，所以gameboy马上卸下身上的背包去接。但由于小径两侧都不能站人，所以他只能在小径上接。由于gameboy平时老呆在房间里玩游戏，虽然在游戏中是个身手敏捷的高手，但在现实中运动神经特别迟钝，每秒种只有在移动不超过一米的范围内接住坠落的馅饼。现在给这条小径如图标上坐标：

为了使问题简化，假设在接下来的一段时间里，馅饼都掉落在0-10这11个位置。开始时gameboy站在5这个位置，因此在第一秒，他只能接到4,5,6这三个位置中期中一个位置上的馅饼。问gameboy最多可能接到多少个馅饼？（假设他的背包可以容纳无穷多个馅饼）

Input

输入数据有多组。每组数据的第一行为以正整数n(0 < n < 100000)，表示有n个馅饼掉在这条小径上。在结下来的n行中，每行有两个整数x,T(0 < T < 100000),表示在第T秒有一个馅饼掉在x点上。同一秒钟在同一点上可能掉下多个馅饼。n=0时输入结束。

Output

每一组输入数据对应一行输出。输出一个整数m，表示gameboy最多可能接到m个馅饼。

提示：本题的输入数据量比较大，建议用scanf读入，用cin可能会超时。

Example Input

Example Output

一开始写的超时代码：

/*
先按照时间顺序优先到后以及落得位置有小到大排好序，然后再依次进行比较
*/
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
struct st
{int p,t;
}a[100005],temp;
int cmp(st a,st b)
{return a.t<b.t;
}
int main()
{int t,q=0;while(scanf("%d",&t),t){q=0;for(int i=0;i<t;i++)scanf("%d%d",&a[i].p,&a[i].t);stable_sort(a,a+t,cmp);int l=0,_max=0,dp[100000]={0},k=0;for(int i=0;i<t-1;++i){for(int j=i+1;j<t&&i!=j;++j){if((a[j].p-a[i].p)>=-1&&(a[j].p-a[i].p)<=1&&(a[j].t-a[i].t)>=-1&&(a[j].t-a[i].t)<=1){k=1;dp[j]=dp[i]+1;}}}for(int i=0;i<t;i++)if(_max<dp[i])_max=dp[i];printf("%d\n",_max+k);}return 0;
}
后来用的是二维数组来写。

设a[i][j]为第i秒的j位置掉下的馅饼数量，f[i][j]为第i秒在j位置接馅饼最多可以接到的最多馅饼数量。由于每秒只能移动一个位置，因此这一状态可能由三种情况达到：
f[i - 1][j - 1]
f[i - 1][j]
f[i - 1][j + 1]
这三种情况中的最大值加上当前位置可以接到的馅饼数即是当前位置可以接到的最大馅饼数量：
f [ i ] [ j ] = max ( f [ i - 1 ] [ j - 1 ] , f [ i - 1 ] [ j ] , f [ i - 1 ] [ j + 1 ] ) + a [ i ] [ j ] ;
可以看出，当前状态与之前同一阶段的多个状态有关，而类似于 Ugly Numbers 等动态规划中的每一阶段都只有一个状态，为区分两者，我将它称为 二维动态规划 。当然，这道题只是二维动态规划中最简单的一种罢了。

#include<stdio.h>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int dp[100005][12];
int main()
{int n,i,j,maxt;int x,t;while(scanf("%d",&n),n){maxt=0;memset(dp,0,sizeof(dp));for(i=0;i<n;i++){scanf("%d%d",&x,&t);dp[t][x]++;if(maxt<t)  maxt=t;}for(i=maxt-1;i>=0;i--){dp[i][0]+=max(dp[i+1][1],dp[i+1][0]);for(j=1;j<11;j++){dp[i][j]+=max(max(dp[i+1][j-1],dp[i+1][j]),dp[i+1][j+1]);}}printf("%d\n",dp[0][5]);}return 0;
}