传送门
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();
}