HDU 5713 & 2016"百度之星" 复赛(Astar Round3)1002 k个联通块

原创 2016年05月30日 20:20:19

题意:众所周知,度度熊喜欢图,尤其是联通的图。
今天,它在图上又玩出了新花样,新高度。有一张无重边的无向图, 求有多少个边集,使得删掉边集里的边后,图里恰好有K个连通块。

思路:首先可以很容易想到状态压缩,但比赛的时候复杂度算错了导致写挂了…….由二项展开可以知道枚举所有子集的所有子集的复杂度应该是O(3n)而不是O(4n)

这道题用状态dp[s][k]表示点集为s,k个联通块的边集数量,那么可以得到一个状态转移方程
dp[s][k]=sigma(dp[s1][k1]+f[s^s1])
其中s1s的一个子集,f[s]表示s这个点集为一个联通块的边集数量,可能重复枚举,所以我们每次只枚举lowbit(s)这个点所在的连通块,那么我们现在的任务就是计算f[s]了。
直接计算f[s]并不好做,我们可以先计算使这个点集不连通的方案数然后用总方案数减去不连通的方案就是f[s]
假设使这个点集不连通的方案数为h[s],那么可以用类似计算dp数组的方法计算,状态转移方程为
h[s]=sigma(f[s1]g[s^s1])
其中g数组是这个边集中总的方案数,也就是说我们每次枚举lowbit(s)所在的连通块,因为总的是不连通的,所以剩下的部分可以任意取边,所以我们就可以得到如上的状态转移方程。
所以我们只需要再预处理每个点集的边的总数然后按照上述算法来计算即可。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
#include<stack>
#include<string>
#include<map>
#include<set>
#include<ctime>
#define eps 1e-6
#define LL long long
#define pii pair<int, int>
//#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
const int mod = 1e9+9;
const int maxn = 15;

int n, m, q, num[1<<maxn], G[maxn][maxn];
LL f[1<<maxn], h[1<<maxn], p[maxn*maxn], dp[1<<maxn][maxn];
void init() {
    memset(G, 0, sizeof(G));
    memset(num, 0, sizeof(num));
    memset(f, 0, sizeof(f));
    memset(h, 0, sizeof(h));
}
inline int lowbit(int x) {
    return x&-x;
}
int cal_edge(int e, int s) {
    int ans = 0;
    for (int i = 0; i < n; i++) 
        if ((1<<i) == e) {
            e = i;
            break;
        }
    for (int i = 0; i < n; i++) 
        if ((1<<i)&s)
            ans += G[i][e];
    return ans;
}

int main()
{
    freopen("input.txt", "r", stdin);
    int T;
    cin >> T;
    int kase = 0;
    p[0] = 1;
    for (int i = 1; i <= 200; i++) {
        p[i] = p[i-1] << 1;
        p[i] %= mod;
    }
    while (T--) {
        scanf("%d%d%d", &n, &m, &q);
        init();
        int maxs = 1 << n;
        for (int i = 1; i <= m; i++) {
            int u, v;
            scanf("%d%d", &u, &v);
            u--;
            v--;
            G[u][v] = G[v][u] = 1;
        }
        num[0] = 0;
        for (int i = 1; i < maxs; i++) {
            int t = lowbit(i);
            num[i] = num[i^t] + cal_edge(t, i);
        } 
        for (int i = 0; i < n; i++) 
            f[1<<i] = p[num[1<<i]];
        for (int i = 1; i < maxs; i++) {
            int t = lowbit(i);
            for (int j = i^t; j; j=(j-1)&(i^t)) {
                h[i] += f[i^j] * p[num[j]];
                h[i] %= mod;
            }
            f[i] = p[num[i]] - h[i];
            f[i] %= mod; 
        }
        for (int i = 1; i < maxs; i++)
            dp[i][1] = f[i];
        for (int i = 1; i < maxs; i++) {
            for (int k = 2; k <= q; k++) {
                dp[i][k] = 0;
                int t = lowbit(i);
                for (int j = i^t; j; j=(j-1)&(i^t)) {
                    dp[i][k] += dp[j][k-1] * f[i^j];
                    dp[i][k] %= mod;
                }
            }
        }
        if (dp[maxs-1][q] < 0)
            dp[maxs-1][q] += mod;
        printf("Case #%d:\n%I64d\n", ++kase, dp[maxs-1][q]);
    }
    return 0;
}















版权声明:本文为Godspeed原创文章,未经博主允许不得转载。

到底是Tier还是Layer?

当Tier遇到Layer看到论坛上有人问到Tier和Layer的区别,这里详细论述一下,大家共同探讨。1。横向与纵向有一种说法是“Layer是水平方向的,tier是垂直方向的”这种说法首先容易给人造成...

"巴卡斯杯" 中国大学生程序设计竞赛 - 女生专场(重现)解题思路

"巴卡斯杯" 中国大学生程序设计竞赛 - 女生专场(重现)解题思路

HDU5713 2016"百度之星" - 复赛(Astar Round3)K个联通块

一道状态压缩动态规划题 题目描述:众所周知,度度熊喜欢图,尤其是联通的图。今天,它在图上又玩出了新花样,新高度。有一张无重边的无向图, 求有多少个边集,使得删掉边集里的边后,图里恰好有K个连通块。 ...

HDU 5715 & 2016"百度之星" 复赛(Astar Round3)1004 XOR游戏(Trie)

题意:众所周知,度度熊喜欢XOR运算(XOR百科)。 今天,它发明了一种XOR新游戏,最开始,它有一个长度为N的数组,度度熊可以任意添加分割线,将数组划分为M段,且每段长度小于等于L。 当然这是个...

2016"百度之星" - 复赛(Astar Round3)1003 1004

1003 拍照Problem Description小明在旅游的路上看到了一条美丽的河,河上有许多船只, 有的船只向左航行,有的船只向右航行。 小明希望拍下这一美丽的风景,并且把尽可能多的船只都完整地...

2016"百度之星" - 复赛(Astar Round3) 题解 (待续)

拍照扫描线 模版题#include #include #include #include #include #include #include #include #include #include u...

HDU 5713 K个联通块 状压dp枚举子集 (2016百度之星复赛)

题意 众所周知,度度熊喜欢图,尤其是联通的图。 今天,它在图上又玩出了新花样,新高度。有一张无重边的无向图, 求有多少个边集,使得删掉边集里的边后,图里恰好有K个连通块。  #inclu...

HDU 5690 All X的多种算法(2016"百度之星" - 初赛(Astar Round2A)1001)

传送门 All XTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total...

2016"百度之星" - 初赛(Astar Round2B)1003 瞬间移动(HDU5698)(逆元+快速幂)

2016”百度之星” - 初赛(Astar Round2B)1003 瞬间移动(HDU5698)(逆元+快速幂)Time Limit: 4000/2000 MS (Java/Others) ...

2016"百度之星" - 资格赛(Astar Round1)-(模拟+线段树+乘法逆元)

Problem A    Accepts: 1351    Submissions: 9951  Time Limit: 2000/1000 MS (Java/Others)    ...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 5713 & 2016"百度之星" 复赛(Astar Round3)1002 k个联通块
举报原因:
原因补充:

(最多只允许输入30个字)