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
",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