#!/usr/local/bin/tcc -run
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define NUM 11 // 母集団数は11
typedef struct _people{
int height;
int weight;
}people; // 身長、体重を有するpeople という名の構造体を定義

 void cal_MAX(people data[]){ // 身長、体重の最大値を計算
 int hMAX = data[0].height; // 身長の最大値を仮定
 int wMAX = data[0].weight; // 体重の最大値を仮定
 int i;// カウンタ
 for(i=0; i<NUM; i++){
   if(hMAX < data[i].height){
      hMAX = data[i].height; // 身長の最大値を更新
   }
   if(wMAX < data[i].weight){
    wMAX = data[i].weight; // 体重の最大値を更新
     }
   }
   printf("%d 人の身長の最大値は%d[cm]、体重の最大値は%d[kg]\n", NUM, hMAX, wMAX); // 身長、体重の最大値を表示
  }

 void cal_min(people data[]){ // 身長、体重の最小値を計算
 int hmin = data[0].height; // 身長の最小値を仮定
 int wmin = data[0].weight; // 体重の最小値を仮定
 int i;// カウンタ
 for(i=0; i<NUM; i++){
   if(hmin > data[i].height){
     hmin = data[i].height; // 身長の最小値を更新
   }
   if(wmin > data[i].weight){
 wmin = data[i].weight; // 体重の最小値を更新
   }
   }
 printf("%d 人の身長の最小値は%d[cm]、体重の最小値は%d[kg]\n", NUM, hmin, wmin);
 // 身長、体重の最小値を表示
  }

void cal_Med(people data[]){
 int hTmp, wTmp, hMed, wMed;
 int i, j;// カウンタ
 for(i=0; i<NUM; i++){ /* バブルソートで昇順に並べ替え*/
    for(j=0; j<NUM -i; j++){
      if(data[j].height > data[j + 1].height){ /* 前の要素が大きかったら*/
 hTmp = data[j].height; /* 並べ替え*/
 data[j].height = data[j + 1].height;
 data[j + 1].height = hTmp;
      }
      if(data[j].weight > data[j + 1].weight){
 wTmp = data[j].weight;
 data[j].weight = data[j + 1].weight;
 data[j + 1].weight = wTmp;
      }
      }
     }
 hMed = data[ NUM / 2 + 1/ 2].height; // 身長の中央値を計算
 wMed = data[ NUM / 2 + 1/ 2].weight; //体重の中央値を計算
 printf("%d 人の身長の中央値は%d[cm]、体重の中央値は%d[kg]\n", NUM, hMed, wMed);
 // 身長、体重の中央値を表示
  }

void cal_Mean_STDEV_CC(people data[]){
 // 身長・体重の平均値、標準偏差、相関係数を計算、表示する関数
 double h_Sum = 0.0, xh_Sum = 0.0, w_Sum = 0.0, xw_Sum = 0.0, xhw_Sum;
 // 身長の和・平方和、体重の和・平方和、身長と体重の積和の変数
 double h_Mean = 0.0, w_Mean = 0.0; // 身長・体重の平均値の変数
 double hDisp = 0.0, wDisp = 0.0, hSTDEV = 0.0, wSTDEV = 0.0;
 // 身長・体重の分散、標準偏差の変数
 double Cov = 0.0, CC = 0.0; // 身長・体重の共分散、相関係数の変数
 int i; // カウンタ
 for(i=0; i <NUM; i++){
 h_Sum += data[i].height;// 和
 w_Sum += data[i].weight;
 xh_Sum += (data[i].height * data[i].height);// 平方和
 xw_Sum += (data[i].weight * data[i].weight);
 xhw_Sum += (data[i].height) * (data[i].weight);// 積和
  }

 h_Mean = h_Sum / NUM;// 平均値
 w_Mean = w_Sum / NUM;

 hDisp = (xh_Sum / NUM) - pow(h_Mean, 2);// 分散
 wDisp = (xw_Sum / NUM) - pow(w_Mean, 2);

 hSTDEV = sqrt(hDisp);// 標準偏差
 wSTDEV = sqrt(wDisp);

 Cov = (xhw_Sum / NUM) - (h_Mean * w_Mean);// 共分散

 CC = Cov / (hSTDEV * wSTDEV );// 相関係数

 printf("%d 人の身長の平均値は%.1f[cm]、体重の平均値は%.1f[kg]\n", NUM, h_Mean, w_Mean);
 // 身長・体重の平均値の表示
   printf("%d 人の身長の標準偏差は%.2f、体重の標準偏差は%.2f\n", NUM, hSTDEV, wSTDEV);
 // 身長・体重の標準偏差の表示
  printf("%d 人の身長と体重の相関係数は%.3f\n", NUM, CC); // 身長と体重の相関係数の表示
}
  int main(void){
  people data[] = { // 各配列に身長と体重を与える
    {168, 60},
    {162, 55},
    {171, 65},
    {159, 60},
    {174, 72},
    {166, 57},
    {175, 61},
    {176, 66},
    {168, 56},
    {182, 69},
    {179, 62}
 };
 cal_MAX(data); // 与えられたデータの最大値を計算、表示
 cal_min(data); // 与えられたデータの最小値を計算、表示
 cal_Mean_STDEV_CC(data); // 与えられたデータの平均値、標準偏差、相関係数を計算、表示

 cal_Med(data); // 与えられたデータの中央値を計算

 return 0;

   }