bzoj3659: Which Dreamed It

题面在这里

题意:

有n个房间,每个房间有若干把钥匙能够打开某个房间的门。
最初你在房间1。

每当你到达一个房间,你可以选择该房间的一把钥匙,前往该钥匙对应的房间,并将该钥匙丢到垃圾桶中。

你希望最终回到房间1,且垃圾桶中有所有的钥匙。
求方案数。两组方案不同,当且仅当使用钥匙的顺序不同。
每把钥匙都是不同的。
房间数小于等于100,钥匙数小于等于200000.

做法:

我做的都是些啥题啊(雾)感觉自己学了一些奇怪的定理就没啥了

普及一些奇怪的定理。。
best theorem(最好定理??(雾
有向图中以i为起点的欧拉回路个数为: i×i=1n(degi1)! 以 i 为 根 的 树 形 图 个 数 × ∏ i = 1 n ( d e g i − 1 ) ! .
matrix tree theorem(矩阵树定理)
以i为根的树形图个数=基尔霍夫矩阵去掉第i行第i列的行列式.

基尔霍夫矩阵:度数矩阵减去邻接矩阵(邻接矩阵要记录度数)

于是 ans=1×deg1 a n s = 以 1 为 起 点 的 欧 拉 回 路 个 数 × d e g 1

至于行列式怎么求,就高斯消元一下好了。

代码:

/*************************************************************
    Problem: bzoj 3659 Which Dreamed It
    User: fengyuan
    Language: C++
    Result: Accepted
    Time: 808 ms
    Memory: 2172 kb
    Submit_Time: 2018-01-24 14:50:07
*************************************************************/

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cctype>
#include<cstdlib>
#include<cmath>
using namespace std;
typedef long long ll;

inline ll read() {
    char ch = getchar(); ll x = 0; int op = 1;
    for(; !isdigit(ch); ch = getchar()) if(ch == '-') op = -1;
    for(; isdigit(ch); ch = getchar()) x = x*10+ch-'0';
    return x*op;
}
inline void write(ll a) {
    if(a < 0) putchar('-'), a = -a;
    if(a >= 10) write(a/10); putchar('0'+a%10);
}

const int N = 110, M = 200010;
const int mod = 1000003;
int n, m;
int a[N][N], e[N][N], fac[M], vis[N], in[N], out[N];

inline int ksm(int x, int p) {
    int ret = 1;
    for(; p; p >>= 1, x = (ll)x*x%mod) if(p&1) ret = (ll)ret*x%mod;
    return ret;
}
inline int inv(int x) { return ksm(x, mod-2); }
inline void prepare() { fac[0] = 1; for(int i = 1; i < M; i ++) fac[i] = (ll)fac[i-1]*i%mod; }
inline void dfs(int u) {
    vis[u] = ++ m;
    for(int i = 1; i <= n; i ++) if(e[u][i] && !vis[i]) dfs(i);
}
inline int det(int n) {
    int ans = 1; bool flag = 0;
    for(int i = 1; i <= n; i ++) for(int j = 1; j <= n; j ++) a[i][j] = (a[i][j]%mod+mod)%mod;
    for(int i = 1; i <= n; i ++) {
        int r = i;
        for(int j = i+1; j <= n; j ++) if(a[j][i]) { r = j; break; }
        if(i != r) { for(int j = i; j <= n; j ++) swap(a[i][j], a[r][j]); flag ^= 1; }
        int ttt = inv(a[i][i]);
        for(int j = i+1; j <= n; j ++) {
            int t = (ll)a[j][i]*ttt%mod;
            for(int k = i; k <= n; k ++) a[j][k] = (a[j][k]+mod-(ll)t*a[i][k]%mod)%mod;
        }
        ans = (ll)ans*a[i][i]%mod;
        if(!ans) return 0;
    }
    if(flag) ans = (mod-ans)%mod;
    return ans;
}
int main() {
    prepare();
    while(1) {
        n = read(); if(!n) break;
        memset(a, 0, sizeof a); memset(e, 0, sizeof e); memset(vis, 0, sizeof vis);
        memset(in, 0, sizeof in); memset(out, 0, sizeof out);
        for(int i = 1; i <= n; i ++) {
            int k = read();
            while(k --) { int x = read(); e[i][x] ++; }
        } m = 0;
        dfs(1); bool flag = 1;
        for(int i = 1; i <= n; i ++)
            if(!vis[i] && e[i]) { flag = 0; break; }
        if(!flag) { puts("0"); continue; }
        for(int i = 1; i <= n; i ++)
            for(int j = 1; j <= n; j ++) if(e[i][j]) {
                int x = vis[i], y = vis[j];//vis[i]是给节点重新编号,剔除不存在的点
                in[y] += e[i][j]; out[x] += e[i][j];
                a[x-1][y-1] = (a[x-1][y-1]-e[i][j]+mod)%mod; a[x-1][x-1] = (a[x-1][x-1]+e[i][j])%mod;
            } flag = 1;
        for(int i = 1; i <= m; i ++) if(in[i] != out[i]) { flag = 0; break; }
        if(!flag) { puts("0"); continue; }
        if(m == 1) { write(fac[e[1][1]]); continue; }
        int ans = (ll)det(m-1)*in[1]%mod;
        for(int i = 1; i <= m; i ++) ans = (ll)ans*fac[in[i]-1]%mod;
        write(ans); puts("");
    } return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值