【DP】钓鱼

【问题描述】  
  
　　约翰是个垂钓谜，星期天他决定外出钓鱼h小时（1≤h≤16），约翰家附近共有n个池塘（2≤n≤25），这些池塘分布在一条直线上，约翰将这些池塘按离家的距离由近到远编上号，依次为L1，L2，…，Ln，约翰家门外就是第一个池塘，所以他到第一个池塘是不用花时间的。 


　　约翰可以任选若干个池塘由近到远地垂钓，并且在每个池塘他都可以呆上任意长的时间，但呆的时间必须为5分钟的倍数（即5分钟为一个单位时间），已知从池塘Li到池塘Li+1要化去约翰ti个单位时间。每个池塘的上鱼率预先也是已知的，池塘Li在第一个单位时间内能钓到的鱼为Fi（0≤Fi≤100），并且每当他在某一个鱼塘呆上一个单位时间后，该鱼塘单位时间内能钓到的鱼将减少一个常数di（0≤di≤100）。 


　　现在请你编一个程序计算约翰最多能钓到多少鱼。
 
    
 【输入格式】  
  
　　第一行为一个整数n，第二行为一个整数h，第三行为n个用空格隔开的整数，表示Fi（i=1，2，…，n），第四行为n个用空格隔开的整数，表示di（i=1，2，…，n），第五行为n-1个用空格隔开的整数，表示ti（i=1，2，…，n-1）
 
    
 【输出格式】  
   
　　一个整数，表示约翰最多能钓到鱼的数量。
 
    
 【输入样例】   
   
2
1
10 1
2 5
2


 
    
 【输出样例】  
   
31


 
    
 【数据范围】  
   
2≤n≤25 , 1≤h≤16 , 0≤ti≤100 , 0≤Fi≤100 , 0≤di≤100

最先想的是二维设状态函数f(i,j)表示在i个时段内,在前j个鱼塘能够钓的最多鱼数。

但是题目中每个鱼塘的每个时段所钓的鱼是根据约翰在这个鱼塘呆的时间所变化的，如果单纯地用两维来描述不能清楚的体现鱼塘中单位分钟能钓的鱼随着钓的时间的变化。

所以应该增加一维

设状态函数f(i,j,k)表示约翰在i个时段内,在第j个鱼塘呆了k个时段的最大鱼数。

考虑第i个时段,显然约翰面临着两个选择

1.继续在此鱼塘钓鱼 f(i+1,j,k+1)(下一个时段依然在这个鱼塘,在这个鱼塘呆的时间为k+1)=max(f(i+1,j,k+1),f(i,j,k)+当前该鱼塘可钓的鱼数p);

2.换鱼塘 f(i+t[j+1],j+1,0)(在下一个鱼塘呆的时间是0)=max(f(i+t[j+1],j+1,0),f(i,j,k))

边界f(0,1,0)=0;(约翰一开始在第一个鱼塘,还没有开始钓鱼)

其余为-oo;(约翰不能从除了第一个鱼塘以外的地方开始)

注意一个小时有12个时段。

#include<cstdio>
#include<cstdlib>
#include<cctype>
#include<cstring>
#include<algorithm>
#define oo 200000010
using namespace std;
const int maxn=30;
int n,h;
int F[maxn],d[maxn],t[maxn];
int f[16*12+5][maxn][16*12+5];
/*
	f(i,j,k)=在i个单位时间内,在第j个鱼塘待了k个单位时间所能够钓到的最大鱼数
	f(i+1,j,k+1)=max(f(i+1,j,k+1),f(i,j,k)+ff[i]);
	f(i+t[j+1],j+1,0)=max(f(i+t[j+1],j+1,0),f(i,j,k));
*/

int ff(int i,int k)//计算第i号鱼塘当前可钓的鱼数 
{
	int t=F[i];
	t=t-k*d[i];
	return t;
}

int main()
{
	//freopen("my.in","r",stdin);
	//freopen("my.out","w",stdout);
	
	scanf("%d%d",&n,&h);
	for(int i=1;i<=n;i++)scanf("%d",&F[i]);
	for(int i=1;i<=n;i++)scanf("%d",&d[i]);
	for(int i=2;i<=n;i++)scanf("%d",&t[i]);
	
	for(int i=0;i<=12*h;i++)
	for(int j=1;j<=n;j++)
	for(int k=0;k<=12*h;k++)
	f[i][j][k]=-oo;
	
	f[0][1][0]=0;
	
	int ans=0;
	for(int i=0;i<=12*h;i++)
	for(int j=1;j<=n;j++)
	for(int k=0;k<=12*h;k++)
	{
		int p=ff(j,k),t1,t2;
		t1=t2=0;
		if(i+1<=12*h && k+1<=12*h)
		{
			t1=max(f[i+1][j][k+1],f[i][j][k]+p);
			f[i+1][j][k+1]=t1;
		}
		if(i+t[j+1]<=12*h && j+1<=n)
		{
			t2=max(f[i+t[j+1]][j+1][0],f[i][j][k]);
			f[i+t[j+1]][j+1][0]=t2;
		}

		ans=max(ans,f[i][j][k]);
		
		//if(i+t[j+1]<=12*h && j+1<=n && i+1<=12*h && k+1<=12*h)ans=max(f[i+1][j][k+1],f[i+t[j+1]][j+1][0]);
	}
	
		
	printf("%d\n",ans);
	return 0;
}