#404784 (C++14) No.941 商とあまり

提出ソース

結果

問題	No.941 商とあまり
ユーザー	pekempey
提出日時	2019-12-04 21:40:51
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	TLE
実行時間	-
コード長	3,407 bytes
コンパイル時間	2,327 ms
コンパイル使用メモリ	189,500 KB
実行使用メモリ	22,252 KB
最終ジャッジ日時	2024-12-14 11:13:58
合計ジャッジ時間	56,485 ms
ジャッジサーバーID （参考情報）	judge1 / judge5

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	TLE	-
testcase_01	AC	4 ms 13,636 KB
testcase_02	AC	7 ms 16,020 KB
testcase_03	AC	4 ms 13,636 KB
testcase_04	AC	3 ms 14,444 KB
testcase_05	AC	3 ms 12,908 KB
testcase_06	AC	3 ms 13,636 KB
testcase_07	TLE	-
testcase_08	TLE	-
testcase_09	TLE	-
testcase_10	AC	35 ms 12,908 KB
testcase_11	AC	8 ms 13,632 KB
testcase_12	AC	8 ms 12,916 KB
testcase_13	AC	3 ms 13,640 KB
testcase_14	AC	7 ms 22,252 KB
testcase_15	AC	10 ms 13,636 KB
testcase_16	AC	12 ms 21,620 KB
testcase_17	AC	3 ms 13,644 KB
testcase_18	AC	3 ms 19,696 KB
testcase_19	AC	3 ms 14,824 KB
testcase_20	AC	3 ms 17,516 KB
testcase_21	AC	3 ms 14,700 KB
testcase_22	AC	3 ms 13,644 KB
testcase_23	AC	3 ms 13,640 KB
testcase_24	AC	2 ms 14,964 KB
testcase_25	AC	3 ms 13,636 KB
testcase_26	AC	11 ms 6,820 KB
testcase_27	AC	21 ms 6,816 KB
testcase_28	AC	8 ms 6,820 KB
testcase_29	AC	11 ms 6,820 KB
testcase_30	AC	11 ms 6,816 KB
testcase_31	AC	8 ms 6,820 KB
testcase_32	AC	27 ms 10,444 KB
testcase_33	AC	3 ms 6,820 KB
testcase_34	AC	3 ms 6,816 KB
testcase_35	AC	3 ms 6,816 KB
testcase_36	AC	9 ms 6,820 KB
testcase_37	AC	93 ms 6,816 KB
testcase_38	TLE	-
testcase_39	AC	3 ms 6,816 KB
testcase_40	AC	3 ms 6,816 KB
testcase_41	AC	3 ms 6,820 KB
testcase_42	AC	3 ms 6,820 KB
testcase_43	AC	4 ms 6,820 KB
testcase_44	TLE	-
testcase_45	AC	12 ms 6,984 KB
testcase_46	AC	11 ms 8,392 KB
testcase_47	AC	6 ms 6,816 KB
testcase_48	AC	3 ms 6,820 KB
testcase_49	AC	15 ms 6,820 KB
testcase_50	TLE	-
testcase_51	AC	889 ms 6,824 KB
testcase_52	TLE	-
testcase_53	AC	790 ms 6,820 KB
testcase_54	AC	12 ms 7,116 KB
testcase_55	AC	12 ms 7,244 KB
testcase_56	AC	14 ms 7,120 KB
testcase_57	AC	3 ms 6,816 KB
testcase_58	AC	3 ms 6,820 KB
testcase_59	TLE	-
testcase_60	TLE	-
testcase_61	TLE	-
testcase_62	AC	3 ms 6,824 KB
testcase_63	AC	3 ms 6,820 KB
testcase_64	AC	3 ms 6,816 KB
testcase_65	AC	4 ms 6,816 KB
testcase_66	TLE	-
testcase_67	AC	3 ms 6,820 KB
testcase_68	TLE	-
testcase_69	AC	4 ms 6,820 KB
testcase_70	AC	4 ms 6,820 KB
testcase_71	TLE	-
testcase_72	TLE	-
testcase_73	AC	3 ms 6,820 KB
testcase_74	AC	4 ms 6,820 KB
testcase_75	AC	4 ms 8,008 KB
testcase_76	AC	3 ms 6,820 KB
testcase_77	AC	4 ms 6,820 KB
testcase_78	AC	3 ms 6,820 KB
testcase_79	AC	3 ms 6,820 KB
testcase_80	AC	3 ms 6,816 KB
testcase_81	AC	4 ms 6,816 KB
testcase_82	AC	4 ms 6,816 KB
testcase_83	AC	3 ms 7,884 KB
testcase_84	AC	4 ms 6,816 KB
testcase_85	AC	4 ms 7,876 KB
testcase_86	AC	3 ms 6,816 KB
testcase_87	AC	3 ms 6,816 KB
testcase_88	AC	4 ms 6,820 KB
testcase_89	AC	4 ms 6,816 KB
testcase_90	AC	3 ms 6,816 KB
testcase_91	AC	4 ms 6,816 KB
testcase_92	AC	4 ms 6,820 KB
testcase_93	AC	3 ms 6,816 KB
testcase_94	AC	3 ms 6,816 KB
testcase_95	AC	3 ms 6,816 KB
testcase_96	AC	3 ms 6,820 KB
testcase_97	AC	3 ms 6,816 KB
testcase_98	AC	3 ms 6,816 KB
testcase_99	AC	3 ms 6,820 KB
testcase_100	AC	3 ms 6,816 KB
testcase_101	AC	3 ms 6,820 KB
testcase_102	AC	3 ms 6,820 KB
testcase_103	AC	3 ms 6,816 KB
testcase_104	AC	2 ms 6,820 KB
testcase_105	AC	2 ms 6,816 KB
testcase_106	AC	2 ms 6,816 KB
testcase_107	AC	3 ms 6,820 KB
testcase_108	AC	3 ms 6,820 KB
testcase_109	WA	-

権限があれば一括ダウンロードができます

ソースコード

raw source code

#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); }

bool can[500001];

set< vector<ll> > st;

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;
}

void dfs2(vector<ll> a) {
  if (a.size() == 1) {
    can[a[0]] = true;
    return;
  }
  sort(range(a));
  if (st.count(a)) return;
  st.insert(a);
  ll p = 1;
  rep(i, a.size()) p *= a[i] + 1;
  can[p - 1] = true;
  rep(i, a.size()) repe(j, i + 1, a.size()) {
    auto tmp = a;
    ll x = a[i];
    ll y = a[j];
    tmp.erase(tmp.begin() + j);
    tmp.erase(tmp.begin() + i);
    ll pp = p;
    pp /= x + 1;
    pp /= y + 1;
    set<ll> nxt;
    for (ll t = y + 2; (x * t + y + 1) * pp - 1 <= X; t++) {
      nxt.insert(x * t + y);
    }
    for (ll t = x + 2; (y * t + x + 1) * pp - 1 <= X; t++) {
      nxt.insert(y * t + x);
    }
    for (ll t : nxt) {
      tmp.push_back(t);
      dfs2(tmp);
      tmp.pop_back();
    }
  }
}

int main() {
  beg = clock();
  cin >> N >> X;
  rep(i, N) cin >> A[i];
  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;
    }
  }
  ll g = 0;
  rep(i, N) g = __gcd(g, A[i]);
  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;
  }
  if (p < 5000) {
    dfs2(vector<ll>(A, A + N));
    string ans(X, '0');
    for (int i = 1; i <= X; i++) {
      if (can[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;
  }
  cerr << "here" << endl;
  string ans(X, '0');
  int cont = 0;
  for (int i = p - 1; i <= X; i++) {
    if ((i - (p - 1)) % g != 0) continue;
    if (cont >= 100 || dfs(i, (1 << N) - 1)) {
      ans[i - 1] = '1';
      cont++;
    } else {
      cont = 0;
    }
  }
  cout << ans << endl;
  cerr << cnt << endl;
}

yukicoder

結果

テストケース

ソースコード