BZOJ3208:花神的秒题计划Ⅰ(记忆化搜索DP)
Description
背景【backboard】:
Memphis等一群蒟蒻出题中,花神凑过来秒题……
描述【discribe】:
花花山峰峦起伏,峰顶常年被雪,Memphis打算帮花花山风景区的人员开发一个滑雪项目。
我们可以把风景区看作一个n*n的地图,每个点有它的初始高度,滑雪只能从高处往低处滑【严格大于】。但是由于地势经常变动【比如雪崩、滑坡】,高度经常变化;同时,政府政策规定对于每个区域都要间歇地进行保护,防止环境破坏。现在,滑雪项目的要求是给出每个n*n个点的初始高度,并给出m个命令,C a b c表示坐标为a,b的点的高度改为c;S a b c d表示左上角为a,b右下角为c,d的矩形地区开始进行保护,即不能继续滑雪;B a b c d表示左上角为a b,右下角为c d的矩形地区取消保护,即可以开始滑雪;Q表示询问现在该风景区可以滑雪的最长路径为多少。对于每个Q要作一次回答。
花神一看,这不是超简单!立刻秒出了标算~
Input
第一行n,第二行开始n*n的地图,意义如上;接下来一个m,然后是m个命令,如上
Output
对于每一个Q输出单独一行的回答
Sample Input
5
1 2 3 4 5
10 9 8 7 6
11 12 13 14 15
20 19 18 17 16
21 22 23 24 25
5
C 1 1 3
Q
S 1 3 5 5
S 3 1 5 5
Q
1 2 3 4 5
10 9 8 7 6
11 12 13 14 15
20 19 18 17 16
21 22 23 24 25
5
C 1 1 3
Q
S 1 3 5 5
S 3 1 5 5
Q
Sample Output
24
3
样例解释:
第一个Q路线为:25->24->23->22….->3->2
第二个Q的路线为:10->9->2
3
样例解释:
第一个Q路线为:25->24->23->22….->3->2
第二个Q的路线为:10->9->2
HINT
100%的数据:1<=n<=700;1<=m<=1000000;其中Q、S、B操作总和<=100;
题中所有数据不超过2*10^9
Solution
暴力可过 心态崩了
Code
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N (709) 5 using namespace std; 6 7 int n,m,x,y,z,c,ans,a[N][N],f[N][N],lim[N][N]; 8 int dx[7]={0,1,-1,0,0},dy[7]={0,0,0,1,-1}; 9 char opt[2]; 10 11 int DP(int x,int y) 12 { 13 if (lim[x][y]) return -1; 14 if (f[x][y]!=-1) return f[x][y]; 15 f[x][y]=1; 16 for (int i=1; i<=4; ++i) 17 { 18 int xx=x+dx[i], yy=y+dy[i]; 19 if (xx<1 || xx>n || yy<1 || yy>n || a[x][y]<=a[xx][yy]) continue; 20 f[x][y]=max(f[x][y],DP(xx,yy)+1); 21 } 22 return f[x][y]; 23 } 24 25 void Update(int x,int y,int xx,int yy,int opt) 26 { 27 for (int i=x; i<=xx; ++i) 28 for (int j=y; j<=yy; ++j) lim[i][j]=opt; 29 } 30 31 inline int read() 32 { 33 int x=0; char ch=getchar(); 34 while (ch<'0' || ch>'9') ch=getchar(); 35 while (ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar(); 36 return x; 37 } 38 39 int main() 40 { 41 n=read(); 42 for (int i=1; i<=n; ++i) 43 for (int j=1; j<=n; ++j) 44 a[i][j]=read(); 45 m=read(); 46 while (m--) 47 { 48 scanf("%s",opt); 49 if (opt[0]!='Q') x=read(),y=read(),z=read(); 50 if (opt[0]=='C') a[x][y]=z; 51 if (opt[0]=='S') c=read(),Update(x,y,z,c,1); 52 if (opt[0]=='B') c=read(),Update(x,y,z,c,0); 53 if (opt[0]=='Q') 54 { 55 memset(f,-1,sizeof(f)); ans=0; 56 for (int i=1; i<=n; ++i) 57 for (int j=1; j<=n; ++j) 58 ans=max(ans,DP(i,j)); 59 printf("%d\n",ans); 60 } 61 } 62 }