题解 排序 [矩阵树]

题解 排序 [矩阵树]

数据范围: $n \leq 100000 , m \leq 50 $

解题思路

把每行每列看成一个点，如果图中$(i,j)$位置为.,则将第$i$行和第$j$行连一条边，那么题目即要求生成树的个数

虽然矩阵树定理是$O(n^3)$的，但这道题的图有些特殊可以优化复杂度。

考虑基尔霍夫矩阵的样子

可以发现在$n\times n$之间只有对角线上有值，其他位置在$n$和$m$里的边的对应位置有值，那么高斯消元可以只消下面的$m$部分，时间复杂度为$O(n{m}^2)$。

代码

/*******************************
Author:galaxy yr
LANG:C++
Created Time:2019年09月27日 星期五 10时11分07秒
*******************************/
#include<iostream>
#include<cstdio>

using namespace std;

const int maxn=1e4+100;
const int mod=998244353;
int n,m,md,ans,a[maxn][maxn];
char ch[maxn][55];
int h[maxn],l[55];

int ksm(int a,int b)
{
    if(a>=mod)a%=mod;
    int res=1;
    while(b)
    {
        if(b&1)res=1ll*res*a%mod;
        a=1ll*a*a%mod;
        b>>=1;
    }
    return res;
}

void solve()
{
    for(int i=1;i<=n+m;i++)
    {
        for(int j=max(n+1,i+1);j<=n+m;j++)
        {
            int tmp=1ll*a[j][i]*ksm(a[i][i],mod-2)%mod;
            a[j][i]=0;
            for(int k=max(i,n+1);k<=n+m;k++)
                a[j][k]=(a[j][k]-1ll*a[i][k]*tmp%mod+mod)%mod;
        }
    }
    ans=1;
    for(int i=1;i<n+m;i++)
    {
        ans=1ll*ans*a[i][i]%mod;
    }
    printf("%d %d\n",n+m-1,ans);
}

int main()
{
    //freopen("copper.in","r",stdin);
    //freopen("copper.out","w",stdout);
    ios::sync_with_stdio(false);
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            cin>>ch[i][j];
            if(ch[i][j]=='.')
                a[i][i]++,a[j+n][j+n]++,a[i][j+n]=a[j+n][i]=(mod-1);
        }
    solve();
    return 0;
}