#include <iostream>
#include <complex>

float pi()
{
  return 3.141593f;
}

int main()
{
  // 実部1.0f、虚部2.0fの複素数オブジェクトを作る
  std::complex<float> c(1.0f, 2.0f);

  // ストリーム出力
  std::cout << "output : " << c << std::endl;

  // 各要素の取得
  float real = c.real(); // 実部
  float imag = c.imag(); // 虚部
  std::cout << "real : " << real << std::endl;
  std::cout << "imag : " << imag << std::endl;

  // 演算
  std::complex<float> a(1.0f, 2.0f);
  std::complex<float> b(2.0f, 3.0f);
  std::cout << "a + b : " << a + b << std::endl;
  std::cout << "a - b : " << a - b << std::endl;
  std::cout << "a * b : " << a * b << std::endl;
  std::cout << "a / b : " << a / b << std::endl;

  // 各複素数の値を取得する
  std::cout << "abs : " << std::abs(c) << std::endl;   // 絶対値
  std::cout << "arg : " << std::arg(c) << std::endl;   // 偏角
  std::cout << "norm : " << std::norm(c) << std::endl; // ノルム
  std::cout << "conj : " << std::conj(c) << std::endl; // 共役複素数
  std::cout << "proj : " << std::proj(c) << std::endl; // リーマン球面への射影
  std::cout << "polar : " << std::polar(1.0f, pi() / 2.0f); // 極座標(絶対値：1.0、偏角：円周率÷2.0)から複素数を作る
}

#include <print>
#include <complex>
#include <vector>
#include <ranges>

int main()
{
  // 構造化束縛で各要素を取得する
  {
    std::complex<float> c{1.0f, 2.0f};
    auto [real, imag] = c;
    std::println("real:{} imag:{}", real, imag);
  }

  // 定数の添字を指定して各要素を取得する
  {
    std::complex<float> c{1.0f, 2.0f};
    float real = std::get<0>(c);
    float imag = std::get<1>(c);
    std::println("real:{} imag:{}", real, imag);
  }

  // complexの配列から実部のみを抽出する
  {
    std::vector<std::complex<float>> v {
      {1.0f, 2.0f},
      {3.0f, 4.0f},
      {5.0f, 6.0f}
    };

    for (float real : v | std::views::elements<0>) {
      std::println("real:{}", real);
    }
  }
}