思路

  • 把每个坐标转化为一个数,如果"1"相邻的话进行合并
  • 因为方案数等于连通块数量的阶乘乘上每一个连通块的大小
  • 把格子从0改成1后,进行除原先的大小乘上合并后的大小

代码

// Problem: 点一成零
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/9981/D
// Memory Limit: 524288 MB
// Time Limit: 4000 ms
// Powered by CP Editor (https://github.com/cpeditor/cpeditor)

#include <bits/stdc++.h>
using namespace std;
#define pb push_back
#define mp(aa,bb) make_pair(aa,bb)
#define _for(i,b) for(int i=(0);i<(b);i++)
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define per(i,b,a) for(int i=(b);i>=(a);i--)
#define mst(abc,bca) memset(abc,bca,sizeof abc)
#define X first
#define Y second
#define lowbit(a) (a&(-a))
#define debug(a) cout<<#a<<":"<<a<<"\n"
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
typedef long double ld;
const int N=505;
const int INF=0x3f3f3f3f;
const int mod=1e9+7;
const double eps=1e-6;
const double PI=acos(-1.0);

int fa[N*N],siz[N*N];
void init(int n){
    for(int i=1;i<=n;i++) 
        fa[i]=i,siz[i]=1;
}
//查询树的根
int find(int x){
    if(x!=fa[x]) fa[x]=find(fa[x]);
    return fa[x];
}
//合并a和b所属的集合   
void unite(int a,int b){
    a=find(a),b=find(b);
    if(a!=b){
        siz[b]+=siz[a];
        fa[a]=fa[b];
    }
}
//判断a和b是否属于同一个集合
bool same(int a,int b){
    return find(a)==find(b);
}

char a[N][N];
int dx[4]={0,0,1,-1},dy[4]={1,-1,0,0};


ll fpow(ll a,ll b){
    ll ans=1%mod;
    while(b){
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}

void solve(){
    int n;cin>>n;
    rep(i,1,n) rep(j,1,n) cin>>a[i][j];
    init(n*n+n);
    rep(i,1,n) rep(j,1,n){
        if(a[i][j]=='1'){
            if(a[i-1][j]=='1') unite(i*n+j,(i-1)*n+j);
            if(a[i][j-1]=='1') unite(i*n+j,i*n+j-1); 
            if(a[i+1][j]=='1') unite(i*n+j,(i+1)*n+j); 
            if(a[i][j+1]=='1') unite(i*n+j,i*n+j+1); 
        }
    }
    ll cnt=0,ans=1;
    rep(i,1,n) rep(j,1,n){
        if(a[i][j]=='1'&&fa[i*n+j]==i*n+j){
            cnt++;
            ans*=siz[i*n+j]%mod,ans%=mod;    
        }
    }
    rep(i,1,cnt){ // cnt!
        ans*=i%mod,ans%=mod;
    }

    int k;cin>>k;
    while(k--){
        int x,y;cin>>x>>y;
        x++,y++;
        if(a[x][y]=='1'){
            cout<<ans<<"\n";
        }
        else{
            a[x][y]='1';
            cnt++;
            ans*=cnt%mod,ans%=mod;
            for(int i=0;i<4;i++){
                int xx=x+dx[i],yy=y+dy[i];
                if(a[xx][yy]=='1'){
                    int f1=find(xx*n+yy),f2=find(x*n+y);
                    if(f1!=f2){
                        ans*=fpow(cnt,mod-2),ans%=mod; //ans/=cnt;
                        ans*=fpow(siz[f1],mod-2),ans%=mod; //ans/=siz[f1];
                        ans*=fpow(siz[f2],mod-2),ans%=mod; //ans/=siz[f2];
                        ans*=(siz[f1]+siz[f2])%mod,ans%=mod;
                        cnt--;
                        unite(f1,f2);
                    }
                }
            }
            cout<<ans<<"\n";
        }
    }
}


int main(){
    ios::sync_with_stdio(0);cin.tie(0);
//    int t;cin>>t;while(t--)
    solve();
    return 0;
}