[BZOJ 3534] 重建

　　題意

　　給定一張 n (1 < n <= 50) 個點 m 條邊的無向圖, 每條邊有一定出現的概率, 問出現的邊恰好形成一個生成樹的概率.

　　實現

　　註意到對於出現概率為 1 的邊, 分母可能為 0 . 我們對 d = 1 變換為 d = 1-EPS 就好了.

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <cstdlib>
 4 #include <cctype>
 5 #include <cmath>
 6 #include <algorithm>
 7
 
 using namespace std;
 8 #define F(i, a, b) for (register int i = (a); i <= (b); i++)
 9 #define db double
10 
11 const db EPS = 1e-9;
12 inline int sign(db x) { return fabs(x)<EPS ? 0 : x<0 ? -1 : 1; }
13 inline int cmp(db x, db y) { return sign(x-y); }
14 
15 const int N = 55;
16 
17 int n; db f[N][N], Div = 1
 
;
18 inline void Init(int x, int y, db w) {
19     f[x][x] += w, f[y][y] += w, f[x][y] -= w, f[y][x] -= w;
20 }
21 db Det(void) {
22     n--;
23     db Ret = 1;
24     F(i, 1, n) {
25         int id = i;
26         while (sign(f[id][i]) == 0) id++;
27         if (id > n) return 0;
28         if (id != i) {
 
29             F(k, 1, n) swap(f[i][k], f[id][k]);
30             Ret = -Ret;
31         }
32         Ret *= f[i][i];
33         F(j, i+1, n) if (sign(f[j][i]) != 0) {
34             db t = f[j][i] / f[i][i];
35             F(k, i, n) f[j][k] -= t * f[i][k];
36         }
37     }
38     return Ret;
39 }
40 
41 int main(void) {
42     #ifndef ONLINE_JUDGE
43         freopen("bzoj3534.in", "r", stdin);
44     #endif
45     
46     scanf("%d", &n);
47     F(i, 1, n) F(j, 1, n) {
48         db p; scanf("%lf", &p);
49         if (i < j) {
50             if (cmp(p, 1) == 0) p -= EPS;
51             Div *= 1-p;
52             Init(i, j, p/(1-p));
53         }
54     }
55     printf("%0.8lf\n", Div * Det());
56     
57     return 0;
58 }

[BZOJ 1016] 最小生成樹計數

　　題意

　　給定一張 n (n <= 16) 個點的圖, 求有多少棵最小生成樹.

　　分析

　　Kruskal 算法告訴我們, 每種邊使用的數量相同, 它們構成的連通性相同.

　　在 Kruskal 的同時, 該權值作用之前的連通塊當做若幹個點, 它們在該權值作用之後形成若幹個連通塊, 對每個連通塊進行生成樹計數.

　　實現

　　使用 vector 進行分類重標號.

Mashmokh‘s Designed Problem　　01生成樹

　　題意

　　給一張 n (1 <= n <= 100000) 個點 m (1 <= m <= 100000) 條邊的無向連通圖, 每條邊的權值可能為 0 或 1 , 問是否存在一種生成樹權值為 x ?

　　分析

　　考慮構建最小生成樹, 最大生成樹, 最小生成樹的權值為 Min, 最大生成樹的權值為 Max , 顯然能構造出來的權值只可能在 [Min, Max] 之間.

　　最小生成樹一定可以通過若幹次消圈得到最大生成樹, 而邊權又是 0 或 1, 所以變化是連續的, 所以 [Min, Max] 之間的所有權值都可以構造出來.

　　綜上, 存在一種生成樹權值為 x 當且僅當 Min <= x 且 x <= Max .

[Hdu 5304]　　基環樹計數

　　題意　　

　　給定一張 n 個點 m 條邊的無向圖, 問其有多少個基環生成樹.

　　n <= 16 , 無重邊, 無自環.

　　分析

　　利用狀壓DP 求每種環的狀態的個數.

　　枚舉環, 縮點, 利用 Matrix-Tree定理 進行計數.

　　實現

　　f[s][i] 表示從 s 狀態中的最小值出發, 終點在 i 的方案數.

　　邊界 f[2 ^ i][i] = 1 .

　　統計答案的時候, 對於 f[s][i] , 若 s 中的個數大於 2 , 且 i 能到達 s 狀態中的最小值, 則 cnt[s] += f[s][i] .

  1 #include <cstdio>
  2 #include <cstring>
  3 #include <cstdlib>
  4 #include <cctype>
  5 #include <algorithm>
  6 using namespace std;
  7 #define F(i, a, b) for (register int i = (a); i <= (b); i++)
  8 #define D(i, a, b) for (register int i = (a); i >= (b); i--)
  9 inline int rd(void) {
 10     int f = 1; char c = getchar(); for (; !isdigit(c); c = getchar()) if (c == ‘-‘) f = -1;
 11     int x = 0; for (; isdigit(c); c = getchar()) x = x*10+c-‘0‘; return x*f;
 12 }
 13 
 14 const int N = 20;
 15 const int S = 70000;
 16 const int MOD = 998244353;
 17 
 18 inline void Add(int &x, int y) { x = (x + y) % MOD; }
 19 inline int Pow(int x, int y) {
 20     int Mul = 1;
 21     for (; y > 0; y >>= 1, x = 1LL * x * x % MOD)
 22         if (y & 1) Mul = 1LL * Mul * x % MOD;
 23     return Mul;
 24 }
 25 inline int Inv(int x) { return Pow(x, MOD-2); }
 26 
 27 int n, m; bool G[N][N];
 28 int Full, cnt[S], f[S][N];
 29 int Lab[N], tot, Mat[N][N], ans;
 30 
 31 void Build(int s) {
 32     tot = 1; F(i, 1, n) Lab[i] = (s >> i-1 & 1 ? 1 : ++tot);
 33     memset(Mat, 0, sizeof Mat);
 34     F(i, 1, n) F(j, i+1, n)
 35         if (G[i][j] && Lab[i] != Lab[j]) {
 36             Mat[Lab[i]][Lab[i]]++, Mat[Lab[j]][Lab[j]]++;
 37             Mat[Lab[i]][Lab[j]]--, Mat[Lab[j]][Lab[i]]--;
 38         }
 39 }
 40 int Gauss(void) {
 41     tot--;
 42     int Mul = 1;
 43     F(i, 1, tot) {
 44         if (!Mat[i][i]) {
 45             F(j, i+1, tot) if (Mat[j][i] != 0) {
 46                 F(k, 1, tot) swap(Mat[i][k], Mat[j][k]);
 47                 Mul = -Mul;
 48                 break;
 49             }
 50         }
 51         if (!Mat[i][i]) return 0;
 52         Mul = 1LL * Mul * Mat[i][i] % MOD;
 53         int I = Inv(Mat[i][i]);
 54         F(j, i+1, tot) if (Mat[j][i] != 0) {
 55             int t = 1LL * I * Mat[j][i] % MOD;
 56             F(k, i, tot)
 57                 Add(Mat[j][k], -1LL * Mat[i][k] * t % MOD);
 58         }
 59     }
 60     return Mul;
 61 }
 62 
 63 int main(void) {
 64     #ifndef ONLINE_JUDGE
 65         freopen("hdu5304.in", "r", stdin);
 66     #endif
 67     
 68     int inv2 = Inv(2);
 69     while (~scanf("%d %d", &n, &m)) {
 70         memset(G, false, sizeof G);
 71         F(i, 1, m) {
 72             int x = rd(), y = rd();
 73             G[x][y] = G[y][x] = true;
 74         }
 75         
 76         Full = (1 << n) - 1, memset(cnt, 0, sizeof cnt), memset(f, 0, sizeof f);
 77         F(i, 1, n) f[1 << i-1][i] = 1;
 78         F(s, 1, Full) {
 79             int Min = 0; D(i, n, 1) if (s >> i-1 & 1) Min = i;
 80             F(i, Min, n) if ((s >> i-1 & 1) && f[s][i] != 0) {
 81                 F(j, Min+1, n) if (!(s >> j-1 & 1) && G[i][j])
 82                     Add(f[s | 1 << j-1][j], f[s][i]);
 83             }
 84         }
 85         F(s, 1, Full) {
 86             int tot = 0; for (int x = s; x > 0 && tot <= 2; x ^= (x & -x), tot++); if (tot <= 2) continue;
 87             int Min = 0; D(i, n, 1) if (s >> i-1 & 1) Min = i;
 88             F(i, Min+1, n) if ((s >> i-1 & 1) && G[i][Min])
 89                 Add(cnt[s], f[s][i]);
 90             cnt[s] = 1LL * cnt[s] * inv2 % MOD;
 91         }
 92         
 93         ans = 0;
 94         F(s, 1, Full) if (cnt[s] != 0) {
 95             Build(s);
 96             Add(ans, 1LL * cnt[s] * Gauss() % MOD);
 97         }
 98         printf("%d\n", (ans + MOD) % MOD);
 99     }
100     
101     return 0;
102 }

生成樹