Submission #309936

#	Time	Username	Problem	Language	Result	Execution time	Memory
309936		phathnv	새로운 문제 (COCI19_akvizna)	C++11	0 / 130	1598 ms	40828 KiB

Main

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

#define X first
#define Y second
#define mp make_pair

typedef long long ll;
typedef pair<int, int> ii;

const int N = 3001;

struct line{
    double a, b;
    int cnt;
    line(double _a = 0, double _b = 0, int _cnt = 0){
        a = _a;
        b = _b;
        cnt = _cnt;
    }
};

struct convexHull{
    vector <line> d;
    int ptr;
    void reset(){
        d.clear();
        ptr = 0;
    }
    double calcY(const line &l, ll x){
        return l.a * x + l.b;
    }
    bool bad(const line &l1, const line &l2, const line &l3){
        return (l2.a - l3.a) * (l2.b - l1.b) >= (l1.a - l2.a) * (l3.b - l2.b);
    }
    void add(const line &newLine){
        while (d.size() > 1){
            if (bad(d[d.size() - 2], d[d.size() - 1], newLine))
                d.pop_back();
            else
                break;
        }
        d.push_back(newLine);
    }
    double getMax(double x){
        assert(!d.empty());
        ptr = min(ptr, (int) d.size() - 1);
        while (ptr < (int) d.size() - 1)
            if (calcY(d[ptr], x) <= calcY(d[ptr + 1], x))
                ptr++;
            else
                break;
        return calcY(d[ptr], x);
    }
    int getCnt(double x){
        assert(!d.empty());
        ptr = min(ptr, (int) d.size() - 1);
        while (ptr < (int) d.size() - 1)
            if (calcY(d[ptr], x) <= calcY(d[ptr + 1], x))
                ptr++;
            else
                break;
        return d[ptr].cnt;
    }
} cvh;

int n, k;
double dp[N][N];

void readInput(){
    cin >> n >> k;
}

void solve(){
    dp[1][0] = 0;
    for(int j = 1; j <= n; j++)
         dp[1][j] = 1;
    for(int i = 2; i <= k; i++){
        cvh.reset();
        for(int j = 1; j <= n; j++){
            cvh.add(line(dp[i - 1][j - 1], -(j - 1)));
            dp[i][j] = cvh.getMax(j) / j + 1;
        }
    }
    cout << dp[k][n];
}

int main(){
    readInput();
    solve();
    return 0;
}

• Subtask 1: 0 / 5
• Subtask 2: 0 / 5
• Subtask 3: 0 / 5
• Subtask 4: 0 / 5
• Subtask 5: 0 / 5
• Subtask 6: 0 / 5
• Subtask 7: 0 / 5
• Subtask 8: 0 / 5
• Subtask 9: 0 / 5
• Subtask 10: 0 / 5
• Subtask 11: 0 / 5
• Subtask 12: 0 / 5
• Subtask 13: 0 / 5
• Subtask 14: 0 / 5
• Subtask 15: 0 / 5
• Subtask 16: 0 / 5
• Subtask 17: 0 / 5
• Subtask 18: 0 / 5
• Subtask 19: 0 / 5
• Subtask 20: 0 / 5
• Subtask 21: 0 / 5
• Subtask 22: 0 / 5
• Subtask 23: 0 / 5
• Subtask 24: 0 / 5
• Subtask 25: 0 / 5
• Subtask 26: 0 / 5