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

原创 2016年05月30日 17:52:26

用dp[i][k] 表示 状态i分成k个联通块的方案数。状态i的二进制第j位代表编号为j的点。

f[i]表示状态i内删边,仍保持一个联通块的方案数。

则有 转移方程 dp[i][k] = sum(dp[i-j][k-1]*f[j]).   j 是i的子状态,但不是i的所有子状态。为什么呢?

假设i 表示 点集 {1,2,3,4} 然后分成 {1,2} {3,4} 如果枚举所有子状态,则在j = {1,2}的时候算了一次这种方案, j = {3,4}的时候又算了一次这种方案。

处理方法是,i状态分成k个联通块,k个联通块中,必有一个包含i中编号最小的点。则我们只枚举包含最小编号点的子状态j。避免了重复


f[i]怎么求?

直接求删边后仍是一个连通块的方案较难。用总方案数减去删边后被分割的方案。

假设g[i]是状态i删边后不是一个联通块的方案数。 则f[i] = 2^cnt[i] - g[i];    (cnt[i]是状态i内的边数,易求)

g[i] = sum(2^cnt[i-j]*f[j]). j是i的子方案(和上面dp方程一样,不是所有子方案,而是包含i中最小点的子方案)

------------------------------------------------------------------------------------

f[i] 其实就是 dp[i][1]

#include<stdio.h>
#include<string.h>
#include<ctype.h>
#include<algorithm>
#include<vector>
using namespace std;
void fre(){freopen("t.txt","r",stdin);}
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b>a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b<a)a = b; }
#define MS(x,y) memset(x,y,sizeof(x))
#define MC(x, y) memcpy(x, y, sizeof(x))
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
const int INF = 0x3f3f3f3f;
const int M = 1000000009;
const double eps = 1e-8;

int n,m,K;
int cnt[20000],e[15][15];
LL dp[20000][15],f[20000],g[20000],bit[150];
int calc(int s,int p)
{
    int ret = 0;
    p = lower_bound(bit,bit+15,p)-bit;
    for(int i = 0; i < n; ++i) if(s&bit[i]) ret += e[i][p];
    return ret;
}
void solve()
{
    int u,v;
    scanf("%d%d%d",&n,&m,&K);
    MS(e,0);MS(g,0);MS(dp,0),MS(f,0);
    while(m--)
    {
        scanf("%d%d",&u,&v);
        e[u-1][v-1]++;
        if(u!=v) e[v-1][u-1]++;
    }

    for(int i = 1; i < bit[n]; ++i) cnt[i] = cnt[i-(i&-i)] + calc(i,i&-i);
    for(int i = 1; i < bit[n]; ++i)
    {
        int j = i;
        do
        {
            j = j-1&i;
            if(j&(i&-i)) g[i] = (g[i] + bit[cnt[i-j]]*f[j])%M;
        }while(i!=j);
        f[i] = (bit[cnt[i]]-g[i])%M;//刚开始没注意到这里会减成负数,于是wa啦
    }
    dp[0][0] = 1;
    for(int i = 1; i < bit[n]; ++i)
    {
        for(int k = 1; k <= K; ++k)
        {
            int j = i;
            do
            {
                j = j-1&i;
                if(j&(i&-i)) dp[i][k] = (dp[i][k] + dp[i-j][k-1]*f[j]%M)%M;
            }while(i!=j);
        }
    }
    printf("%I64d\n",(dp[bit[n]-1][K]+M)%M);
}
int main()
{
    //fre();
    bit[0] = 1; for(int i = 1; i <= 115; ++i) bit[i] = (bit[i-1]<<1)%M;
    int T,cas = 0;
    scanf("%d",&T);
    while(T--)
        printf("Case #%d:\n",++cas) ,solve();
    return 0;
}


版权声明:转载注明出处。有任何疑问可留言。

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

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

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

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

2017百度之星 复赛

作为一只菜鸡,这可能是我OI生涯中最后一次打百度之星了,希望明年这个时候能不退役(虽然不大可能)开始贴题解吧。会持续更新。1001 Arithmetic of Bomb 煞笔模拟题,1A#inclu...
  • qq_35866453
  • qq_35866453
  • 2017年08月18日 21:36
  • 395

HDU 5713 K个联通块(状压DP)

Description众所周知,度度熊喜欢图,尤其是联通的图。今天,它在图上又玩出了新花样,新高度。有一张无重边的无向图, 求有多少个边集,使得删掉边集里的边后,图里恰好有K个连通块。 Input第一...
  • V5ZSQ
  • V5ZSQ
  • 2017年08月29日 09:28
  • 209

百度之星2014-2016资格赛试题及解答

百度之星2014-2016资格赛试题、解答及代码实现
  • u014800380
  • u014800380
  • 2017年07月05日 11:34
  • 544

2017百度之星复赛1001 Arithmetic of Bomb

算是想法题吧  写的有点乱 #include #include #include #include #include #include #include #include #include #in...
  • Frozensmile
  • Frozensmile
  • 2017年08月18日 17:04
  • 133

2016"百度之星" - 复赛(Astar Round3)1001 D++游戏 hdu 5712

这道题一看当时一猜不是dp就是搜 但是无论是搜还是dp都要解决一个问题就是 一个序列被删除后剩下的两部分是连在一起的,这个时候怎么搞 当时我想的是枚举一开始的等差数列那么再去求等差数列的相邻两项的值就...
  • meixiuxiudd
  • meixiuxiudd
  • 2016年05月30日 13:47
  • 698

百度之星2017初赛题解(A)

T1: 简单数论,问满足(a0+a1*B+...+an*B^n)=a0+a1+...+an(mod P)的P的个数 即满足P|(B-1)a1+(B^2-1)a2+(B^3-1)a3+...,即P|...
  • yfzcsc
  • yfzcsc
  • 2017年08月13日 21:14
  • 465

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

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

2016"百度之星" - 初赛(Astar Round2A)解题报告

2016"百度之星" - 初赛(Astar Round2A)解题报告
  • queuelovestack
  • queuelovestack
  • 2016年05月21日 19:33
  • 6027
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:hdu 5713 K个联通块 2016百度之星复赛1002 DP
举报原因:
原因补充:

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