【bzoj4894】天赋

Portal-->bzoj4894

Solution

　　这题的话其实，一句话题意就是求。。外向树（方向是根往叶子）。。

　　然后关于有向图的生成树计数的话，求外向树就是将度数矩阵改成入度，内向树就是改成出度

​　　然后其他的一样的

　　

　　代码大概长这个样子

#include<iostream>

#include<cstdio>

#include<cstring>

using namespace std;

const int N=310,MOD=1e9+7;

int a[N][N],mp[N][N],d[N][N];

int n,m;

int solve(int n);

int main(){

#ifndef ONLINE_JUDGE

    freopen("a.in","r",stdin);

#endif

    char ch;

    scanf("%d\n",&n);

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

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

            scanf("%c",&ch);

            mp[i][j]=ch-'0';

            if (mp[i][j]) ++d[j][j];

        }

        scanf("\n");

    }

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

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

            a[i][j]=(d[i][j]-mp[i][j]+MOD)%MOD;

    printf("%d\n",solve(n));

}

int solve(int n){

    int id,mark=1;

    int tmp;

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

        id=i;

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

            if (a[j][i]){id=j;break;}

        if (id!=i){

            mark=-mark;

            for (int j=2;j<=n;++j) swap(a[id][j],a[i][j]);

        }

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

            while (a[j][i]){

                tmp=a[j][i]/a[i][i];

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

                    a[j][k]=(1LL*a[j][k]+MOD-1LL*tmp*a[i][k]%MOD)%MOD;

                if (a[j][i]==0) break;

                mark=-mark;

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

                    swap(a[j][k],a[i][k]);

            }

        }

    }

    int ret=mark;

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

        ret=1LL*ret*a[i][i]%MOD;

    return (ret+MOD)%MOD;

}