結果

問題 No.137 貯金箱の焦り
コンテスト
ユーザー Min_25
提出日時 2016-07-10 12:27:01
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 41 ms / 5,000 ms
コード長 1,660 bytes
コンパイル時間 1,304 ms
コンパイル使用メモリ 94,864 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-07 18:17:10
合計ジャッジ時間 2,561 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize ("O3")
#pragma GCC target ("avx") // yukicoder
// #pragma GCC target ("sse4.2") // SPOJ, codechef

#include <cstdio>
#include <cassert>
#include <cmath>
#include <cstring>

#include <algorithm>
#include <iostream>
#include <vector>
#include <map>
#include <set>
#include <functional>
#include <tuple>

#define _rep(_1, _2, _3, _4, name, ...) name
#define rep2(i, n) rep3(i, 0, n)
#define rep3(i, a, b) rep4(i, a, b, 1)
#define rep4(i, a, b, c) for (int i = int(a); i < int(b); i += int(c))
#define rep(...) _rep(__VA_ARGS__, rep4, rep3, rep2, _)(__VA_ARGS__)

using namespace std;

using i64 = long long;
using u8 = unsigned char;
using u32 = unsigned;
using u64 = unsigned long long;
using f80 = long double; 

u32 A[100];
u32 poly[50010];

const u32 MOD = 1234567891;

u32 add_mod(u32 a, u32 b) {
  return u32(a += b) >= MOD ? a - MOD : a;
}

// overlap も考慮してベクトル化される．
void vec_add(int n, u32* a, u32* b) {
  for (int i = n - 1; i >= 0; --i) {
    a[i] = add_mod(a[i], b[i]);
  }
}

void solve() {
  i64 N, M;
  while (~scanf("%lld %lld", &N, &M)) {
    i64 s = 0;
    rep(i, N) scanf("%u", &A[i]), s += A[i];
    sort(A, A + N);
    int t = 0;
    poly[t++] = 1;
    while (M) {
      fill(poly + t, poly + s + t, 0);
      rep(i, N) {
        vec_add(t, poly + A[i], poly);
        t += A[i];
      }
      int p = M & 1;
      M >>= 1;
      t = (t + 1 - p) >> 1;
      rep(i, t) poly[i] = poly[2 * i + p];
    }
    printf("%u\n", poly[0]);
  }
}

int main() {
  clock_t beg = clock();
  solve();
  clock_t end = clock();
  fprintf(stderr, "%.3f sec\n", double(end - beg) / CLOCKS_PER_SEC);
  return 0;
}
0