結果

問題 No.68 よくある棒を切る問題 (2)
ユーザー LayCurseLayCurse
提出日時 2014-11-13 02:12:59
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 258 ms / 5,000 ms
コード長 2,151 bytes
コンパイル時間 1,180 ms
コンパイル使用メモリ 160,376 KB
実行使用メモリ 10,368 KB
最終ジャッジ日時 2024-06-10 02:37:07
合計ジャッジ時間 5,746 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 258 ms
10,340 KB
testcase_01 AC 252 ms
10,240 KB
testcase_02 AC 249 ms
10,240 KB
testcase_03 AC 255 ms
10,368 KB
testcase_04 AC 247 ms
10,112 KB
testcase_05 AC 241 ms
10,368 KB
testcase_06 AC 235 ms
9,728 KB
testcase_07 AC 240 ms
10,112 KB
testcase_08 AC 247 ms
9,600 KB
testcase_09 AC 248 ms
9,856 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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[500001];
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);

  assert(1 <= N && N <= 200000);
  rep(i,N) assert(1 <= A[i] && A[i] <= 1000000000);
  assert(1 <= Q && Q <= 100000);
  rep(i,Q) assert(1 <= K[i] && K[i] <= 500000);

  hp.malloc(N);
  hp.init(N);
  rep(i,N) hp.change(i, -A[i]), use[i] = 1;
  REP(i,1,500001){
    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;
}
0