[JZOJ3809] 设备塔

Description

为了封印辉之环，古代塞姆利亚大陆的人民在异空间中建造了一座设备塔。
简单的说，这座设备塔是一个漂浮在异空间中的圆柱体，圆柱体两头的圆是计算核心，而侧面则是
传输信息所用的数据通道，划分成N *m 个区块。
然而，随着工作的继续进行，他们希望把侧面的一部分区块也改造成其他模块。然而，任何时候都
必须保证存在一条数据通道，能从圆柱体的一端通向另一端。
由于无法使用辉之环掌控下的计算系统，他们寻求你的帮助来解决这个问题。他们将逐个输入想要
改造的区域，而你则执行所有可行的改造并忽略可能导致数据中断的改造（按顺序）。
输出可行改造数。
• 对于分值为30 的子任务1，保证N;M <=100;K<= 5000
• 对于分值为30 的子任务2，保证N;M <=3000;K <= 5000
• 对于分值为40 的子任务3，保证N;M <=3000;K <= 300000。

Solution

容易想到，只有构成了一种情况路才会被封住，就是一条连着的线（八方向）连到左右边界。并且因为这是环，所以还要左右相接。

明显可以复制一遍。
然后对于所有操作，目前是合法的，加上这个点以及它复制后的对应点如果要隔断是怎么样的？

如果我们将相连的块用一个并查集维护，那么因为复制前后是完全相同的，所以如果这两个对应点加上后在同一个集合内，显然是不合法的。

为什么？

显然有一条路径从原点穿过复制的边界到对应点，那么将它分成两部分，复制后的部分在复制前的部分一定存在，因此构成了一个环将路封住。

所以我们只需判断原点和对应点八方向有没有在同一集合就行了。

Code

#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define N 3005
#define M N*N*2
#define INF 1802201963
using namespace std;
int f[M],n,m,k,fx[8][2]={{1,-1},{0,-1},{-1,-1},{1,1},{0,1},{-1,1},{1,0},{-1,0}},l1[8],l2[8];
bool bz[N][2*N];
int g(int x,int y)
{
    return (x-1)*2*m+y;
}
int getf(int k)
{
    if(f[k]==0||f[k]==k) return k;
    f[k]=getf(f[k]);
    return f[k];
}
void hb(int x,int y)
{
    if(x==-1||y==-1) return;
    int fx=getf(x),fy=getf(y);
    if(fx==fy) return;
    f[fy]=fx;
}
int main()
{
    cin>>n>>m>>k;
    int s=0;
    memset(f,0,sizeof(f));
    fo(i,1,k)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        bool b1=0;
        if(bz[x][y]==1) continue;
        fo(j,0,7)
        {
            int x2=x+fx[j][0],y2=y+fx[j][1];
            if(y2==0) y2=2*m;
            if(x2<=n&&x2>0&&y2<=2*m&&bz[x2][y2]) l1[j]=getf(g(x2,y2));
            else l1[j]=-1;
        }
        fo(j,0,7)
        {
            int x3=x+fx[j][0],y3=y+fx[j][1]+m;
            if(y3==2*m+1) y3=1;
            if(x3<=n&&x3>0&&y3<=2*m&&y3>0&&bz[x3][y3]) l2[j]=getf(g(x3,y3));
            else l2[j]=-1;
        }
        fo(j,0,7)
            fo(p,0,7)
                if(l1[j]>0&&l2[p]>0&&l1[j]==l2[p])
                {
                    b1=1;
                    break;
                }
        if(!b1) 
        {
            bz[x][y]=bz[x][y+m]=1;
            fo(j,0,7) hb(l1[j],g(x,y)),hb(l2[j],g(x,y+m));
            s++;
        }
    }
    cout<<s;
}