[HDOJ5890]Eighty seven（暴力，dp，bitset）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5890

题意：50个数，10W个询问，每次问删掉第i,j,k个数后，是否存在一种选10个数和为87的方案，只需要输出 ’Yes’ 或者 ’No’

题解：暴力：不同的询问大概2W个，每个暴力bitset DP，抠一抠能卡着过。优化1：先求出一组解，如果询问和解没交就是’Yes’，否则暴力，不同的询问大概1W个；优化2：先预处理出所有询问的答案，能方便的复用之前的DP数组，不用每次从头开始重新求。

学了一波bitset的用法，好机智，巧妙地记下了选i个数的所有可能的和的情况。

dp(i,j)表示选择i个数，j为和的种类。

dp(j)|=dp(j-1)<<a(i)。

一开始想直接背包做，发现存不下那么大的结果。

#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
using namespace std;
#define fr first
#define sc second
#define cl clear
#define BUG puts("here!!!")
#define W(a) while(a--)
#define pb(a) push_back(a)
#define Rint(a) scanf("%d", &a)
#define Rll(a) scanf("%I64d", &a)
#define Rs(a) scanf("%s", a)
#define Cin(a) cin >> a
#define FRead() freopen("in", "r", stdin)
#define FWrite() freopen("out", "w", stdout)
#define Rep(i, len) for(int i = 0; i < (len); i++)
#define For(i, a, len) for(int i = (a); i < (len); i++)
#define Cls(a) memset((a), 0, sizeof(a))
#define Clr(a, x) memset((a), (x), sizeof(a))
#define Full(a) memset((a), 0x7f7f7f, sizeof(a))
#define lrt rt << 1
#define rrt rt << 1 | 1
#define pi 3.14159265359
#define RT return
#define lowbit(x) x & (-x)
#define onecnt(x) __builtin_popcount(x)
typedef long long LL;
typedef long double LD;
typedef unsigned long long ULL;
typedef pair<int, int> pii;
typedef pair<string, int> psi;
typedef pair<LL, LL> pll;
typedef map<string, int> msi;
typedef vector<int> vi;
typedef vector<LL> vl;
typedef vector<vl> vvl;
typedef vector<bool> vb;

const int maxn = 55;
int n;
int a[maxn];
bool sta[maxn][maxn][maxn];
bitset<88> dp[11];

bool ok(int x, int y, int z) {
  Rep(i, 11) dp[i].reset();
  dp[0][0] = 1;
  For(i, 1, n+1) {
    if(i == x || i == y || i == z || a[i] > 87) continue;
    for(int j = 10; j >= 1; j--) {
      dp[j] |= dp[j-1] << a[i];
    }
  }
  return dp[10][87];
}

signed main() {
    //FRead();
  int T, q, x[5];
  Rint(T);
  W(T) {
    Rint(n);
    For(i, 1, n+1) Rint(a[i]);
    For(i, 1, n+1) {
      For(j, i, n+1) {
        For(k, j, n+1) {
          sta[i][j][k] = ok(i, j, k);
        }
      }
    }
    Rint(q);
    W(q) {
      Rep(i, 3) Rint(x[i]);
      sort(x, x+3);
      if(sta[x[0]][x[1]][x[2]]) puts("Yes");
      else puts("No");
    }
  }
    RT 0;
}