結果

問題 No.1413 Dynamic Sushi
ユーザー evimaevima
提出日時 2021-02-26 16:43:47
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 544 ms / 4,000 ms
コード長 5,452 bytes
コンパイル時間 5,035 ms
コンパイル使用メモリ 268,692 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-10 13:06:14
合計ジャッジ時間 15,732 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 472 ms
6,940 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 478 ms
6,940 KB
testcase_05 AC 479 ms
6,940 KB
testcase_06 AC 544 ms
6,940 KB
testcase_07 AC 502 ms
6,940 KB
testcase_08 AC 486 ms
6,940 KB
testcase_09 AC 518 ms
6,944 KB
testcase_10 AC 493 ms
6,944 KB
testcase_11 AC 510 ms
6,944 KB
testcase_12 AC 516 ms
6,940 KB
testcase_13 AC 526 ms
6,940 KB
testcase_14 AC 538 ms
6,940 KB
testcase_15 AC 522 ms
6,940 KB
testcase_16 AC 521 ms
6,940 KB
testcase_17 AC 526 ms
6,944 KB
testcase_18 AC 531 ms
6,944 KB
testcase_19 AC 516 ms
6,944 KB
testcase_20 AC 229 ms
6,940 KB
testcase_21 AC 4 ms
6,944 KB
testcase_22 AC 539 ms
6,940 KB
testcase_23 AC 534 ms
6,940 KB
testcase_24 AC 510 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

// Enjoy your stay. Code by evima

#include <atcoder/convolution>
#include <atcoder/dsu>
#include <atcoder/fenwicktree>
#include <atcoder/lazysegtree>
#include <atcoder/math>
#include <atcoder/maxflow>
#include <atcoder/mincostflow>
#include <atcoder/modint>
#include <atcoder/scc>
#include <atcoder/segtree>
#include <atcoder/string>
#include <atcoder/twosat>
using namespace atcoder;
using mint = modint998244353;// modint1000000007;// modint;

#include <bits/stdc++.h>
using namespace std;

using ld = long double;
using ll = long long;
using lll = __int128;
using vl = vector<ll>;
using vs = vector<string>;
using LOOPVAR_TYPE = ll;

#define all(x) (x).begin(), (x).end()
#define sq(x) ((x) * (x))
#define sz(x) ll((x).size())
#define GET_MACRO(_1, _2, _3, NAME, ...) NAME
#define rep1(i, n) rep2(i, 0, n)
#define rep2(i, a, b) for(LOOPVAR_TYPE i = LOOPVAR_TYPE(a); i < LOOPVAR_TYPE(b); i++)
#define rep(...) GET_MACRO(__VA_ARGS__, rep2, rep1)(__VA_ARGS__)

#define eb emplace_back
#define fir first
#define sec second
#define mp make_pair
#define mt make_tuple

const ld EPS = 1e-9;
const ld PI = 3.14159265358979323846L;
const ll INF = 1070000000LL;
const ll MOD = 998244353LL;// 1000000007LL;

void fast_io() { cin.tie(NULL); ios_base::sync_with_stdio(false); }
ll lin() { ll x; cin >> x; return x; }
string input() { string s; cin >> s; return s; }
vl vin(int n) { vector<ll> v(n); rep(i, n) cin >> v[i]; return v; }
template<typename T> bool chmin(T& a, const T& b) { return (b < a) ? (a = b, true) : false; }
template<typename T> bool chmax(T& a, const T& b) { return (a < b) ? (a = b, true) : false; }
template<typename T> map<T, ll> freq(const vector<T>& v){ map<T, ll> ret; for(T x: v) ++ret[x]; return ret; }
template<typename T> vector<T> reversed(vector<T> v) { reverse(all(v)); return v; }
template<typename T> vector<T> sorted(vector<T> v) { sort(all(v)); return v; }
template<typename T> vector<T> sub(const vector<T>& v, int from, int to) { vector<T> ret; copy(&v[from], &v[to], back_inserter(ret)); return ret; }
template<typename T> string str(const T& x) { stringstream ss; ss << x; return ss.str(); }
template<typename Tuple, size_t... I> array<int, sizeof...(I)> str_tuple_impl(stringstream& ss, Tuple&& t, index_sequence<I...>){ return {{ (void( ss << str(get<I>(t)) << " " ), 0)... }}; }
template<typename Tuple> string str_tuple(Tuple&& t){ stringstream ss; str_tuple_impl(ss, forward<Tuple>(t), make_index_sequence<tuple_size<decay_t<Tuple>>::value>{}); string s = ss.str(); if(s.size() > 0) s.pop_back(); return s; }
template<typename T, typename U> string str(const pair<T, U>& p) { return str_tuple(p); }
template<typename... T> string str(const tuple<T...>& t) { return str_tuple(t); }
template<typename T> string str(const vector<T>& v) { stringstream ss; rep(i, sz(v)) ss << v[i] << (i < sz(v) - 1 ? " " : ""); return ss.str(); }
template<typename T> void print1(T&& x, const string& end) { cout << str(x) << end; }
void print() { print1("", "\n"); }
template<typename T, typename... U> void print(T&& head, U&&... tail) { print1(head, " "); print(forward<U>(tail)...); }
template<typename T> void eprint1(T&& x, const string& end) { cerr << str(x) << end; }
void eprint() { eprint1("", "\n"); }
template<typename T, typename... U> void eprint(T&& head, U&&... tail) { eprint1(head, " "); eprint(forward<U>(tail)...); }
template<typename... T> void quit(T&&... x) { print(forward<T>(x)...); exit(0); }
template<typename T, typename U> void yn(bool cnd, T&& yes = "Yes", U&& no = "No") { if(cnd) print(yes); else print(no); }
void zip(vl& v) { vl w = v; sort(all(w)); int n = unique(all(w)) - w.begin(); for(ll& x: v) x = lower_bound(&w[0], &w[n], x) - &w[0]; }
vl zipped(vl v) { zip(v); return v; }
map<ll, ll> zipmap(vl v) { map<ll, ll> ret; sort(all(v)); v.erase(unique(all(v)), v.end()); rep(i, sz(v)) ret[v[i]] = i; return ret; }

const int MN = 12, MX = 100;

int N;
double W, X[MN + 1], Y[MN + 1], R[MN + 1], V[MN + 1], A[MN + 1];
double dp[1 << MN][MN + 1];

double earliest(int from, int to, double t){
    double lo = 0, hi = 6 * MX / W;
    double sa = (V[from] * t + A[from]) * PI / 180;
    double sx = X[from] + R[from] * cos(sa);
    double sy = Y[from] + R[from] * sin(sa);
    rep(_, 60){
        double mi = (lo + hi) / 2;
        double a = (V[to] * (t + mi) + A[to]) * PI / 180;
        double tx = X[to] + R[to] * cos(a);
        double ty = Y[to] + R[to] * sin(a);
        double dx = tx - sx;
        double dy = ty - sy;
        double d = W * mi;
        (dx * dx + dy * dy < d * d ? hi : lo) = mi;
    }
    return t + lo;
}

void solveOne(){
    cin >> N >> W;
    rep(i, N){
        cin >> X[i] >> Y[i] >> R[i] >> V[i] >> A[i];
    }
    rep(mask, 1 << N){
        rep(i, N + 1){
            dp[mask][i] = mask == 0 && i == N ? 0 : INF;
        }
    }
    rep(mask, (1 << N) - 1){
        rep(i, N + 1){
            if(mask == 0 && i == N || mask >> i & 1){
                rep(j, N){
                    if(!(mask >> j & 1)){
                        chmin(dp[mask | 1 << j][j], earliest(i, j, dp[mask][i]));
                    }
                }
            }
        }
    }
    cout << *min_element(dp[(1 << N) - 1], dp[(1 << N) - 1] + N + 1) << endl;
}

int main(){
    fast_io();
    cout << setprecision(15);
    
    int num_tc = 1;
    // cin >> num_tc;
    
    
    
    rep(tc, 1, num_tc + 1){
        // cout << "Case #" << tc << ": " ;// << endl;
        solveOne();
    }
}
0