double napier(int n) {
int box[CNTS];
for(int i = 0; i < CNTS; i ++) { // 箱を1で初期化
box[i] = 1;
}
int cnt = 0;
int j;
for(int i = 0; i < CNTS; i ++) {
j = (int) CNTS * genrand_res53(); // 入る箱の番号
box[j] = 0;
}
for(int i = 0; i < CNTS; i ++) { // 残った箱をカウント
cnt += box[i];
}
double napier = (double) CNTS / cnt;
return napier;
}