結果

問題 No.186 中華風 (Easy)
コンテスト
ユーザー te-sh
提出日時 2016-09-13 16:48:40
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,126 bytes
コンパイル時間 1,748 ms
コンパイル使用メモリ 171,900 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-12 04:20:05
合計ジャッジ時間 2,627 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.array, std.container, std.range, std.bitmanip;
import std.numeric, std.math, std.bigint, std.random, core.bitop;
import std.string, std.regex, std.conv, std.stdio, std.typecons;

T exEuclid(T)(T a, T b, ref T x, ref T y)
{
  auto g = a;
  x = 1;
  y = 0;
  if (b != 0) {
    g = exEuclid(b, a % b, y, x);
    y -= a / b * x;
  }
  return g;
}

void main()
{
  auto rd1 = readln.split.map!(a => BigInt(a));
  auto x1 = rd1[0], y1 = rd1[1];
  auto rd2 = readln.split.map!(a => BigInt(a));
  auto x2 = rd2[0], y2 = rd2[1];
  auto rd3 = readln.split.map!(a => BigInt(a));
  auto x3 = rd3[0], y3 = rd3[1];

  BigInt a, b, g, l, r;

  g = exEuclid(y1, y2, a, b);

  if (g > 1 && x1 % g != x2 % g) {
    writeln(-1);
    return;
  }

  l = y1 / g * y2;
  r = y1 / g * x2 * a + y2 / g * x1 * b;
  if (r < 0) r += ((-r / l) + 1) * l;
  r %= l;

  BigInt x4 = r, y4 = l;

  g = exEuclid(y3, y4, a, b);

  if (g > 1 && x3 % g != x4 % g) {
    writeln(-1);
    return;
  }

  l = y3 / g * y4;
  r = y3 / g * x4 * a + y4 / g * x3 * b;
  if (r < 0) r += ((-r / l) + 1) * l;
  r %= l;

  writeln(r == 0 ? l : r);
}
0