結果

問題 No.2955 Pizza Delivery Plan
ユーザー GOTKAKOGOTKAKO
提出日時 2024-11-08 22:07:32
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 125 ms / 2,000 ms
コード長 2,638 bytes
コンパイル時間 2,685 ms
コンパイル使用メモリ 212,080 KB
実行使用メモリ 38,144 KB
最終ジャッジ日時 2024-11-08 22:07:40
合計ジャッジ時間 5,490 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 39 ms
16,512 KB
testcase_09 AC 39 ms
16,768 KB
testcase_10 AC 45 ms
16,768 KB
testcase_11 AC 50 ms
16,768 KB
testcase_12 AC 58 ms
20,224 KB
testcase_13 AC 63 ms
20,224 KB
testcase_14 AC 71 ms
23,808 KB
testcase_15 AC 75 ms
23,936 KB
testcase_16 AC 82 ms
27,392 KB
testcase_17 AC 86 ms
27,392 KB
testcase_18 AC 95 ms
30,848 KB
testcase_19 AC 100 ms
30,976 KB
testcase_20 AC 107 ms
34,560 KB
testcase_21 AC 114 ms
34,560 KB
testcase_22 AC 121 ms
38,016 KB
testcase_23 AC 121 ms
38,144 KB
testcase_24 AC 125 ms
38,016 KB
testcase_25 AC 58 ms
20,224 KB
testcase_26 AC 50 ms
16,512 KB
testcase_27 AC 45 ms
16,640 KB
testcase_28 AC 95 ms
30,848 KB
testcase_29 AC 100 ms
30,976 KB
testcase_30 AC 124 ms
38,144 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

template<typename T>
struct Point{
    public:
    T x,y;
    Point() : x(0),y(0) {}
    Point(T a,T b) : x(a),y(b) {}
    Point &operator=(const Point &b) = default;
    Point operator+(const Point &b){return Point(x+b.x,y+b.y);}
    Point operator-(const Point &b){return Point(x-b.x,y-b.y);}
    Point operator+=(const Point &b){return *this=*this+b;}
    Point operator-=(const Point &b){return *this=*this-b;}
    bool operator==(const Point &b){return x==b.x && y==b.y;}
    bool operator!=(const Point &b){return x!=b.x || y!=b.y;}
    bool operator<(const Point &b){
        if(x == b.x) return y < b.y;
        else return x < b.x; 
    }
    bool operator<=(const Point &b){return (*this) < b || (*this) == b;}
    bool operator>(const Point &b){return !((*this) <= b);}
    bool operator>=(const Point &b){return !((*this) < b);}

    friend T inner(const Point a,const Point b){return a.x*b.x+a.y*b.y;}
    friend T cross(const Point a,const Point b){return a.x*b.y-a.y*b.x;}
    friend T twodist(const Point a,const Point b){return sqrtl((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
};
template<typename T>
istream &operator>>(istream &is,Point<T> &a){
    is >> a.x >> a.y;
    return is;
}
template<typename T>
ostream &operator<<(ostream &os,Point<T> &a){
    os << a.x << " " << a.y;
    return os;
}

int main(){
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int N,K; cin >> N >> K;
    vector<Point<double>> P(N);
    for(auto &p : P) cin >> p;

    int n2 = 1<<N;
    vector<vector<vector<double>>> dp(n2,vector<vector<double>>(N,vector<double>(K,1e100)));
    for(int i=0; i<N; i++) dp.at(1<<i).at(i).at(K-1) = twodist({0.0,0.0},P.at(i));
    for(int i=1; i<n2; i++){
        for(int k=0; k<N; k++){
            if(!(i&(1<<k))) continue;
            int left = i-(1<<k);
            for(int l=0; l<N; l++){
                if(!(i&(1<<l))) continue;

                double mind = dp.at(left).at(l).at(0);    
                for(int m=0; m<K-1; m++){
                    dp.at(i).at(k).at(m) = min(dp.at(i).at(k).at(m),dp.at(left).at(l).at(m+1)+twodist(P.at(k),P.at(l)));
                    mind = min(mind,dp.at(left).at(l).at(m+1));
                }
                dp.at(i).at(k).at(K-1) = min(dp.at(i).at(k).at(K-1),mind+twodist({0,0},P.at(l))+twodist({0,0},P.at(k)));
            }
        }
    }

    double answer = 1e100;
    for(int i=0; i<N; i++){
        double nowd = twodist({0,0},P.at(i));
        for(int k=0; k<K; k++) answer = min(answer,dp.at(n2-1).at(i).at(k)+nowd);
    }
    cout << fixed << setprecision(20) << answer << endl;
}
0