個人用メモ帳

2018年12月6日

C++ の動作確認をしてみた（３８）

C++の練習を兼ねて, AtCoder Beginner Contest 039 の問題C（ピアニスト高橋君) ～問題D (画像処理高橋君) を解いてみた.

■感想.
1. 問題D は, impossible の判定に, 画像復元 → 画像収縮を行った上で, 入力情報と比較する必要があることに気付いたため, 解答に辿り着くことが出来たように思う.
考え方のベースが, AtCoder Beginner Contest 075 (B – Minesweeper) に近いように思った.
但し, 結局ゴリゴリ書く羽目になって, 非常にきつかった.
2. 解説と, ほぼ同じ方針だったので, 及第点は, 取れたように思う.

本家のサイトABC039解説をご覧下さい.

■C++版プログラム(問題C／AC版).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {
    
    // 1. 入力情報取得.
    string S;
    cin >> S;
    
    // 2. 与えられた文字列データから, 一オクターブ の 開始位置 を取得.
    // 白黒白黒白白黒白黒白黒白 ⇔ "ド ド# レ レ# ミ ファ ファ# ソ ソ# ラ ラ# シ"
    // なので, "WBWBWWBWBWBW" を 一オクターブ と見て考える.
    string ans = "";
    if(S == "WBWBWWBWBWBWWBWBWWBW") ans = "Do";
    if(S == "WBWWBWBWBWWBWBWWBWBW") ans = "Re";
    if(S == "WWBWBWBWWBWBWWBWBWBW") ans = "Mi";
    if(S == "WBWBWBWWBWBWWBWBWBWW") ans = "Fa";
    if(S == "WBWBWWBWBWWBWBWBWWBW") ans = "So";
    if(S == "WBWWBWBWWBWBWBWWBWBW") ans = "La";
    if(S == "WWBWBWWBWBWBWWBWBWWB") ans = "Si";
    
    // 3. 後処理.
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

string S;

cin >> S;

// 2. 与えられた文字列データから, 一オクターブの開始位置を取得.

// 白黒白黒白白黒白黒白黒白 ⇔ "ドド# レレ# ミファファ# ソソ# ララ# シ"

// なので, "WBWBWWBWBWBW" を一オクターブと見て考える.

string ans = "";

if(S == "WBWBWWBWBWBWWBWBWWBW") ans = "Do";

if(S == "WBWWBWBWBWWBWBWWBWBW") ans = "Re";

if(S == "WWBWBWBWWBWBWWBWBWBW") ans = "Mi";

if(S == "WBWBWBWWBWBWWBWBWBWW") ans = "Fa";

if(S == "WBWBWWBWBWWBWBWBWWBW") ans = "So";

if(S == "WBWWBWBWWBWBWBWWBWBW") ans = "La";

if(S == "WWBWBWWBWBWBWWBWBWWB") ans = "Si";

// 3. 後処理.

cout << ans << endl;

return 0;

}

■C++版プログラム(問題D／AC版).

// AtCoder Beginner Contest 075 (B - Minesweeper) を 参照.
#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

// 指定地点(i, j) の周囲に配置されている'#', '$' の個数(上下左右斜め方向)をカウント.
// @param i: (行方向)i番目.
// @param j: (列方向)j番目.
// @param a[]: マス目.
// @return: '#', '$' の個数.
int countSharpAndDollar(int i, int j, string a[]) {
    const char SHARP_SIGN  = '#';
    const char DOLLAR_SIGN = '$';
    int ret = 0;
    if(a[i - 1][j - 1] == SHARP_SIGN || a[i - 1][j - 1] == DOLLAR_SIGN) ret++;    // 左上.
    if(a[i - 1][j]     == SHARP_SIGN || a[i - 1][j]     == DOLLAR_SIGN) ret++;    // 上.
    if(a[i - 1][j + 1] == SHARP_SIGN || a[i - 1][j + 1] == DOLLAR_SIGN) ret++;    // 右上.
    if(a[i][j - 1]     == SHARP_SIGN || a[i][j - 1]     == DOLLAR_SIGN) ret++;    // 左.
    if(a[i][j]         == SHARP_SIGN || a[i][j]         == DOLLAR_SIGN) ret++;    // 中央.
    if(a[i][j + 1]     == SHARP_SIGN || a[i][j + 1]     == DOLLAR_SIGN) ret++;    // 右.
    if(a[i + 1][j - 1] == SHARP_SIGN || a[i + 1][j - 1] == DOLLAR_SIGN) ret++;    // 左下.
    if(a[i + 1][j]     == SHARP_SIGN || a[i + 1][j]     == DOLLAR_SIGN) ret++;    // 下.
    if(a[i + 1][j + 1] == SHARP_SIGN || a[i + 1][j + 1] == DOLLAR_SIGN) ret++;    // 右下.
    return ret;
}

// 指定地点(i, j) の周囲に配置されている'#' の個数(上下左右斜め方向)をカウント.
// @param i: (行方向)i番目.
// @param j: (列方向)j番目.
// @param a[]: マス目.
// @return: '#' の個数.
int countSharp(int i, int j, string a[]) {
    const char SHARP_SIGN  = '#';
    int ret = 0;
    if(a[i - 1][j - 1] == SHARP_SIGN) ret++;    // 左上.
    if(a[i - 1][j]     == SHARP_SIGN) ret++;    // 上.
    if(a[i - 1][j + 1] == SHARP_SIGN) ret++;    // 右上.
    if(a[i][j - 1]     == SHARP_SIGN) ret++;    // 左.
    if(a[i][j + 1]     == SHARP_SIGN) ret++;    // 右.
    if(a[i + 1][j - 1] == SHARP_SIGN) ret++;    // 左下.
    if(a[i + 1][j]     == SHARP_SIGN) ret++;    // 下.
    if(a[i + 1][j + 1] == SHARP_SIGN) ret++;    // 右下.
    return ret;
}

int main() {
    
    // 1. 入力情報取得.
    int H, W;
    cin >> H >> W;
    // 調査対象のマス目の外側を, '$' で囲む.
    // 例)
    // [入力値]
    // .....
    // .#.#.
    // .....
    // 
    // ↓
    // 
    // ['$'で囲んだ状態]
    // $$$$$$$
    // $.....$
    // $.#.#.$
    // $.....$
    // $$$$$$$
    
    const char DOLLAR_SIGN = '$';
    const char SHARP_SIGN = '#';
    const char A_SIGN = 'A';
    
    string a;
    FOR(i, 0, W + 2) a += DOLLAR_SIGN;
    
    string I[H + 2] = {};
    I[0] = a;
    I[H + 1] = a;
    FOR(i, 1, H + 1) {
        string s;
        cin >> s;
        I[i] = (DOLLAR_SIGN + s + DOLLAR_SIGN);
    }
    
    // 2. マス目を一つ一つ調べ, '#', '$' の個数を数える.
    // 出力用のマス目は, シャープ記号で初期化しておく.
    string b;
    FOR(i, 0, W) b += SHARP_SIGN;
    string O[H] = {};
    FOR(i, 0, H) O[i] = b;
    
    // 3. 上下左右斜め の 八方向でカウント.
    FOR(i, 1, H + 1) FOR(j, 1, W + 1) if(countSharpAndDollar(i, j, I) != 9) O[i - 1][j - 1] = '.';
    
    // 4. 出力用のマス目から, 再度収縮を行って, チェック用のマス目を作成し, 
    // 入力データと比較.
    string CHECK[H + 2] = {};
    CHECK[0] = a;
    CHECK[H + 1] = a;
    // 4-1. チェック用のマス目作成.
    FOR(i, 1, H + 1) CHECK[i] = (DOLLAR_SIGN + O[i - 1] + DOLLAR_SIGN);
    FOR(i, 1, H + 1) {
        FOR(j, 1, W + 1) {
            // cout << CHECK[i][j] << " " << countSharp(i, j, CHECK) << endl;
            if(CHECK[i][j] == SHARP_SIGN || CHECK[i][j] == A_SIGN)   continue;
            if(countSharp(i, j, CHECK) > 0 && CHECK[i][j] != A_SIGN) CHECK[i][j] = A_SIGN;
        }
    }
    // 4-2. 入力時のマス目(I) と チェック用のマス目(CHECK)について,
    // I について '#' の個数, CHECK について, '#' と 'A' の 合計個数を比較して, 
    // 同じだったら, "possible", 異なったら, "impossible" と判定.
    bool ok = true;
    int ci = 0;
    int cc = 0;
    FOR(i, 1, H + 1) FOR(j, 1, W + 1) if(I[i][j] == SHARP_SIGN) ci++;
    FOR(i, 1, H + 1) FOR(j, 1, W + 1) if(CHECK[i][j] == SHARP_SIGN || CHECK[i][j] == A_SIGN) cc++;
    if(ci != cc) ok = false;
    // FOR(i, 0, H + 2) cout << CHECK[i] << endl;
    // FOR(i, 0, H) cout << O[i] << endl;
    
    // 5. 出力 ～ 後処理.
    if(ok) {
        cout << "possible" << endl;
        FOR(i, 0, H) cout << O[i] << endl;
    } else {
        cout << "impossible" << endl;
    }
    return 0;
    
}

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

// AtCoder Beginner Contest 075 (B - Minesweeper) を参照.

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

// 指定地点(i, j) の周囲に配置されている'#', '$' の個数(上下左右斜め方向)をカウント.

// @param i: (行方向)i番目.

// @param j: (列方向)j番目.

// @param a[]: マス目.

// @return: '#', '$' の個数.

int countSharpAndDollar(int i, int j, string a[]) {

const char SHARP_SIGN = '#';

const char DOLLAR_SIGN = '$';

int ret = 0;

if(a[i - 1][j - 1] == SHARP_SIGN || a[i - 1][j - 1] == DOLLAR_SIGN) ret++; // 左上.

if(a[i - 1][j] == SHARP_SIGN || a[i - 1][j] == DOLLAR_SIGN) ret++; // 上.

if(a[i - 1][j + 1] == SHARP_SIGN || a[i - 1][j + 1] == DOLLAR_SIGN) ret++; // 右上.

if(a[i][j - 1] == SHARP_SIGN || a[i][j - 1] == DOLLAR_SIGN) ret++; // 左.

if(a[i][j] == SHARP_SIGN || a[i][j] == DOLLAR_SIGN) ret++; // 中央.

if(a[i][j + 1] == SHARP_SIGN || a[i][j + 1] == DOLLAR_SIGN) ret++; // 右.

if(a[i + 1][j - 1] == SHARP_SIGN || a[i + 1][j - 1] == DOLLAR_SIGN) ret++; // 左下.

if(a[i + 1][j] == SHARP_SIGN || a[i + 1][j] == DOLLAR_SIGN) ret++; // 下.

if(a[i + 1][j + 1] == SHARP_SIGN || a[i + 1][j + 1] == DOLLAR_SIGN) ret++; // 右下.

return ret;

}

// 指定地点(i, j) の周囲に配置されている'#' の個数(上下左右斜め方向)をカウント.

// @param i: (行方向)i番目.

// @param j: (列方向)j番目.

// @param a[]: マス目.

// @return: '#' の個数.

int countSharp(int i, int j, string a[]) {

const char SHARP_SIGN = '#';

int ret = 0;

if(a[i - 1][j - 1] == SHARP_SIGN) ret++; // 左上.

if(a[i - 1][j] == SHARP_SIGN) ret++; // 上.

if(a[i - 1][j + 1] == SHARP_SIGN) ret++; // 右上.

if(a[i][j - 1] == SHARP_SIGN) ret++; // 左.

if(a[i][j + 1] == SHARP_SIGN) ret++; // 右.

if(a[i + 1][j - 1] == SHARP_SIGN) ret++; // 左下.

if(a[i + 1][j] == SHARP_SIGN) ret++; // 下.

if(a[i + 1][j + 1] == SHARP_SIGN) ret++; // 右下.

return ret;

}

int main() {

// 1. 入力情報取得.

int H, W;

cin >> H >> W;

// 調査対象のマス目の外側を, '$' で囲む.

// 例)

// [入力値]

// .....

// .#.#.

// .....

// ↓

// ['$'で囲んだ状態]

// $$$$$$$

// $.....$

// $.#.#.$

// $.....$

// $$$$$$$

const char DOLLAR_SIGN = '$';

const char SHARP_SIGN = '#';

const char A_SIGN = 'A';

string a;

FOR(i, 0, W + 2) a += DOLLAR_SIGN;

string I[H + 2] = {};

I[0] = a;

I[H + 1] = a;

FOR(i, 1, H + 1) {

string s;

cin >> s;

I[i] = (DOLLAR_SIGN + s + DOLLAR_SIGN);

}

// 2. マス目を一つ一つ調べ, '#', '$' の個数を数える.

// 出力用のマス目は, シャープ記号で初期化しておく.

string b;

FOR(i, 0, W) b += SHARP_SIGN;

string O[H] = {};

FOR(i, 0, H) O[i] = b;

// 3. 上下左右斜めの八方向でカウント.

FOR(i, 1, H + 1) FOR(j, 1, W + 1) if(countSharpAndDollar(i, j, I) != 9) O[i - 1][j - 1] = '.';

// 4. 出力用のマス目から, 再度収縮を行って, チェック用のマス目を作成し,

// 入力データと比較.

string CHECK[H + 2] = {};

CHECK[0] = a;

CHECK[H + 1] = a;

// 4-1. チェック用のマス目作成.

FOR(i, 1, H + 1) CHECK[i] = (DOLLAR_SIGN + O[i - 1] + DOLLAR_SIGN);

FOR(i, 1, H + 1) {

FOR(j, 1, W + 1) {

// cout << CHECK[i][j] << " " << countSharp(i, j, CHECK) << endl;

if(CHECK[i][j] == SHARP_SIGN || CHECK[i][j] == A_SIGN) continue;

if(countSharp(i, j, CHECK) > 0 && CHECK[i][j] != A_SIGN) CHECK[i][j] = A_SIGN;

}

// 4-2. 入力時のマス目(I) とチェック用のマス目(CHECK)について,

// I について '#' の個数, CHECK について, '#' と 'A' の合計個数を比較して,

// 同じだったら, "possible", 異なったら, "impossible" と判定.

bool ok = true;

int ci = 0;

int cc = 0;

FOR(i, 1, H + 1) FOR(j, 1, W + 1) if(I[i][j] == SHARP_SIGN) ci++;

FOR(i, 1, H + 1) FOR(j, 1, W + 1) if(CHECK[i][j] == SHARP_SIGN || CHECK[i][j] == A_SIGN) cc++;

if(ci != cc) ok = false;

// FOR(i, 0, H + 2) cout << CHECK[i] << endl;

// FOR(i, 0, H) cout << O[i] << endl;

// 5. 出力～後処理.

if(ok) {

cout << "possible" << endl;

FOR(i, 0, H) cout << O[i] << endl;

} else {

cout << "impossible" << endl;

}

return 0;

}

■参照サイト
AtCoder Beginner Contest 039

2018年12月4日2018年12月4日

Docker の動作確認をしてみた（２）

環境は、Window 10 Home (64bit) 上で行った。

※1. 基本的に, 参照URL① の Tutorial に従って確認した(※ほぼ, Tutorial通りのため, 手順は割愛している).
※2. [Dockerfile]ファイルには, FROM python:3.7-slim を指定した.
※3. ブラウザ確認時の IPアドレスが, PC(※ここでは, 192.168.56.1)でなくて,
docker instance の IPアドレス(※ここでは, 192.168.99.100) となることに注意が必要とのこと.

■結果.

1. Run the app
practice $ docker run -p 4000:80 friendlyhello
 * Serving Flask app "app" (lazy loading)
 * Environment: production
   WARNING: Do not use the development server in a production environment.
   Use a production WSGI server instead.
 * Debug mode: off
 * Running on http://0.0.0.0:80/ (Press CTRL+C to quit)

2. check web browser.
①対象URL: http://localhost:4000/
②対象URL: http://192.168.56.1:4000/
-> このサイトにアクセスできません との error メッセージが表示された.

3. check curl command.
practice $ curl http://localhost:4000
curl: (7) Failed to connect to localhost port 4000: Connection refused

4. 参照URL② を 参考に, docker instance の IPアドレス を チェック.
practice $ docker-machine ip default
192.168.99.100
practice $ curl http://192.168.99.100:4000
<h3>Hello World!</h3><b>Hostname:</b> 4cb6191f7590<br/><b>Visits:</b> <i>cannot connect to Redis, counter disabled</i>practice $

1. Run the app

practice $ docker run -p 4000:80 friendlyhello

* Serving Flask app "app" (lazy loading)

* Environment: production

WARNING: Do not use the development server in a production environment.

Use a production WSGI server instead.

* Debug mode: off

* Running on http://0.0.0.0:80/ (Press CTRL+C to quit)

2. check web browser.

①対象URL: http://localhost:4000/

②対象URL: http://192.168.56.1:4000/

-> このサイトにアクセスできませんとの error メッセージが表示された.

3. check curl command.

practice $ curl http://localhost:4000

curl: (7) Failed to connect to localhost port 4000: Connection refused

4. 参照URL② を参考に, docker instance の IPアドレスをチェック.

practice $ docker-machine ip default

192.168.99.100

practice $ curl http://192.168.99.100:4000

<h3>Hello World!</h3><b>Hostname:</b> 4cb6191f7590<br/><b>Visits:</b> <i>cannot connect to Redis, counter disabled</i>practice $

5. check web browser again.
対象URL: http://192.168.99.100:4000/
-> Hello World! を含むメッセージが表示された.

5. check web browser again.

対象URL: http://192.168.99.100:4000/

-> Hello World! を含むメッセージが表示された.

6. stop docker container.
practice $ docker container stop friendlyhello
Error response from daemon: No such container: friendlyhello
practice $ docker container stop 4cb6191f7590
4cb6191f7590

7. check web browser.
対象URL: http://192.168.99.100:4000/
-> このサイトにアクセスできません との error メッセージが表示されたので, docker container の停止を確認できた.

6. stop docker container.

practice $ docker container stop friendlyhello

Error response from daemon: No such container: friendlyhello

practice $ docker container stop 4cb6191f7590

4cb6191f7590

7. check web browser.

対象URL: http://192.168.99.100:4000/

-> このサイトにアクセスできませんとの error メッセージが表示されたので, docker container の停止を確認できた.

■参照サイト
【参照URL①】Get Started, Part 2: Containers.
【参照URL②】Docker – Failed to connect to localhost port 4000: Connection refused.

2018年12月3日2018年12月3日

C++ の動作確認をしてみた（３７）

C++の練習を兼ねて, AtCoder Beginner Contest 114 の問題C（755) ～問題D (756) を解いてみた.

■感想.
1. 問題C, 問題D とも, 解答に至るまでの条件が複雑だったため, 結局, ゴリゴリ書く羽目になって, 非常にきつかった.
2. 但し, 解説と, ほぼ同じ方針だったので, 及第点は, 取れたように思う.

本家のサイトABC 114 解説をご覧下さい.

■C++版プログラム(問題C／AC版).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

// 文字'3', '5', '7' を使って, "3" ～ "777777777" までの数を生成.
// @param v:   前回生成した文字列の配列.
// @param l:   再帰処理の深さ(終了条件).
// @param d:   追加対象の文字列の長さ.
// @param ret: 生成した文字列の配列を返却.
// ex. l == 1 の場合, "33" ～ "777" の 36個の数字が生成される.
vector<string> make753(vector<string> v, int l, int d) {
    string base = "357";
    vector<string> ret = v;
    // Time Limit Exceeded が, テストケース b04 で発生したので, 修正.
    // -> N == 1 の コードテスト で, 
    // terminate called after throwing an instance of 'std::bad_alloc' what():  std::bad_alloc
    // を確認したので, l== 0 -> l <= 0 に修正.
    if(l <= 0) return ret;
    else       for(auto &p : v) FOR(i, 0, 3) if(p.size() == d) ret.push_back(p + base[i]);
    return make753(ret, --l, ++d);
}

int main() {
    
    // 1. 入力情報取得.
    string N;
    cin >> N;
    
    // 2. 各桁について, '3', '5', '7' と比較して, 場合の数を計算.
    int ans = 0;
    int l = N.size();
    
    // イメージとして, "3" ～ "777777777" までの数について,
    // N よりも小さい七五三数を抽出出来れば, 解答となるはず.
    // 但し, '3', '5', '7' が, それぞれ, 1回以上出現している必要がある.
    // i は, 比較対象とする七五三数(もどき, 333 なども比較していくため)の桁数とする.
    vector<string> v = {"33", "35", "37", "53", "55", "57", "73", "75", "77"};
    vector<string> cur = make753(v, l - 2, 2);
    
    // 3. 七五三数をカウントアップ.
    // for(auto &p : cur) cout << p << endl;
    for(auto &p : cur) {
        map<char, int> m;
        m['3'] = 0;
        m['5'] = 0;
        m['7'] = 0;
        FOR(i, 0, p.size()) m[p[i]]++;
        if(m['3'] > 0 && m['5'] > 0 && m['7'] > 0) {
            // 文字列長 が N よりも短い or 文字列の順序が, N より先の場合, カウントアップ.
            if(p.size() < N.size() || p <= N) {
                ans++;
                // cout << p << endl;
            }
        }
    }
    
    // 4. 後処理.
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

// 文字'3', '5', '7' を使って, "3" ～ "777777777" までの数を生成.

// @param v: 前回生成した文字列の配列.

// @param l: 再帰処理の深さ(終了条件).

// @param d: 追加対象の文字列の長さ.

// @param ret: 生成した文字列の配列を返却.

// ex. l == 1 の場合, "33" ～ "777" の 36個の数字が生成される.

vector<string> make753(vector<string> v, int l, int d) {

string base = "357";

vector<string> ret = v;

// Time Limit Exceeded が, テストケース b04 で発生したので, 修正.

// -> N == 1 のコードテストで,

// terminate called after throwing an instance of 'std::bad_alloc' what(): std::bad_alloc

// を確認したので, l== 0 -> l <= 0 に修正.

if(l <= 0) return ret;

else for(auto &p : v) FOR(i, 0, 3) if(p.size() == d) ret.push_back(p + base[i]);

return make753(ret, --l, ++d);

}

int main() {

// 1. 入力情報取得.

string N;

cin >> N;

// 2. 各桁について, '3', '5', '7' と比較して, 場合の数を計算.

int ans = 0;

int l = N.size();

// イメージとして, "3" ～ "777777777" までの数について,

// N よりも小さい七五三数を抽出出来れば, 解答となるはず.

// 但し, '3', '5', '7' が, それぞれ, 1回以上出現している必要がある.

// i は, 比較対象とする七五三数(もどき, 333 なども比較していくため)の桁数とする.

vector<string> v = {"33", "35", "37", "53", "55", "57", "73", "75", "77"};

vector<string> cur = make753(v, l - 2, 2);

// 3. 七五三数をカウントアップ.

// for(auto &p : cur) cout << p << endl;

for(auto &p : cur) {

map<char, int> m;

m['3'] = 0;

m['5'] = 0;

m['7'] = 0;

FOR(i, 0, p.size()) m[p[i]]++;

if(m['3'] > 0 && m['5'] > 0 && m['7'] > 0) {

// 文字列長が N よりも短い or 文字列の順序が, N より先の場合, カウントアップ.

if(p.size() < N.size() || p <= N) {

ans++;

// cout << p << endl;

}

// 4. 後処理.

cout << ans << endl;

return 0;

}

■C++版プログラム(問題D／AC版).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

// a を b で何回割り切れるかをカウント.
// @param a: 割られる数.
// @param b: 割る数.
// @param c: 割り切れた回数.
// @return: 割り切れる回数を返却.
int countDivide(int a, int b, int c){
    if(a % b == 0) return countDivide(a / b, b, ++c);
    return c;
}

int main() {
    
    // 1. 入力情報取得.
    int N;
    cin >> N;
    
    // 2. 約数をカウント.
    int prime[25] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97};
    
    // 各素数が何回出現するかをカウント.
    map<int, int> m;
    FOR(i, 1, N + 1) {
        int n = i;
        FOR(j, 0, 25) {
            if(i < prime[j]) break;
            int c = countDivide(i, prime[j], 0);
            // cout << i << " " << j << " " << c << endl;
            if(c != 0) m[prime[j]] += c;
        }
    }
    // for(auto &p : m) cout << p.first << " " << p.second << endl;
    
    // 3. 出現した素数の数に応じて, 分類.
    // 74回以上        現れた素数の種類数 … c74.
    // 24回以上74回未満現れた素数の種類数 … c14.
    // 14回以上24回未満現れた素数の種類数 … c14.
    //  4回以上14回未満現れた素数の種類数 …  c4.
    //  2回以上 4回未満現れた素数の種類数 …  c2.
    // を計算.
    int c2  = 0;
    int c4  = 0;
    int c14 = 0;
    int c24 = 0;
    int c74 = 0;
    for(auto &p : m) {
        int t = p.second;
        if(t >= 74)           c74++;
        if(t >= 24 && t < 74) c24++;
        if(t >= 14 && t < 24) c14++;
        if(t >= 4  && t < 14) c4++;
        if(t >= 2  && t <  4) c2++;
    }
    
    // 4. 七五数を計算.
    int ans = 0;
    int total = c2 + c4 + c14 + c24 + c74;
    
    // 4-1. c74 だけで七五数を構築.
    // 75 * 1 = 75.
    ans += c74;
    
    // 4-2. c24, と もう一つ で七五数を構築など.
    // 25 * 3 = 75.
    ans += (total - c2 - c4 - c14) * (total - 1);
    
    // 4-3. c14, と もう一つ で七五数を構築など.
    // 15 * 5 = 75.
    ans += (total - c2 - c4) * (total - c2 - 1);
    
    // 4-4. c2, と もう二つ, もしくは, c2 以外 で七五数を構築.
    // 3 * 5 * 5 = 75.
    ans += (total - c2) * (total - c2 - 1) * c2 / 2;
    ans += (total - c2) * (total - c2 - 1) * (total - c2 - 2) / 2;

    // 5. 出力 ～ 後処理.
    cout << ans << endl;
    return 0;

}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

// a を b で何回割り切れるかをカウント.

// @param a: 割られる数.

// @param b: 割る数.

// @param c: 割り切れた回数.

// @return: 割り切れる回数を返却.

int countDivide(int a, int b, int c){

if(a % b == 0) return countDivide(a / b, b, ++c);

return c;

}

int main() {

// 1. 入力情報取得.

int N;

cin >> N;

// 2. 約数をカウント.

int prime[25] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97};

// 各素数が何回出現するかをカウント.

map<int, int> m;

FOR(i, 1, N + 1) {

int n = i;

FOR(j, 0, 25) {

if(i < prime[j]) break;

int c = countDivide(i, prime[j], 0);

// cout << i << " " << j << " " << c << endl;

if(c != 0) m[prime[j]] += c;

}

// for(auto &p : m) cout << p.first << " " << p.second << endl;

// 3. 出現した素数の数に応じて, 分類.

// 74回以上現れた素数の種類数 … c74.

// 24回以上74回未満現れた素数の種類数 … c14.

// 14回以上24回未満現れた素数の種類数 … c14.

// 4回以上14回未満現れた素数の種類数 … c4.

// 2回以上 4回未満現れた素数の種類数 … c2.

// を計算.

int c2 = 0;

int c4 = 0;

int c14 = 0;

int c24 = 0;

int c74 = 0;

for(auto &p : m) {

int t = p.second;

if(t >= 74) c74++;

if(t >= 24 && t < 74) c24++;

if(t >= 14 && t < 24) c14++;

if(t >= 4 && t < 14) c4++;

if(t >= 2 && t < 4) c2++;

}

// 4. 七五数を計算.

int ans = 0;

int total = c2 + c4 + c14 + c24 + c74;

// 4-1. c74 だけで七五数を構築.

// 75 * 1 = 75.

ans += c74;

// 4-2. c24, ともう一つで七五数を構築など.

// 25 * 3 = 75.

ans += (total - c2 - c4 - c14) * (total - 1);

// 4-3. c14, ともう一つで七五数を構築など.

// 15 * 5 = 75.

ans += (total - c2 - c4) * (total - c2 - 1);

// 4-4. c2, ともう二つ, もしくは, c2 以外で七五数を構築.

// 3 * 5 * 5 = 75.

ans += (total - c2) * (total - c2 - 1) * c2 / 2;

ans += (total - c2) * (total - c2 - 1) * (total - c2 - 2) / 2;

// 5. 出力～後処理.

cout << ans << endl;

return 0;

}

■参照サイト
AtCoder Beginner Contest 114

2018年12月2日2018年12月2日

Docker の動作確認をしてみた（１）

環境は、Window 10 Home (64bit) 上で行った。

※1. install時の注意点として, [Docker Quickstart Terminal]ショートカットのリンク先情報(下記, 参照)に, Git install が, 前提だったので, 再度, Git installする形で, 対応した.
※2. install手順は, 他サイトに, 詳細があるので, 割愛した.

[リンク先情報]
"C:\Program Files\Git\bin\bash.exe" --login -i "C:\Program Files\Docker Toolbox\start.sh"

1 2	[リンク先情報] "C:\Program Files\Git\bin\bash.exe" --login -i "C:\Program Files\Docker Toolbox\start.sh"

■結果
1. Docker Quickstart Terminal の起動.

2. Docker の動作確認.

> docker run hello-world

1	> docker run hello-world

■参照サイト
【参照URL①】Install Docker on Windows 10 Home.
【参照URL②】Docker Toolbox overview.
【参照URL③】Docker quick start treminal fails to start on Windows 10 with VBoxManage.exe error #473.
【参照URL④】Cannot boot Windows guest in VirtualBox without kernel module error.
【参照URL⑤】How do I figure out the display name of a Windows service from its “short” name?.
【参照URL⑥】“Docker Quick Start Terminal” shortcut, an error occurs. #335.

2018年11月30日

C++ の動作確認をしてみた（３６）

C++の練習を兼ねて, AtCoder Beginner Contest 045 の問題C（たくさんの数式 / Many Formulas) ～問題D (すぬけ君の塗り絵 / Snuke’s Coloring) を解いてみた.

■感想.
1. とりあえず, 解説見ずに解けたので良かったと思う.
2. 問題C は, 方針を決めることが出来たが, 結構ゴリゴリ書く羽目になってしまったように思う.
3. 問題D は, 計算量を減らすのに, 黒いマスに関係する3行3列の連続するマス目に限定させる方法に, 辿り着いたので解答出来たと思う.
4. 解説と, ほぼ同じ方針だったので, 及第点は, 取れたと思う.

本家のサイトABC 045 / ARC 061 解説をご覧下さい.

■C++版プログラム(問題C／AC版).

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {
    
    // 1. 入力情報取得.
    const int a[10] = {1, 2, 4, 8, 16, 32, 64, 128, 256, 512};
    string S;
    cin >> S;

    // 2. 考えられる数式を組み立てて, 総和を計算.
    // ex.
    // i == 25 -> 00011001 と読み替えて, 1 を '+' と見る.
    // この場合, 123456789 -> 1234 + 5 + 678 + 9 となる.
    // -> 上記の場合, 数式自体は, 2の8乗 == 256通り(00000000 ～ 11111111) と理解される.
    int l = S.size();
    int limit = a[l - 1];
    LL ans = 0;
    FOR(i, 0, limit) {
        // 2-1. i を 2進数表記に変換.
        stringstream ss;
        ss << bitset<9>(i);
        string s = ss.str();
        
        // 2-2. 2進数表記の桁数を, 文字列Sの長さ - 1 に揃える.
        // ex. 
        // l == 3 の場合, 000000011 -> 11 のように変換.
        s = s.substr(10 - l, l - 1);
        // cout << s << endl;
        
        // 2-3. 数式を組み立て, 保存.
        LL f = S[0] - '0';
        vector<LL> v;
        FOR(j, 0, s.size()) {
            if(s[j] == '0') f = 10 * f + S[j + 1] - '0';
            if(s[j] == '1') {
                v.push_back(f);
                f = S[j + 1] - '0';
            }
        }
        v.push_back(f);
        
        // 2-4. 数式の和を計算.
        LL t = 0;
        for(auto &p : v) t += p;
        // cout << t << endl;
        ans += t;
    }
    
    // 3. 出力 ～ 後処理.
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

const int a[10] = {1, 2, 4, 8, 16, 32, 64, 128, 256, 512};

string S;

cin >> S;

// 2. 考えられる数式を組み立てて, 総和を計算.

// ex.

// i == 25 -> 00011001 と読み替えて, 1 を '+' と見る.

// この場合, 123456789 -> 1234 + 5 + 678 + 9 となる.

// -> 上記の場合, 数式自体は, 2の8乗 == 256通り(00000000 ～ 11111111) と理解される.

int l = S.size();

int limit = a[l - 1];

LL ans = 0;

FOR(i, 0, limit) {

// 2-1. i を 2進数表記に変換.

stringstream ss;

ss << bitset<9>(i);

string s = ss.str();

// 2-2. 2進数表記の桁数を, 文字列Sの長さ - 1 に揃える.

// ex.

// l == 3 の場合, 000000011 -> 11 のように変換.

s = s.substr(10 - l, l - 1);

// cout << s << endl;

// 2-3. 数式を組み立て, 保存.

LL f = S[0] - '0';

vector<LL> v;

FOR(j, 0, s.size()) {

if(s[j] == '0') f = 10 * f + S[j + 1] - '0';

if(s[j] == '1') {

v.push_back(f);

f = S[j + 1] - '0';

}

v.push_back(f);

// 2-4. 数式の和を計算.

LL t = 0;

for(auto &p : v) t += p;

// cout << t << endl;

ans += t;

}

// 3. 出力～後処理.

cout << ans << endl;

return 0;

}

■C++版プログラム(問題D／AC版).

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

    // 1. 入力情報取得.
    LL H, W, N;
    cin >> H >> W >> N;
    // 3行3列の連続するマス目(以下, 領域X と略記)の中心座標を, 一次元に変換後, 
    // カウントアップし, その情報を保存.
    map<LL, LL> m;
    
    // 2. 黒いマスの座標(a, b) について, 
    // どの領域X から カウントアップされるかを確認.
    FOR(i, 0, N) {
        LL a, b;
        cin >> a >> b;
        FOR(r, -1, 2) {
            FOR(c, -1, 2) {
                // 黒いマスの座標(a, b) に対応する, 領域X の 中心座標(ca, cb) を 取得.
                LL ca = a + r;
                LL cb = b + c;
                // 中心座標(ca, cb) を 1次元に変換.
                if(ca >= 2 && cb >= 2 && ca <= H - 1 && cb <= W - 1) {
                    // 対象となる 領域X について, 中心座標を, 一次元に変換後, カウントアップ.
                    LL cIndex = cb + (ca - 1) * W;
                    m[cIndex]++;
                }
            }
        }
    }
    
    // 3. 上記で収集した, 領域X の 中心座標 及び, カウントアップに関する情報を集計.
    // -> ans[i] に, 黒いマスが, i個あるものを保管.
    LL ans[10] = {};
    LL total = 0;
    for(auto &p : m) {
        total++;
        ans[p.second]++;
    }
    ans[0] = (H - 2) * (W - 2) - total;

    // 4. 出力 ～ 後処理.
    FOR(i, 0, 10) cout << ans[i] << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

LL H, W, N;

cin >> H >> W >> N;

// 3行3列の連続するマス目(以下, 領域X と略記)の中心座標を, 一次元に変換後,

// カウントアップし, その情報を保存.

map<LL, LL> m;

// 2. 黒いマスの座標(a, b) について,

// どの領域X からカウントアップされるかを確認.

FOR(i, 0, N) {

LL a, b;

cin >> a >> b;

FOR(r, -1, 2) {

FOR(c, -1, 2) {

// 黒いマスの座標(a, b) に対応する, 領域X の中心座標(ca, cb) を取得.

LL ca = a + r;

LL cb = b + c;

// 中心座標(ca, cb) を 1次元に変換.

if(ca >= 2 && cb >= 2 && ca <= H - 1 && cb <= W - 1) {

// 対象となる領域X について, 中心座標を, 一次元に変換後, カウントアップ.

LL cIndex = cb + (ca - 1) * W;

m[cIndex]++;

}

// 3. 上記で収集した, 領域X の中心座標及び, カウントアップに関する情報を集計.

// -> ans[i] に, 黒いマスが, i個あるものを保管.

LL ans[10] = {};

LL total = 0;

for(auto &p : m) {

total++;

ans[p.second]++;

}

ans[0] = (H - 2) * (W - 2) - total;

// 4. 出力～後処理.

FOR(i, 0, 10) cout << ans[i] << endl;

return 0;

}

■参照サイト
AtCoder Beginner Contest 045

2018年11月29日

C++ の動作確認をしてみた（３５）

C++の練習を兼ねて, AtCoder Beginner Contest 043 の問題D（アンバランス / Unbalanced) を解いてみた.

■感想.
1. とりあえず, 解説見ずに解けたので良かったと思う.
2. 問題D は, 過半数の解釈に時間かかったように思う, ソース上のコメント通りであるが,
2パターンに集約できることに気付いたので, 解答に到達できた.
※なお, WA版は, パターン① しか考慮出来てなかったことが原因と思われる.

本家のサイトABC043 / ARC059 解説をご覧下さい.

■C++版プログラム(問題D／WA版, 1_08, 1_09, 1_10, 2_18, 2_19, 2_20 で, Wrong Answer).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

    // 1. 入力情報取得.
    string S;
    cin >> S;
    
    // 2. アンバランスな (連続する) 部分文字列が存在するか判定.
    // 1_08, 1_09, 1_10, 2_18, 2_19, 2_20 で, error.
    // -> 以下の反例がある.
    // 文字列: 'aba' は, アンバランスな筈.
    int l = S.size();
    char bef = S[0];
    int ans[2] = {-1, -1};
    FOR(i, 1, l) {
        char cur = S[i];
        if(bef == cur && ans[0] == -1) {
            ans[0] = i;
            ans[1] = i + 1;
            break;
        }
        bef = cur;
    }

    // 3. 出力 ～ 後処理.
    cout << ans[0] << " " << ans[1] << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

string S;

cin >> S;

// 2. アンバランスな (連続する) 部分文字列が存在するか判定.

// 1_08, 1_09, 1_10, 2_18, 2_19, 2_20 で, error.

// -> 以下の反例がある.

// 文字列: 'aba' は, アンバランスな筈.

int l = S.size();

char bef = S[0];

int ans[2] = {-1, -1};

FOR(i, 1, l) {

char cur = S[i];

if(bef == cur && ans[0] == -1) {

ans[0] = i;

ans[1] = i + 1;

break;

}

bef = cur;

}

// 3. 出力～後処理.

cout << ans[0] << " " << ans[1] << endl;

return 0;

}

■C++版プログラム(問題D／AC版).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

    // 1. 入力情報取得.
    string S;
    cin >> S;
    
    // 2. アンバランスな (連続する) 部分文字列が存在するか判定.
    // 1_08, 1_09, 1_10, 2_18, 2_19, 2_20 で, error.
    // -> 以下の反例がある.
    // 文字列: 'aba' は, アンバランスなはず.
    // 
    // ex.
    // アンバランスな文字列とは？
    // パターン①: 'aa'    -> 'a' が,  100% なので, 過半数と分かる.
    // パターン②: 'aba'   -> 'a' が, 約66% なので, 過半数と分かる.
    // パターン③: 'abbaa' -> 'a' が,   60% なので, 過半数と分かる.
    // -> 但し, パターン③ は, 文字列 'aa' を抽出できるので, パターン① に集約されるはず.
    int l = S.size();
    char bef_2 = S[0]; // 2文字前.
    char bef_1 = S[1]; // 1文字前.
    int ans[2] = {-1, -1};
    
    // 2-1. 1, 2文字目をチェックしておく.
    if(S[0] == S[1]) {
        ans[0] = 1;
        ans[1] = 2;
    }
    
    // 2-2. 3文字目以降もチェックする.
    FOR(i, 2, l) {
        char cur = S[i];
        if(ans[0] != -1) break;
        // 文字が連続しているかチェック(パターン①).
        if(bef_1 == cur && ans[0] == -1) {
            ans[0] = i;
            ans[1] = i + 1;
            break;
        }
        // 文字が1文字空けて再度出現しているかチェック(パターン②).
        if(bef_2 == cur && ans[0] == -1) {
            ans[0] = i - 1;
            ans[1] = i + 1;
            break;
        }
        bef_2 = bef_1;
        bef_1 = cur;
    }

    // 3. 出力 ～ 後処理.
    cout << ans[0] << " " << ans[1] << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

string S;

cin >> S;

// 2. アンバランスな (連続する) 部分文字列が存在するか判定.

// 1_08, 1_09, 1_10, 2_18, 2_19, 2_20 で, error.

// -> 以下の反例がある.

// 文字列: 'aba' は, アンバランスなはず.

// ex.

// アンバランスな文字列とは？

// パターン①: 'aa' -> 'a' が, 100% なので, 過半数と分かる.

// パターン②: 'aba' -> 'a' が, 約66% なので, 過半数と分かる.

// パターン③: 'abbaa' -> 'a' が, 60% なので, 過半数と分かる.

// -> 但し, パターン③ は, 文字列 'aa' を抽出できるので, パターン① に集約されるはず.

int l = S.size();

char bef_2 = S[0]; // 2文字前.

char bef_1 = S[1]; // 1文字前.

int ans[2] = {-1, -1};

// 2-1. 1, 2文字目をチェックしておく.

if(S[0] == S[1]) {

ans[0] = 1;

ans[1] = 2;

}

// 2-2. 3文字目以降もチェックする.

FOR(i, 2, l) {

char cur = S[i];

if(ans[0] != -1) break;

// 文字が連続しているかチェック(パターン①).

if(bef_1 == cur && ans[0] == -1) {

ans[0] = i;

ans[1] = i + 1;

break;

}

// 文字が1文字空けて再度出現しているかチェック(パターン②).

if(bef_2 == cur && ans[0] == -1) {

ans[0] = i - 1;

ans[1] = i + 1;

break;

}

bef_2 = bef_1;

bef_1 = cur;

}

// 3. 出力～後処理.

cout << ans[0] << " " << ans[1] << endl;

return 0;

}

■参照サイト
AtCoder Beginner Contest 043

2018年11月27日2018年11月27日

PyTorch の動作確認をしてみた（２５）

1. 環境は、Window 10 Home (64bit) 上で行った。

2. Anaconda3 (64bit) – Spyder上で、動作確認を行った。

3. python のバージョンは、python 3.7.0 である。

4. pytorch のバージョンは、pytorch 0.4.1 である。

5. Flask のバージョンは、Flask 1.0.2 である。

6. GPU は, NVIDIA社の GeForce GTX 1050 である。

7. CPU は, Intel社の Core(TM) i7-7700HQ である。

今回確認した内容は、現場で使える! PyTorch開発入門深層学習モデルの作成とアプリケーションへの実装 (AI & TECHNOLOGY) の 7.2 Flaskを用いたWebAPI化 (P.175 – P.182) である。

※1. プログラムの詳細は, 書籍を参考(P.175 – P.182)にして下さい.

※2. Unix環境で無いため, 原書のような Gunicorn は使用せず, 代案となるが,
参照URL① にある HTTPクライアント(Advanced REST client)を使う形で, 動作確認を行った.
原書の確認方法から, かなり逸脱してしまったが, 学習済みモデルに, テスト画像を読み込ませて,
予測した情報をレスポンスさせる形で確認できたので, とりあえず良しとした.

※3. 下記, ソースのコメントにも記載したが, “loadするモデル” に注意が必要である,
具体的には, resnet18 以外のネットワークを使った学習モデルを保存（ここでは, taco_and_burrito_09.pth）した場合,
Missing key(s) in state_dict: に関する RuntimeError が生じる.
これを回避するため, resnet18 のネットワークを改変せず, そのまま taco画像, burrito画像を訓練したモデル
（ここでは, resnet18_19.pth）を保存し, このモデルを loadする形で, 動作確認を行った.

[モデル／load error]

[モデル／再保存]

[モデル／load ok]

※4. 下記, ソースのコメントにも記載したが, 最初, 変数 img の取得が出来なかったので, 参照URL② を参考に, Debug Mode(SET FLASK_ENV=development を指定する形)で, 動作確認を行った.

[KeyError: ‘img’]

[Advanced REST client／setting 1]

[Advanced REST client／setting 2]

■フォルダ構成.

■対象プログラム(原書を一部改変).

# -*- coding: utf-8 -*-
# 1. library import.
from flask import Flask, request, jsonify
from PIL import Image

# 2. create an application.
def create_app(classifier):
～(略)～
    def predict():
        # get handler of received file.
        # Debug Mode.
        # http://flask.pocoo.org/docs/1.0/quickstart/#debug-mode
        # SET FLASK_ENV=development
        # -> 以下の error を抽出できた.
        # raise exceptions.BadRequestKeyError(key)
        # werkzeug.exceptions.HTTPException.wrap.<locals>.newcls: 
        # 400 Bad Request: KeyError: 'img'
        # print("request:" + str(request.files)) # ImmutableMultiDict([])
        img_file = request.files["img"]
        # check if the file is empty.
        print("img_file.filename: " + str(img_file.filename) + " => OK!")
～(略)～

# -*- coding: utf-8 -*-

# 1. library import.

from flask import Flask, request, jsonify

from PIL import Image

# 2. create an application.

def create_app(classifier):

～(略)～

def predict():

# get handler of received file.

# Debug Mode.

# http://flask.pocoo.org/docs/1.0/quickstart/#debug-mode

# SET FLASK_ENV=development

# -> 以下の error を抽出できた.

# raise exceptions.BadRequestKeyError(key)

# werkzeug.exceptions.HTTPException.wrap.<locals>.newcls:

# 400 Bad Request: KeyError: 'img'

# print("request:" + str(request.files)) # ImmutableMultiDict([])

img_file = request.files["img"]

# check if the file is empty.

print("img_file.filename: " + str(img_file.filename) + " => OK!")

～(略)～

# -*- coding: utf-8 -*-
# 1. library import.
from torch import nn
from torchvision import transforms, models

# 2. create network.
def create_network():
～(略)～

# 3. declare a class.
class Classifier(object):
～(略)～

# -*- coding: utf-8 -*-

# 1. library import.

from torch import nn

from torchvision import transforms, models

# 2. create network.

def create_network():

～(略)～

# 3. declare a class.

class Classifier(object):

～(略)～

# -*- coding: utf-8 -*-
# 1. library import.
import os, torch
from flaskTest import create_app
from flaskTest.classifier import Classifier

# 2. set path.
folder_path = os.path.expanduser('~')
folder_path = folder_path + '\\.spyder-py3\\pytorch\\flaskTest\\'
# pth_file = folder_path + "taco_and_burrito_09.pth"
pth_file = folder_path + "resnet18_19.pth"

# 3. load parameters.
# RuntimeError: Error(s) in loading state_dict for ResNet:
# Missing key(s) in state_dict: "conv1.weight", "bn1.weight", ...
# "0.11.running_var", "0.11.num_batches_tracked", "1.weight", "1.bias".
# -> このエラーは, load する pthファイルの不具合と推測される.
# -> resnet18 以外のネットワークを使った学習モデルを保存したことに起因する, パラメータ不具合と推測.
# -> resnet18 のネットワークを改変しないで, そのまま taco画像, burrito画像を訓練したモデルを保存して, 再確認.
～(略)～

# -*- coding: utf-8 -*-

# 1. library import.

import os, torch

from flaskTest import create_app

from flaskTest.classifier import Classifier

# 2. set path.

folder_path = os.path.expanduser('~')

folder_path = folder_path + '\\.spyder-py3\\pytorch\\flaskTest\\'

# pth_file = folder_path + "taco_and_burrito_09.pth"

pth_file = folder_path + "resnet18_19.pth"

# 3. load parameters.

# RuntimeError: Error(s) in loading state_dict for ResNet:

# Missing key(s) in state_dict: "conv1.weight", "bn1.weight", ...

# "0.11.running_var", "0.11.num_batches_tracked", "1.weight", "1.bias".

# -> このエラーは, load する pthファイルの不具合と推測される.

# -> resnet18 以外のネットワークを使った学習モデルを保存したことに起因する, パラメータ不具合と推測.

# -> resnet18 のネットワークを改変しないで, そのまま taco画像, burrito画像を訓練したモデルを保存して, 再確認.

～(略)～

■実行結果.
1. 訓練データ（taco）
画像：taco_000.jpg
予測：taco
⇒ 正解

2. 訓練データ（burrito）
画像：burrito_000.jpg
予測：burrito
⇒ 正解

3. テストデータ（taco）
画像：taco_380.jpg
予測：taco
⇒ 正解

4. テストデータ（burrito）
画像：burrito_383.jpg
予測：taco
⇒ 不正解

5. テストデータ（burrito）
画像：burrito_374.jpg
予測：burrito
⇒ 正解

■参照サイト
【参照URL①】Advanced REST client
【参照URL②】Debug Mode.

■参考書籍
現場で使える! PyTorch開発入門深層学習モデルの作成とアプリケーションへの実装 (AI & TECHNOLOGY)

2018年11月26日

PyTorch の動作確認をしてみた（２４）

1. 環境は、Window 10 Home (64bit) 上で行った。

2. Anaconda3 (64bit) – Spyder上で、動作確認を行った。

3. python のバージョンは、python 3.7.0 である。

4. pytorch のバージョンは、pytorch 0.4.1 である。

5. Flask のバージョンは、Flask 1.0.2 である。

6. GPU は, NVIDIA社の GeForce GTX 1050 である。

7. CPU は, Intel社の Core(TM) i7-7700HQ である。

今回確認した内容は、現場で使える! PyTorch開発入門深層学習モデルの作成とアプリケーションへの実装 (AI & TECHNOLOGY) の 7.2 Flaskを用いたWebAPI化 (P.174 – P.174) である。

※1. プログラムの詳細は, 書籍を参考(P.174 – P.174)にして下さい.
※2. 今回は, PyTorch というよりは, 寧ろ Flask の動作確認となったように思う.

■Flaskの動作確認.

1. [Anaconda Navigator] から [Environments] - [tf_env_3.7] の 右三角ボタンを押下し, [Open Terminal]押下.

(tf_env_3.7) C:XXXXXXXXXXXXXXXXXXXX>prompt $g 

> conda install -c anaconda flask
～(略)～
Downloading and Extracting Packages
zeromq-4.2.5         | 9.5 MB    | ################################################################ | 100%
flask-1.0.2          | 137 KB    | ################################################################ | 100%
lzo-2.10             | 154 KB    | ################################################################ | 100%
bzip2-1.0.6          | 145 KB    | ################################################################ | 100%
libtiff-4.0.9        | 829 KB    | ################################################################ | 100%
hdf5-1.10.2          | 34.8 MB   | ################################################################ | 100%
libsodium-1.0.16     | 585 KB    | ################################################################ | 100%
sqlite-3.25.3        | 937 KB    | ################################################################ | 100%
yaml-0.1.7           | 103 KB    | ################################################################ | 100%
libxslt-1.1.32       | 458 KB    | ################################################################ | 100%
snappy-1.1.7         | 82 KB     | ################################################################ | 100%
icu-58.2             | 21.9 MB   | ################################################################ | 100%
curl-7.61.0          | 120 KB    | ################################################################ | 100%
zlib-1.2.11          | 128 KB    | ################################################################ | 100%
libpng-1.6.35        | 589 KB    | ################################################################ | 100%
jpeg-9b              | 313 KB    | ################################################################ | 100%
certifi-2018.10.15   | 138 KB    | ################################################################ | 100%
tk-8.6.8             | 3.8 MB    | ################################################################ | 100%
Preparing transaction: done
Verifying transaction: done
Executing transaction: done

> exit

1. [Anaconda Navigator] から [Environments] - [tf_env_3.7] の右三角ボタンを押下し, [Open Terminal]押下.

(tf_env_3.7) C:XXXXXXXXXXXXXXXXXXXX>prompt $g

> conda install -c anaconda flask

～(略)～

Downloading and Extracting Packages

zeromq-4.2.5 | 9.5 MB | ################################################################ | 100%

flask-1.0.2 | 137 KB | ################################################################ | 100%

lzo-2.10 | 154 KB | ################################################################ | 100%

bzip2-1.0.6 | 145 KB | ################################################################ | 100%

libtiff-4.0.9 | 829 KB | ################################################################ | 100%

hdf5-1.10.2 | 34.8 MB | ################################################################ | 100%

libsodium-1.0.16 | 585 KB | ################################################################ | 100%

sqlite-3.25.3 | 937 KB | ################################################################ | 100%

yaml-0.1.7 | 103 KB | ################################################################ | 100%

libxslt-1.1.32 | 458 KB | ################################################################ | 100%

snappy-1.1.7 | 82 KB | ################################################################ | 100%

icu-58.2 | 21.9 MB | ################################################################ | 100%

curl-7.61.0 | 120 KB | ################################################################ | 100%

zlib-1.2.11 | 128 KB | ################################################################ | 100%

libpng-1.6.35 | 589 KB | ################################################################ | 100%

jpeg-9b | 313 KB | ################################################################ | 100%

certifi-2018.10.15 | 138 KB | ################################################################ | 100%

tk-8.6.8 | 3.8 MB | ################################################################ | 100%

Preparing transaction: done

Verifying transaction: done

Executing transaction: done

> exit

2. 以下のソースを記載し, 保存.
[.spyder-py3]フォルダ内 の PyTorch_7-1_01.py を 編集(※UTF-8保存).

[編集内容]
# -*- coding: utf-8 -*-
# 1. library import.
import flask
from flask import Flask

# 2. set app.
app = Flask(__name__)

@app.route('/')
def hello():
    return str('Hello, Flask ver:') + str(flask.__version__) + str('!')

2. 以下のソースを記載し, 保存.

[.spyder-py3]フォルダ内の PyTorch_7-1_01.py を編集(※UTF-8保存).

[編集内容]

# -*- coding: utf-8 -*-

# 1. library import.

import flask

from flask import Flask

# 2. set app.

app = Flask(__name__)

@app.route('/')

def hello():

return str('Hello, Flask ver:') + str(flask.__version__) + str('!')

3. コマンドプロンプトから, 下記コマンドを実行.
※参照サイト③ を参照した.

> cd
C:XXXXXXXXXXXXXXXXXXXX\.spyder-py3

> FLASK_APP=PyTorch_7-1_01.py flask run
'FLASK_APP' は、内部コマンドまたは外部コマンド、
操作可能なプログラムまたはバッチ ファイルとして認識されていません。

> SET FLASK_APP=PyTorch_7-1_01.py

> flask run
 * Serving Flask app "PyTorch_7-1_01.py"
 * Environment: production
   WARNING: Do not use the development server in a production environment.
   Use a production WSGI server instead.
 * Debug mode: off
 * Running on http://127.0.0.1:5000/ (Press CTRL+C to quit)

3. コマンドプロンプトから, 下記コマンドを実行.

※参照サイト③ を参照した.

> cd

C:XXXXXXXXXXXXXXXXXXXX\.spyder-py3

> FLASK_APP=PyTorch_7-1_01.py flask run

'FLASK_APP' は、内部コマンドまたは外部コマンド、

操作可能なプログラムまたはバッチファイルとして認識されていません。

> SET FLASK_APP=PyTorch_7-1_01.py

> flask run

* Serving Flask app "PyTorch_7-1_01.py"

* Environment: production

WARNING: Do not use the development server in a production environment.

Use a production WSGI server instead.

* Debug mode: off

* Running on http://127.0.0.1:5000/ (Press CTRL+C to quit)

4. サンプルページの表示.
対象URL: http://127.0.0.1:5000/

1 2	4. サンプルページの表示. 対象URL: http://127.0.0.1:5000/

■参照サイト
【参照URL①】flask 1.0.2
【参照URL②】Flask 1.0.2
【参照URL③】Env. Variables not set while running Minimal Flask application

■参考書籍
現場で使える! PyTorch開発入門深層学習モデルの作成とアプリケーションへの実装 (AI & TECHNOLOGY)

2018年11月24日

C++ の動作確認をしてみた（３４）

C++の練習を兼ねて, AtCoder Beginner Contest 047 の問題B（すぬけ君の塗り絵 2 イージー / Snuke’s Coloring 2 (ABC Edit)) ～問題D（高橋君と見えざる手 / An Invisible Hand) を解いてみた.

■感想.
1. とりあえず, 解説見ずに解けたので良かったと思う.
2. 問題B は, 最初, 方針が全く見えなかったが, 数列 a の値を確認する度に, 黒く塗りつぶす x, y座標の最小値, 最大値が変化することに気付いたので, 解答に到達できた(※個人的には, 非常に面白い問題と感じた).
2. 問題D は, 数列A の各要素の後方について, 最大値の抽出が必要があることに気付けた(※下記, ソースでは, 数列B として抽出している)ため, 解答に到達できた.

本家のサイトARC #063 / ABC #047 解説をご覧下さい.

■C++版プログラム(問題B／AC版).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {
    
    // 1. 入力情報取得.
    int W, H, N;
    cin >> W >> H >> N;

    // 2. 長方形を黒く塗りつぶしていく.
    int xMin = 0;
    int xMax = W;
    int yMin = 0;
    int yMax = H;
    
    FOR(i, 0, N) {
        int x, y, a;
        cin >> x >> y >> a;
        if(a == 1) xMin = max(x, xMin);
        if(a == 2) xMax = min(x, xMax);
        if(a == 3) yMin = max(y, yMin);
        if(a == 4) yMax = min(y, yMax);
    }

    // 3. 後処理.
    int ans = (xMax - xMin) * (yMax - yMin);
    if(xMax <= xMin || yMax <= yMin) ans = 0;
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

int W, H, N;

cin >> W >> H >> N;

// 2. 長方形を黒く塗りつぶしていく.

int xMin = 0;

int xMax = W;

int yMin = 0;

int yMax = H;

FOR(i, 0, N) {

int x, y, a;

cin >> x >> y >> a;

if(a == 1) xMin = max(x, xMin);

if(a == 2) xMax = min(x, xMax);

if(a == 3) yMin = max(y, yMin);

if(a == 4) yMax = min(y, yMax);

}

// 3. 後処理.

int ans = (xMax - xMin) * (yMax - yMin);

if(xMax <= xMin || yMax <= yMin) ans = 0;

cout << ans << endl;

return 0;

}

■C++版プログラム(問題C／AC版).

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {
    
    // 1. 入力情報取得.
    string S;
    cin >> S;
    
    // 2. 必要な操作回数をカウントする.
    int ans = 0;
    char bef = S[0];    // 前回出現した文字
    FOR(i, 0, S.size()) {
        // 2-1. 現在の石をチェック.
        char cur = S[i];
        
        // 2-2. 前回と異なる石の色であれば, カウントアップ.
        if(cur != bef) ans++;
        
        // 2-3. 前回出現した文字を更新.
        bef = cur;
    }
    
    // 3. 出力 ～ 後処理.
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

string S;

cin >> S;

// 2. 必要な操作回数をカウントする.

int ans = 0;

char bef = S[0]; // 前回出現した文字

FOR(i, 0, S.size()) {

// 2-1. 現在の石をチェック.

char cur = S[i];

// 2-2. 前回と異なる石の色であれば, カウントアップ.

if(cur != bef) ans++;

// 2-3. 前回出現した文字を更新.

bef = cur;

}

// 3. 出力～後処理.

cout << ans << endl;

return 0;

}

■C++版プログラム(問題D／AC版).

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(LL i = (a); i < (b); ++i)
#define FORR(i, a, b) for(LL i = (a); i >= (b); --i)

int main() {

    // 1. 入力情報取得.
    LL N, T;
    cin >> N >> T;
    LL A[N] = {};
    FOR(i, 0, N) cin >> A[i];
    
    // 2. 高橋君の収益を少なくとも, 1円下げるために必要な合計コストの最小値を計算.
    // i < j において, A[i] < A[j] となっている組み合わせに限定したとき, 
    // A[j] - A[i] が最大となるものの個数をカウントしていく問題と読み替えることができる.
    // ex.
    // 10 100
    // 7 10 4 5 9 3 6 8 2 1
    // -> A[2] = 4, A[4] = 9 と A[5] = 3, A[7] = 8 の 場合に, 最大5 となることが分かる.
    // 
    // 右側から最大値を順次確認.
    // A =  7 10 4 5 9 3 6 8 2 1
    // B = 10  9 9 9 8 8 8 2 1 1
    // -> 例えば, 4より右側の最大値は 9, 6より右側の最大値は, 8 などと読み取ることが出来る.
    LL B[N] = {};
    B[N - 1] = A[N - 1];
    FORR(i, N - 2, 0) {
        if(A[i + 1] > B[i + 1]) B[i] = A[i + 1];
        else                    B[i] = B[i + 1];
    }
    // FOR(i, 0, N) cout << B[i] << endl;
    
    // 3. 数列 A, B の差分を保管.
    // 上記の例で確認したように, A[2] = 4, A[4] = 9 と A[5] = 3, A[7] = 8 の 場合に, 最大5 となるが,
    // 最大値を取るパターンが, 合計2回であることを抽出する必要がある.
    map<LL, LL> m;
    FOR(i, 0, N) {
        LL diff = max(0LL, B[i] - A[i]);
        m[diff]++;
    }

    // 4. 出力 ～ 後処理.
    // Display the last element in map : map iterator « map multimap « C++ Tutorial.
    // http://www.java2s.com/Tutorial/Cpp/0460__map-multimap/Displaythelastelementinmap.htm
    auto itr = m.end();
    --itr;
    LL ans = itr->second;
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(LL i = (a); i < (b); ++i)

#define FORR(i, a, b) for(LL i = (a); i >= (b); --i)

int main() {

// 1. 入力情報取得.

LL N, T;

cin >> N >> T;

LL A[N] = {};

FOR(i, 0, N) cin >> A[i];

// 2. 高橋君の収益を少なくとも, 1円下げるために必要な合計コストの最小値を計算.

// i < j において, A[i] < A[j] となっている組み合わせに限定したとき,

// A[j] - A[i] が最大となるものの個数をカウントしていく問題と読み替えることができる.

// ex.

// 10 100

// 7 10 4 5 9 3 6 8 2 1

// -> A[2] = 4, A[4] = 9 と A[5] = 3, A[7] = 8 の場合に, 最大5 となることが分かる.

// 右側から最大値を順次確認.

// A = 7 10 4 5 9 3 6 8 2 1

// B = 10 9 9 9 8 8 8 2 1 1

// -> 例えば, 4より右側の最大値は 9, 6より右側の最大値は, 8 などと読み取ることが出来る.

LL B[N] = {};

B[N - 1] = A[N - 1];

FORR(i, N - 2, 0) {

if(A[i + 1] > B[i + 1]) B[i] = A[i + 1];

else B[i] = B[i + 1];

}

// FOR(i, 0, N) cout << B[i] << endl;

// 3. 数列 A, B の差分を保管.

// 上記の例で確認したように, A[2] = 4, A[4] = 9 と A[5] = 3, A[7] = 8 の場合に, 最大5 となるが,

// 最大値を取るパターンが, 合計2回であることを抽出する必要がある.

map<LL, LL> m;

FOR(i, 0, N) {

LL diff = max(0LL, B[i] - A[i]);

m[diff]++;

}

// 4. 出力～後処理.

// Display the last element in map : map iterator « map multimap « C++ Tutorial.

// http://www.java2s.com/Tutorial/Cpp/0460__map-multimap/Displaythelastelementinmap.htm

auto itr = m.end();

--itr;

LL ans = itr->second;

cout << ans << endl;

return 0;

}

■参照サイト
AtCoder Beginner Contest 047

2018年11月22日

C++ の動作確認をしてみた（３３）

C++の練習を兼ねて, AtCoder Beginner Contest 048 の問題C（Boxes and Candies), 問題D（An Ordinary Game) を解いてみた.

■感想.
1. とりあえず, 解説見ずに解けたので良かったと思う.
2. 問題D は, 方針が見えなかったが, 文字列の両端に着目してみたところ, 規則性を抽出できたので, 解答に辿り着いた.
※解説みたら, だいたい同じ方針だったので, 着眼点については, 及第点に到達できたと思う.

本家のサイトABC 048 / ARC 064 解説をご覧下さい.

■C++版プログラム(問題C／AC版).

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(LL i = (a); i < (b); ++i)

int main() {
    
    // 1. 入力情報取得.
    LL N, X;
    cin >> N >> X;
    LL A[N] = {};
    FOR(i, 0, N) cin >> A[i];
    
    // 2. 必要な操作回数をカウントする.
    LL ans = 0;
    FOR(i, 0, N - 1) {
        // 2-1. キャンディを食べる個数.
        LL eat = max(0LL, A[i] + A[i + 1] - X);
        
        // 2-2. i番目の箱について, キャンディの個数を更新.
        A[i] -= max(0LL, eat - A[i + 1]);
        
        // 2-3. (i + 1)番目の箱について, キャンディの個数を更新.
        A[i + 1] = max(0LL, A[i + 1] - eat);
        
        // 2-4. 食べたキャンディの個数をカウント.
        ans += eat;
    }
    
    // 3. 出力 ～ 後処理.
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

#define FOR(i, a, b) for(LL i = (a); i < (b); ++i)

int main() {

// 1. 入力情報取得.

LL N, X;

cin >> N >> X;

LL A[N] = {};

FOR(i, 0, N) cin >> A[i];

// 2. 必要な操作回数をカウントする.

LL ans = 0;

FOR(i, 0, N - 1) {

// 2-1. キャンディを食べる個数.

LL eat = max(0LL, A[i] + A[i + 1] - X);

// 2-2. i番目の箱について, キャンディの個数を更新.

A[i] -= max(0LL, eat - A[i + 1]);

// 2-3. (i + 1)番目の箱について, キャンディの個数を更新.

A[i + 1] = max(0LL, A[i + 1] - eat);

// 2-4. 食べたキャンディの個数をカウント.

ans += eat;

}

// 3. 出力～後処理.

cout << ans << endl;

return 0;

}

■C++版プログラム(問題D／AC版).

#include <bits/stdc++.h>

using namespace std;

int main() {

    // 1. 入力情報取得.
    string S;
    cin >> S;
    
    // 2. 文字列の両端を確認する.
    string ans;
    int l = S.size();
    bool same = (S[0] == S[l - 1]);
    
    // 2-1. 両端が等しい場合.
    if(same) {
        if(l % 2 == 0) ans = "First";
        else           ans = "Second";
    }

    // 2-2. 両端が異なる場合.
    if(!same) {
        if(l % 2 == 0) ans = "Second";
        else           ans = "First";
    }

    // 3. 出力 ～ 後処理.
    cout << ans << endl;
    return 0;
    
}

#include <bits/stdc++.h>

using namespace std;

int main() {

// 1. 入力情報取得.

string S;

cin >> S;

// 2. 文字列の両端を確認する.

string ans;

int l = S.size();

bool same = (S[0] == S[l - 1]);

// 2-1. 両端が等しい場合.

if(same) {

if(l % 2 == 0) ans = "First";

else ans = "Second";

}

// 2-2. 両端が異なる場合.

if(!same) {

if(l % 2 == 0) ans = "Second";

else ans = "First";

}

// 3. 出力～後処理.

cout << ans << endl;

return 0;

}

■参照サイト
AtCoder Beginner Contest 048

C++の練習を兼ねて, AtCoder Beginner Contest 039 の 問題C（ピアニスト高橋君) ～ 問題D (画像処理高橋君) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 114 の 問題C（755) ～ 問題D (756) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 045 の 問題C（たくさんの数式 / Many Formulas) ～ 問題D (すぬけ君の塗り絵 / Snuke’s Coloring) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 043 の 問題D（アンバランス / Unbalanced) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 047 の 問題B（すぬけ君の塗り絵 2 イージー / Snuke’s Coloring 2 (ABC Edit)) ～ 問題D（高橋君と見えざる手 / An Invisible Hand) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 048 の 問題C（Boxes and Candies), 問題D（An Ordinary Game) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 039 の問題C（ピアニスト高橋君) ～問題D (画像処理高橋君) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 114 の問題C（755) ～問題D (756) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 045 の問題C（たくさんの数式 / Many Formulas) ～問題D (すぬけ君の塗り絵 / Snuke’s Coloring) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 043 の問題D（アンバランス / Unbalanced) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 047 の問題B（すぬけ君の塗り絵 2 イージー / Snuke’s Coloring 2 (ABC Edit)) ～問題D（高橋君と見えざる手 / An Invisible Hand) を解いてみた.

C++の練習を兼ねて, AtCoder Beginner Contest 048 の問題C（Boxes and Candies), 問題D（An Ordinary Game) を解いてみた.