bzoj3659: Which Dreamed It

題面在這裏

題意:

有n個房間,每個房間有若干把鑰匙能夠打開某個房間的門。
最初你在房間1。

每當你到達一個房間,你可以選擇該房間的一把鑰匙,前往該鑰匙對應的房間,並將該鑰匙丟到垃圾桶中。

你希望最終回到房間1,且垃圾桶中有所有的鑰匙。
求方案數。兩組方案不同,當且僅當使用鑰匙的順序不同。
每把鑰匙都是不同的。
房間數小於等於100,鑰匙數小於等於200000.

做法:

我做的都是些啥題啊(霧)感覺自己學了一些奇怪的定理就沒啥了

普及一些奇怪的定理。。
best theorem(最好定理??(霧
有向圖中以i爲起點的歐拉回路個數爲:i×i=1n(degi1)! .
matrix tree theorem(矩陣樹定理)
以i爲根的樹形圖個數=基爾霍夫矩陣去掉第i行第i列的行列式.

基爾霍夫矩陣:度數矩陣減去鄰接矩陣(鄰接矩陣要記錄度數)

於是ans=1×deg1

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

代碼:

/*************************************************************
    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;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章