[bzoj5133][矩阵树]白金元首与独舞

传送门
http://www.lydsy.com/JudgeOnline/problem.php?id=5133

sol:

考虑建一个汇点，那么这道题求的就是有向图的树的个数。但是n^4的矩阵太大了。注意到原本题目中已经给出的点对角线上消出来的肯定是1（可以使用显然证明法，给出的点显然不会有贡献）。所以已经确定的点我们就可以不考虑了。对于一个不确定的点，他4个方向走出去，不是走到另一个不确定的点就是走到汇点。这样用关键点建出来的矩阵就可以接受了。

非常的难受啊，晚上脑子好像不怎么转的。高消写错了好多地方。。wor
全网最长代码系列，相信你们也不是来看代码的。

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<cmath>
using namespace std;
typedef long long ll;
typedef double s64;
int n,m;
inline int read()
{
    char c;
    int res,flag=0;
    while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
    res=c-'0';
    while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
    return flag?-res:res;
}
const int N=310;
char sr[N][N];
int tot,b[N][N],c[N][N];
int to[N][N];
const int dx[4]={0,0,1,-1};
const int dy[4]={1,-1,0,0};
const int pyz=1e9+7;
int ans;
inline int ksm(int s,int t)
{
    int res=1;
    while(t)
    {
        if(t&1) res=(ll)res*s%pyz;
        s=(ll)s*s%pyz;
        t>>=1;
    }
    return res;
}
inline void guass()
{
/*  for(int i=1;i<=tot;++i){
    for(int j=1;j<=tot;++j)
    printf("%d ",c[i][j]);
printf("\n");}*/

    ans=1;
    for(int i=1;i<=tot;++i)
    {
        int row;
        for(row=i;row<=tot;++row) if(c[row][i]) break;
        if(row>tot) continue;
        if(row!=i)
        {
            for(int j=1;j<=tot;++j) swap(c[row][j],c[i][j]);
            ans*=-1;
        }
        int inv=ksm(c[i][i],pyz-2);
        for(int j=i+1;j<=tot;++j)
        if(c[j][i])
        {
            int tmp=(ll)c[j][i]*inv%pyz;
            for(int k=1;k<=tot;++k)
            c[j][k]=(c[j][k]+pyz-(ll)c[i][k]*tmp%pyz)%pyz;
        }
    }
    ans=(ans+pyz)%pyz;
    for(int i=1;i<=tot;++i)
    ans=(ll)ans*c[i][i]%pyz;
}
inline int d(int i,int j)
{
    if(sr[i][j]=='R') return 0;
    if(sr[i][j]=='L') return 1;
    if(sr[i][j]=='D') return 2;
    if(sr[i][j]=='U') return 3;
}
inline int get(int i,int j,int k)
{
    int x,y;
    if(to[i][j]) return to[i][j];
    x=i+dx[k];
    y=j+dy[k];
    if(1<=x&&x<=n&&1<=y&&y<=m)
    {
        if(sr[x][y]=='.') return b[x][y];
        return get(x,y,d(x,y));
    }
    else return tot+1;
}
int father[N*N];
inline int find(int x)
{
    return father[x]==x?x:father[x]=find(father[x]);
}
inline void unit(int x,int y)
{
    int fx=find(x),fy=find(y);
    if(fx!=fy) father[fy]=fx;
}
inline void solve()
{
    n=read();
    m=read();
    tot=0;
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    memset(to,0,sizeof(to));
    for(int i=1;i<=n;++i)
    scanf("%s",sr[i]+1);
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    {
        b[i][j]=++tot;
        father[tot]=tot;
    }
    father[tot+1]=tot+1;
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    if(sr[i][j]!='.')
    {
        int k=d(i,j);
        int x,y;
        x=i+dx[k];
        y=j+dy[k];
        if(1<=x&&x<=n&&1<=y&&y<=m&&sr[x][y]!='.')
        unit(b[i][j],b[x][y]);
        else unit(b[i][j],tot+1);
    }
    int flag=0;
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    if(sr[i][j]!='.')
    if(find(b[i][j])==b[i][j]) ++flag;
    if(find(tot+1)==tot+1) ++flag;
    if(flag!=1)
    {
        printf("0\n");
        return;
    }

    tot=0;
    memset(b,0,sizeof(b));
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    if(sr[i][j]=='.') b[i][j]=++tot;
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    {
        if(sr[i][j]=='.')
        for(int k=0;k<=3;++k)
        {
            int t=get(i,j,k);
            c[t][b[i][j]]--;
            c[b[i][j]][b[i][j]]++;
        }
        else to[i][j]=get(i,j,d(i,j));
    }
    for(int i=1;i<=tot;++i)
    for(int j=1;j<=tot;++j)
    c[i][j]=(c[i][j]%pyz+pyz)%pyz;
    guass();
    printf("%d\n",ans);
}
int main()
{
//  freopen("5133.in","r",stdin);
//  freopen(".out","w",stdout);
    int T=read();
    while(T--) solve();
}