結果

問題 No.186 中華風 (Easy)
コンテスト
ユーザー moti
提出日時 2015-11-29 23:16:57
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,180 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 535 ms
コンパイル使用メモリ 63,804 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-14 05:19:44
合計ジャッジ時間 1,455 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 8 WA * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)

typedef unsigned long long ull;

ull extgcd(ull a, ull b, ull& x, ull& y) {
  ull d = a;
  if(b != 0) {
    d = extgcd(b, a % b, y, x);
    y -= (a / b) * x;
  }
  else {
    x = 1; y = 0;
  }
  return d;
}

ull mod_inverse(ull a, ull m) { ull x, y; extgcd(a, m, x, y); return (m + x % m) % m; }

namespace math {
template<class T>
pair<T, T> linear_congruence(vector<T> const& A, vector<T> const& B, vector<T> const& M) {
  T x = 0, m = 1;
  rep(i, A.size()) {
    T a = A[i] * m, b = B[i] - A[i] * x, d = __gcd(M[i], a);
    if(b % d != 0) { return {0, -1}; }
    T t = b / d * mod_inverse(a / d, M[i] / d) % (M[i] / d);
    x += m * t;
    m *= M[i] / d;
  }
  return {x % m, m};
}
}

int main() {

  vector<ull> A(3), B(3),M(3);
  rep(i, 3) {
    cin >> B[i] >> M[i];
    A[i] = 1;
  }
  
  auto ans = math::linear_congruence(A, B, M);
  if(ans.second == -1) {
    cout << -1 << endl;
    exit(0);
  }

  if(ans.first == 0) {
    cout << ans.second << endl;
  }
  else {
    cout << ans.first << endl;
  }

  return 0;
}
0