| | | @@ 変更前 3行目〜 / 変更後 3行目〜 @@ |
| 3 | 3 | | |
| 4 | 4 | | /** |
| 5 | 5 | | * ローレル指数。 |
| 6 | | − | * v1: 身長・体重を float に読み、体重 ÷ (身長 ÷ 100) ** 3 × 10 を計算して、区分は丸める前の指数を < で比べ、 |
| 7 | | − | * 指数とふつうの体重の目安を sprintf('%.1f') で書く素直な書き方。 |
| 6 | + | * v2: 身長・体重の文字列を float に読まず、小数点を外した数字の列(10^桁数 倍した整数)として10進のまま読み、 |
| 7 | + | * 指数を「分子 ÷ 分母」の整数の組で正確に持つ。区分は分子と「境目 × 分母」を比べ、 |
| 8 | + | * 表示は小数第1位に「0.5 ちょうどは大きい方へ」四捨五入する。 |
| 9 | + | * v1 の float では 160cm・59.392kg の指数 145 ちょうどが 144.99999999999997 になって「ふつう(標準)」に落ち、 |
| 10 | + | * 160cm・12.8kg の 31.25 が 31.249999999999996 になって sprintf('%.1f') で 31.2 に落ちた。 |
| 11 | + | * 身長の3乗に 10 の累乗を掛けると int の上限(約 9.2 × 10^18)を超えうるので、整数は数字の列のまま筆算する |
| 12 | + | * (bcmath・gmp は入っていない環境があるので使わない)。 |
| 8 | 13 | | */ |
| 14 | + | const DECIMAL = '/^[0-9]+(\.[0-9]+)?$/D'; // D: $ を末尾の改行の前で一致させない |
| 15 | + | const BANDS = [['100', 'やせ'], ['115', 'やせ気味'], ['145', 'ふつう(標準)'], ['160', '太り気味']]; |
| 9 | 16 | | |
| 10 | | − | function roehrerCategory(float $index): string |
| 17 | + | /** 先頭の 0 を外す(0 は "0") */ |
| 18 | + | function digitsNorm(string $a): string |
| 11 | 19 | | { |
| 12 | | − | if ($index < 100) return 'やせ'; |
| 13 | | − | if ($index < 115) return 'やせ気味'; |
| 14 | | − | if ($index < 145) return 'ふつう(標準)'; |
| 15 | | − | if ($index < 160) return '太り気味'; |
| 16 | | − | return '肥満'; |
| 20 | + | $a = ltrim($a, '0'); |
| 21 | + | return $a === '' ? '0' : $a; |
| 22 | + | } |
| 23 | + | |
| 24 | + | /** 0 以上の整数(数字の列)の大小。-1・0・1 */ |
| 25 | + | function digitsCmp(string $a, string $b): int |
| 26 | + | { |
| 27 | + | return strlen($a) <=> strlen($b) ?: strcmp($a, $b) <=> 0; |
| 28 | + | } |
| 29 | + | |
| 30 | + | function digitsAdd(string $a, string $b): string |
| 31 | + | { |
| 32 | + | $n = max(strlen($a), strlen($b)); |
| 33 | + | $a = str_pad($a, $n, '0', STR_PAD_LEFT); |
| 34 | + | $b = str_pad($b, $n, '0', STR_PAD_LEFT); |
| 35 | + | $out = ''; |
| 36 | + | $carry = 0; |
| 37 | + | for ($i = $n - 1; $i >= 0; $i--) { |
| 38 | + | $s = (int) $a[$i] + (int) $b[$i] + $carry; |
| 39 | + | $out = ($s % 10) . $out; |
| 40 | + | $carry = intdiv($s, 10); |
| 41 | + | } |
| 42 | + | return digitsNorm(($carry ? '1' : '') . $out); |
| 43 | + | } |
| 44 | + | |
| 45 | + | /** $a - $b($a >= $b のとき) */ |
| 46 | + | function digitsSub(string $a, string $b): string |
| 47 | + | { |
| 48 | + | $b = str_pad($b, strlen($a), '0', STR_PAD_LEFT); |
| 49 | + | $out = ''; |
| 50 | + | $borrow = 0; |
| 51 | + | for ($i = strlen($a) - 1; $i >= 0; $i--) { |
| 52 | + | $d = (int) $a[$i] - (int) $b[$i] - $borrow; |
| 53 | + | $borrow = $d < 0 ? 1 : 0; |
| 54 | + | $out = ($d + 10 * $borrow) . $out; |
| 55 | + | } |
| 56 | + | return digitsNorm($out); |
| 57 | + | } |
| 58 | + | |
| 59 | + | function digitsMul(string $a, string $b): string |
| 60 | + | { |
| 61 | + | $la = strlen($a); |
| 62 | + | $lb = strlen($b); |
| 63 | + | $acc = array_fill(0, $la + $lb, 0); |
| 64 | + | for ($i = $la - 1; $i >= 0; $i--) { |
| 65 | + | for ($j = $lb - 1; $j >= 0; $j--) { |
| 66 | + | $acc[$i + $j + 1] += (int) $a[$i] * (int) $b[$j]; |
| 67 | + | } |
| 68 | + | } |
| 69 | + | for ($k = $la + $lb - 1; $k > 0; $k--) { |
| 70 | + | $acc[$k - 1] += intdiv($acc[$k], 10); |
| 71 | + | $acc[$k] %= 10; |
| 72 | + | } |
| 73 | + | return digitsNorm(implode('', $acc)); |
| 74 | + | } |
| 75 | + | |
| 76 | + | /** 筆算の割り算。[商, 余り] */ |
| 77 | + | function digitsDivMod(string $a, string $b): array |
| 78 | + | { |
| 79 | + | $q = ''; |
| 80 | + | $r = '0'; |
| 81 | + | for ($i = 0, $n = strlen($a); $i < $n; $i++) { |
| 82 | + | $r = digitsNorm($r . $a[$i]); |
| 83 | + | $d = 0; |
| 84 | + | while (digitsCmp($r, $b) >= 0) { |
| 85 | + | $r = digitsSub($r, $b); |
| 86 | + | $d++; |
| 87 | + | } |
| 88 | + | $q .= $d; |
| 89 | + | } |
| 90 | + | return [digitsNorm($q), $r]; |
| 91 | + | } |
| 92 | + | |
| 93 | + | /** 0 以上の $num ÷ $den を小数第1位に四捨五入(0.5 ちょうどは大きい方へ)して "31.3" の形にする */ |
| 94 | + | function fmt1(string $num, string $den): string |
| 95 | + | { |
| 96 | + | [$q, $r] = digitsDivMod($num . '0', $den); // 10 倍して割る |
| 97 | + | if (digitsCmp(digitsAdd($r, $r), $den) >= 0) { |
| 98 | + | $q = digitsAdd($q, '1'); |
| 99 | + | } |
| 100 | + | $q = str_pad($q, 2, '0', STR_PAD_LEFT); |
| 101 | + | return substr($q, 0, -1) . '.' . substr($q, -1); |
| 102 | + | } |
| 103 | + | |
| 104 | + | /** "21.97" → ['2197', 2](値 = 数字の列 ÷ 10^桁数)。読めなければ null */ |
| 105 | + | function readDecimal($s): ?array |
| 106 | + | { |
| 107 | + | if (!is_string($s) || !preg_match(DECIMAL, $s)) { |
| 108 | + | return null; |
| 109 | + | } |
| 110 | + | [$int, $frac] = explode('.', $s . '.'); |
| 111 | + | return [digitsNorm($int . $frac), strlen($frac)]; |
| 17 | 112 | | } |
| 18 | 113 | | |
| 19 | 114 | | /** @param array{height_cm: string, weight_kg: string} $in */ |
| 20 | 115 | | function roehrer(array $in): array |
| 21 | 116 | | { |
| 22 | | − | if (!is_numeric($in['height_cm']) || !is_numeric($in['weight_kg'])) { |
| 117 | + | $h = readDecimal($in['height_cm'] ?? null); |
| 118 | + | $w = readDecimal($in['weight_kg'] ?? null); |
| 119 | + | if ($h === null || $w === null || $h[0] === '0' || $w[0] === '0') { |
| 23 | 120 | | return ['error' => 'invalid_input']; |
| 24 | 121 | | } |
| 25 | | − | $height = (float) $in['height_cm']; |
| 26 | | − | $weight = (float) $in['weight_kg']; |
| 27 | | − | if ($height <= 0 || $weight <= 0) { |
| 28 | | − | return ['error' => 'invalid_input']; |
| 122 | + | // 身長(m)³ = h³ ÷ 10^(3 × 桁数 + 6) |
| 123 | + | $cubeNum = digitsMul(digitsMul($h[0], $h[0]), $h[0]); |
| 124 | + | $cubeDen = '1' . str_repeat('0', 3 * $h[1] + 6); |
| 125 | + | // 指数 = 体重 ÷ 身長(m)³ × 10 = (w × cubeDen × 10) ÷ (10^桁数 × cubeNum) |
| 126 | + | $num = $w[0] . substr($cubeDen, 1) . '0'; |
| 127 | + | $den = $cubeNum . str_repeat('0', $w[1]); |
| 128 | + | $category = '肥満'; |
| 129 | + | foreach (BANDS as [$bound, $name]) { |
| 130 | + | if (digitsCmp($num, digitsMul($bound, $den)) < 0) { |
| 131 | + | $category = $name; |
| 132 | + | break; |
| 133 | + | } |
| 29 | 134 | | } |
| 30 | | − | $hm = $height / 100; |
| 31 | | − | $index = $weight / $hm ** 3 * 10; |
| 32 | 135 | | return [ |
| 33 | | − | 'index' => sprintf('%.1f', $index), |
| 34 | | − | 'category' => roehrerCategory($index), |
| 35 | | − | 'normal_low' => sprintf('%.1f', 115 * $hm ** 3 / 10), |
| 36 | | − | 'normal_high' => sprintf('%.1f', 145 * $hm ** 3 / 10), |
| 136 | + | 'index' => fmt1($num, $den), |
| 137 | + | 'category' => $category, |
| 138 | + | 'normal_low' => fmt1(digitsMul('115', $cubeNum), $cubeDen . '0'), |
| 139 | + | 'normal_high' => fmt1(digitsMul('145', $cubeNum), $cubeDen . '0'), |
| 37 | 140 | | ]; |
| 38 | 141 | | } |