結果

問題 No.941 商とあまり
コンテスト
ユーザー pekempey
提出日時 2019-12-04 21:09:19
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 2,284 bytes
コンパイル時間 2,167 ms
コンパイル使用メモリ 175,140 KB
実行使用メモリ 13,768 KB
最終ジャッジ日時 2024-12-14 10:14:46
合計ジャッジ時間 77,532 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample RE * 6
other AC * 71 RE * 17 TLE * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
using ll = long long;
 
#define rep(i, n)      for (int i = 0; i < (n); i++)
#define repr(i, n)     for (int i = (n) - 1; i >= 0; i--)
#define repe(i, l, r)  for (int i = (l); i < (r); i++)
#define reper(i, l, r) for (int i = (r) - 1; i >= (l); i--)
#define repi(i, l, r)  for (int i = (l); i <= (r); i++)
#define repir(i, l, r) for (int i = (r); i >= (l); i--)
#define range(a) a.begin(), a.end()
void initio() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); }

bitset<500001> dp[1 << 11];
bitset<500001> checked[1 << 11];
int N;
ll A[500010];
ll prod[1 << 11];
int root[1 << 11];
int cnt;
ll X;

clock_t beg;

bool dfs(int x, int u) {
  u = root[u];
  if (checked[u][x]) return dp[u][x];
  cnt++;
  checked[u][x] = true;
  if (prod[u] - 1 > x) return dp[u][x] = false;
  if ((u & u - 1) == 0) return dp[u][x] = prod[u] - 1 == x;
  rep(i, N) if (u >> i & 1) {
    if ((x - (prod[u] - 1)) % A[i] == 0) return dp[u][x] = true;
  }
  for (int i = (u - 1) & u; i != 0; i = (i - 1) & u) {
    if (root[i] != i) continue;
    int j = u ^ i;
    for (int k = prod[j]; k <= x && prod[i] - 1 <= x / k; k++) if (x % k != 0) {
      if (dfs(x % k, j) && dfs(x / k, i)) {
        return dp[u][x] = true;
      }
    }
  }
  return dp[u][x] = false;
}

int main() {
  beg = clock();
  cin >> N >> X;
  rep(i, N) cin >> A[i];
  if (N <= 3) {
    abort();
  }
  sort(A, A + N);
  ll p = 1;
  rep(i, N) {
    p *= A[i] + 1;
    if (p - 1 > X) {
      cout << string(X, '0') << endl;
      return 0;
    }
  }
  if (count(A, A + N, 1) > 0 && N >= 2) {
    string ans(X, '0');
    for (int i = p - 1; i <= X; i++) {
      ans[i - 1] = '1';
    }
    cout << ans << endl;
    return 0;
  }
  rep(i, 1 << N) prod[i] = 1;
  rep(i, 1 << N) rep(j, N) if (i >> j & 1) {
    prod[i] *= A[j] + 1;
  }
  rep(i, 1 << N) {
    vector<ll> a;
    rep(j, N) if (i >> j & 1) a.push_back(A[j]);
    int mask = 0;
    int k = 0;
    rep(j, N) {
      if (k < a.size() && a[k] == A[j]) {
        mask |= 1 << j;
        k++;
      }
    }
    root[i] = mask;
  }
  string ans(X, '0');
  for (int i = p - 1; i <= X; i++) {
    if (dfs(i, (1 << N) - 1)) {
      ans[i - 1] = '1';
    }
  }
  cout << ans << endl;
  cerr << cnt << endl;
}
0