東京工業大学プログラミングコンテスト2015 D: 文字列と素数

競プロに逃げすぎてレポートが進みません。

問題概要

原文 → D: 文字列と素数 - 東京工業大学プログラミングコンテスト2015 | AtCoder

文字列 が与えられる。各文字を、同じ文字は同じ数字・違う文字は違う数字になるように '1', '3', '5', '7', '9' のいずれかに変換する。文字列を素数に変換できるならばその素数を出力し、できないならば -1 を出力せよ。

解説

map と permutation を使って全探索するだけです。next_permutation とかそんなに使わないので備忘録的に書きました。

文字列を反転させておくことで計算が楽になります。文字列から数字にしたりとか、桁数がやたら多かったりとかするときはこれで幸せになることが多い気がします。

ソースコード

bool isprime(int n) {
    if(n == 1) return false;
    for(int i=2; i*i<=n; i++) {
        if(!(n%i)) return false;
    }
    return true;
}

signed main() {
    string s; cin >> s;
    map<char, int> m;
    int cnt = 0;
    reverse(s.begin(), s.end());
    rep(i,0,s.length()) {
        if(!m.count(s[i])) {
            m[s[i]] = cnt++;
        }
    }

    if(cnt > 5) {
        cout << -1 << endl;
        return 0;
    }
    int a[] = {1, 3, 5, 7, 9};
    do {
        int d = 1;
        int num = 0;
        rep(i,0,s.length()) {
            num += d * a[ m[s[i]] ];
            d *= 10;
        }
        if(isprime(num)) {
            cout << num << endl;
            return 0;
        }
    }while(next_permutation(a, a+5));
    cout << -1 << endl;
    return 0;
}

AOJ 0037: Path on a Grid

バグが取れなくて相当苦労した。

問題概要

原文 → 格子状の経路 | Aizu Online Judge

格子の各辺が壁であるかどうかの情報が与えられる。壁に右手をついたまま1周するときの経路を出力せよ。

解説

前回の記事 と同様の方法で解いてみました。

今見ている方向を とおくと、

• 右手に壁がない場合、 の方向に向き直す
• 進行方向に壁がある場合、 の方向に向き直す

というように動けばよいです。終了条件にも気をつけましょう。

ソースコード

#define S 6
// right, down, left, up
int dx[] = {0, 1, 0, -1};
int dy[] = {1, 0, -1, 0};

int main() {
    bool board[S][S][4] = {};
    rep(i,0,9) {
        int l, n;
        if(i % 2 == 0) l = 4;
        else l = 5;
        rep(j,0,l) {
            scanf("%1d", &n);
            if(n == 1) {
                if(i % 2 == 0) {
                    board[i/2][j+1][1] = true;
                    board[i/2+1][j+1][3] = true;
                }
                else {
                    board[i/2+1][j][0] = true;
                    board[i/2+1][j+1][2] = true;
                }
            }
        }
    }

    string d = "RDLU";
    string ret;
    int dir = 0;
    int nx = 0, ny = 1;
    for(int i=0; i<150; i++) {
        if(nx == 1 && ny == 0 && dir == 2) break;
        if(nx == 0 && ny == 0 && dir == 3) break;
        int x, y;
        x = nx + dx[dir], y = ny + dy[dir];
        if(board[nx][ny][(dir + 1) % 4] == 0) {
            dir = (dir + 1) % 4;
            nx = nx + dx[dir], ny = ny + dy[dir];
        }
        else if(board[nx][ny][dir] == 1) {
            ret += d[dir];
            dir = (dir - 1 + 4) % 4;
        }
        else {
            ret += d[dir];
            nx = x, ny = y;
        }
    }
    cout << ret << endl;
    return 0;
}

自分の実装の方法が悪いのか、バグ直しに相当時間がかかってしまった・・・。他の実装方法も考えたい。まずは enum 使って書きなおしたほうがいいよなあ。

ICPC 国内予選 2010 B: 迷図と命ず

あけましておめでとうございます。私はひたすら AOJ-ICPC を埋めています。

実装が面倒だなあこれ・・・と思ったのでほぼ自分用記事を書きます。

問題概要

原文 → Amazing Mazes | Aizu Online Judge

縦 、横 の迷路が与えられる。迷路の壁の情報が与えられるので、スタートからゴールまでの最短距離を求めよ。ゴールに辿りつけない迷路の場合は を出力せよ。

解説

解法としてはもちろん幅優先探索なんですが、今回は迷路の「壁」が与えられるため少し面倒です。いつもの方針と少し違う配列を持つ必要があります。

board[i][j][k] := 座標(i, j) において方向 k に動けるか という bool 配列を用意します。これさえ用意できればあとは普通に探索するだけです。が、バグが怖いですね。

ソースコード

int dx[]={1, -1, 0, 0};
int dy[]={0, 0, 1, -1};

signed main() {
    int h, w;
    while(cin >> w >> h, h || w) {
        bool board[40][40][4];
        rep(i,0,h) rep(j,0,w) rep(k,0,4) board[i][j][k] = true;

        rep(i,0,h) board[i][0][3]   = false; // left
        rep(i,0,h) board[i][w-1][2] = false; // right
        rep(i,0,w) board[0][i][1]   = false; // up
        rep(i,0,w) board[h-1][i][0] = false; // down

        int inp;
        rep(i,0,2*h-1) {
            int l;
            if(i % 2 == 0) l = w-1;
            else l = w;
            rep(j,0,l) {
                cin >> inp;
                if(inp == 1) {
                    if(i % 2 == 0) {
                        board[i/2][j][2] = false;
                        board[i/2][j+1][3] = false;
                    }
                    else {
                        board[i/2][j][0] = false;
                        board[i/2+1][j][1] = false;
                    }
                }
            }
        }

        int dist[40][40];
        rep(i,0,h) rep(j,0,w) dist[i][j] = INF;
        dist[0][0] = 1;
        queue<pii> q;
        q.push(pii(0, 0));

        while(!q.empty()) {
            pii t = q.front(); q.pop();
            rep(i,0,4) {
                if(!board[t.fr][t.sc][i]) continue;
                int x = t.fr + dx[i];
                int y = t.sc + dy[i];
                if(dist[x][y] <= dist[t.fr][t.sc] + 1) continue;
                dist[x][y] = dist[t.fr][t.sc] + 1;
                q.push(pii(x, y));
            }
        }

        if(dist[h-1][w-1] == INF) cout << 0 << endl;
        else cout << dist[h-1][w-1] << endl;
    }
    return 0;
}

なんか壁与えられる問題他にもあった気がする。