BZOJ 4031: [HEOI2015]小Z的房间 [矩阵树定理 行列式取模]

http://www.lydsy.com/JudgeOnline/problem.php?id=4031

---

裸题........

问题在于模数是$10^9$

我们发现消元的目的是让一个地方为0

辗转相除法也可以做到这一点

只不过取模用减整除来代替就好了

注意本题需要分配$id$，因为柱子不能算

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=,MOD=1e9;
inline int read(){
    char c=getchar();int x=;
    while(c<''||c>''){c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x;
}
int n,m,a[N][N];
char s[][];
void Gauss(int n){
    int s=;
    for(int i=;i<=n;i++){//printf("i %d\n",i);
        int r=i;
        for(;r<=n;r++) if(a[r][i]) break;
        if(r==n+){puts("");return;}
        if(r!=i){
            s^=;
            for(int k=i;k<=n;k++) swap(a[i][k],a[r][k]);
        }
        for(int j=i+;j<=n;j++)
            while(a[j][i]){//printf("j %d %d\n",j,a[j][i]);
                ll t=a[j][i]/a[i][i];
                for(int k=i;k<=n;k++) a[j][k]=(a[j][k]-t*a[i][k]%MOD+MOD)%MOD;
                if(a[j][i]==) break;
                s^=;
                for(int k=i;k<=n;k++) swap(a[i][k],a[j][k]);
            }
    }
    ll ans=;
    for(int i=;i<=n;i++) ans=ans*a[i][i]%MOD;
    if(s) ans=(-ans+MOD)%MOD;
    printf("%lld",ans);
}
int id[N][N],tot;
void buildEquation(){
    for(int i=;i<=m;i++) for(int j=;j<=n;j++) if(s[i][j]=='.') id[i][j]=++tot;
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++) if(id[i][j]){
            int u=id[i][j],v;
            if(i!=&&s[i-][j]=='.'){
                v=id[i-][j];
                a[u][u]++;a[v][v]++;
                a[u][v]--;a[v][u]--;
            }
            if(j!=&&s[i][j-]=='.'){
                v=id[i][j-];
                a[u][u]++;a[v][v]++;
                a[u][v]--;a[v][u]--;
            }
        }
    for(int i=;i<=m*n;i++) for(int j=;j<=m*n;j++) a[i][j]=(a[i][j]+MOD)%MOD;//printf("%d%c",a[i][j],j==n?'\n':' ');
}
int main(){
    freopen("in","r",stdin);
    m=read();n=read();
    for(int i=;i<=m;i++) scanf("%s",s[i]+);
    buildEquation();//puts("hi");
    Gauss(tot-);
}