【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;
}
posted @ 2018-06-25 21:41  yoyoball  阅读(297)  评论(0编辑  收藏  举报