結果

問題 No.1794 Continued Fraction
ユーザー tnakao0123tnakao0123
提出日時 2021-12-24 16:37:54
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 980 bytes
コンパイル時間 866 ms
コンパイル使用メモリ 42,092 KB
実行使用メモリ 5,156 KB
最終ジャッジ日時 2023-10-19 19:00:37
合計ジャッジ時間 9,556 ms
ジャッジサーバーID
(参考情報)
judge14 / judge15
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
権限があれば一括ダウンロードができます

ソースコード

diff #

/* -*- coding: utf-8 -*-
 *
 * 1794.cc:  No.1794 Continued Fraction - yukicoder
 */

#include<cstdio>
#include<algorithm>
 
using namespace std;

/* constant */

const int MAX_L = 62;

/* typedef */

/* global variables */

int as[MAX_L];

/* subroutines */

template <typename T>
T gcd(T m, T n) {  // m > 0, n > 0
  if (m < n) swap(m, n);
  while (n > 0) {
    T r = m % n;
    m = n;
    n = r;
  }
  return m;
}

void reduce(int &n, int &m) {
  if (n > 0 && m > 0) {
    int g = gcd(n, m);
    n /= g, m /= g;
  }
}

/* main */

int main() {
  /*
    Let xi=ai+1/x{i+1}, xL=aL
    -> n/m = x0 = a0+1/x1 -> 1/x1 = n/m-a0 = (n-m*a0)/m
    -> x1 = m/(n-m*a0)
   */

  int n, m;
  scanf("%d%d", &n, &m);
  reduce(n, m);

  int l = 0;
  while (n > 0) {
    int d = n / m;
    n -= m * d;
    as[l++] = d;
    reduce(n, m);
    swap(n, m);
  }

  while (l > 0 && as[l - 1] == 0) l--;

  for (int i = 0; i < l; i++)
    printf("%d%c", as[i], (i + 1 < l) ? ' ' : '\n');
  return 0;
}
0