結果

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

ソースコード

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 long long ll;

ll gcd(ll a, ll b) {
  if(b == 0) { return a; }
  return gcd(b, a%b);
}

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

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

pair<ll, ll> linear_congruence(const vector<ll>& B, const vector<ll>& M)
{
  ll x = 0, m = 1;
  rep(i, B.size()) {
    ll a = m, b = B[i] - x, d = gcd(M[i],a);
    if(b % d != 0) { return std::make_pair(0,-1); }
    ll t = b / d * mod_inverse(a / d, M[i] / d) % (M[i] / d);
    x = x + m * t;
    m *= M[i] / d;
  }
  return std::make_pair(x % m, m);
}

int main() {

  vector<ll> B(3),M(3);
  rep(i, 3) {
    cin >> B[i] >> M[i];
  }
  
  auto ans = linear_congruence(B, M);
  if(ans.second == -1) {
    cout << -1 << endl;
  }
  else {
    if(ans.first == 0) {
      cout << ans.second << endl;
    }
    else {
      cout << ans.first << endl;
    }
  }


  return 0;
}
0