結果
問題 | No.68 よくある棒を切る問題 (2) |
ユーザー |
![]() |
提出日時 | 2014-11-13 01:49:12 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 445 ms / 5,000 ms |
コード長 | 1,981 bytes |
コンパイル時間 | 1,624 ms |
コンパイル使用メモリ | 160,356 KB |
実行使用メモリ | 14,480 KB |
最終ジャッジ日時 | 2024-12-31 10:06:17 |
合計ジャッジ時間 | 9,269 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define REP(i,a,b) for(i=a;i<b;i++)#define rep(i,n) REP(i,0,n)#define mygc(c) (c)=getchar_unlocked()#define mypc(c) putchar_unlocked(c)void reader(int *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}template <class T>struct heapEx {int *hp, *place, size; T *val;void malloc(int N){hp=(int*)std::malloc(N*sizeof(int));place=(int*)std::malloc(N*sizeof(int));val=(T*)std::malloc(N*sizeof(T));}void free(){free(hp);free(place);free(val);}void* malloc(int N, void *workMemory){hp=(int*)workMemory;workMemory=(void*)(hp+N);place=(int*)workMemory;workMemory=(void*)(place+N);val=(T*)workMemory;workMemory=(void*)(val+N);return workMemory;}void init(int N){int i;size=0;rep(i,N)place[i]=-1;}void up(int n){int m;while(n){m=(n-1)/2;if(val[hp[m]]<=val[hp[n]])break;swap(hp[m],hp[n]);swap(place[hp[m]],place[hp[n]]);n=m;}}void down(int n){int m;for(;;){m=2*n+1;if(m>=size)break;if(m+1<size&&val[hp[m]]>val[hp[m+1]])m++;if(val[hp[m]]>=val[hp[n]])break;swap(hp[m],hp[n]);swap(place[hp[m]],place[hp[n]]);n=m;}}void change(int n, T v){T f=val[n];val[n]=v;if(place[n]==-1){place[n] = size;hp[size++] = n;up(place[n]);}else{if(f < v) down(place[n]); else if(f >v) up(place[n]);}}int pop(void){int res=hp[0];place[res]=-1;size--;if(size)hp[0]=hp[size],place[hp[0]]=0,down(0);return res;}};int N, Q;int A[200000], use[200000];double res[1000001];int K[100000];int main(){int i, j, k;heapEx<double> hp;reader(&N);rep(i,N) reader(A+i);reader(&Q);rep(i,Q) reader(K+i);hp.malloc(N);hp.init(N);rep(i,N) hp.change(i, -A[i]), use[i] = 1;REP(i,1,1000001){k = hp.hp[0];res[i] = -hp.val[k];use[k]++;hp.change(k, -A[k]/(double)use[k]);}rep(i,Q) printf("%.15f\n",res[K[i]]);return 0;}