結果

問題 No.2955 Pizza Delivery Plan
ユーザー ATMATM
提出日時 2024-11-08 23:00:30
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 234 ms / 2,000 ms
コード長 2,950 bytes
コンパイル時間 1,831 ms
コンパイル使用メモリ 177,136 KB
実行使用メモリ 40,448 KB
最終ジャッジ日時 2024-11-08 23:00:42
合計ジャッジ時間 6,303 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 21 ms
8,960 KB
testcase_09 AC 22 ms
8,960 KB
testcase_10 AC 45 ms
11,264 KB
testcase_11 AC 69 ms
13,824 KB
testcase_12 AC 94 ms
16,128 KB
testcase_13 AC 118 ms
18,816 KB
testcase_14 AC 141 ms
20,992 KB
testcase_15 AC 165 ms
23,424 KB
testcase_16 AC 173 ms
25,856 KB
testcase_17 AC 184 ms
28,416 KB
testcase_18 AC 190 ms
30,720 KB
testcase_19 AC 201 ms
33,280 KB
testcase_20 AC 203 ms
35,456 KB
testcase_21 AC 208 ms
38,144 KB
testcase_22 AC 234 ms
40,320 KB
testcase_23 AC 209 ms
40,448 KB
testcase_24 AC 211 ms
40,448 KB
testcase_25 AC 92 ms
16,128 KB
testcase_26 AC 71 ms
13,952 KB
testcase_27 AC 44 ms
11,392 KB
testcase_28 AC 197 ms
30,720 KB
testcase_29 AC 201 ms
33,152 KB
testcase_30 AC 211 ms
40,448 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#define CPP_STR(x) CPP_STR_I(x)
#define CPP_CAT(x, y) CPP_CAT_I(x, y)
#define CPP_STR_I(args...) #args
#define CPP_CAT_I(x, y) x##y

#define ASSERT(expr...) assert((expr))

using i8 = int8_t;
using u8 = uint8_t;
using i16 = int16_t;
using u16 = uint16_t;
using i32 = int32_t;
using u32 = uint32_t;
using i64 = int64_t;
using u64 = uint64_t;

using f32 = float;
using f64 = double;
// }}}

constexpr i64 INF = 1'010'000'000'000'000'017LL;

constexpr i64 MOD = 998244353LL;

constexpr f64 EPS = 1e-12;

constexpr f64 PI = 3.14159265358979323846;

#define M5 100007
#define M9 1000000000

#define F first
#define S second

// util {{{
#define FOR(i, start, end) for (i64 i = (start), CPP_CAT(i, xxxx_end) = (end); i < CPP_CAT(i, xxxx_end); ++i)
#define REP(i, n) FOR(i, 0, n)

#define all(x) (x).begin(), (x).end()
#define ll long long int
#define VI vector<ll>
#define VVI vector<VI>

#define ISD true
#define debug(x) \
    if (ISD)     \
    cout << #x << ": " << x << endl

template <typename T, typename U, typename Comp = less<>>
bool chmax(T &xmax, const U &x, Comp comp = {})
{
    if (comp(xmax, x))
    {
        xmax = x;
        return true;
    }
    return false;
}

template <typename T, typename U, typename Comp = less<>>
bool chmin(T &xmin, const U &x, Comp comp = {})
{
    if (comp(x, xmin))
    {
        xmin = x;
        return true;
    }
    return false;
}

int main()
{
    ll N, K;
    cin >> N >> K;
    vector<ll> X(N + 1), Y(N + 1);
    REP(i, N)
    {
        cin >> X[i + 1] >> Y[i + 1];
    }

    auto dis = [&](ll x1, ll y1, ll x2, ll y2)
    {
        return sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
    };
    auto calc = [&](int i, int j, bool is0)
    {
        if (is0)
        {
            return dis(0, 0, X[i], Y[i]) + dis(0, 0, X[j], Y[j]);
        }
        else
        {
            return dis(X[i], Y[i], X[j], Y[j]);
        }
    };
    double IN = 100000000000000;
    vector<vector<vector<double>>> dp((1 << N), vector<vector<double>>(K + 1, vector<double>(N + 1, IN)));
    dp[0][K][0] = 0;
    REP(bit, (1 << N))
    {
        REP(k, K + 1)
        {
            REP(i, N + 1)
            {
                if (dp[bit][k][i] == IN)
                {
                    continue;
                }
                if (k > 0)
                {
                    REP(j, N)
                    {
                        chmin(dp[bit | (1 << j)][k - 1][j + 1], dp[bit][k][i] + calc(i, j + 1, 0));
                    }
                }
                REP(j, N)
                {
                    chmin(dp[bit | (1 << j)][K - 1][j + 1], dp[bit][k][i] + calc(i, j + 1, 1));
                }
            }
        }
    }
    double ans = IN;
    REP(k, K + 1)
    {
        REP(i, N + 1)
        {
            chmin(ans, dp[((1 << N) - 1)][k][i] + calc(0, i, 0));
        }
    }
    cout << fixed << setprecision(10) << ans << endl;
}
0