【51nod 1634】【HDU5729】刚体图/Rigid Frameworks（二分图、DP）

最新推荐文章于 2022-04-24 23:34:47 发布

Hany01

最新推荐文章于 2022-04-24 23:34:47 发布

阅读量346

点赞数 1

分类专栏： 51nod DP 二分图 HDU

本文链接：https://blog.csdn.net/hhaannyyii/article/details/79509211

版权

DP 同时被 3 个专栏收录

83 篇文章 0 订阅

订阅专栏

HDU

13 篇文章 0 订阅

订阅专栏

51nod

7 篇文章 0 订阅

订阅专栏

Description

link1
link2

Solution

2.1

发现加一条对角线可以使该行、列的边互相垂直，将每一行、列看成二分图的点，从而可以把题目转化成求二分图的连通图方案。

2.2

设 $dp_{i,j}$ 为左边有 $i$ 个点、右边有 $j$ 个点的连通图的方案，转移方程：

d p n, m = \sum i = 1 n \sum j = 0 m (n - 1 i - 1) (m j) d p i, j 3 (n - i) (m - j)

$dp_{n,m}=\sum_{i=1}^{n}\sum_{j=0}^{m}\binom{n-1}{i-1}\binom{m}{j}dp_{i,j}3^{(n-i)(m-j)}$
解释一下，我们这里是通过左边节点

1 1 $1$ 所在的联通块大小判断图是否联通的。
其中

i

$i$ 枚举左边和

1 1 $1$ 相联通的点数（包括

1

$1$ 本身），同理

j j $j$ 枚举右边，剩下的边随便乱放，即

3

$3$ 的幂次。

Source

/************************************************
 * Au: Hany01
 * Date: Mar 10th, 2018
 * Prob: [51nod 1634] 刚体图
 * Email: hany01@foxmail.com
************************************************/

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define fir first
#define sec second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000007)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia

template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

inline int read()
{
    register int _, __; register char c_;
    for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
    for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

inline void File()
{
#ifdef hany01
    freopen("51nod1634.in", "r", stdin);
    freopen("51nod1634.out", "w", stdout);
#endif
}

int C[11][11], dp[11][11], Pow3[101];

inline void Init()
{
    Pow3[0] = 1;
    For(i, 1, 100) Pow3[i] = Pow3[i - 1] * 3ll % Mod;

    For(i, 0, 10) {
        C[i][0] = C[i][i] = 1;
        For(j, 1, i - 1) C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % Mod;
    }
}

inline void Solve()
{
    dp[0][1] = dp[1][0] = 1;
    For(n, 1, 10) For(m, 1, 10) {
        dp[n][m] = Pow3[n * m];
        For(i, 1, n) For(j, 0, m) if (i != n || j != m)
            (dp[n][m] -= (LL)dp[i][j] * (LL)C[n - 1][i - 1] % Mod * (LL)C[m][j] % Mod * (LL)Pow3[(n - i) * (m - j)] % Mod) %= Mod;
        if (dp[n][m] < 0) dp[n][m] += Mod;
    }
}

int main()
{

    File();

    Init();

    Solve();

    static int n, m;
    while (~scanf("%d%d", &n, &m)) printf("%d\n", dp[n][m]);

    return 0;
}
//笙歌散后酒初醒，深院月斜人静。
//    -- 司马光《西江月·宝髻松松挽就》