// heap-ki
#include <stdio.h>
#include <stdlib.h>
#define height 4
#define MAX (1<<height) //ビットシフト演算 2^height と同じ
int t[MAX+1]; //配列外アクセス防止のためのダミーで+1
int sz = 0;
void swap(int *x, int *y){
int tmp = *x;
*x = *y;
*y = tmp;
}
void initTree(int n){
int i;
for(i=0;i<MAX;i++){
t[i] = -1;
}
}
void printA(){
int i;
for(i
=1;i
<=sz
;i
++) printf("%d ",t
[i
]); }
void printT(int i){
int x = i;
while(x/2!=0){
x/=2;
}
}
int goP(int i){
if(i/2 == 0) return 0;
else return i/2;
}
int goL(int i){
if(2*i >= MAX) return 0;
else return 2*i;
}
int goR(int i){
if(2*i+1 >= MAX) return 0;
else return 2*i+1;
}
void preOrder(int i){
if(t[i] == -1) return;
printT(i);
preOrder(goL(i));
preOrder(goR(i));
}
void inOrder(int i){
if(t[i] == -1) return;
inOrder(goL(i));
printT(i);
inOrder(goR(i));
}
void postOrder(int i){
if(t[i] == -1) return;
postOrder(goL(i));
postOrder(goR(i));
printT(i);
}
void insBT(int x){
int k,i = 1;
for(k=0;k<height;k++){
if(t[i]==-1){
t[i] = x;
sz++;
return;
}
if(x < t[i]) i = goL(i);
else i = goR(i);
}
printf("Error : too high -> %d\n",x
); }
//先頭の要素を取り出す
//ダウンヒープ
int popHeap(){
// 取り出して
// 末尾の要素を根に持ってきて
// 子のうちの大きい方と比較
// 子の方が大きかったら交換
// をHeapが完成するまで繰り返す
// 事前学習資料のP.56~60
int a,i,j,k=1;
a=t[1];
t[1]=-1;
swap(&t[1],&t[sz]);
sz--;
while(1){
i=k*2;
j=k*2+1;
if(t[i]>t[j] && t[i]>t[k]){
swap(&t[i],&t[k]);
k=i;
}
else if(t[i]<t[j] && t[j]>t[k]){
swap(&t[j],&t[k]);
k=i;
}
else if(t[k]>t[i] && t[k]>t[j]){
break;
}
else if(i>sz || j>sz){
break;
}
}
return a;
}
//末尾に要素を追加する
//アップヒープ
void pushHeap(int x){
// 末尾に追加
// 親と比較して、子の方が大きかったら交換
// をHeapが完成するまで繰り返す
int i;
sz++;
i = sz;
t[i] = x;
while(i>1&&t[i]>t[i/2]){
swap(&t[i],&t[i/2]);
i=i/2;
}
}
int main(void){
int i,n;
// 木の初期化
initTree(n);
for(i=0;i<n;i++){
}
sz = n;
// 中間順で表示
inOrder(1);
// pop
int a=popHeap();
inOrder(1);
a = popHeap();
inOrder(1);
pushHeap(100);
inOrder(1);
pushHeap(30);
inOrder(1);
return 0;
}
Ly8gaGVhcC1raQogCiNpbmNsdWRlIDxzdGRpby5oPgojaW5jbHVkZSA8c3RkbGliLmg+CiAKI2RlZmluZSBoZWlnaHQgNAojZGVmaW5lIE1BWCAoMTw8aGVpZ2h0KSAgLy/jg5Pjg4Pjg4jjgrfjg5Xjg4jmvJTnrpcgMl5oZWlnaHQg44Go5ZCM44GYCiAKaW50IHRbTUFYKzFdOyAvL+mFjeWIl+WkluOCouOCr+OCu+OCuemYsuatouOBruOBn+OCgeOBruODgOODn+ODvOOBp++8i++8kQppbnQgc3ogPSAwOwogCnZvaWQgc3dhcChpbnQgKngsIGludCAqeSl7CiAgICBpbnQgdG1wID0gKng7CiAgICAqeCA9ICp5OwogICAgKnkgPSB0bXA7Cn0KIAp2b2lkIGluaXRUcmVlKGludCBuKXsKICAgIGludCBpOwogICAgZm9yKGk9MDtpPE1BWDtpKyspewogICAgICAgIHRbaV0gPSAtMTsKICAgIH0KfQogCnZvaWQgcHJpbnRBKCl7CiAgICBpbnQgaTsKICAgIGZvcihpPTE7aTw9c3o7aSsrKSBwcmludGYoIiVkICIsdFtpXSk7CiAgICBwcmludGYoIlxuIik7Cn0KIAp2b2lkIHByaW50VChpbnQgaSl7CiAgICBpbnQgeCA9IGk7CiAgICB3aGlsZSh4LzIhPTApewogICAgICAgIHByaW50ZigiICAiKTsKICAgICAgICB4Lz0yOwogICAgfQogICAgcHJpbnRmKCIlZFxuIix0W2ldKTsKfQogCmludCBnb1AoaW50IGkpewogICAgaWYoaS8yID09IDApIHJldHVybiAwOwogICAgZWxzZSByZXR1cm4gaS8yOwp9CiAKaW50IGdvTChpbnQgaSl7CiAgICBpZigyKmkgPj0gTUFYKSByZXR1cm4gMDsKICAgIGVsc2UgcmV0dXJuIDIqaTsKfQogCmludCBnb1IoaW50IGkpewogICAgaWYoMippKzEgPj0gTUFYKSByZXR1cm4gMDsKICAgIGVsc2UgcmV0dXJuIDIqaSsxOwp9CiAKdm9pZCBwcmVPcmRlcihpbnQgaSl7CiAgICBpZih0W2ldID09IC0xKSByZXR1cm47CiAgICBwcmludFQoaSk7CiAgICBwcmVPcmRlcihnb0woaSkpOwogICAgcHJlT3JkZXIoZ29SKGkpKTsKfQogCnZvaWQgaW5PcmRlcihpbnQgaSl7CiAgICBpZih0W2ldID09IC0xKSByZXR1cm47CiAgICBpbk9yZGVyKGdvTChpKSk7CiAgICBwcmludFQoaSk7CiAgICBpbk9yZGVyKGdvUihpKSk7Cn0KIAp2b2lkIHBvc3RPcmRlcihpbnQgaSl7CiAgICBpZih0W2ldID09IC0xKSByZXR1cm47CiAgICBwb3N0T3JkZXIoZ29MKGkpKTsKICAgIHBvc3RPcmRlcihnb1IoaSkpOwogICAgcHJpbnRUKGkpOwp9CiAKdm9pZCBpbnNCVChpbnQgeCl7CiAgICBpbnQgayxpID0gMTsKICAgIGZvcihrPTA7azxoZWlnaHQ7aysrKXsKICAgICAgICBpZih0W2ldPT0tMSl7CiAgICAgICAgICAgIHRbaV0gPSB4OwogICAgICAgICAgICBzeisrOwogICAgICAgICAgICByZXR1cm47CiAgICAgICAgfQogICAgICAgIGlmKHggPCB0W2ldKSBpID0gZ29MKGkpOwogICAgICAgIGVsc2UgaSA9IGdvUihpKTsKICAgIH0KICAgIHByaW50ZigiRXJyb3IgOiB0b28gaGlnaCAtPiAlZFxuIix4KTsKfQogCi8v5YWI6aCt44Gu6KaB57Sg44KS5Y+W44KK5Ye644GZCi8v44OA44Km44Oz44OS44O844OXCmludCBwb3BIZWFwKCl7Ci8vIOWPluOCiuWHuuOBl+OBpgovLyDmnKvlsL7jga7opoHntKDjgpLmoLnjgavmjIHjgaPjgabjgY3jgaYKLy8g5a2Q44Gu44GG44Gh44Gu5aSn44GN44GE5pa544Go5q+U6LyDCi8vIOWtkOOBruaWueOBjOWkp+OBjeOBi+OBo+OBn+OCieS6pOaPmwovLyDjgpJIZWFw44GM5a6M5oiQ44GZ44KL44G+44Gn57mw44KK6L+U44GZCi8vIOS6i+WJjeWtpue/kuizh+aWmeOBrlAuNTbvvZ42MAogICAgaW50IGEsaSxqLGs9MTsKCiAgICBhPXRbMV07Cgl0WzFdPS0xOwoJCglzd2FwKCZ0WzFdLCZ0W3N6XSk7Cglzei0tOwoKCXdoaWxlKDEpewoJCWk9ayoyOwoJCWo9ayoyKzE7CgkJCgkJaWYodFtpXT50W2pdICYmIHRbaV0+dFtrXSl7CgkJCXN3YXAoJnRbaV0sJnRba10pOwoJCQlrPWk7CgkJfQoJCWVsc2UgaWYodFtpXTx0W2pdICYmIHRbal0+dFtrXSl7CgkJCXN3YXAoJnRbal0sJnRba10pOwoJCQlrPWk7CgkJfQoJCWVsc2UgaWYodFtrXT50W2ldICYmIHRba10+dFtqXSl7CgkJCWJyZWFrOwoJCX0KCQllbHNlIGlmKGk+c3ogfHwgaj5zeil7CgkJCWJyZWFrOwoJCX0KCX0KCQoJcmV0dXJuIGE7CiAgIAogCn0KIAovL+acq+WwvuOBq+imgee0oOOCkui/veWKoOOBmeOCiwovL+OCouODg+ODl+ODkuODvOODlwp2b2lkIHB1c2hIZWFwKGludCB4KXsKLy8g5pyr5bC+44Gr6L+95YqgCi8vIOimquOBqOavlOi8g+OBl+OBpuOAgeWtkOOBruaWueOBjOWkp+OBjeOBi+OBo+OBn+OCieS6pOaPmwovLyDjgpJIZWFw44GM5a6M5oiQ44GZ44KL44G+44Gn57mw44KK6L+U44GZCglpbnQgaTsKICAgIHN6Kys7CiAgICBpID0gc3o7CiAgICB0W2ldID0geDsKCiAgICB3aGlsZShpPjEmJnRbaV0+dFtpLzJdKXsKICAgICAgICBzd2FwKCZ0W2ldLCZ0W2kvMl0pOwogICAgICAgIGk9aS8yOwogICAgfQp9CiAKaW50IG1haW4odm9pZCl7CiAgICBpbnQgaSxuOwogICAgc2NhbmYoIiVkIiwmbik7CiAgICAvLyDmnKjjga7liJ3mnJ/ljJYKICAgIGluaXRUcmVlKG4pOwogICAgZm9yKGk9MDtpPG47aSsrKXsKICAgICAgICBzY2FuZigiJWQiLCZ0W2krMV0pOwogICAgfQogICAgc3ogPSBuOwogICAgLy8g5Lit6ZaT6aCG44Gn6KGo56S6Cglpbk9yZGVyKDEpOwogCgkvLyBwb3AKCWludCBhPXBvcEhlYXAoKTsKCXByaW50ZigicG9wIDogJWRcbiIsYSk7Cglpbk9yZGVyKDEpOwogCglhID0gcG9wSGVhcCgpOwoJcHJpbnRmKCJwb3AgOiAlZFxuIixhKTsKIAoJaW5PcmRlcigxKTsKIAoJcHJpbnRmKCJwdXNoIDogMTAwXG4iKTsKCXB1c2hIZWFwKDEwMCk7Cglpbk9yZGVyKDEpOwogCglwcmludGYoInB1c2ggOiAzMFxuIik7CglwdXNoSGVhcCgzMCk7Cglpbk9yZGVyKDEpOwogICAgcmV0dXJuIDA7Cn0KIA==