当前位置: 代码迷 >> 综合 >> 麻烦的DP-BZOJ-1605-[Usaco2008 Open]Crisis on the Farm 牧场危机
  详细解决方案

麻烦的DP-BZOJ-1605-[Usaco2008 Open]Crisis on the Farm 牧场危机

热度:67   发布时间:2023-12-20 21:20:18.0

Description
约翰和他的奶牛组建了一只乐队“后街奶牛”,现在他们正在牧场里排练.奶牛们分成一堆一堆,共N(1≤N≤1000)堆.每一堆里,30只奶牛一只踩在另一只的背上,叠成一座牛塔.牧场里还有M(1≤M≤1000)个高高的草垛. 作为出色的指挥家,约翰可以通过口哨指挥奶牛们移动.他的口哨有四个音,分别能使所有的牛塔向东南西北四个方向移动一格.
每一次,当一个牛塔到达了一个草垛所在的格子,牛塔最上方的奶牛就会跳到草垛上,而且不再下来,而其他奶牛仍然呈塔状站在草垛所在的格子里.当牛塔只剩一只奶牛时,这只奶牛也会跳到草垛上. 突然,约翰大惊失色:原来邻家的奶缸爆炸了!滚滚而下的牛奶正朝着约翰的牧场冲来,不久就要将牧场淹没.约翰必须马上行动,用口哨声挽救奶牛们的生命.他要指挥奶牛尽量多地跳上草垛,草垛上的奶牛将不会被淹死. 约翰还有K次吹口哨的机会.那他最多还能救多少奶牛呢?请计算最多能挽救的奶牛数,以及达到这个数目约翰需要吹的口哨调子序列.序列用E,W,S,N表示东西南北.如果有多种序列能达到
要求,输出作为字符串最小的.

Input
第1行输入三个整数N,M,K,之后N行每行输入一对整数(Xi,Yi)表示一座牛塔所在的位置,1<=K<=30
之后M行每行输入一对整数(Xi,Yi)表示一个草垛所在的位置.1≤Xi≤1000;1≤Yi≤1000.

Output
第1行输出最多能挽救的奶牛数.第2行输出口哨调子序列.

Sample Input
3 6 3

3 4

6 2

5 7

8 2

9 2

6 4

5 4

6 7

8 7

Sample Output
Use the ‘east’ whistle three times, at which point the milk floods

the area. Each haystack ends up saving 1 cow.

6

EEE


题解:
这道题真的是超级麻烦,思想上非常好想,就是一个DP加方案的输出,但是这题无论是预处理还是最后方案输出都比较麻烦,WA了很久都是方案输出的问题。
首先要注意的是K是在30以内的,所以不用去考虑那个牛塔没牛的情况,那么我们先预处理一下在移动delta(x),delta(y)的情况下,能够有多少的牛到草棚上,记为get[x][y]。
然后dp数组为3维的,分别为x,y,k,也就是说dp[i][j][t]是在第t次移动后,总共移动为i,j的最大权值(即最大获救牛数量)。所以可以很轻松的写出状态转移方程:dp[i][j][t]=max(dp[i+1][j][t-1],dp[i+1][j][t-1],dp[i][j-1][t-1],dp[i][j+1][t-1])+get[i][j]。那么就能轻松得到最大权值了。
接下来是讨厌的方案输出,由于非常容易有多个解,此题要求输出字典序最小的方案,那么我们可以知道东南西北的英文缩写的字典序为E,N,S,W。由于我们要倒着来求方案,所以也按着W,S,N,E的顺序来规划局部方案,这样最后得到的就是字典序最小的了。
有些地方没做处理,最后跑了48ms,还算勉强。


#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <vector>
#include <algorithm>
#include <cmath>
#include <queue>
#include <set>
using namespace std;
int n,m,k;
typedef struct node
{int x,y;
} node;
node cow[1005];
bool grass[1005][1005];
int  get[105][105];
int dp[105][105][35];
char Turn(int x,int y)
{if(x==-1)return 'W';if(x==1)return 'E';if(y==-1)return 'S';return 'N';
}
int moved[4][2]= {
   {
   1,0},{
   0,1},{
   0,-1},{-1,0}};
bool check(int x,int y)
{if(x<0 || x>1000 || y<0 || y>1000)return 0;return 1;
}
char pre[105][105][35];
char nn[4]= {
   'E','N','S','W'};
int main()
{int x,y;cin >> n >> m >> k;int fix=k+1;for(int i=0; i<n; i++)scanf("%d%d",&cow[i].x,&cow[i].y);for(int i=0; i<m; i++){scanf("%d%d",&x,&y);grass[x][y]=1;}for(int i=1; i<=2*k+1; i++)for(int j=1; j<=2*k+1; j++){if(abs(i-fix)+abs(j-fix)>k)continue;for(int z=0; z<n; z++){x=cow[z].x+i-fix,y=cow[z].y+j-fix;if(!check(x,y))continue;if(grass[x][y])get[i][j]++;}}int out=0;for(int t=1; t<=k; t++)for(int i=1; i<=2*k+1; i++)for(int j=1; j<=2*k+1; j++){int tmp=abs(i-fix)+abs(j-fix);if(tmp>t)continue;if(tmp==t || (t-tmp)%2==0){dp[i][j][t]=max(max(dp[i-1][j][t-1],dp[i+1][j][t-1]),max(dp[i][j-1][t-1],dp[i][j+1][t-1]))+get[i][j];if(t==k && dp[i][j][t]>=out)out=dp[i][j][t];}}for(int i=1; i<=2*k+1; i++)for(int j=1; j<=2*k+1; j++)if(dp[i][j][k]==out)pre[i][j][k]='F';cout << out << endl;for(int t=k-1; t>=0; t--)for(int i=1; i<=2*k+1; i++)for(int j=1; j<=2*k+1; j++){int tmp=abs(i-fix)+abs(j-fix);if(tmp>t)continue;for(int z=3; z>=0; z--)if(dp[i][j][t]+get[i+moved[z][0]][j+moved[z][1]]==dp[i+moved[z][0]][j+moved[z][1]][t+1] && pre[i+moved[z][0]][j+moved[z][1]][t+1]!=0)pre[i][j][t]=nn[z];}int nowx=fix,nowy=fix;for(int i=0; i<k; i++){printf("%c",pre[nowx][nowy][i]);switch(pre[nowx][nowy][i]){case 'E':nowx++;break;case 'N':nowy++;break;case 'S':nowy--;break;case 'W':nowx--;break;}}return 0;
}
  相关解决方案