NOIP 2010 提高组 复赛 第二题 乌龟棋 tortoise AC代码(多种状态的动态规划dp)+20分代码(深搜dfs)

NOIP 2010 提高组 复赛 第二题  乌龟棋  tortoise  AC代码(多种状态的动态规划dp)+20分代码(深搜dfs)

总目录详见：NOIP 提高组 复赛 试题 目录 信奥 历年

在线测评地址：https://www.luogu.com.cn/problem/P1541

1.AC代码(多种状态的动态规划dp)

部分内容(有改动)摘自https://www.luogu.com.cn/problem/solution/P1541

1.f[a][b][c][d]:表示你出了a张爬行牌1，b张爬行牌2，c张爬行牌3，d张爬行牌4时的得分

2.cnt[x]:表示牌x一共有多少张

题干中说如何出牌,那我们就不妨DP一下每一种牌的出牌张数

初始化：

f[0][0][0][0]=s[1];

 显然,乌龟开始时在s[1],题中说乌龟棋子自动获得起点格子的分数,故未出牌时（f[0][0][0][0]）分数为s[1]

之后边输入边存每一种牌的张数（输入数据第3行:M个整数，b1b2……bM，表示M张爬行卡片上的数字,故卡1~卡4张数一定）:

for(int i=1;i<=m;i++){
	scanf("%d",&x);
	cnt[x]++;
}

之后便可以开始DP了:

起始状态f[0][0][0][0]=s[1]，即不出任何爬行卡;之后对于每一张卡片，我都可以选择放与不放， 设当前放的卡1数量为a，卡2数量为b，卡3数量为c，卡4数量为d（以下出现a~d均为这个意思），则对于卡1:

比较卡1的放与不放,只需决策卡1的放与不放，即取f[a][b][c][d],f[a-1][b][c][d]+num[r]的最大值。又由于a有一定数量,所以我们可以得出关于a的转移方程:

	f[a][b][c][d]=max(f[a][b][c][d],f[a-1][b][c][d]+s[r])

其中r=1+a+b2+c3+d4(至于r在a+b2+c3+d4加1原因见后)

DP 数量a：

for(int a=0;a<=g[1];a++){
	int r=1+a+b*2+c*3+d*4;
	if(a!=0) f[a][b][c][d]=max(f[a][b][c][d],F[a-1][b][c][d]+num[1+a+b*2+c*3+d*4])
}

这不就是个“物品占的空间”为1，“价值”为s[r]的多重背包嘛！！

至于这个(a!=0)，显然，你要是调用f[a-1][b][c][d]，肯定得保证a-1>=0吧。a显然作为卡1个数不可能<0,故取(a!=0)即可

根据多维背包的思想，背包DP几个“价值”（即爬行牌种类）开几维即可,故 转移方程为:

		f[a][b][c][d]=max(f[a-1][b][c][d],f[a][b-1][c][d],f[a][b][c-1][d],f[a][b][c][d-1])+s[1+1*a+2*b+3*c+4*d]

最后DP出来的f[cnt[1]][cnt[2]][cnt[3]][cnt[4]]即为答案。

DP代码如下：

	for(int a=0;a<=cnt[1];a++)
		for(int b=0;b<=cnt[2];b++)
			for(int c=0;c<=cnt[3];c++)
				for(int d=0;d<=cnt[4];d++){
					int r=1+a+b*2+c*3+d*4;//千万千万别忘了加一,因为乌龟从s[1]出发,设前进i步，则到达s[i+1],我就是因为这调了一个小时死活没找出毛病 
					if(a!=0)	f[a][b][c][d]=max(f[a][b][c][d],f[a-1][b][c][d]+s[r]); //a!=0原因见上
					if(b!=0)    f[a][b][c][d]=max(f[a][b][c][d],f[a][b-1][c][d]+s[r]);
					if(c!=0)    f[a][b][c][d]=max(f[a][b][c][d],f[a][b][c-1][d]+s[r]);
					if(d!=0)	f[a][b][c][d]=max(f[a][b][c][d],f[a][b][c][d-1]+s[r]);
				}	

在for循环中将f[a][b][c][d]与f[a-1][b][c][d]+s[1+1a+2b+3c+4d],f[a][b-1][c][d]+s[1+1a+2b+3c+4d],f[a][b][c-1][d]+s[1+1a+2b+3c+4d],f[a][b][c][d-1]+s[1+1a+2b+3c+4d]逐个比较，实现了转移方程的更新(c++中max函数貌似只能比较两个数)
当然，我们知道总牌数n和每种卡的张数，则实际写代码时我们完全可以考虑减一维,可我至今还是没想起来。。。。。。（我还是个蒟蒻）​

AC代码如下：

#include <bits/stdc++.h>
using namespace std;
int s[355],cnt[5],n,m,f[45][45][45][45];
int main(){
	int i,x,a,b,c,d;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)scanf("%d",&s[i]);
	for(i=1;i<=m;i++){
		scanf("%d",&x);
		cnt[x]++;
	}
	f[0][0][0][0]=s[1];
	for(a=0;a<=cnt[1];a++)
		for(b=0;b<=cnt[2];b++)
			for(c=0;c<=cnt[3];c++)
				for(d=0;d<=cnt[4];d++){
					if(a)f[a][b][c][d]=max(f[a][b][c][d],f[a-1][b][c][d]+s[1+a*1+b*2+c*3+d*4]);
					if(b)f[a][b][c][d]=max(f[a][b][c][d],f[a][b-1][c][d]+s[1+a*1+b*2+c*3+d*4]);
					if(c)f[a][b][c][d]=max(f[a][b][c][d],f[a][b][c-1][d]+s[1+a*1+b*2+c*3+d*4]);
					if(d)f[a][b][c][d]=max(f[a][b][c][d],f[a][b][c][d-1]+s[1+a*1+b*2+c*3+d*4]);
				}
	printf("%d\n",f[cnt[1]][cnt[2]][cnt[3]][cnt[4]]);
	return 0;
}

2.20分代码(深搜dfs)

12!=479001600=5*10^8

超时是必然，试试看，能得几分。

样例模拟如下：

输入：
9 5
6 10 14 2 8 8 18 5 17
1 3 1 2 1
输出：
73

存储牌的数组的位置1 2 3 4 5
存储牌的数组的数值1 3 1 2 1

关于牌的位置的选择，对应5的排列情况
1 2 3 4 5
1 2 3 5 4
1 2 5 3 4
1 2 5 4 3
......

这正是深搜(dfs)所擅长的。


小明使用爬行卡片顺序为1,1,3,1,2得到的分数为6+10+14+8+18+17=73。注意，由于起点是1，所以自动获得第1格的分数6。

对应牌的位置的选择
1 3 2 5 4

即
b[1]=1 b[3]=1 b[2]=3 b[5]=1 b[4]=2

20分代码(深搜dfs) 代码如下：

#include <bits/stdc++.h>
using namespace std;
int n,m,mx;
int a[35],b[15],c[15],vis[15];
void dfs(int step){
	int i,j,now,val;
	if(step==m+1){
		now=1;
		val=a[1];
		for(j=1;j<=m;j++){
			now+=c[j];
			val+=a[now];
		}
		mx=max(mx,val);
		return ;
	}
	for(i=1;i<=m;i++)
		if(!vis[i]){
			vis[i]=1;
			c[step]=b[i];
			dfs(step+1);
			vis[i]=0;
		}
}
int main(){
	int i;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)scanf("%d",&a[i]);
	for(i=1;i<=m;i++)scanf("%d",&b[i]);
	dfs(1);
	printf("%d\n",mx);
	return 0;
}