数学 - JavaScript - パスカルの三角形と剰余 - パスカルの三角形に色を塗る - 剰余と座標のリスト、描画ソフトへの対応、連想リストとしてのデータ

素数夜曲

開発環境

• macOS Sierra - Apple (OS)
• Emacs (Text Editor)
• JavaScript (プログラミング言語)
• D3.js (JavaScript Library)
• Safari(Web browser)
• 参考書籍
  • JavaScript 第6版 (David Flanagan(著)、村上 列(翻訳)、オライリージャパン)
  • JavaScriptリファレンス 第6版(David Flanagan(著)、木下 哲也(翻訳)、オライリージャパン)
  • インタラクティブ・データビジュアライゼーション(Scott Murray (著)、長尾 高弘 (翻訳)、オライリージャパン)

素数夜曲 (吉田 武 (著)、東海大学出版会)の付録D(女王陛下のLISP)、D.6(パスカルの三角形と剰余)、D.6.2(パスカルの三角形に色を塗る)、剰余と座標のリスト、D.6.3(描画ソフトへの対応)、連想リストとしてのデータを JavaScript で取り組んでみる。

コード(Emacs)

表示が違うの、計算による誤差の発生(Math.round で丸めない場合よりは綺麗になった)が原因か、コードに間違いがあるのか。。

HTML5

<div id="graph0"></div>
<label for="row0">段数: </label>
<input id="row0" type="number" min="0" step="1" value="64">
<label for="m0">除数: </label>
<input id="m0" type="number" min="1" step="1" value="2">
<label for="r0">点の半径: </label>
<input id="r0" type="number" min="1" step="1" value="2">

<script type="text/javascript" src="https://cdnjs.cloudflare.com/ajax/libs/d3/4.2.6/d3.min.js" integrity="sha256-5idA201uSwHAROtCops7codXJ0vja+6wbBrZdQ6ETQc=" crossorigin="anonymous"></script>

<script src="sample5.js"></script>

JavaScript

{
    let div_graph = document.querySelector('#graph0'),
        input_row = document.querySelector('#row0'),
        input_m = document.querySelector('#m0'),
        input_r = document.querySelector('#r0'),
        inputs = [input_row, input_m, input_r],
        width = 600,
        height = 600,
        padding = 20,
        svg;

    let iota = (start, end, step=1) => {
        let result = [];
        
        for (let i = start; i <= end; i += step) {
            result.push(i);
        }
        return result;
    };
    let fact = (n) => {
        let result = 1;
        for (let i = 2; i <= n; i += 1) {
            result *= i;
        }
        return result;
    };
    let bc = (n, r) => {
        return Math.round(fact(n) / (fact(r) * fact(n - r)));
    };
    let row = (n) => {
        return iota(0, n).map((x) => bc(n, x));
    };
    let row_mod = (n, m) => {
        return row(n).map((x) => Math.round(x % m));
    };
    let tri_n = (n) => {
        let y = -n;
        
        return iota(y, n, 2).map((x) => [x, y]);
    };
    let mod_xy = (n, m) => {
        let a = tri_n(n),
            b = row_mod(n, m),
            result = [];

        for (let i = 0; i < a.length; i += 1) {
            let x = a[i].slice();
            x.push(b[i]);
            result.push(x);
        }
        return result;
    };
    let plot_data = (n, m) => {
        return iota(0, n).map((x) => mod_xy(x, m))
            .reduce((x, y) => x.concat(y), []);
    };
    let extreme = (lst) => {
        let a = [];

        for (let i = 0; i < lst[0].length; i += 1) {
            let b = [];
            for (let j = 0; j < lst.length; j += 1) {
                b.push(lst[j][i]);
            }
            a.push(b);
        }
        return a.map((x) => [Math.min.apply(null, x), Math.max.apply(null, x)]);
    };

    let key1 = (obj, lst) => {
        let loop = (lst, tmp) => {
            if (lst.length === 0) {
                return tmp;
            }
            if (lst[0][2] === obj) {
                tmp.unshift(lst[0]);
                return loop(lst.slice(1), tmp);
            }
            return loop(lst.slice(1), tmp);
        };
        return loop(lst, []);
    }
    let build = (n, m) => {
        let lst = plot_data(n, m),
            min_v = extreme(lst)[2][0],
            max_v = extreme(lst)[2][1];
        
        let loop = (k, tmp) => {
            if (max_v < k) {
                tmp.reverse();
                return tmp;
            }
            return loop(k + 1, key1(k, lst).concat(tmp));
        }
        return loop(min_v, []);
    };
    let draw = () => {
        let n = parseInt(input_row.value, 10) - 1,
            m = parseInt(input_m.value, 10),
            r = parseInt(input_r.value, 10),
            dataset = build(n, m),
            xscale,
            yscale;

        xscale = d3.scaleLinear()
            .domain([-n, n])
            .range([padding, width - padding]);
        yscale = d3.scaleLinear()
            .domain([0, -n])
            .range([padding, height * Math.sqrt(3) / 2 - padding]);

        div_graph.innerHTML = '';
        svg = d3.select('#graph0')
            .append('svg')
            .attr('width', width)
            .attr('height', height);

        svg.selectAll('circle')
            .data(dataset)
            .enter()
            .append('circle')
            .attr('cx', (d) => xscale(d[0]))
            .attr('cy', (d) => yscale(d[1]))
            .attr('r', r)
            .attr('fill', (d, i) => {
                let n = d[2] % 3;

                if (n === 0) {
                    return 'red';
                }
                if (n == 1) {
                    return 'green';
                }
                return 'blue';
            });
    };

    inputs.forEach((input) => input.onchange = draw);
    draw();
}

段数: 除数: 点の半径: