首页 > 其他 > 详细

[BZOJ4894]天赋

时间:2018-04-21 10:58:49      阅读:222      评论:0      收藏:0      [点我收藏+]

bzoj

sol

矩阵树定理裸题
其实可以不用辗转相除直接乘逆元就好了。然而我还是写的辗转相除

code

#include<cstdio>
#include<algorithm>
using namespace std;
int gi()
{
    int x=0,w=1;char ch=getchar();
    while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
    if (ch=='-') w=0,ch=getchar();
    while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return w?x:-x;
}
const int N = 305;
const int mod = 1e9+7;
int n,a[N][N],ans=1;char s[N][N];
int main()
{
    n=gi();
    for (int i=1;i<=n;++i)
    {
        scanf("%s",s[i]+1);
        for (int j=1;j<=n;++j)
            if (s[i][j]=='1') a[i][j]--,a[j][j]++;
    }
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j)
            a[i][j]=(a[i][j]+mod)%mod;
    for (int i=2;i<=n;++i)
    {
        for (int j=i+1;j<=n;++j)
            while (a[j][i])
            {
                int t=a[i][i]/a[j][i];
                for (int k=i;k<=n;++k) a[i][k]=(a[i][k]-1ll*t*a[j][k]%mod+mod)%mod,swap(a[i][k],a[j][k]);
                ans=(mod-ans)%mod;
            }
        ans=1ll*ans*a[i][i]%mod;
    }
    printf("%d\n",ans);return 0;
}

[BZOJ4894]天赋

原文:https://www.cnblogs.com/zhoushuyu/p/8899076.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!