SRM667 div1 easy OrderOfOperations

考えたこと

• bitDPをしてくださいみたいな雰囲気をしている
• 制約的に考えると dp[S] = (集合Sのセルを実行しているときの最小コスト) みたいな持ち方になりそう
• n個全部を実行したことについての判定をどうするかちょっと悩んだけど命令s[i]の和集合が全て実行されていれば命令をコスト0で実行できるのでどうでもいい
• O(n2^m)で投げたら通る

かなりTLぎりぎりで怖い
計算量5*10^7くらいな気がしたけどpopcountの定数倍が重い…？

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef vector<int> VI;
typedef vector<VI> VVI;
typedef vector<ll> VL;
typedef vector<VL> VVL;
typedef pair<int, int> PII;

#define FOR(i, a, n) for (ll i = (ll)a; i < (ll)n; ++i)
#define REP(i, n) FOR(i, 0, n)
#define ALL(x) x.begin(), x.end()
#define IN(a, b, x) (a<=x&&x<b)
#define MP make_pair
#define PB push_back
const int INF = (1LL<<30);
const ll LLINF = (1LL<<60);
const double PI = 3.14159265359;
const double EPS = 1e-12;
const int MOD = 1000000007;
//#define int ll

template <typename T> T &chmin(T &a, const T &b) { return a = min(a, b); }
template <typename T> T &chmax(T &a, const T &b) { return a = max(a, b); }

int dx[] = {0, 1, 0, -1}, dy[] = {1, 0, -1, 0};

int dp[1<<20];
class OrderOfOperations {
   public:
   int minTime(vector <string> s)
  {
    int n = s.size(), m = s[0].size();
    VI v(n);
    int all = 0;
    REP(i, n) {
      REP(j, m) {
        if(s[i][j] == '1') v[i] |= 1<<j;
      }
      all |= v[i];
    }

    REP(i, 1<<m) dp[i] = INF;
    dp[0] = 0;

    REP(i, 1<<m) {
      REP(j, n) {
        // iで0でv[j]で1
        int tmp = __builtin_popcount(i|v[j]) - __builtin_popcount(i);
        chmin(dp[i | v[j]], dp[i] + tmp*tmp);
      }
    }

    return dp[all];
  }
};