BZOJ 3208: 花神的秒题计划Ⅰ

3208: 花神的秒题计划Ⅰ

Time Limit: 16 Sec  Memory Limit: 128 MB
Submit: 695  Solved: 474
[Submit][Status][Discuss]

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

[Submit][Status][Discuss]

看到数据范围都不相信自己眼睛，额，花神系列什么时候出了这么一道水题？

完完全全的暴力，每次暴力修改高度，暴力标记保护，询问就BFS一遍算答案。

#include <cstdio>
#include <cstring>

const int siz = 705;

int n, m;
int h[siz][siz];
int v[siz][siz];
int c[siz][siz];
int f[siz][siz];

const int mv[4][2] =
{
    {0, +1},
    {0, -1},
    {-1, 0},
    {+1, 0}
};

inline int path(void)
{
    static int que[siz * siz][2], hd, tl;
    
    memset(f, 0, sizeof(f));
    memset(c, 0, sizeof(c));
    
    hd = 0, tl = 0;
    
    int ans = 0;
    
    for (int i = 1; i <= n; ++i)
        for (int j = 1; j <= n; ++j)
            if (!v[i][j])
            {
                for (int k = 0; k < 4; ++k)
                {
                    int x = i + mv[k][0];
                    int y = j + mv[k][1];
                    
                    if (x < 1 || x > n)continue;
                    if (y < 1 || y > n)continue;
                    
                    if (!v[x][y] && h[x][y] > h[i][j])
                        ++c[i][j];
                }
                
                if (!c[i][j])
                {
                    que[tl][0] = i;
                    que[tl][1] = j;
                    f[i][j] = 1;
                    ++tl;
                }
            }
    
    while (hd != tl)
    {
        int x = que[hd][0];
        int y = que[hd][1];
        ++hd;
        
        if (ans < f[x][y])
            ans = f[x][y];
        
        for (int k = 0; k < 4; ++k)
        {
            int i = x + mv[k][0];
            int j = y + mv[k][1];
            
            if (i < 1 || i > n)continue;
            if (j < 1 || j > n)continue;
            
            if (!v[i][j] && h[x][y] > h[i][j])
            {
                if (f[i][j] < f[x][y] + 1)
                    f[i][j] = f[x][y] + 1;
                
                if (--c[i][j] == 0)
                {
                    que[tl][0] = i;
                    que[tl][1] = j;
                    ++tl;
                }
            }
        }
    }
    
    return ans;
}

signed main(void)
{
    scanf("%d", &n);
    
    for (int i = 1; i <= n; ++i)
        for (int j = 1; j <= n; ++j)
            scanf("%d", &h[i][j]);
            
    scanf("%d", &m);
    
    while (m--)
    {
        static int a, b, c, d;
        
        static char s[5];
        
        scanf("%s", s);
        
        switch(s[0])
        {
            case 'Q':
                printf("%d
", path()); 
                break;
            case 'C':
                scanf("%d%d%d", &a, &b, &c);
                h[a][b] = c;
                break;
            case 'S':
                scanf("%d%d%d%d", &a, &b, &c, &d);
                for (int i = a; i <= c; ++i)
                    for (int j = b; j <= d; ++j)
                        v[i][j] = true;
                break;
            case 'B':
                scanf("%d%d%d%d", &a, &b, &c, &d);
                for (int i = a; i <= c; ++i)
                    for (int j = b; j <= d; ++j)
                        v[i][j] = false;
                break;
        }
    }
}

@Author: YouSiki