QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#383248 | #5537. Storing Eggs | kevinyang# | TL | 64ms | 187384kb | C++20 | 12.3kb | 2024-04-09 07:18:49 | 2024-04-09 07:18:49 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")
/* Macros {{{ */
/* A lot of this is from some of Benq's submissions
[https://codeforces.com/profile/Benq]
Ugly af to the eyes, but with vim fold its barable
Hopefully c++20 concepts can make all this stuff must cleaner */
/* Basics {{{ */
using ll = long long;
using ld = long double;
using str = string;
using pi = pair<int, int>;
using pll = pair<ll, ll>;
using pld = pair<ld, ld>;
#define mp make_pair
#define fi first
#define se second
#define arr array
#define ve vector
using vi = vector<int>;
using vll = vector<ll>;
using vld = vector<ld>;
using vpi = vector<pi>;
using vpll = vector<pll>;
using vpld = vector<pld>;
using vvi = vector<vi>;
using vvll = vector<vll>;
using vvld = vector<vld>;
using vvpi = vector<vpi>;
using vvpll = vector<vpll>;
using vvpld = vector<vpld>;
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define sz size()
#define rsz(a) resize(a)
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define For(i, a, b) for (int i = a; i < b; ++i)
#define Rof(i, a, b) for (int i = (b)-1; i >= (a); --i)
#define rep(a) For(_, 0, a)
#define each(a, x) for (auto &a : x)
#define reach(a, x) for (auto a = x.rbegin(); a != x.rend(); ++a)
template <typename T, typename U>
inline void cmin(T &x, U y) {
if (y < x) x = y;
}
template <typename T, typename U>
inline void cmax(T &x, U y) {
if (x < y) x = y;
}
/*}}}*/
/* IO {{{ */
/* Template Macros {{{ */
#define tcT template <class T
#define tcTU tcT, class U
#define tcTUU tcT, class... U
/*}}}*/
inline namespace Helpers { /*{{{*/
tcT, class = void > struct is_iterable : false_type {};
tcT > struct is_iterable<
T, void_t<decltype(begin(declval<T>())), decltype(end(declval<T>()))>>
: true_type {};
tcT > constexpr bool is_iterable_v = is_iterable<T>::value;
tcT, class = void > struct is_readable : false_type {};
tcT > struct is_readable<T, typename std::enable_if_t<is_same_v<
decltype(cin >> declval<T &>()), istream &>>>
: true_type {};
tcT > constexpr bool is_readable_v = is_readable<T>::value;
tcT, class = void > struct is_printable : false_type {};
tcT > struct is_printable<T, typename std::enable_if_t<is_same_v<
decltype(cout << declval<T>()), ostream &>>>
: true_type {};
tcT > constexpr bool is_printable_v = is_printable<T>::value;
} /* namespace Helpers */
/*}}}*/
inline namespace Input { /*{{{*/
tcT > constexpr bool needs_input_v = !is_readable_v<T> && is_iterable_v<T>;
tcTUU > void re(T &t, U &...u);
tcTU > void re(pair<T, U> &p); /* pairs */
/* re: read{{{ */
tcT > typename enable_if<is_readable_v<T>, void>::type re(T &x) {
cin >> x;
} /* default */
tcT > typename enable_if<needs_input_v<T>, void>::type re(
T &i); // vectors, arrays, etc...
tcTU > void re(pair<T, U> &p) { re(p.fi, p.se); } // pairs
tcT > typename enable_if<needs_input_v<T>, void>::type re(T &i) {
each(x, i) re(x);
}
tcTUU > void re(T &t, U &...u) {
re(t);
re(u...);
} /* read multiple}}} */
/* rv: resize and read vectors{{{ */
void rv(size_t) {}
tcTUU > void rv(size_t N, ve<T> &t, U &...u);
template <class... U>
void rv(size_t, size_t N2, U &...u);
tcTUU > void rv(size_t N, ve<T> &t, U &...u) {
t.rsz(N);
re(t);
rv(N, u...);
}
template <class... U>
void rv(size_t, size_t N2, U &...u) {
rv(N2, u...);
} /*}}}*/
/* dumb shortcuts to read in ints{{{ */
void decrement() {} /* subtract one from each */
tcTUU > void decrement(T &t, U &...u) {
--t;
decrement(u...);
}
#define ints(...) \
int __VA_ARGS__; \
re(__VA_ARGS__);
#define int1(...) \
ints(__VA_ARGS__); \
decrement(__VA_ARGS__); /*}}}*/
} /* namespace Input */
/*}}}*/
inline namespace ToString { /*{{{*/
tcT > constexpr bool needs_output_v = !is_printable_v<T> && is_iterable_v<T>;
/* ts: string representation to print */
tcT > typename enable_if<is_printable_v<T>, str>::type ts(T v) {
stringstream ss;
ss << fixed << setprecision(15) << v;
return ss.str();
} /* default */
tcT > str bit_vec(T t) { /* bit vector to string */
str res = "{";
For(i, 0, t.sz) res += ts(t[i]);
res += "}";
return res;
}
str ts(ve<bool> v) { return bit_vec(v); }
template <size_t SZ>
str ts(bitset<SZ> b) {
return bit_vec(b);
} /* bit vector */
tcTU > str ts(pair<T, U> p); /* pairs */
tcT > typename enable_if<needs_output_v<T>, str>::type ts(
T v); /* vectors, arrays */
tcTU > str ts(pair<T, U> p) { return "(" + ts(p.fi) + ", " + ts(p.se) + ")"; }
tcT > typename enable_if<is_iterable_v<T>, str>::type ts_sep(T v, str sep) {
/* convert container to string w/ separator sep */
bool fst = 1;
str res = "";
for (const auto &x : v) {
if (!fst) res += sep;
fst = 0;
res += ts(x);
}
return res;
}
tcT > typename enable_if<needs_output_v<T>, str>::type ts(T v) {
return "{" + ts_sep(v, ", ") + "}";
}
/* for nested DS */
template <int, class T>
typename enable_if<!needs_output_v<T>, ve<str>>::type ts_lev(const T &v) {
return {ts(v)};
}
template <int lev, class T>
typename enable_if<needs_output_v<T>, ve<str>>::type ts_lev(const T &v) {
if (lev == 0 || !v.sz) return {ts(v)};
ve<str> res;
for (const auto &t : v) {
if (res.sz) res.back() += ",";
ve<str> tmp = ts_lev<lev - 1>(t);
res.insert(end(res), all(tmp));
}
For(i, 0, res.sz) {
str bef = " ";
if (i == 0) bef = "{";
res[i] = bef + res[i];
}
res.back() += "}";
return res;
}
} /* namespace ToString */
/*}}}*/
inline namespace Output { /*{{{*/
template <class T>
void pr_sep(ostream &os, str, const T &t) {
os << ts(t);
}
template <class T, class... U>
void pr_sep(ostream &os, str sep, const T &t, const U &...u) {
pr_sep(os, sep, t);
os << sep;
pr_sep(os, sep, u...);
}
/* print w/ no spaces */
template <class... T>
void pr(const T &...t) {
pr_sep(cout, "", t...);
}
/* print w/ spaces, end with newline */
void ps() { cout << "\n"; }
template <class... T>
void ps(const T &...t) {
pr_sep(cout, " ", t...);
ps();
}
/* debug to cerr */
template <class... T>
void dbg_out(const T &...t) {
pr_sep(cerr, " | ", t...);
cerr << endl;
}
void loc_info(int line, str names) {
cerr << "Line(" << line << ") -> [" << names << "]: ";
}
template <int lev, class T>
void dbgl_out(const T &t) {
cerr << "\n\n" << ts_sep(ts_lev<lev>(t), "\n") << "\n" << endl;
}
} /* namespace Output */
/*}}}}}}}}}*/
int n, k;
ve<str> carton;
int dp[101][301][101][64] = { };
int iterative_transition() {
for(int ck=0; ck<=k; ++ck) for(int x=0; x<n; ++x) for(int s=0; s<(1<<6); ++s)
dp[n][ck][x][s] = -1*(ck>0) + (INT_MAX/4)*(ck==0);
for(int i=n-1; i>=0; --i) {
for(int ck=0; ck<=k; ++ck) {
for(int x=0; x<n; ++x) {
for(int s=0; s<64; ++s) {
int s0 = s>>4, s1 = (s&15)>>2, s2 = s&3;
int &res = dp[i][ck][x][s];
res = max(1, dp[i+1][ck][x][s]);
for(int c0=0; c0<=carton[0][i]; ++c0) {
for(int c1=0; c1<=carton[1][i]; ++c1) {
for(int c2=0; c2<=carton[2][i]; ++c2) {
if(c0 == 0 && c1 == 0 && c2 == 0 || ((ck-c0-c1-c2) < 0)) continue;
int opt = INT_MAX/4;
// if(c1 & (c0 | c2)) cmin(opt, 1);
if(c1 & (c0 | c2)) continue;
if(c0 & c2) cmin(opt, 4);
if(s0 != 0) {
int dx = (i - (x-1+(s0&1)));
int dx2 = dx*dx;
if(c0) cmin(opt, dx2);
else if(c1) cmin(opt, dx2 + 1);
else if(c2) cmin(opt, dx2 + 4);
}
if(s1 != 0) {
int dx = (i - (x-1+(s1&1)));
int dx2 = dx*dx;
if(c1) cmin(opt, dx2);
else if(c0|c2) cmin(opt, dx2+1);
}
if(s2 != 0) {
int dx = (i - (x-1+(s2&1)));
int dx2 = dx*dx;
if(c2) cmin(opt, dx2);
else if(c1) cmin(opt, dx2 + 1);
else if(c0) cmin(opt, dx2 + 4);
}
if(opt < res) continue;
int ns;
if(x == i-1) ns = ((((s0<<1)|c0)&3)<<4) | ((((s1<<1)|c1)&3)<<2) | ((((s2<<1)|c2)&3));
else ns = (c0<<4) | (c1<<2) | c2;
cmin(opt, dp[i+1][ck-c0-c1-c2][i][ns]);
cmax(res, opt);
}
}
}
}
}
}
}
return dp[0][k][0][0];
}
inline int transition(int i, int k, int x, int s0, int s1, int s2) {
if(k <= 0) return INT_MAX/4;
if(i >= n) return -1;
int s = (s0<<4) | (s1<<2) | s2;
if(dp[i][k][x][s] != 0) return dp[i][k][x][s];
int &res = dp[i][k][x][s];
res = max(1, transition(i+1, k, x, s0, s1, s2));
for(int c0=0; c0<=carton[0][i]; ++c0) {
for(int c1=0; c1<=carton[1][i]; ++c1) {
for(int c2=0; c2<=carton[2][i]; ++c2) {
if(c0 == 0 && c1 == 0 && c2 == 0) continue;
int opt = INT_MAX/4;
// if(c1 & (c0 | c2)) cmin(opt, 1);
if(c1 & (c0 | c2)) continue;
if(c0 & c2) cmin(opt, 4);
if(s0 != 0) {
int dx = (i - (x - __builtin_ctz(s0)));
int dx2 = dx*dx;
if(c0) cmin(opt, dx2);
else if(c1) cmin(opt, dx2 + 1);
else if(c2) cmin(opt, dx2 + 4);
}
if(s1 != 0) {
int dx = (i - (x - __builtin_ctz(s1)));
int dx2 = dx*dx;
if(c1) cmin(opt, dx2);
else if(c0|c2) cmin(opt, dx2+1);
}
if(s2 != 0) {
int dx = (i - (x - __builtin_ctz(s2)));
int dx2 = dx*dx;
if(c2) cmin(opt, dx2);
else if(c1) cmin(opt, dx2 + 1);
else if(c0) cmin(opt, dx2 + 4);
}
if(opt < res) continue;
if(x == i-1) cmin(opt, transition(i+1, k-c0-c1-c2, i, ((s0<<1)|c0)&3, ((s1<<1)|c1)&3, ((s2<<1)|c2)&3));
else cmin(opt, transition(i+1, k-c0-c1-c2, i, c0, c1, c2));
cmax(res, opt);
}
}
}
return res;
}
void solve() {
re(n, k), rv(3, carton);
for(int r=0; r<3; ++r) for(int c=0; c<n; ++c) carton[r][c] = (carton[r][c] == '.');
int tot_dots = 0;
for(int r=0; r<3; ++r) for(int c=0; c<n; ++c) tot_dots += carton[r][c];
if(tot_dots < k) { ps(-1); return; }
// ps(sqrt(transition(0, k, 0, 0, 0, 0)));
ps(sqrt(iterative_transition()));
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(NULL);
/* cout << fixed << setprecision(6); */
int t = 1;
// cin >> t;
for (int i = 0; i < t; i++) solve();
return 0;
// you should actually read the stuff at the bottom
}
/* stuff you should look for
* int overflow, array bounds
* special cases (n=1?)
* do smth instead of nothing and stay organized
* WRITE STUFF DOWN
* DON'T GET STUCK ON ONE APPROACH
*/
詳細信息
Test #1:
score: 100
Accepted
time: 0ms
memory: 14232kb
input:
5 2 #.... ..... ....#
output:
4.472135954999580
result:
ok found '4.4721360', expected '4.4721360', error '0.0000000'
Test #2:
score: 0
Accepted
time: 2ms
memory: 16240kb
input:
5 6 ##.## ##### .....
output:
1.000000000000000
result:
ok found '1.0000000', expected '1.0000000', error '0.0000000'
Test #3:
score: 0
Accepted
time: 0ms
memory: 10024kb
input:
3 4 ..# ... ...
output:
1.414213562373095
result:
ok found '1.4142136', expected '1.4142140', error '0.0000003'
Test #4:
score: 0
Accepted
time: 0ms
memory: 3616kb
input:
2 6 .. .# ..
output:
-1
result:
ok found '-1.0000000', expected '-1.0000000', error '-0.0000000'
Test #5:
score: 0
Accepted
time: 0ms
memory: 5888kb
input:
1 2 . . .
output:
2.000000000000000
result:
ok found '2.0000000', expected '2.0000000', error '0.0000000'
Test #6:
score: 0
Accepted
time: 32ms
memory: 187384kb
input:
100 2 .................................................................................................... .................................................................................................... ...............................................................................................
output:
99.020199959402220
result:
ok found '99.0202000', expected '99.0202000', error '0.0000000'
Test #7:
score: 0
Accepted
time: 64ms
memory: 177960kb
input:
100 3 .................................................................................................... .................................................................................................... ...............................................................................................
output:
49.040799340956916
result:
ok found '49.0407993', expected '49.0407990', error '0.0000000'
Test #8:
score: -100
Time Limit Exceeded
input:
100 100 .................................................................................................... .................................................................................................... .............................................................................................