免费馅饼(HDU1176) 背包问题

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可能会超时。

Sample Input

6
5 1
4 1
6 1
7 2
7 2
8 3
0

Sample Output

4

思路：

这是一个背包问题，在第i秒第j个位置能接到最多的饼数=第i-1秒的第j-1,j,j+1的位置中所能接到的最多的位置加上第i秒第j个位置要下落的饼数。
用 dp[time][position] 数组来保存状态,用 cnt[time][position] 数组来记录在time秒position位置要下落的饼数。

那么状态转移方程为:
dp[time][pos]=max(dp[time-1][pos-1],dp[time-1][pos],dp[time-1][pos+1])+cnt[time][pos];

不过这道题由于初始位置在5，且每次只能移动一格，所以前5秒的状态转移方程要单独计算，并且0和10的边界位置的状态转移方程也要单独计算。

下面是AC代码：

#include <bits/stdc++.h>
#define mem(a,n) memset(a,n,sizeof(a))
#define rd1(x) scanf("%d",&x)
#define rd2(x,y) scanf("%d%d",&x,&y)
#define prt1(x) printf("%d\n",x)
typedef unsigned long long ull;
using namespace std;

int cnt[100005][12],dp[100005][12];

int main()
{
    int n;
    while(rd1(n)!=EOF&&n)
    {
        int time,pos,maxtime=1,ans=0;
        mem(cnt,0);
        mem(dp,0);
        for(int i=0;i<n;i++)
        {
            rd2(pos,time);
            //记录最长的时间
            if(time>maxtime) maxtime=time;
            cnt[time][pos]++;
        }
        for(int i=1;i<=maxtime;i++)
        {
            //由于初始位置在5，要对时间进行分类
            if(i>5)
            {
                dp[i][0]=max(dp[i-1][0],dp[i-1][1])+cnt[i][0];
                dp[i][10]=max(dp[i-1][10],dp[i-1][9])+cnt[i][10];
                for(int j=1;j<10;j++)
                {
                    dp[i][j]=max(dp[i-1][j+1],max(dp[i-1][j-1],dp[i-1][j]))+cnt[i][j];
                }
            }
            else
            {
                //第i秒时左右各到达最远为5-i，5+i处
                //左右边界各两个需要单独计算
                dp[i][5-i]=dp[i-1][5-i+1]+cnt[i][5-i];//左1
                dp[i][5-i+1]=max(dp[i-1][5-i+1],dp[i-1][5-i+2])+cnt[i][5-i+1];//左2
                dp[i][5+i-1]=max(dp[i-1][5+i-1],dp[i-1][5+i-2])+cnt[i][5+i-1];//右1
                dp[i][5+i]=dp[i-1][5+i-1]+cnt[i][5+i];//右2
                for(int j=5-i+2;j<5+i-1;j++)
                {
                    dp[i][j]=max(dp[i-1][j+1],max(dp[i-1][j-1],dp[i-1][j]))+cnt[i][j];
                }
            }
        }
        for(int i=0;i<11;i++)
        {
            if(dp[maxtime][i]>ans) ans=dp[maxtime][i];
        }
        prt1(ans);
    }
    return 0;
}