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

原创 2016年05月31日 22:15:34

题意


众所周知,度度熊喜欢图,尤其是联通的图。

今天,它在图上又玩出了新花样,新高度。有一张无重边的无向图, 求有多少个边集,使得删掉边集里的边后,图里恰好有K个连通块。 

推荐题目:Codeforces #332 div 2 E Nuts (599E) 状压+树形+计数


#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
#include <vector>
#include <iostream>
#include <map>
#define pii pair<int,int>
#define mp(a,b) make_pair<a,b>
#define xx first
#define yy second
#define mem(a) memset( a, 0, sizeof(a) )
using namespace std;
typedef long long ll;

ll dp[20][1<<15], f[1<<15], g[1<<15], f2[1000];
int fa[20], vis[1<<15];
int w[1<<15], n, m, k;
const ll mod = 1e9+9;

/*
dp[i][st]存的是现在有i个连通块,加入的点的情况是st
f[st]存的是这个连通块有多少个作为连通块的形态数(就是这个连通块包含的边怎么个删法使他还是个连通块,有多少种情况)
g[st]辅助数组
抓住  第一个点一定属于一个连通分量(所以通俗来说枚举子集时只需要枚举一半即可)

*/
int root( int x )
{
    if( fa[x] == x ) return x;
    else return root( fa[x] );
}

void merge( int a, int b )
{
    if( root(a) == root(b) ) return ;
    int pa, pb;
    pa = root(a);
    pb = root(b);
    fa[pa] = pb;
    return ;
}

int lowbit( int x )
{
   return x&-x;
}


void solve()
{
    int i, j, st;
    cin >> n >> m >> k;
    mem(f);mem(g);mem(vis);
    for( i = 0; i < n; i ++ )fa[i] = i;
    for( i = 1; i <= m; i ++ ){
        int x, y;
        scanf("%d %d", &x, &y);
        x --; y -- ;
        for( j = 0; j < (1<<n); j ++ ){
            if( ( j&(1<<x) ) && ( j&(1<<y) ) ){
                f[j] ++;
            }
        }
        merge( x, y );
    }// 用并查集判连通
    for( st = 1; st < (1<<n); st ++ ){
        int p = lowbit(st), to = -1, fail = 0;
        for( i = 0; i < n; i ++ ){
            if( st & ( 1 << i ) ){
                if( to == -1 ){
                    to = root(i);
                }
                else{
                    if( root(i) != to ){
                        fail = 1;
                        break;
                    }
                }
            }
        }
        if( fail ){
            continue;
        }
        else vis[st] = 1;
        f[st] = f2[f[st]];
        g[st] = f[st];
        for( j = st; j ; j = (j-1) & st ){
            if( j == st ) continue;
            if( !vis[j] ) continue;
            if( j & p ){
                f[st] -= ((ll)f[j]*g[j^st])%mod;
                f[st] %= mod;
            }
        }
    }// vis数组判这个子集本身是否连通,不联通f g = 0(可以不用判断)
    memset( dp, 0, sizeof(dp) );
    dp[0][0] = 1;
    for( i = 1; i <= k; i ++ ){
        for( st = 1; st < (1<<n); st ++ ){
            int p = lowbit(st), sub;
            for( sub = st; sub; sub = ( sub - 1 )&st ){
                if( !vis[sub] ) continue;
                if( p & sub ){
                    dp[i][st] += ( (ll)dp[i-1][sub^st]*f[sub] )%mod;
                    dp[i][st] %= mod;
                }
            }
        }
    }
    ll ans = (dp[k][(1<<n)-1]+mod)%mod;
    printf("%lld\n", ans);
}

void init()
{
    f2[0] = 1;
    for( ll i = 1; i <= 999; i ++ ){
        f2[i] = f2[i-1]*2%mod;
    }
}

int main()
{
    int cas = 1, T;
    init();
    cin >> T;
    while( T -- ){
        printf("Case #%d:\n", cas ++);
        solve();
    }
}


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

相关文章推荐

UVA 11825 - Hackers' Crackdown 状态压缩 dp 枚举子集

11825 - Hackers' Crackdown 状态压缩 dp 枚举子集 ACM 题目地址:11825 - Hackers' Crackdown 题意:  有一个由编号0~n-1的n台计算...
  • hcbbt
  • hcbbt
  • 2014-06-27 21:35
  • 1465

HDU Untitled(状压DP OR dfs枚举子集)

Untitled Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total...

我是如何成为一名python大咖的?

人生苦短,都说必须python,那么我分享下我是如何从小白成为Python资深开发者的吧。2014年我大学刚毕业..

hdu-5823-color II-状压DP-枚举子集的子集

http://acm.hdu.edu.cn/showproblem.php?pid=5823 给你一个无向图,然后要枚举所有点的子集的最小染色方案数,然后输出 ∑ans[i]*233^i m...

UVA 11825 Hackers' Crackdown dp+位运算状态压缩(子集枚举)

题意:       有N台电脑(编号为0,1,2......n-1),每台电脑有N(N Input There will be multiple test cases in the inpu...

【CDOJ 1355】柱爷与三叉戟不得不说的故事 【状压DP+子集枚举】

柱爷与三叉戟不得不说的故事 Time Limit: 500/500MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Othe...

二进制法枚举子集

二进制法枚举子集:  n个数有2^n个子集,每个子集对应一个二进数,每位二进制对应一个数。二进制的位权为0表示子集不包含那个数,二进制的位权为1表示子集包含那个数。 案例  2  1 ...

hdu 5713 K个联通块 2016百度之星复赛1002 DP

用dp[i][j] 表示 状态i分成k个联通块的方案数。状态i的二进制第j位代表编号为j的点。 f[i]表示状态i内删边,仍保持一个联通块的方案数。 则有 转移方程 dp[i][k] = sum(dp...

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

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

HDU 5714 拍照 优先队列前后扫 (2016百度之星复赛)

题意: 小明在旅游的路上看到了一条美丽的河,河上有许多船只,有的船只向左航行,有的船只向右航行。小明希望拍下这一美丽的风景,并且把尽可能多的船只都完整地拍到一张照片中。 小明位于河的边...

2016百度之星复赛 1003 拍照 扫描线

原题见1004有n条船平行于x轴,已知初始时刻的位置,向平行于x轴正向或反向走。所有船速度大小一致。现在海哥站在x轴上,要给船拍照,角度为朝向+y方向的固定90度,可以在任意时刻、任意位置拍,问何时拍...
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

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