「GXOI / GZOI2019」宝牌一大堆 （DP）

题意

​​LOJ传送门​​

题解

可以发现「七对子」 和 「国士无双」直接暴力就行了。

唯一的就是剩下的"3*4+2"。

考试的时候写了个爆搜剪枝，开了O2有50pts。写的时候发现可以DP，但是没写。

然后下来写了发现就4个转移。。。

用

$dp[i][j][k][a][b][c]$

表示当考虑前

$i$

张牌，有

$j$

个雀头，

$k$

个面子，

$i-2$

用了

$a$

张，

$i-1$

用了

$b$

张，

$i$

用了

$c$

张时，前

$i-3$

张牌的最大分数。

注意这里是"前

$i-3$

张牌的最大分数"，也就是没有计算

$a,b,c$

这最后

$3$

种牌的贡献。

有个性质就是杠子一定不会选，因为选杠子一定劣于面子。

因为

$C(4,4)=1<C(4,3)=4$

，就算这张牌是宝牌，杠子的贡献是

$C(4,4)*2^4<C(4,3)*2^3$

所以我们没有考虑杠子。

上代码

CODE

luogu上开了O2才过

#include <bits/stdc++.h>
using namespace std;
inline int getid(char *s) {
  if(isdigit(s[0])) {
    int re = s[0]-'0';
    if(s[1] == 'm') return re;
    if(s[1] == 'p') return 9+re;
    if(s[1] == 's') return 18+re;
  }
  if(s[0] == 'E') return 28;
  if(s[0] == 'S') return 29;
  if(s[0] == 'W') return 30;
  if(s[0] == 'N') return 31;
  if(s[0] == 'Z') return 32;
  if(s[0] == 'B') return 33;
  if(s[0] == 'F') return 34;
}
typedef long long LL;
int cnt[50], bel[50];
LL ans, c[5][5], pw6[10];
bool bao[50];

int arr13[14] = { 0, 1, 9, 10, 18, 19, 27, 28, 29, 30, 31, 32, 33, 34 }; //指定的13张牌的编号
void solve13() {
  bool flg = 1;
  for(int i = 1; i <= 13; ++i) flg &= bool(cnt[arr13[i]]);
  if(flg) {
    LL now = 13;
    for(int i = 1; i <= 13; ++i)
      now = now * c[cnt[arr13[i]]][1] * (bao[arr13[i]] ? 2 : 1);
    for(int i = 1; i <= 13; ++i)
      ans = max(ans, now / c[cnt[arr13[i]]][1] * c[cnt[arr13[i]]][2] * (bao[arr13[i]] ? 2 : 1));
  }
}

void dfs(int i, int j, LL now) { //暴力枚举七对子
  if(j == 7) { ans = max(now, ans); return; }
  if(i > 34) return;
  LL tmp = now*pw6[7-j]*(1ll<<((7-j)<<1)); //剪枝，预估之后的最大值，如果小于等于当前答案就不用了搜下去了
  if(tmp <= ans) return;
  if(cnt[i] >= 2) {
    cnt[i] -= 2;
    dfs(i+1, j+1, now*c[cnt[i]+2][2]*(bao[i]?4:1));
    cnt[i] += 2;
  }
  if(tmp <= ans) return;
  dfs(i+1, j, now);
}
void solve7() { dfs(1, 0, 7); }

LL dp[36][2][5][5][5][5];
inline void upd(LL &x, LL y) { x = y > x ? y : x; }
inline int calc(int i, int j) { return i > 0 ? c[cnt[i]][j] * (bao[i] ? (1<<j) : 1) : 1; } //计算i选j张的贡献
void solve() {
  memset(dp, 0, sizeof dp);
  dp[1][0][0][0][0][0] = 1;
  for(int i = 1; i <= 34; ++i)
    for(int j = 0; j < 2; ++j)
      for(int k = 0; k <= 4; ++k)
        for(int a = 0; a <= 4; ++a)
          for(int b = 0; b <= 4; ++b)
            for(int c = 0; c <= 4; ++c) {
              LL val = dp[i][j][k][a][b][c];
              if(val) {
                upd(dp[i+1][j][k][b][c][0], val * calc(i-2, a)); 
                if(!j && cnt[i]-c >= 2)
                  upd(dp[i][j+1][k][a][b][c+2], val);
                if(k<4 && cnt[i]-c >= 3)
                  upd(dp[i][j][k+1][a][b][c+3], val);
                if(k<4 && bel[i] && i>2 && bel[i]==bel[i-2] && cnt[i]-c >= 1 && cnt[i-1]-b >= 1 && cnt[i-2]-a >= 1)
                  upd(dp[i][j][k+1][a+1][b+1][c+1], val);
              }
          }
  for(int a = 0; a <= 4; ++a)
    for(int b = 0; b <= 4; ++b)
      for(int c = 0; c <= 4; ++c) //最后答案算上a,b,c
        upd(ans, dp[34][1][4][a][b][c] * calc(32, a) * calc(33, b) * calc(34, c)); 
}

int main () {
  //freopen("doraippai.in", "r", stdin);
//  freopen("doraippai.out", "w", stdout);
  for(int i = 1; i <= 9; ++i) bel[i] = 1;
  for(int i = 10; i <= 18; ++i) bel[i] = 2;
  for(int i = 19; i <= 27; ++i) bel[i] = 3;
  c[0][0] = 1;
  for(int i = 1; i <= 4; ++i) {
    c[i][0] = c[i][i] = 1;
    for(int j = 1; j < i; ++j)
      c[i][j] = c[i-1][j-1] + c[i-1][j];
  }
  pw6[0] = 1;
  for(int i = 1; i <= 7; ++i) pw6[i] = 6ll * pw6[i-1];
  int T; scanf("%d", &T); while(T--) {
    for(int i = 1; i <= 34; ++i)
      cnt[i] = 4, bao[i] = 0;
    char s[5];
    while(scanf("%s", s), s[0] != '0') {
      int num = getid(s);
      --cnt[num];
    }
    while(scanf("%s", s), s[0] != '0') {
      int num = getid(s);
      bao[num] = 1;
    }
    ans = 0;
    solve13();
    solve7();
    solve();
    printf("%lld\n", ans);
  }
}