BZOJ4031: [HEOI2015]小Z的房间

BZOJ4031

任意两点若都不是柱子，视为可连边。
然后矩阵树定理就好。

【代码】

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#define N 200005
#define Mod 1000000000
#define INF 200001
using namespace std;
typedef long long ll;

ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m,Num;ll ans; 
char ch[10][10];
ll a[100][100];
int id[10][10];
int dx[]={0,1};
int dy[]={1,0};

void debug()
{
    for(int i=1;i<Num;i++,puts(""))
    for(int j=1;j<Num;j++) printf("%lld ",a[i][j]);
    printf("\n");
}

void Add(int x,int y)
{
    a[x][x]++,a[y][y]++;
    a[x][y]--,a[y][x]--;
}

ll Gauss_Elimination()
{
    ll rtn=1;
    for(int i=1;i<Num;i++)
    {
        for(int j=i+1;j<Num;j++)
        {
            while(a[j][i]) 
            {
                ll t=a[i][i]/a[j][i];
                for(int k=i;k<Num;k++)
                    a[i][k]=(a[i][k]-t*a[j][k])%Mod;
                rtn=-rtn;
                for(int k=i;k<Num;k++) swap(a[i][k],a[j][k]);
            }
        }
        rtn=rtn*a[i][i]%Mod;
        if(rtn==0) return 0;
    }
    return (rtn+Mod)%Mod;   
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)
    {
        scanf("%s",ch[i]+1);
        for(int j=1;j<=m;j++) if(ch[i][j]=='.') id[i][j]=++Num;
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) if(ch[i][j]=='.')
        {
            for(int k=0;k<2;k++)
            {
                int xx=i+dx[k],yy=j+dy[k];
                if(!xx||!yy||xx>n||yy>m||ch[xx][yy]=='*') continue;
                Add(id[i][j],id[xx][yy]);
            }
        }
    printf("%lld\n",Gauss_Elimination());
    return 0;
}