結果

問題 No.1794 Continued Fraction
ユーザー tnakao0123tnakao0123
提出日時 2021-12-24 16:41:08
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 989 bytes
コンパイル時間 512 ms
コンパイル使用メモリ 41,780 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-19 19:04:27
合計ジャッジ時間 1,939 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ（β）

テストケース

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

ソースコード

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 && m > 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