Yukicoder No.54 Happy Hallowe'en

,

http://yukicoder.me/problems/no/54

ひどい嘘だなあと思ったので後から落とされて通知が飛ぶのを回避するために想定解も提出しておいたが、嘘解法の方が面白いのでこちらを貼っておく。

solution

嘘解法。使用順序を乱択で決定しbitsetで加速したDPで殴る。 サンプル$2$にもあるような以下のような入力には弱いが、これは$T_i$で整列すれば倒せるのでやっておく。

10000
1 1
1 2
1 3
1 4
1 5
1 6
1 7
1 8
...

乱択は完全に嘘。 ほとんどの部分でサンプル$2$のようになっているがサンプル$4$のように$T_i$の貪欲では倒せない部分が混じっているケースがあると死ぬだろうが、これは作れるはず。

implementation

http://yukicoder.me/submissions/142689

#include <iostream>
#include <vector>
#include <algorithm>
#include <numeric>
#include <bitset>
#include <memory>
#include <chrono>
#define repeat(i,n) for (int i = 0; (i) < int(n); ++(i))
#define whole(f,x,...) ([&](decltype((x)) whole) { return (f)(begin(whole), end(whole), ## __VA_ARGS__); })(x)
using namespace std;
template <class T> void setmax(T & a, T const & b) { if (a < b) a = b; }
template <class T> void setmin(T & a, T const & b) { if (b < a) a = b; }
int main() {
    chrono::high_resolution_clock::time_point clock_begin = chrono::high_resolution_clock::now();
    int n; cin >> n;
    vector<int> v(n), t(n); repeat (i,n) cin >> v[i] >> t[i];
    int ans_upper = 0;
    repeat (i,n) setmax(ans_upper, t[i]-1 + v[i]);
    setmin(ans_upper, whole(accumulate, v, 0));
    int ans = 0;
    auto dp = make_unique<bitset<2*10000+1> >();
    auto it = make_unique<bitset<2*10000+1> >();
    vector<int> ix(n); whole(iota, ix, 0);
    whole(iota, ix, 0);
    whole(sort, ix, [&](int i, int j) { return t[i] < t[j]; });
    while (true) {
        chrono::high_resolution_clock::time_point clock_end = chrono::high_resolution_clock::now();
        if (chrono::duration_cast<chrono::milliseconds>(clock_end - clock_begin).count() >= 4900) break;
        dp->reset();
        (*dp)[0] = true;
        for (int i : ix) {
            it->set();
            *it <<= t[i];
            it->flip();
            *it &= *dp;
            *it <<= v[i];
            *dp |= *it;
        }
        for (int i = ans_upper; i >= ans; -- i) {
            if ((*dp)[i]) setmax(ans, i);
        }
        whole(random_shuffle, ix);
    }
    cout << ans << endl;
    return 0;
}