1. 程式人生 > >HDU-6321___Dynamic Graph Matching——解題報告 狀壓DP

HDU-6321___Dynamic Graph Matching——解題報告 狀壓DP

題目大意:

  給出n個點和m個操作,每次操作可以可以在兩個點之間連線或者刪除這兩個點之間的連線,問匹配數為1、2...n/2的邊的數量,意思就是求出互不相交的j條邊的組數,j為1~n/2。

解題思路:

  用圖的思想來理解點,首先最多有10個點,每個點都有被佔用和未被佔用兩種狀態,所以一共有2^10種狀態。

  每次增加一條邊,則在0~1023種狀態中找,a和b在二進位制位數下都是0的狀態,然後更新這個狀態的a和b相應位數都變為1之後的新狀態

程式碼思路:

  1、用num[a] 表示1<<(a-1),用cnt[i] 表示二進位制下i的1的數目

  2、dp[ i+num[a]+num[b] ]直接加上當前dp[ i ]的值並取模

  3、同理,若是減邊,則是dp[ i ]減去dp[ i-num[a]-num[b] ]的值並取模

核心:狀態轉換為二進位制,並不斷地用原狀態來更新下一個狀態,記錄答案,注意取模時為負數!

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll mod = 1e9+7;
ll dp[1025], cnt[1025];    //記錄每一種情況 
ll num[15];     
ll ans[15];    //計算成功匹配邊數的數目 
int n;

void add(int a, int b)
{
    for(int i=0; i<(1<<n); i++)
    {
        if((i&num[a])==0 && (i&num[b])==0)    //a和b兩個點都沒佔用 
        {    
            if(cnt[i]%2==0)    //記錄該種情況加邊後的dp 並更新答案 
            {
                dp[i+num[a]+num[b]] = (dp[i+num[a]+num[b]] + dp[i])%mod;
                ans[cnt[i]/2+1] = (ans[cnt[i]/2+1]+dp[i])%mod;
            }
        }
    }
}

void sub(int a, int b)
{
    for(int i=0; i<(1<<n); i++)    //a和b兩個點都佔用  
    {
        if((i&num[a])>0 && (i&num[b])>0)
        {
            if(cnt[i]%2==0)    //記錄該種情況減邊後的dp 並更新答案 
            {
                dp[i] = (dp[i] - dp[i-num[a]-num[b]]+mod)%mod;
                ans[cnt[i]/2] = (ans[cnt[i]/2] - dp[i-num[a]-num[b]]+mod)%mod;
            }
        }
    }
}

int main()
{
    int caz;
    scanf("%d", &caz);
    for(int i=0; i<11; i++) 
        num[i]=1<<i;
    for(int i=1;i<=1023;i++)
        cnt[i]=__builtin_popcount(i);
    
    while(caz--)
    {
        int m, a, b;
        scanf("%d %d", &n, &m);
        memset(dp, 0, sizeof(dp));
        memset(ans, 0, sizeof(ans));
        dp[0]=1; //最初開始加邊時從0開始加 
        char ch[10];
        for(int i=0; i<m; i++)
        {
            scanf("%s%d%d", ch, &a, &b);
            if(ch[0]=='+') add(a-1, b-1);
            else sub(a-1, b-1);
            for(int i=1; i<=n/2; i++)
            { 
                printf("%lld", ans[i]);
                if(i==n/2) printf("\n");
                else printf(" ");
             } 
        }
    }
    return 0;
}