#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
); }
//先頭の要素を取り出す
//ダウンヒープ
void pushHeap(int x){
sz++;
t[sz] = x;
int i = sz;
int p = goP(i);
while(p != 0 && t[p] < t[i]){
swap(&t[p], &t[i]);
i = p;
p = goP(i);
}
}
//末尾に要素を追加する
//アップヒープ
int popHeap(){
if(sz == 0) return -1;
int ret = t[1];
t[1] = t[sz];
t[sz] = -1;
sz--;
int i = 1;
while(1){
int L = goL(i);
int R = goR(i);
if(L == 0 || t[L] == -1) break;
int larger = L;
if(R != 0 && t[R] != -1 && t[R] > t[L]){
larger = R;
}
if(t[larger] > t[i]){
swap(&t[larger], &t[i]);
i = larger;
} else {
break;
}
}
return ret;
}
int main(void){
int i,x,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;
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxzdGRsaWIuaD4KIAojZGVmaW5lIGhlaWdodCA0CiNkZWZpbmUgTUFYICgxPDxoZWlnaHQpICAvL+ODk+ODg+ODiOOCt+ODleODiOa8lOeulyAyXmhlaWdodCDjgajlkIzjgZgKIAppbnQgdFtNQVgrMV07IC8v6YWN5YiX5aSW44Ki44Kv44K744K56Ziy5q2i44Gu44Gf44KB44Gu44OA44Of44O844Gn77yL77yRCmludCBzeiA9IDA7CiAKdm9pZCBzd2FwKGludCAqeCwgaW50ICp5KXsKICAgIGludCB0bXAgPSAqeDsKICAgICp4ID0gKnk7CiAgICAqeSA9IHRtcDsKfQogCnZvaWQgaW5pdFRyZWUoaW50IG4pewogICAgaW50IGk7CiAgICBmb3IoaT0wO2k8TUFYO2krKyl7CiAgICAgICAgdFtpXSA9IC0xOwogICAgfQp9CiAKdm9pZCBwcmludEEoKXsKICAgIGludCBpOwogICAgZm9yKGk9MTtpPD1zejtpKyspIHByaW50ZigiJWQgIix0W2ldKTsKICAgIHByaW50ZigiXG4iKTsKfQogCnZvaWQgcHJpbnRUKGludCBpKXsKICAgIGludCB4ID0gaTsKICAgIHdoaWxlKHgvMiE9MCl7CiAgICAgICAgcHJpbnRmKCIgICIpOwogICAgICAgIHgvPTI7CiAgICB9CiAgICBwcmludGYoIiVkXG4iLHRbaV0pOwp9CiAKaW50IGdvUChpbnQgaSl7CiAgICBpZihpLzIgPT0gMCkgcmV0dXJuIDA7CiAgICBlbHNlIHJldHVybiBpLzI7Cn0KIAppbnQgZ29MKGludCBpKXsKICAgIGlmKDIqaSA+PSBNQVgpIHJldHVybiAwOwogICAgZWxzZSByZXR1cm4gMippOwp9CiAKaW50IGdvUihpbnQgaSl7CiAgICBpZigyKmkrMSA+PSBNQVgpIHJldHVybiAwOwogICAgZWxzZSByZXR1cm4gMippKzE7Cn0KIAp2b2lkIHByZU9yZGVyKGludCBpKXsKICAgIGlmKHRbaV0gPT0gLTEpIHJldHVybjsKICAgIHByaW50VChpKTsKICAgIHByZU9yZGVyKGdvTChpKSk7CiAgICBwcmVPcmRlcihnb1IoaSkpOwp9CiAKdm9pZCBpbk9yZGVyKGludCBpKXsKICAgIGlmKHRbaV0gPT0gLTEpIHJldHVybjsKICAgIGluT3JkZXIoZ29MKGkpKTsKICAgIHByaW50VChpKTsKICAgIGluT3JkZXIoZ29SKGkpKTsKfQogCnZvaWQgcG9zdE9yZGVyKGludCBpKXsKICAgIGlmKHRbaV0gPT0gLTEpIHJldHVybjsKICAgIHBvc3RPcmRlcihnb0woaSkpOwogICAgcG9zdE9yZGVyKGdvUihpKSk7CiAgICBwcmludFQoaSk7Cn0KIAp2b2lkIGluc0JUKGludCB4KXsKICAgIGludCBrLGkgPSAxOwogICAgZm9yKGs9MDtrPGhlaWdodDtrKyspewogICAgICAgIGlmKHRbaV09PS0xKXsKICAgICAgICAgICAgdFtpXSA9IHg7CiAgICAgICAgICAgIHN6Kys7CiAgICAgICAgICAgIHJldHVybjsKICAgICAgICB9CiAgICAgICAgaWYoeCA8IHRbaV0pIGkgPSBnb0woaSk7CiAgICAgICAgZWxzZSBpID0gZ29SKGkpOwogICAgfQogICAgcHJpbnRmKCJFcnJvciA6IHRvbyBoaWdoIC0+ICVkXG4iLHgpOwp9CiAKLy/lhYjpoK3jga7opoHntKDjgpLlj5bjgorlh7rjgZkKLy/jg4Djgqbjg7Pjg5Ljg7zjg5cKdm9pZCBwdXNoSGVhcChpbnQgeCl7CiAgICBzeisrOwogICAgdFtzel0gPSB4OwoKICAgIGludCBpID0gc3o7CiAgICBpbnQgcCA9IGdvUChpKTsKCiAgICB3aGlsZShwICE9IDAgJiYgdFtwXSA8IHRbaV0pewogICAgICAgIHN3YXAoJnRbcF0sICZ0W2ldKTsKICAgICAgICBpID0gcDsKICAgICAgICBwID0gZ29QKGkpOwogICAgfQp9CgogCi8v5pyr5bC+44Gr6KaB57Sg44KS6L+95Yqg44GZ44KLCi8v44Ki44OD44OX44OS44O844OXCmludCBwb3BIZWFwKCl7CiAgICBpZihzeiA9PSAwKSByZXR1cm4gLTE7CgogICAgaW50IHJldCA9IHRbMV07CgogICAgdFsxXSA9IHRbc3pdOwogICAgdFtzel0gPSAtMTsKICAgIHN6LS07CgogICAgaW50IGkgPSAxOwoKICAgIHdoaWxlKDEpewogICAgICAgIGludCBMID0gZ29MKGkpOwogICAgICAgIGludCBSID0gZ29SKGkpOwoKICAgICAgICBpZihMID09IDAgfHwgdFtMXSA9PSAtMSkgYnJlYWs7CgogICAgICAgIGludCBsYXJnZXIgPSBMOwogICAgICAgIGlmKFIgIT0gMCAmJiB0W1JdICE9IC0xICYmIHRbUl0gPiB0W0xdKXsKICAgICAgICAgICAgbGFyZ2VyID0gUjsKICAgICAgICB9CgogICAgICAgIGlmKHRbbGFyZ2VyXSA+IHRbaV0pewogICAgICAgICAgICBzd2FwKCZ0W2xhcmdlcl0sICZ0W2ldKTsKICAgICAgICAgICAgaSA9IGxhcmdlcjsKICAgICAgICB9IGVsc2UgewogICAgICAgICAgICBicmVhazsKICAgICAgICB9CiAgICB9CgogICAgcmV0dXJuIHJldDsKfQogCmludCBtYWluKHZvaWQpewogICAgaW50IGkseCxuOwogICAgc2NhbmYoIiVkIiwmbik7CiAgICAvLyDmnKjjga7liJ3mnJ/ljJYKICAgIGluaXRUcmVlKG4pOwogICAgZm9yKGk9MDtpPG47aSsrKXsKICAgICAgICBzY2FuZigiJWQiLCZ0W2krMV0pOwogICAgfQogICAgc3ogPSBuOwogICAgLy8g5Lit6ZaT6aCG44Gn6KGo56S6Cglpbk9yZGVyKDEpOwogCgkvLyBwb3AKCWludCBhPXBvcEhlYXAoKTsKCXByaW50ZigicG9wIDogJWRcbiIsYSk7Cglpbk9yZGVyKDEpOwogCglhID0gcG9wSGVhcCgpOwoJcHJpbnRmKCJwb3AgOiAlZFxuIixhKTsKIAoJaW5PcmRlcigxKTsKIAoJcHJpbnRmKCJwdXNoIDogMTAwXG4iKTsKCXB1c2hIZWFwKDEwMCk7Cglpbk9yZGVyKDEpOwogCglwcmludGYoInB1c2ggOiAzMFxuIik7CglwdXNoSGVhcCgzMCk7Cglpbk9yZGVyKDEpOwogICAgcmV0dXJuIDA7Cn0KIA==