#pragma once

#include <vector>


struct Image {
  std::vector<unsigned char> img;
  int width;
  int height;
};
void resize(Image* img, int w, int h) {
  img->width = w;
  img->height = h;
  img->img.resize(w * h);
}

//! @brief 画像のピクセルを着色する
//! @param [in,out] img 画像
//! @param [in] x 座標
//! @param [in] y 座標
//! @return なし
void plot(Image& img, int x, int y) {
  img.img[y * img.width + x] = 255;
}


enum class CPOS {
  left, right, top, bottom
};

//! @brief 円の指定した1/4を描画する
//! @param [out] img 画像
//! @param [in] cx 中心座標
//! @param [in] cy 中心座標
//! @param [in] r  半径
//! @param [in] bx 円の四分割したエリアの指定
//! @param [in] by 円の四分割したエリアの指定
//! @return なし
void arc4(Image& img, const int cx, const int cy, const int r, const CPOS bx, const CPOS by) {
  int x, y;//描画位置
  const int rr = r * r;//半径の二乗
  int dx, dy;
  int limx;
  int limy;

  if (bx == CPOS::right && by == CPOS::bottom) {
    // □□
    // □■
    //初期位置
    x = r;
    y = 0;
    dx = -1;
    dy = +1;
    limx = 0;
    limy = r;
  }
  else if (bx == CPOS::left && by == CPOS::bottom) {
    // □□
    // ■□
    x = 0;
    y = r;
    dx = -1;
    dy = -1;
    limx = -r;
    limy = 0;
  }
  else if (bx == CPOS::left && by == CPOS::top) {
    // ■□
    // □□
    x = -r;
    y = 0;
    dx = +1;
    dy = -1;
    limx = 0;
    limy = -r;
  }
  else if (bx == CPOS::right && by == CPOS::top) {
    // □■
    // □□
    x = 0;
    y = -r;
    dx = +1;
    dy = +1;
    limx = r;
    limy = 0;
  }

  //最初の一点
  //plot(img,x + cx, y + cy);

  while ((x != limx) || (y != limy)) {

    /////////////////////////
    // 次に描画するポイントを決定

    //次に描画するのは以下のいずれか。
    int xy1 = pow(x + dx, 2) + pow(y, 2);
    int xy2 = pow(x + dx, 2) + pow(y + dy, 2);
    int xy3 = pow(x, 2) + pow(y + dy, 2);

    // x*x + y*y == r*rなので
    // 上記各xyのうち最もr*rとの差が少ないx*x+y*yを選択
    xy1 = abs(xy1 - rr);
    xy2 = abs(xy2 - rr);
    xy3 = abs(xy3 - rr);

    if (xy1 < xy2) {
      if (xy1 < xy3) {
        // xy1
        x = x + dx;
        y = y;
      }
      else {
        // xy2
        x = x + dx;
        y = y + dy;
      }
    }
    else {
      if (xy2 < xy3) {
        // xy2
        x = x + dx;
        y = y + dy;
      }
      else {
        // xy3
        x = x;
        y = y + dy;
      }
    }

    /////////////////////////
    // 描画
    if (
      (img, x + cx < 0) ||
      (img, x + cx >= img.width) ||
      (img, y + cy < 0) ||
      (img, y + cy >= img.height)
      ) {
      continue;
    }
    plot(img, x + cx, y + cy);


  }

}