結果

問題 No.576 E869120 and Rings
ユーザー parukiparuki
提出日時 2017-10-16 10:36:53
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
MLE  
実行時間 -
コード長 2,920 bytes
コンパイル時間 1,830 ms
コンパイル使用メモリ 172,884 KB
実行使用メモリ 227,420 KB
最終ジャッジ日時 2024-04-28 23:22:32
合計ジャッジ時間 7,401 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 MLE -
testcase_01 TLE -
testcase_02 TLE -
testcase_03 TLE -
testcase_04 TLE -
testcase_05 TLE -
testcase_06 TLE -
testcase_07 TLE -
testcase_08 TLE -
testcase_09 TLE -
testcase_10 TLE -
testcase_11 TLE -
testcase_12 TLE -
testcase_13 TLE -
testcase_14 TLE -
testcase_15 TLE -
testcase_16 TLE -
testcase_17 TLE -
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
testcase_21 TLE -
testcase_22 MLE -
testcase_23 MLE -
testcase_24 MLE -
testcase_25 MLE -
testcase_26 MLE -
testcase_27 MLE -
testcase_28 MLE -
testcase_29 MLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i))
#define rep(i,j) FOR(i,0,j)
#define each(x,y) for(auto &(x):(y))
#define mp make_pair
#define mt make_tuple
#define all(x) (x).begin(),(x).end()
#define debug(x) cout<<#x<<": "<<(x)<<endl
#define smax(x,y) (x)=max((x),(y))
#define smin(x,y) (x)=min((x),(y))
#define MEM(x,y) memset((x),(y),sizeof (x))
#define sz(x) (int)(x).size()
#define pb push_back
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<ll> vll;

vector<int> a;
template<class Val, class Cmp>
class SparseTableRMQ {
public:
    SparseTableRMQ() {
    }
    SparseTableRMQ(Val *ar, int n) : n(n), data(n) {
        for (int i = 0; i < n; ++i) data[i] = ar[i];
        table = buildRMQ();
    }

    int queryPos(int l, int r) {
        int k = lgInt(r - l - 1);
        int lid = table[l + n*k], rid = table[r + n*k - (1 << k)];
        return (cmp(data[lid], data[rid])) ? lid : rid;
    }

    Val queryValue(int l, int r) {
        return data[queryPos(l, r)];
    }
private:
    int n;
    vector<Val> data;
    vector<int> table;
    Cmp cmp;
    int lgInt(int v) {
        static const unsigned b[] = { 0x2, 0xC, 0xF0, 0xFF00, 0xFFFF0000 };
        static const unsigned S[] = { 1, 2, 4, 8, 16 };
        int res = 0;
        for (int i = 4; i >= 0; i--) {
            if (v & b[i]) {
                v >>= S[i];
                res |= S[i];
            }
        }
        return res;
    }

    vector<int> buildRMQ() {
        int logn = 1;
        for (int k = 1; k < n; k *= 2) ++logn;
        // vector<int> a(n*logn);
        for (int i = 0; i < n; ++i)a[i] = i;
        int cur = 0;
        for (int k = 1; k < n; k *= 2) {
            copy(a.begin() + cur, a.begin() + cur + n, a.begin() + cur + n);
            cur += n;
            for (int i = 0; i < n - k; ++i) {
                int idl = a[cur + i], idr = a[cur + i + k];
                if (cmp(data[idl], data[idr])) {
                    a[cur + i] = idl;
                } else {
                    a[cur + i] = idr;
                }
            }
        }
        return a;
    }
};

int N, K, color[500005];

double dat[500005*2];

int check(double x) {
    rep(i, N) {
        dat[i] = dat[i + N] = color[i] - x;
    }
    rep(i, 2 * N)dat[i + 1] += dat[i];
    SparseTableRMQ<double, greater<double> > rmq(dat, N * 2);

    double pre = 0;
    rep(i, N) {
        double ma = rmq.queryValue(i + K - 1, i + N);
        if (ma - pre >= 0)return 1;
        pre = dat[i];
    }
    return 0;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    string S;
    cin >> N >> K >> S;
    rep(i, N)color[i] = S[i] - '0';
    a = vector<int>(500000*50);
    double l = 0, r = 1, x;
    rep(ITER, 20) {
        x = (l + r)*0.5;
        (check(x) ? l : r) = x;
    }
    cout << fixed << setprecision(20);
    cout << r << endl;
}
0