BZOJ4894:天赋(矩阵树定理)

Description

小明有许多潜在的天赋，他希望学习这些天赋来变得更强。正如许多游戏中一样，小明也有n种潜在的天赋，但有一些天赋必须是要有前置天赋才能够学习得到的。

也就是说，有一些天赋必须是要在学习了另一个天赋的条件下才能学习的。比如，要想学会"开炮"，必须先学会"开枪"。

一项天赋可能有多个前置天赋，但只需习得其中一个就可以学习这一项天赋。

上帝不想为难小明，于是小明天生就已经习得了1号天赋-----"打架"。于是小明想知道学习完这n种天赋的方案数，答案对1,000,000,007取模。

Input

第一行一个整数n。

接下来是一个n*n的01矩阵，第i行第j列为1表示习得天赋j的一个前置天赋为i。

数据保证第一列和主对角线全为0。

n<=300

Output

第一行一个整数，问题所求的方案数。

Sample Input

8
01111111
00101001
01010111
01001111
01110101
01110011
01111100
01110110

Sample Output

72373

Solution

终于来填矩阵树的坑了……

其实也没啥难的，就是高斯消元解个行列式就完了……（行列式的性质可以看这里）

这个题其实是让你求以$1$为根的外向树生成树个数。

无向图生成树：度数矩阵-邻接矩阵

有向图外向生成树：入度矩阵-邻接矩阵

有向图内向生成树：出度矩阵-邻接矩阵

把这个矩阵删掉一行一列然后求出的行列式的值就是生成树个数。

注意有向图的生成树需要删除根所在的行列来计算行列式。

Code

 #include<iostream>

 #include<cstdio>

 #include<cstdlib>

 #define N (309)

 #define LL long long

 #define MOD (1000000007)

 using namespace std;

 LL n,f[N][N];

 char s[N];

 LL inv(LL a)

 {

     LL ans=,b=MOD-;

     while (b)

     {

         if (b&) ans=ans*a%MOD;

         a=a*a%MOD; b>>=;

     }

     return ans;

 }

 void Gauss(LL n)

 {

     int w=,ans=;

     for (int i=; i<=n; ++i)

     {

         int num=i;

         for (int j=i; j<=n; ++j)

             if (abs(f[j][i])>abs(f[num][i])) num=j;

         if (num!=i) swap(f[num],f[i]), w=-w;

         for (int j=i+; j<=n; ++j)

         {

             int t=f[j][i]*inv(f[i][i])%MOD;

             for (int k=i; k<=n; ++k)

                 f[j][k]=(f[j][k]-t*f[i][k])%MOD;

         }

     }

     for (int i=; i<=n; ++i)

         ans=ans*f[i][i]%MOD;

     ans=(ans*w%MOD+MOD)%MOD;

     printf("%lld\n",ans);

 }

 int main()

 {

     scanf("%lld",&n);

     for (int i=; i<n; ++i)

     {

         scanf("%s",s);

         for (int j=; j<n; ++j)

             f[i][j]-=s[j]-'';

     }

     for (int i=; i<n; ++i)

         for (int j=; j<n; ++j)

             if (i!=j && f[i][j]==-) f[j][j]++;

     Gauss(n-);

 }