結果

問題 No.576 E869120 and Rings
ユーザー fumofumofunifumofumofuni
提出日時 2022-10-16 22:16:53
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 654 ms / 1,500 ms
コード長 2,857 bytes
コンパイル時間 2,175 ms
コンパイル使用メモリ 209,800 KB
実行使用メモリ 18,652 KB
最終ジャッジ日時 2023-09-09 23:32:44
合計ジャッジ時間 16,138 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 549 ms
18,396 KB
testcase_01 AC 551 ms
18,652 KB
testcase_02 AC 654 ms
18,540 KB
testcase_03 AC 643 ms
18,500 KB
testcase_04 AC 557 ms
18,568 KB
testcase_05 AC 575 ms
18,224 KB
testcase_06 AC 573 ms
18,484 KB
testcase_07 AC 577 ms
18,480 KB
testcase_08 AC 600 ms
18,260 KB
testcase_09 AC 522 ms
18,640 KB
testcase_10 AC 505 ms
18,556 KB
testcase_11 AC 492 ms
18,340 KB
testcase_12 AC 503 ms
18,232 KB
testcase_13 AC 371 ms
16,212 KB
testcase_14 AC 517 ms
16,476 KB
testcase_15 AC 162 ms
15,996 KB
testcase_16 AC 295 ms
16,056 KB
testcase_17 AC 401 ms
16,352 KB
testcase_18 AC 509 ms
18,324 KB
testcase_19 AC 542 ms
18,256 KB
testcase_20 AC 504 ms
18,488 KB
testcase_21 AC 546 ms
18,564 KB
testcase_22 AC 1 ms
4,376 KB
testcase_23 AC 2 ms
4,380 KB
testcase_24 AC 1 ms
4,376 KB
testcase_25 AC 2 ms
4,376 KB
testcase_26 AC 1 ms
4,380 KB
testcase_27 AC 1 ms
4,376 KB
testcase_28 AC 2 ms
4,380 KB
testcase_29 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
//#pragma GCC optimize("Ofast")

#define rep(i,n) for(ll i=0;i<n;i++)
#define repl(i,l,r) for(ll i=(l);i<(r);i++)
#define per(i,n) for(ll i=(n)-1;i>=0;i--)
#define perl(i,r,l) for(ll i=r-1;i>=l;i--)
#define fi first
#define se second
#define pb push_back
#define ins insert
#define pqueue(x) priority_queue<x,vector<x>,greater<x>>
#define all(x) (x).begin(),(x).end()
#define CST(x) cout<<fixed<<setprecision(x)
#define vtpl(x,y,z) vector<tuple<x,y,z>>
#define rev(x) reverse(x);
using ll=long long;
using vl=vector<ll>;
using vvl=vector<vector<ll>>;
using pl=pair<ll,ll>;
using vpl=vector<pl>;
using vvpl=vector<vpl>;
const ll MOD=1000000007;
const ll MOD9=998244353;
const int inf=1e9+10;
const ll INF=4e18;
const ll dy[9]={1,0,-1,0,1,1,-1,-1,0};
const ll dx[9]={0,1,0,-1,1,-1,1,-1,0};
template<class T> inline bool chmin(T& a, T b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template<class T> inline bool chmax(T& a, T b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
template<class T>
class foldable_queue{
    function<T(T,T)> fx;
    using P=pair<T,T>;
    stack<P> sf,sb;
    T ex;
    void reload(){
        while(!sb.empty()){
            if(sf.empty())sf.push(make_pair(sb.top().first,sb.top().first));
            else sf.push(make_pair(sb.top().first,fx(sb.top().first,sf.top().second)));
            sb.pop();
        }
    }
    public:
    foldable_queue(function<T(T,T)> _fx,T _ex){
        fx=_fx;ex=_ex;
    }
    int size(){
        return sf.size()+sb.size();
    }
    T front(){
        if(sf.empty())reload();
        return sf.top().first;
    }
    T fold(){
        T ans=(sf.empty()? ex:sf.top().second);
        ans=fx(ans,(sb.empty()? ex:sb.top().second));
        return ans;
    }
    void push(T x){
        if(sb.empty())sb.push(make_pair(x,x));
        else sb.push(make_pair(x,fx(sb.top().second,x)));
    }
    void pop(){
        if(sf.empty())reload();
        sf.pop();
    }
};

int main(){
    CST(10);
    ll n,k;cin >> n >> k;
    string s;cin >> s;s+=s;
    ll m=s.size();
    vector<double> cum(m+1);
    rep(i,m){
        cum[i+1]=cum[i]+(s[i]-'0');
    }
    double ok=0,ng=1.00;
    rep(_,30){
        double mid=(ok+ng)/2;
        foldable_queue<double> que([](double a,double b){return max(a,b);},-inf);
        bool y=false;
        repl(i,k,n+1){
            //cout << cum[i]-mid*i << endl;
            que.push(cum[i]-mid*i);
        }
        rep(i,n){
            double now=cum[i]-mid*i;
            //cout << now << " " << que.fold() << endl;
            if(que.fold()>=now){
                y=true;break;
            }
            que.pop();
            if(i<n-1)que.push(cum[n+i+1]-mid*(n+i+1));
        }
        if(y)ok=mid;
        else ng=mid;
    }
    cout << ok << endl;
}
0