1. 程式人生 > >生成樹

生成樹

swap get tor 生成樹計數 int sed 最小值 數量 基環樹

[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 }
View Code

[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 }
View Code

生成樹