結果

問題 No.124 門松列(3)
コンテスト
ユーザー te-sh
提出日時 2017-02-17 16:12:14
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 45 ms / 5,000 ms
コード長 2,548 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 750 ms
コンパイル使用メモリ 109,452 KB
実行使用メモリ 123,084 KB
最終ジャッジ日時 2024-06-12 07:06:27
合計ジャッジ時間 2,271 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.container; // SList, DList, BinaryHeap

void main()
{
  auto rd = readln.split.to!(int[]), w = rd[0], h = rd[1];
  auto mij = Matrix!int(h.iota.map!(_ => readln.split.to!(int[])).array);

  auto isKadomatsu(point p1, point p2, point p3) {
    auto a1 = mij[p1], a2 = mij[p2], a3 = mij[p3];
    return a1 != a3 && (a2 < a1 && a2 < a3 || a2 > a1 && a2 > a3);
  }

  auto calc() {
    auto vm = new VisitedManager(w, h);
    auto q = DList!qitem(qitem(point(0, 1), point(0, 0), 1),
                         qitem(point(1, 0), point(0, 0), 1));
    while (!q.empty) {
      auto qi = q.front; q.removeFront;

      foreach (sib; sibPoints) {
        auto np = qi.curr + sib;
        if (!mij.validIndex(np) || vm.visited(np, qi.curr) || !isKadomatsu(np, qi.curr, qi.prev)) continue;
        if (np == point(w - 1, h - 1)) return qi.len + 1;
        vm.visit(np, qi.curr);
        q.insertBack(qitem(np, qi.curr, qi.len + 1));
      }
    }
    return -1;
  }

  writeln(calc);
}

struct qitem {
  point curr, prev;
  int len;
}

class VisitedManager {
  bool[][] buf;
  int w, h;

  this(int w, int h) {
    this.w = w;
    this.h = h;
    this.buf = new bool[][](w * h, w * h);
  }

  int idx(point p) {
    return p.y * this.w + p.x;
  }

  void visit(point curr, point prev) {
    this.buf[idx(curr)][idx(prev)] = true;
  }

  bool visited(point curr, point prev) {
    return this.buf[idx(curr)][idx(prev)];
  }
}

struct Point(T) {
  T x, y;

  point opBinary(string op)(point rhs) {
    static if (op == "+") return point(x + rhs.x, y + rhs.y);
  }
}

alias Point!int point;

const auto sibPoints = [point(-1, 0), point(0, -1), point(1, 0), point(0, 1)];

struct Matrix(T)
{
  import std.algorithm, std.conv, std.range, std.traits, std.typecons;

  T[][] m;
  size_t rows, cols;

  mixin Proxy!m;

  this(size_t r, size_t c) { rows = r; cols = c; m = new T[][](rows, cols); }
  this(T[][] s) { rows = s.length; cols = s[0].length; m = s; }

  auto opIndex(U)(U p) { static if (is(U == Point!V, V)) return m[p.y][p.x]; else return m[p]; }
  auto opIndex(size_t y, size_t x) { return m[y][x]; }
  static if (isAssignable!T) {
    auto opIndexAssign(U)(T v, Point!U p) { return m[p.y][p.x] = v; }
    auto opIndexAssign(T v, size_t y, size_t x) { return m[y][x] = v; }
  }
  auto validIndex(U)(Point!U p) { return p.x >= 0 && p.x < cols && p.y >= 0 && p.y < rows; }
  auto points(U)() { return rows.to!U.iota.map!(y => cols.to!U.iota.map!(x => Point!U(x, y))).joiner; }
}
0