This documentation is automatically generated by online-judge-tools/verification-helper
View the Project on GitHub beet-aizu/library
// verification-helper: PROBLEM https://yukicoder.me/problems/1923 // verification-helper: ERROR 1e-6 #include<bits/stdc++.h> using namespace std; #define call_from_test #include "../../io/precision.cpp" #include "../../datastructure/rangeslide.cpp" #undef call_from_test signed main(){ cin.tie(0); ios::sync_with_stdio(0); int n,k; cin>>n>>k; string s; cin>>s; vector<int> a(n); for(int i=0;i<n;i++) a[i]=s[i]-'0'; using D = double; auto check= [&](D x)->int{ vector<D> b(n); for(int i=0;i<n;i++) b[i]=a[i]-x; vector<D> sm(n*2+1,0); for(int i=0;i<n*2;i++) sm[i+1]=sm[i]+b[i%n]; auto cmp=[](D a,D b){return a<b;}; RangeSlide<D, decltype(cmp)> rs(sm,cmp); for(int i=0;i<n;i++) rs.add_range(i,i+n-k+1); auto res=rs.build(); for(int i=0;i<n;i++) if(sm[n+i]>sm[res[i]]) return 1; return 0; }; D L=0,R=1; for(int k=0;k<20;k++){ D M=(L+R)/2; if(check(M)) L=M; else R=M; } cout<<L<<endl; return 0; }
#line 1 "test/yukicoder/1923.test.cpp" // verification-helper: PROBLEM https://yukicoder.me/problems/1923 // verification-helper: ERROR 1e-6 #include<bits/stdc++.h> using namespace std; #define call_from_test #line 1 "io/precision.cpp" #line 3 "io/precision.cpp" using namespace std; #endif //BEGIN CUT HERE struct Precision{ Precision(){ cout<<fixed<<setprecision(12); } }precision_beet; //END CUT HERE #ifndef call_from_test signed main(){ return 0; } #endif #line 1 "datastructure/rangeslide.cpp" #line 3 "datastructure/rangeslide.cpp" using namespace std; #endif //BEGIN CUT HERE template<typename T, typename F> struct RangeSlide{ vector<size_t> ls,rs; vector<T> vs; F cmp; RangeSlide(vector<T> vs,F cmp):vs(vs),cmp(cmp){} void add_range(size_t l,size_t r){ ls.emplace_back(l); rs.emplace_back(r); } vector<size_t> build(){ deque<size_t> deq; vector<size_t> res; for(size_t i=0,l=0,r=0;i<ls.size();i++){ if(r<=ls[i]){ deq.clear(); l=r=ls[i]; } while(r<rs[i]){ while(!deq.empty() and !cmp(vs[deq.back()],vs[r])) deq.pop_back(); deq.emplace_back(r++); } while(l<ls[i]){ if(deq.front()==l++) deq.pop_front(); } res.emplace_back(deq.front()); } return res; } }; //END CUT HERE #ifndef call_from_test //INSERT ABOVE HERE signed main(){ return 0; } #endif #line 10 "test/yukicoder/1923.test.cpp" #undef call_from_test signed main(){ cin.tie(0); ios::sync_with_stdio(0); int n,k; cin>>n>>k; string s; cin>>s; vector<int> a(n); for(int i=0;i<n;i++) a[i]=s[i]-'0'; using D = double; auto check= [&](D x)->int{ vector<D> b(n); for(int i=0;i<n;i++) b[i]=a[i]-x; vector<D> sm(n*2+1,0); for(int i=0;i<n*2;i++) sm[i+1]=sm[i]+b[i%n]; auto cmp=[](D a,D b){return a<b;}; RangeSlide<D, decltype(cmp)> rs(sm,cmp); for(int i=0;i<n;i++) rs.add_range(i,i+n-k+1); auto res=rs.build(); for(int i=0;i<n;i++) if(sm[n+i]>sm[res[i]]) return 1; return 0; }; D L=0,R=1; for(int k=0;k<20;k++){ D M=(L+R)/2; if(check(M)) L=M; else R=M; } cout<<L<<endl; return 0; }