首页 > 代码库 > bzoj3208: 花神的秒题计划Ⅰ

bzoj3208: 花神的秒题计划Ⅰ

3208: 花神的秒题计划Ⅰ

Time Limit: 16 Sec  Memory Limit: 128 MB
Submit: 323  Solved: 211
[Submit][Status]

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

Sample Output

24
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

Source

原创 Memphis

题解:
一看询问这么少肯定暴力搞。。。原来做顺治滑雪的时候写的是记忆化,现在发现暴力建图跑拓扑序求最长链貌似更直观。
代码:
  1 #include<cstdio>  2   3 #include<cstdlib>  4   5 #include<cmath>  6   7 #include<cstring>  8   9 #include<algorithm> 10  11 #include<iostream> 12  13 #include<vector> 14  15 #include<map> 16  17 #include<set> 18  19 #include<queue> 20  21 #include<string> 22  23 #define inf 1000000000 24  25 #define maxn 1000000+5 26  27 #define maxm 705 28  29 #define eps 1e-10 30  31 #define ll long long 32  33 #define pa pair<int,int> 34  35 #define for0(i,n) for(int i=0;i<=(n);i++) 36  37 #define for1(i,n) for(int i=1;i<=(n);i++) 38  39 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 40  41 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 42  43 #define mod 1000000007 44 #define num(i,j) ((i-1)*n+j) 45  46 using namespace std; 47  48 inline int read() 49  50 { 51  52     int x=0,f=1;char ch=getchar(); 53  54     while(ch<0||ch>9){if(ch==-)f=-1;ch=getchar();} 55  56     while(ch>=0&&ch<=9){x=10*x+ch-0;ch=getchar();} 57  58     return x*f; 59  60 } 61 struct edge{int go,next;}e[4*maxn]; 62 int n,tot,a[maxm][maxm],head[maxn],inp[maxn],d[maxn]; 63 queue<int>q; 64 bool flag[maxn]; 65 const int dx[4]={0,0,1,-1}; 66 const int dy[4]={1,-1,0,0}; 67 inline void insert(int x,int y) 68 { 69     e[++tot].go=y;e[tot].next=head[x];head[x]=tot; 70 } 71 void work() 72 { 73     memset(d,0,sizeof(d)); 74     memset(inp,0,sizeof(inp)); 75     memset(head,0,sizeof(head));tot=0; 76     for1(i,n) 77      for1(j,n) 78       if(!flag[num(i,j)]) 79        for0(k,3) 80         { 81           int ii=i+dx[k],jj=j+dy[k]; 82           if(ii<1||ii>n||jj<1||jj>n||flag[num(ii,jj)]||a[ii][jj]>=a[i][j])continue; 83           insert(num(i,j),num(ii,jj));inp[num(ii,jj)]++; 84         } 85     for1(i,n)for1(j,n)if(!flag[num(i,j)]&&!inp[num(i,j)])q.push(num(i,j)),d[num(i,j)]=1; 86     while(!q.empty()) 87     { 88         int x=q.front();q.pop(); 89         for(int i=head[x];i;i=e[i].next) 90         { 91             int y=e[i].go; 92             inp[y]--; 93             d[y]=max(d[y],d[x]+1); 94             if(!inp[y])q.push(y); 95         } 96     } 97     int ans=0; 98     for1(i,n*n)ans=max(ans,d[i]); 99     printf("%d\n",ans);100 }101 102 int main()103 104 {105 106     freopen("input.txt","r",stdin);107 108     freopen("output.txt","w",stdout);109 110     n=read();111     for1(i,n)for1(j,n)a[i][j]=read();112     int m=read();char ch[4];113     while(m--)114     {115         scanf("%s",ch);116         if(ch[0]==Q)work();117         else if(ch[0]==C){int x=read(),y=read();a[x][y]=read();}118         else 119             {120                 int x1=read(),y1=read(),x2=read(),y2=read();121                 for2(i,x1,x2)122                  for2(j,y1,y2)123                   flag[num(i,j)]=ch[0]==S?1:0;124             }125     }126 127     return 0;128 129 } 
View Code

 

bzoj3208: 花神的秒题计划Ⅰ