Topcoder SRM 663 Div2 Hard: CheeseRolling（状压DP）

标签：

{"NN",
 "YN"}

Returns: {0, 2 }

{"NYNY",
 "NNYN",
 "YNNY",
 "NYNN"}

Returns: {8, 0, 16, 0 }

大致题意：

<=16个人举办如图所示的淘汰赛，已给出每两个人之间对决的输赢结果矩阵

求对于每个人来说，有多少种安排位置方案使这个人最后获胜

思路：

状压DP，还是有点复杂的状压

对于第x个人最终获胜有多少种方式可以表示为dp[x][(1<<n)-1]

然后递归向下，枚举x能战胜的对手y

然后把集合S分为两半，一半包含x，一半包含y

则有：dp[x][S] += dp[x][S1]*dp[y][S2]*2(相对位置可以变)

枚举半个集合显然就是枚举子集，所以总的复杂度是O(n*3^n)

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <iostream>
#include <cstring>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <string>
#include <vector>
#include <cstdio>
#include <ctime>
#include <bitset>
#include <algorithm>
#define SZ(x) ((int)(x).size())
#define ALL(v) (v).begin(), (v).end()
#define foreach(i, v) for (__typeof((v).begin()) i = (v).begin(); i != (v).end(); ++ i)
#define reveach(i, v) for (__typeof((v).rbegin()) i = (v).rbegin(); i != (v).rend(); ++ i)
#define REP(i,n) for ( int i=1; i<=int(n); i++ )
#define rep(i,n) for ( int i=0; i<int(n); i++ )
using namespace std;
typedef long long ll;
#define X first
#define Y second
typedef pair<int,int> pii;

template <class T>
inline bool RD(T &ret) {
    char c; int sgn;
    if (c = getchar(), c == EOF) return 0;
    while (c != '-' && (c<'0' || c>'9')) c = getchar();
    sgn = (c == '-') ? -1 : 1;
    ret = (c == '-') ? 0 : (c - '0');
    while (c = getchar(), c >= '0'&&c <= '9') ret = ret * 10 + (c - '0');
    ret *= sgn;
    return 1;
}
template <class T>
inline void PT(T x) {
    if (x < 0) {
        putchar('-');

        x = -x;
    }
    if (x > 9) pt(x / 10);
    putchar(x % 10 + '0');
}

ll dp[20][1<<16];
int bit[1<<16];

int cal(int x){
        int ans = 0;
        while(x) ans += (x&1), x >>= 1;
        return ans;
}
int n;
vector<string> mp;
ll DP(int u,int state){
        ll &ans = dp[u][state];
        if( ans != -1 ) return ans;
        if( bit[state] == 2 ){
                rep(v,n)
                        if( (1<<v)&state && mp[u][v] == 'Y')
                               return ans = 2;
                return ans = 0;
        }
        ans = 0;
        int c = (bit[state]-2)/2;
        rep(v,n){
                if( (1<<v)&state && mp[u][v] == 'Y' ){
                        int all = state & (~(1<<u)) & (~(1<<v));
                        for( int tmp = state; tmp ; tmp = (tmp-1)&all){
                                if( bit[tmp] != c) continue;
                                ans += 2*DP(u, ((1<<u)|tmp) )*DP(v, ((1<<v)|all)&(~tmp) );
                        }
                }
        }
        return ans;
}
class CheeseRolling {
public:
   vector<long long> waysToWin( vector <string> wins ) {
        memset(dp,-1,sizeof(dp));
        foreach(it,wins) mp.push_back(*it);
        n = SZ(wins);
        for(int i = 0; i < (1<<16);i++) bit[i] = cal(i);
        vector<ll> ans;
        rep(i,n) ans.push_back(DP(i,(1<<n)-1));
        return ans;
   }
};

Topcoder SRM 663 Div2 Hard: CheeseRolling（状压DP）

标签：

原文地址：http://blog.csdn.net/kalilili/article/details/47360187